面板空间计量之Stata应用

面板空间计量之Stata应用：学习笔记【同舟共济】更新于2016年4月20日说明目前，在空间计量方面，Stata官方命令语句数量有限且较为零散，尚未形成系统的空间计量工具包。因此，个人建议空间计量的初学者转向Matlab软件，JamesP.LeSage、J.P.Elhorst、DonaldJ.Lacombe等学者所开发的空间计量工具包，其功能相对更加完善，操作起来也比较方便。本人已经习惯了使用stata，初次自学空间计量方面的操作，参考help文件及相关文献，在学习过程中做了简要总结，仅供初学者交流学习。其中若有不当之处，敬请批评指正，谢谢！E-mail:[e-mail]【Stata】AbdElmessihShehata（Econpapers）URL:（Econpapers）URL:（Econpapers）URL：：（Econpapers）URL：目录第一章Stata空间计量命令语句安装1第二章中国31省市自治区（不含港澳台、附属岛屿）shp制作3第三章Stata空间权重制作8第四章Stata空间相关性检验27第五章Stata空间面板数据回归39面板空间计量之Stata应用：学习笔记1第一章Stata空间计量命令包安装更新于2016-03-151.空间计量-Stata命令包Archiveofuser-writtenStatapackagesURL:图1Stata用户自拟命令语句列表另外，在IDEAS（URL:）中可以查询相关命令，顺便推荐几个论坛，大家可以经常逛逛：Stata官方论坛URL:论坛URL::安装单击图1左侧红色框内命令名称，即可下载对应的压缩包，安装过程参考非官方命令手动安装说明（URL：）；单击图1右侧蓝色框内的各命令所对应的描述性语句，即可看到该命令的详细说明及应用举例。或者直接在stata输入命令语句：finditspatial，即可安装相关命令。若stata软件安装在系统C盘，请以管理员身份运行软件，否则可能会出现无法写入的情况。3.注意需要说明的是，Archiveofuser-writtenStatapackages中所示的stata语句均为非官方命令，大家在使用过程中，请仔细阅读对应的帮助文件等资料，有能力的同学研习一下ado文件的源代码，务必保证正确使用命令语句。另外，大家最好Google一下其它相关资料看看，如意大利学者FedericoBelotti等开发的“xsmle”命令（图2），一定要先搞清楚基本原理与适用范围，才可以使用，切记！第一章Stata空间计量命令包安装2图2xsmle命令资料查找面板空间计量之Stata应用：学习笔记3第二章中国31省市地区（不含港澳台、附属岛屿）shp制作—Delete、Merge、Dissovle、Eliminate等命令更新于2016-03-18本文以国界与省界地理信息数据包为例，我们主要使用“bou2_4p.shp”，其简要的说明如下：“bou2_4p.shp”中主要包含以下变量：AREA（面积）；PERIMETER（周长）；BOU2_4M_和BOU2_4M_ID（BOU代表边界的意思，数字1~4代表国家、省、市、县的4级行政划分；4M代表地图比例是400万分之一，ID是其对应的编号）；ADCODE93和ADCODE99（行政区域代码）；NAME（行政区域名称）。需要说明的是，在“bou2_4p.shp”中，每一个省市自治区及其附属岛屿都是用一个多边形来表示的，每个多边形都对应唯一的ID，编号分别从1到925（由于部分省份有很多附属岛屿）。一般情况下，我们在省级层面的空间计量实证研究中，仅涉及中国剔除港澳台地区的31个省市自治区的面板数据（有时还需剔除西藏、新疆等地区），因此在具体应用之前，我们需要使用ArcGIS软件对“bou2_4p.shp”进行精简处理，制作中国31省市地区（不含港澳台、附属岛屿）的shp文件，主要处理对象有：河北省：行政区内陆板块分为两部分，且有附属岛屿，数据/图斑数量为8，如图1所示；图1河北省-ArcMap窗口示意图沿海省份：有很多附属岛屿，上海市数据/图斑数量为9，如图2所示；第二章中国31省市地区（不含港澳台、附属岛屿）shp制作4图2上海市-ArcMap窗口示意图2.1ArcGIS软件图斑处理命令此处不再详细阐述，命令说明及操作方法请参考：2.1.1ArcGIS合并多边形的方法URL:=space&uid=74956&do=blog&id=3870132.1.2ArcGIS中各种合并要素的异同点分析URL:问题：shp文件、属性查询、图斑合并URL:使用之二——DataManagermentToolsURL:消除零碎图斑（Eliminate）URL:=a1oIUZPNwsobmoueWeikYhPxTRbH21kARWlfkfBYdkBcGQhCd8B54TrdPh86-AHa_aI2d9OmhOunzouEOKvuBtbPh-XpZE3elpPGrRnq2Hy关于ArcGIS软件的更多操作，还可参考：2.1.6Arcgis初学者使用教程（人大经济论坛）URL:学习资源帖（人大经济论坛）URL:地理信息系统空间分析实验教程（人大经济论坛）URL:面板空间计量之Stata应用：学习笔记52.2ArcGIS软件：中国31省市地区（不含港澳台、附属岛屿）制作说明2.2.1Prov_Prepare说明：在Bou2_4p.shp的基础上，删除（Delete）香港、澳门和台湾地区，属性表814条数据。图3Prov_Prepare--ArcMap窗口示意图2.2.2Prov_Delete说明：在Prov_Prepare.shp的基础上，继续删除（Delete）各行政区的附属岛屿，属性表31条数据，仅保留面积最大图斑的属性值。如图4所示。注意：河北省内陆板块仅保留了面积大的图斑，红圈内面积较小的内陆板块在地图中消失。图4Prov_Delete--ArcMap窗口示意图第二章中国31省市地区（不含港澳台、附属岛屿）shp制作62.2.3Prov_Merge说明：在Prov_Prepare.shp的基础上，设置“小从属大”的原则，把各行政区与其附属岛屿合并（Merge），属性表含31条数据，Area、Lengh等变量值为面积最大图斑的属性值，如图5所示。注意：与Prov_Delete不同的是，河北省内陆小板块及附属岛屿、东南沿海省市的附属岛屿依旧呈现在地图中。图5Prov_Merge--ArcMap窗口示意图2.2.4Prov_Dissovle说明：在Prov_Prepare.shp的基础上，选择Name、ACODE99为关键指标，将31个行政区与其附属岛屿融合（Dissolve），属性表中31条数据。需要特别强调的是：Shape_Area、Shape_Lengh变量值为各行政区所有图斑Area、Lengh值之和，如图6所示。注意：与Prov_Delete不同的是，河北省内陆小板块及附属岛屿、东南沿海省市的附属岛屿依旧呈现在地图中。图6Prov_Dissolve--ArcMap窗口示意图HebeiShape_Area:19.636559Shape_Length:47.404183面板空间计量之Stata应用：学习笔记72.2.5Prov_Eliminate（未成功）说明：在Prov_Prepare.shp的基础上，设置面积阙值（AREA0.15），清除（Eliminate）各行政区附属岛屿等零碎的小图斑，属性表中含812条数据。注意：按照面积阙值条件，仅成功清除掉2个零碎图斑，其余781个零碎小图斑依旧在属性表及地图中，暂时未弄明白失败原因，还请大家留言指导，谢谢。图7Prov_Dissolve--ArcMap窗口示意图2.3遗留问题2.3.1Eliminate命令原理及操作失败的原因2.3.2哪种版本的shp文件时标准的问题描述：基于以上不同版本（Prov_Delete、Prov_Merge和Prov_Dissolve）的中国31个省市自治区shp文件，生成空间权重矩阵可能会有所差异，这是否会对实证结果产生影响？如果存在显著的影响，在省级层面的空间计量实证研究中，应该选择哪个shp文件呢？第三章Stata空间权重矩阵制作8第三章Stata空间权重矩阵制作更新于2016-04-133.1空间权重矩阵理论空间计量首要解决关键问题就是空间权重的选择和设置，已有文献中经常提及的空间权重形式有：邻接权重矩阵（0-1）、反距离权重矩阵、经济权重矩阵、社交网络权重矩阵以及经济距离等嵌套矩阵，其详细的定义及空间权重矩阵选择等问题，请参考相关文献，例如：王守坤.空间计量模型中权重矩阵的类型与扩展逻辑:述评与展望[C]//中国空间经济学年会.2012.DubinR.Spatialweights[J].TheSageHandbookofSpatialAnalysis,2008,1:125-158.AnselinL.Spatialeconometrics:methodsandmodels[M].SpringerScience&BusinessMedia,2013.AnselinL,BeraAK.Spatialdependenceinlinearregressionmodelswithanintroductiontospatialeconometrics[J].StatisticsTextbooksandMonographs,1998,155:237-290.3.2Stata空间权重命令集合Spmat--Createandmanagespatial-weightingmatrixobjec