自行整理的空间面板计量分析步骤

上传人:共*** 文档编号:87741105 上传时间:2019-04-11 格式:PDF 页数:16 大小:1.54MB
返回 下载 相关 举报
自行整理的空间面板计量分析步骤_第1页
第1页 / 共16页
自行整理的空间面板计量分析步骤_第2页
第2页 / 共16页
自行整理的空间面板计量分析步骤_第3页
第3页 / 共16页
自行整理的空间面板计量分析步骤_第4页
第4页 / 共16页
自行整理的空间面板计量分析步骤_第5页
第5页 / 共16页
点击查看更多>>
资源描述

《自行整理的空间面板计量分析步骤》由会员分享,可在线阅读,更多相关《自行整理的空间面板计量分析步骤(16页珍藏版)》请在金锄头文库上搜索。

1、空间面板计量分析步骤 空间空间面板面板计量计量分析步骤分析步骤 (由于网上的资料操作性不高,作者自行整理空间面板计量分 析步骤,集中与操作,涉及 ArcGIS10.2、Stata13、 Geoda1.6.7、Matlab2010a 不讨论理论,由于作者对计量的认识 不一定充分,步骤可能有误或者未尽完善,希望能够对此问题 有一个正确的系统认识,现展示操作步骤,请各位大神指导。 如有错误或简化的步骤,望留言和联系,谢谢。其他同行如转 载请注明出处,作者 qq:157821891) 一一、总思路总思路 (一一)确定确定有否有否必要使用空间计量,即必要使用空间计量,即使用使用一般面板一般面板计量计量模

2、模 型还是型还是空空间空空间面板计量模型面板计量模型(空间空间自相关分析自相关分析) (二二) 确定确定空间空间滞后模型还是空间误差模型, 还是其他模型滞后模型还是空间误差模型, 还是其他模型, 还是不需要用到空间模型还是不需要用到空间模型(LM) (三三)确定确定随机随机效应还是固定效应效应还是固定效应(Haueman) (四四)确定确定个体个体固定固定、时间固定还是混合固定、时间固定还是混合固定(LR) 空间面板计量分析步骤 一一、ArcGIS 10.2 部分部分 1、收集数据收集数据,形成,形成 excel 总总表表 t 为年份,BSM 为分析地区的唯一 ID; 2、打开打开 ArcMa

3、p,加载,加载“XZQ”图层图层 3、新建新建字段字段 BSM(字段(字段类型为类型为 long) ,按照按照 excel 的的 BSM 标注标注每每 个地区的唯一个地区的唯一 ID 4、运用图层导出,分别导出运用图层导出,分别导出 n 个个“XZQ”图层图层,n 为年份数,如果分为年份数,如果分 析的时间跨度为析的时间跨度为 10 年年,则,则 n=10,分别命名为分别命名为 XZQ2007、 XZQ2008、XZQ2009 5、逐个逐个右键右键“XZQ”图层图层,按照按照年份年份运用运用 join 连接连接上述上述 excel 表,连表,连 接字段为接字段为 BSM 6、连接后,连接后,分

4、别分别右键导出图层,重新命名为右键导出图层,重新命名为“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”,表示分析单元的图斑相互邻接 (通常为面状) ,如果分析的图斑不相互邻接

5、(点状) ,则用 “Rook”; (3)生成的权重为对应 SHP 文件的 GAL 文件,例如 XZQLJ2007.gal 4、进行进行自相关检验自相关检验 (1)全域自相关检验 空间面板计量分析步骤 点击面板上对应图标展开菜单,点击 UniVariate Morans I,弹 出面板,选择需要进行检验的因变量,本例为 LnZons,得到该年份 的散点图和 Morans I 值。查看了相关论文,一般为对全局 Morans I 值进行逐年的趋势分析,得到其趋势,此处不详细说明。 0.25 0.27 0.29 0.31 0.33 0.35 200720082009201020112012201320

6、1420152016 Morans I of lnZCONSMorans I of lnCXCONS 空间面板计量分析步骤 (2)局部自相关检验 点击面板上对应图标展开菜单,点击 UniVariate Local Morans I,同样选择待分析的变量,结果有三种结果图式,可根据相关理论 选择。 空间面板计量分析步骤 (3)回归分析。考虑 GeoDA 主要为对空间截面数据的处理, 如果是空间面板,Geoda 无法实现,将结合 Stata 和 Matlab 进行。 三、三、STATA 部分部分 1、安装一系列的分析插件、安装一系列的分析插件 (1)ssc install spmap/shp2dt

7、a/mif2dta/spweight/spwmatrix (2)search speightxt,点击安装 speightxt 的插件 (3)安装空间回归插件 1)findit spatreg 2)search spregdhp 3)search xsmle(空间面板模型) 空间面板计量分析步骤 2、使用使用 spwmatrix 命令导入命令导入 gal 文件,语法:文件,语法: spwmatrix import using C:UsersAdministratorDesktoptestZ2J2007.gal, wname(Z2J2007) 此阶段生成的矩阵主要供 Matlab 使用,不进行标

8、准化 3、使用、使用 stata 的命令的命令 putexcel 将矩阵导出到将矩阵导出到 excel,语法:,语法: putexcel A1=matrix(Z2J2007) using JUZHEN ID 1 2 3 4 5 6 7 8 9 1 0 1 1 1 1 0 0 0 0 2 1 0 0 1 0 0 0 0 0 3 1 0 0 0 1 0 1 0 1 4 1 1 0 0 1 0 0 1 0 5 1 0 1 1 0 0 1 1 0 6 0 0 0 0 0 0 0 1 0 7 0 0 1 0 1 0 0 0 1 8 0 0 0 1 1 1 0 0 0 9 0 0 1 0 0 0 1 0

9、0 四四、Matlab2010a 部分部分 (一)(一)LM 检验检验 1、下载下载空间计量工具箱空间计量工具箱 jplv7 新建一个工作文件夹,例如在桌面建立一个 TEST2,把 jplv7 所有的内容解压到 TEST 中 空间面板计量分析步骤 2、导入权重导入权重变量和分析变量变量和分析变量 (1)新建变量,命名为 A 把所有变量(首列为因变量,其余为自变量) *注意数据的顺序 *不把变量名称导入 空间面板计量分析步骤 3、保存工作保存工作空间空间 File-Save Workspace As,命名为 TEST2.mat,路径指向 TEST2 文 件夹下 4、设置工作设置工作路径路径 (1

10、)点击左下角的软件开始菜单,点击路径 path (2) 点击 Add with Subsfolders, 把路径指向于文件夹 TEST2 (如 果后面因为 matlab 搜索不到工具而报错,则可尝试一下把所有工具 的文件夹分别指向路径一次,便于 matlab 能够搜索) 空间面板计量分析步骤 把 TEST2 下所有解压出来的 jplv7 文件夹都设置一次路径, 注意 是所有文件夹都要加入 Matlab search path 中, 不然后面的分析会找不 到对应的分析工具而报错。 空间面板计量分析步骤 5、Hausman 检验检验 把数据先视作一般面板数据进行 Hausman 检验,决定是使用随

11、 机效应模型还是固定效应模型,后面看到 LR 检验的工具箱也有空 间面板的 Hausman 检验,但只有对 SLM 模型和 SAC 模型的随机效 应和固定效应的区分,其他模型方法不是太清楚,希望有大神补 全。检验的结果一般为固定效应,所以下面 LM 检验时也是围绕固 定效应,对于随机效应这方面也是并不清晰。 6、空间空间 LM 检验检验 (1) 用 matlab 把 TEST2spatialpaneldemoLMsarsem_panel.m 文 件复制到 TEST2 的目录下,即与 TEST2.mat 相同的目录(主要是为 了 不 影 响 源 文 件 , 以 便 改 错 后 可 以 修 复 )

12、。 打 开 demoLMsarsem_panel.m 空间面板计量分析步骤 (2)加载文件 自己新增变量和权重矩阵后,这两 行是没有用的, 不过可以在对应的的这两个文件分别输入变量和权重, 但我觉得比较不方便,就把他们删掉了,改增加 load TEST.mat (3)修改时间变量、指向因变量和自变量 *T 为总的年数, N 为总的个体数 (区域数) , 对应载入的面板 数据 *修改 y 和 x 所在的列号,例如,如果因变量在第一列,则 y 后面的数字输入为 1;自变量在 2-8 列,则数字对应输入列号,用逗 号分隔。 空间面板计量分析步骤 (4)LM 检验 *分别修改下面几种方法的变量,第一项

13、为因变量,第二项后 面均为自变量,与对应的表格的列相对应 (第一种第一种是最小二乘法估计,含是最小二乘法估计,含 LM 检验检验) (第第二二种种是是对于对于空间空间固定固定滞后模型和空间滞后模型和空间固定固定误差模型误差模型的的空间空间 LM 检验检验和和空间空间 robust 检验,将会根据生成相关值的情况决定是使用空检验,将会根据生成相关值的情况决定是使用空 间滞后模型和空间误差模型,具体可参考相关论文。间滞后模型和空间误差模型,具体可参考相关论文。默认默认情况是个情况是个 空间面板计量分析步骤 体的固定效应,如果要进行时间的固定效应,则把第一行体的固定效应,如果要进行时间的固定效应,则

14、把第一行 model=1 改为改为 model=2,也是修改变量名称即可也是修改变量名称即可) (第第三三种种是是对于对于空间空间混合混合固定固定滞后模型和空间滞后模型和空间混合混合固定固定误差模型误差模型的的 空间空间 LM 检验检验和和空间空间 robust 检验,检验,修改变量名称即可修改变量名称即可) *分别比较各种方法的统计量和显著性,决定使用空间滞后模 型还是空间误差模型 7、LR 检验检验 (1) 用 matlab 把 TEST2spatialpaneldemopanelscompare.m 文件 复制到 TEST2 的目录下。打开 demopanelscompare.m (2)

15、加载文件,加载变量,指向自变量和因变量,此处不再叙 述 空间面板计量分析步骤 (3)LR 检验 以空间固定效应如 LM 检验,需先修改因变量和自变量的名称,此处 要注意 info.model=1 的语句,值为 1 时,代表个体固定效应,值为 2 时,代表时间固定效应,值为 3 时,代表双固定效应,比较统计量和 显著性可以决定使用个体固定、时间固定还是混合固定。 8、对应对应 LM 检验检验中的分析结果,例如分析数据确定了中的分析结果,例如分析数据确定了使用使用空间空间+固固 定定+混合混合+滞后模型,则在滞后模型,则在 LM 检验检验中针对对应的分析结果进行中针对对应的分析结果进行解释解释 即可。即可。

展开阅读全文
相关资源
相关搜索

当前位置:首页 > 大杂烩/其它

电脑版 |金锄头文库版权所有
经营许可证:蜀ICP备13022795号 | 川公网安备 51140202000112号