面板空间计量之 Stata 应用:学习笔记
【同 舟 共 济】
更新于 2016 年 4 月 14 日
说 明
目前,在空间计量方面,Stata 官方命令语句数量有限且较为零散,尚未形成系统的空
间计量工具包。因此,个人建议空间计量的初学者转向 Matlab 软件,James P. LeSage、J. P.
Elhorst、Donald J. Lacombe 等学者所开发的空间计量工具包,其功能相对更加完善,操作起
来也比较方便。本人已经习惯了使用 stata,初次自学空间计量方面的操作,参考 help 文件
及相关文献,在学习过程中做了简要总结,仅供初学者交流学习。其中若有不当之处,敬请
批评指正,谢谢!
E-mail: ares0825@163.com
【Stata】
Abd Elmessih Shehata (Econpapers)
URL: http://econpapers.repec.org/RAS/psh494.htm
Federico Belotti (Econpapers)
URL: http://econpapers.repec.org/RAS/pbe427.htm
P. Wilner Jeanty (Econpapers)
URL:http://econpapers.repec.org/RAS/pje95.htm
Maurizio Pisati
URL:http://www.iq.harvard.edu/people/maurizio-pisati
Yihua Yu (Econpapers)
URL:http://econpapers.repec.org/RAS/pyu79.htm
【Matlab】
James P. LeSage - jplv7
URL: http://www.spatial-econometrics.com/
J. P. Elhorst - Matlab Software to estimate spatial panels
URL: http://www.regroningen.nl/elhorst/software.shtml
Donald J. Lacombe's - MATLAB Code for Spatial Econometrics
URL: http://community.wvu.edu/~djlacombe/matlab.html
Kelley Pace & James P. LeSage - Spatial Statistics Toolbox
URL:http://www.spatial-statistics.com/
目 录
第一章 Stata 空间计量命令语句安装 1
第二章 中国 31 省市自治区(不含港澳台、附属岛屿)shp 制作 3
第三章 Stata 空间权重制作 8
第四章 Stata 空间相关性检验 27
第五章 Stata 空间面板数据回归 39
面板空间计量之 Stata 应用:学习笔记
第一章 Stata 空间计量命令包安装
更新于 2016-03-15
1. 空间计量-Stata 命令包
Archive of user-written Stata packages
URL: http://www.haghish.com/statistics/stata-blog/stata-programming/ssc_stata_package_list.php
另外,在 IDEAS(URL: https://ideas.repec.org/)中可以查询相关命令,顺便推荐几个论坛,大家可以经常逛
图 1 Stata 用户自拟命令语句列表
逛:
Stata 官方论坛 URL: http://www.statalist.org/
UCLA-Idre 论坛 URL: http://www.ats.ucla.edu/stat/stata/
Stata Daily URL: http://statadaily.com/index/
2. 安装
单击图 1 左侧红色框内命令名称,即可下载对应的压缩包,安装过程参考非官方命令手动安装说明(URL:
http://bbs.pinggu.org/thread-2420580-1-1.html);单击图 1 右侧蓝色框内的各命令所对应的描述性语句,即可看到
该命令的详细说明及应用举例。或者直接在 stata 输入命令语句:findit spatial,即可安装相关命令。若 stata 软件
安装在系统 C 盘,请以管理员身份运行软件,否则可能会出现无法写入的情况。
3. 注意
需要说明的是,Archive of user-written Stata packages 中所示的 stata 语句均为非官方命令,大家在使用过程
中,请仔细阅读对应的帮助文件等资料,有能力的同学研习一下 ado 文件的源代码,务必保证正确使用命令语句。
另外,大家最好 Google 一下其它相关资料看看,如意大利学者 Federico Belotti 等开发的“xsmle”命令(图 2),一
定要先搞清楚基本原理与适用范围,才可以使用,切记!
1
第一章 Stata 空间计量命令包安装
图 2 xsmle 命令资料查找
2
面板空间计量之 Stata 应用:学习笔记
第二章 中国 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 所示;
3
第二章 中国 31 省市地区(不含港澳台、附属岛屿)shp 制作
图 2 上海市- ArcMap 窗口示意图
2.1 ArcGIS 软件图斑处理命令
此处不再详细阐述,命令说明及操作方法请参考:
2.1.1 ArcGIS 合并多边形的方法
URL: http://blog.sciencenet.cn/home.php?mod=space&uid=74956&do=blog&id=387013
2.1.2 ArcGIS 中各种合并要素的异同点分析
URL: http://blog.163.com/zwx_gis/blog/static/32434435201192310545948/
2.1.3 ArcGIS 问题:shp 文件、属性查询、图斑合并
URL:http://blog.sina.com.cn/s/blog_76cf3614010176wo.html
2.1.4 ArcToolbox 使用之二——Data Managerment Tools
URL: http://blog.163.com/lele..tiantian/blog/static/60396264201161810164156/
2.1.5 ArcGIS 消除零碎图斑(Eliminate)
URL:http://wenku.baidu.com/link?url=a1oIUZPNwsobmoueWeikYhPxTRbH21kARWlfkfBYdkBcGQhCd8B54TrdP
h86-AHa_aI2d9OmhOunzouEOKvuBtbPh-XpZE3elpPGrRnq2Hy
关于 ArcGIS 软件的更多操作,还可参考:
2.1.6 Arcgis 初学者使用教程(人大经济论坛)
URL: http://bbs.pinggu.org/thread-3058027-1-1.html
2.1.7 Arcgis 学习资源帖(人大经济论坛)
URL: http://bbs.pinggu.org/thread-1422034-1-1.html
2.1.8 ARCGIS 地理信息系统空间分析实验教程(人大经济论坛)
URL: http://bbs.pinggu.org/thread-1084442-1-1.html
4
面板空间计量之 Stata 应用:学习笔记
2.2 ArcGIS 软件:中国 31 省市地区(不含港澳台、附属岛屿)制作说明
2.2.1 Prov_Prepare
说明:在 Bou2_4p.shp 的基础上,删除(Delete)香港、澳门和台湾地区,属性表 814 条数据。
图 3 Prov_Prepare --ArcMap 窗口示意图
2.2.2 Prov_Delete
说明:在 Prov_Prepare.shp 的基础上,继续删除(Delete)各行政区的附属岛屿,属性表 31 条数据,仅保留面积最
大图斑的属性值。如图 4 所示。注意:河北省内陆板块仅保留了面积大的图斑,红圈内面积较小的内陆板块在地图
中消失。
图 4 Prov_Delete --ArcMap 窗口示意图
5