空间面板数据分析——R的splm包（任建辉，暨南大学）ThesplmpackageprovidesmethodsforfittingspatialpaneldatabymaximumlikelihoodandGM.安装R软件及其编辑器Rstudio网址：wwr-project.orgwwrstudio/下载好Rstudio以后，操作都可以Rstudio中完成了，包括命令的编写、命令运行、图形展示，最方便的要数查看数据了。R界面Rstudio界面，形如matlab下面进入正题，了解splm包中的数据、命令及结果展示。所有命令都写在编辑窗口（studio左上区域），可以单独的运行每行命令，也可选取一段一起执行，点run按钮。1、首先，安装splm包并导入，命令如下:intall.packages(“splm”)，选择最近的下载点library(splm)>library(splm)载入需要的程辑包：MASS载入需要的程辑包：nlme载入需要的程辑包：spdep载入需要的程辑包：sp载入需要的程辑包：Matrix载入需要的程辑包：plm载入需要的程辑包：bdsmatrix载入程辑包：‘bdsmatrix’下列对象被屏蔽了from‘package:base’:backsolve载入需要的程辑包：Formula载入需要的程辑包：sandwich载入需要的程辑包：zoo载入程辑包：‘zoo’下列对象被屏蔽了from‘package:base’:as.Date,as.Date.numeric载入需要的程辑包：spam载入需要的程辑包：gridSpamversion0.40-0(2013-09-11)isloaded.Type'help(Spam)'or'demo(spam)'forashortintroductionandoverviewofthispackage.Helpforindividualfunctionsisalsoobtainedbyaddingthesuffix'.spam'tothefunctionname,e.g.'help(chol.spam)'.载入程辑包：‘spam’下列对象被屏蔽了from‘package:bdsmatrix’:backsolve下列对象被屏蔽了from‘package:base’:backsolve,forwardsolve载入需要的程辑包：ibdreg载入需要的程辑包：car载入需要的程辑包：lmtest载入需要的程辑包：Ecdat载入程辑包：‘Ecdat’下列对象被屏蔽了from‘package:car’:Mroz下列对象被屏蔽了from‘package:nlme’:Gasoline下列对象被屏蔽了from‘package:MASS’:SP500下列对象被屏蔽了from‘package:datasets’:Orange载入需要的程辑包：maxLik载入需要的程辑包：miscToolsPleasecitethe'maxLik'packageas:Henningsen,ArneandToomet,Ott(2011).maxLik:ApackageformaximumlikelihoodestimationinR.ComputationalStatistics26(3),443-458.DOI10.1007/s00180-010-0217-1.Ifyouhavequestions,suggestions,orcommentsregardingthe'maxLik'package,pleaseuseaforumor'tracker'atmaxLik'sR-Forgesite:s://r-forge.r-project.org/projects/maxlik/Warningmessage:程辑包‘Matrix’是用R版本3.0.3来建造的注意：在导入splm时，如果发现还有其他配套的包没有安装，需要先安装。2、接着，查看数据及结构，命令如下：data(Produc,package=”Ecdat”)View(Produc)3、引入空间权重矩阵（spatialweightsmatrix），命令如下data(usaww)Views(usaww)4、空间面板数据的广义矩估计，命令spgmGM<-spgm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp,data=Produc,listw=usaww,moments=”fullweights”,spatial.error=TRUE)summary(GM)>GM<-spgm(log(gsp)~log(pcap)+log(pc)+log(emp)+unemp,data=Produc,+listw=usaww,moments="fullweights",spatial.error=TRUE)>summary(GM)SpatialpanelfixedeffectsGMmodelCall:spgm(formula=log(gsp)~log(pcap)+log(pc)+log(emp)+unemp,data=Produc,listw=usaww,spatial.error=TRUE,moments="fullweights")Residuals:Min.1stQu.Median3rdQu.Max.-0.14000-0.01950-0.003160.015300.16800Estimatedspatialcoefficient,variancecomponentsandtheta:Estimaterho0.3277625sigma^2_v0.0012179Coefficients:EstimateStd.Errort-valuePr(>|t|)log(pcap)-0.00224350.0262646-0.08540.9319295log(pc)0.24149790.023582610.2405<2.2e-16***log(emp)0.78132760.028385527.5256<2.2e-16***unemp-0.00360260.0010094-3.56910.0003582***---Signif.codes:0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’15、空间面板极大似然估计，命令spmlfm<-log(gsp)~log(pcap)+log(pc)+log(emp)+unemp##fixedeffectspanelwithspatialerrorsFespaterr<-spml(fm,...