已买,里边的代码如下:
install.packages("splm", repos = "http://R-Forge.R-project.org")
library(splm)
data(Produc, package = "Ecdat")
> Produc <- Produc[Produc$year %in% 1970:1974, ]
data(usaww)
fm <- log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp
res <- spfeml(fm, data = Produc, listw = mat2listw(usaww), effects = "spfe", method = "eigen")
summary(res)
|