第19章 空间面板分析

library(splm)
library(plm)

fm <- gdp ~ pop

fe <- plm(fm, data = hunan_data, effects = "individual",
          model = "within", index = c("code", "year"))

fe_lag <- spml(fm, data = hunan_data, listw = hunan_listw, 
               model = "within", effect = "individual", lag = T, 
               spatial.error = "none", index = c("code", "year"))

fe_err <- update(fe_lag, model = "random")