空間計量模型R操作代碼

shengv-readOGR(".","shengditu")OGRdatasourcewithdriver:esriShapefileSource:".",layer:"shengditu"with30featuresand7fieldsFeaturetype:wkbPolygonwith2dimensionssheng_nbv-poly2nb(sheng)sheng_nblistv-nb2listw(sheng_nb)Yv-log(sheng$lNCOME)x1v-log(sheng$CAiZHENG)x2<-log(sheng$YONGDlAN)M0=lm(sheng$lNCOME~sheng$CAlZHENG+sheng$YONGDlAN)M0Call:lm(formula=sheng$lNCOME~sheng$CAlZHENG+sheng$YONGDlAN)Coefficients:(Intercept)sheng$CAlZHENGsheng$YONGDlAN2242.77621.94090.8433>bptest(M0)studentizedBreusch-Pagantestdata:M0BP=13.3794,df=2,p-value=0.001244Ml=lm(Y~Xl+x2)MlCall:lm(formula=Y~X1+X2)Coefficients:(Intercept)X1X25.90630.14100.2403>bptest(M1)studentizedBreusch-Pagantestdata:MlBP=1.0727,df=2,p-value=0.5849空間同步自回歸誤差模型M5SEARv-errorsarlm(form=Y~x1+x2,listw=sheng_nblist,etype="error")>M5SEARCall:errorsarlm(formula=Y~X1+X2,listw=sheng_nblist,etype="error")Type:errorCoefficients:lambda(Intercept)X1X20.68037056.24397140.16982980.1567685Loglikelihood:16.28129空間杜賓誤差模型M6SDEMv-errorsarlm(formula=Y~x1+x2,listw=sheng_nblist,etype="emixed")>M6SDEMCall:errorsarlm(formula=Y~X1+X2,listw=sheng_nblist,etype="emixed")Type:errorCoefficients:lambda(Intercept)X1x2lag.X1lag.X20.3839413295.6064438950.1749807530.1252941500.1271285720.005312641Loglikelihood:18.43699空間自回歸滯后模型M7SAR<-lagsarlm(Y~x1+x2,listw=sheng_nblist,type="lag")M7SARCall:lagsarlm(formula=Y~X1+X2,listw=sheng_nblist,type="lag")Type:lagCoefficients:rho(Intercept)X1X20.49771382.51428940.15982730.1208989Loglikelihood:19.10202空間杜賓模型M8SDMv-lagsarlm(Y~x1+x2,list=sheng_nblist,type="mixed")M8SDMCall:lagsarlm(formula=Y~X1+X2,listw=sheng_nblist,type="mixed")Type:mixedCoefficients:rho(Intercept)x1x2lag.x10.447347782.941688440.157339440.125535960.03462418lag.x2-0.05020476Loglikelihood:19.35204一般空間模型SACSARARM9SACv-sacsarlm(Y~x1+x2,listw=sheng_nblist,type="sac")>M9SACCall:sacsarlm(formula=Y~x1+x2,listw=sheng_nblist,type="sac")Type:sacCoefficients:rholambda(Intercept)x1x20.6237129-0.38677861.68204350.14026770.1090523Loglikelihood:19.50472M9SARMAv-sacsarlm(Y~x1+x2,listw=sheng_nblist,type="sacmixed")>M9SARMACall:sacsarlm(formula=Y~X1+X2,listw=sheng_nblist,type="sacmixed")Type:sacmixedCoefficients:rholambda(Interc