install.package("splm")
library(splm)
data1<-read.csv(file="C:\\Users\\ASUS\\Desktop\\data.spital.csv")
#载入空间权重矩阵
tine2015<-read.csv("wl.txt",header = TRUE,skip = 1)
attach(tine2015)
wlgwt<-read.gwt2nb("wlgwt.GWT",region.id=ID)
######空间自回归模型#####
#个体固定效应,model那就是填是固定效应还是随机效应,lag表示是否引入因变量滞后项,error那就是是否考虑误差项自回归
splm.result0<-spml(log(gop)~log(tine)+log(hwzz)+log(gktt)+log(jttz)+log(jtcy)+log(gdp),data1,indx=c("country","year"),
listw=wldw,model="within",lag=TRUE,spatial.error="none")
summary(splm.result0)
#随机效应
splm.result1<-spml(log(gop)~log(tine)+log(hwzz)+log(gktt)+log(jttz)+log(jtcy)+log(gdp),data1,indx=c("country","year"),
listw=wldw,model="random",lag=TRUE,spatial.error="none")
summary(splm.result1)
#Hausman检验,原假设是接受随机效应模型
test.hausman<-sphtest(splm.result1,splm.result0)
test.hausman
|