自行整理的空间面板计量分析步骤.pdf
文本预览下载声明
空间面板计量分析步骤
空间面板计量分析步骤
(由于网上的资料操作性不高,作者自行整理空间面板计量分
析步骤,集中与操作,涉及ArcGIS10.2 、Stata13 、
Geoda1.6.7 、Matlab2010a 不讨论理论,由于作者对计量的认识
不一定充分,步骤可能有误或者未尽完善,希望能够对此问题
有一个正确的系统认识,现展示操作步骤,请各位大神指导。
如有错误或简化的步骤,望留言和联系,谢谢。其他同行如转
载请注明出处,作者qq :157821891)
一、总思路
(一)确定有否必要使用空间计量,即使用一般面板计量模
型还是空空间面板计量模型 (空间自相关分析)
(二)确定空间滞后模型还是空间误差模型,还是其他模型,
还是不需要用到空间模型 (LM)
(三)确定随机效应还是固定效应 (Haueman)
(四)确定个体固定、时间固定还是混合固定 (LR)
空间面板计量分析步骤
一、ArcGIS 10.2 部分
1、收集数据,形成excel 总表
t 为年份,BSM 为分析地区的唯一ID ;
2 、打开ArcMap,加载“XZQ”图层
3 、新建字段BSM (字段类型为long),按照excel 的BSM 标注每
个地区的唯一ID
4 、运用图层导出,分别导出n 个“XZQ”图层,n 为年份数,如果分
析的时间跨度为10 年,则n=10 ,分别命名为XZQ2007 、
XZQ2008 、XZQ2009……
5 、逐个右键“XZQ”图层,按照年份运用join 连接上述excel 表,连
接字段为BSM
6 、连接后,分别右键导出图层,重新命名为“XZQLJ2007”、
“XZQLJ2008” 、“XZQLJ2009”……
空间面板计量分析步骤
二、GEODA1.6.7 部分
1、安装GEODA1.6.7
2 、导入每一年的XZQLJ ,以XZQLJ2007 为例
空间面板计量分析步骤
3 、Tools-Weights-Create,创建空间权重矩阵
(1)Weights File ID Variable 选择唯一ID 的BSM ;
(2 )连接方式采取“Queen”,表示分析单元的图斑相互邻接
(通常为面状),如果分析的图斑不相互邻接(点状),则用
“Rook” ;
(3 )生成的权重为对应SHP 文件的GAL 文件,例如
XZQLJ2007.gal
4 、进行自相关检验
(1)全域自相关检验
空间面板计量分析步骤
点击面板上对应图标展开菜单,点击UniVariate Moran’s I ,弹
出面板,选择需要进行检验的因变量,本例为LnZons ,得到该年份
的散点图和Moran’s I 值。查看了相关论文,一般为对全局Moran’s I
值进行逐年的趋势分析,得到其趋势,此处不详细说明。
Morans I of lnZCONS Morans I of lnCXCONS
0.35
0.33
0.31
0.29
0.27
0.25
2007 2008 2009 2010 2011 2012 2013 2014 2015 2016
空间面板计量分析步骤
(2 )局部自相关检验
点击面板上对应图标展开菜单,点击UniVariate Local Moran’s
I ,同样选择待分析的变量,结果有三种结果图式,可根据相关理论
选择。
空间面板计量分析步骤
(3 )回归分析。考虑GeoDA 主要为对空间截面数据的处理,
如果是空间面板,Geoda 无法实现,将结合Stata 和Matlab 进行。
三、STATA 部分
1、安装一
显示全部