由于增加一个维度进行因果识别的难度加大,空间计量经济学经历短暂的停滞之后有所起色,近期关于空间计量经济学的理论建模有所发展。
1.用ArcMap切割地图shp相关数据:
首先需分割出个包含空间单元的底图框架做各种空间权重矩阵的基础,这个还是用到ArcGIS的ArcMap,下面介绍一种目前最快的切割shp相关文件的操作:
首先打开中国地级地区的shp地图,起动Edit,打开属性表:

在属性数据表选中需要删除的地区:

关闭属性表再删除需要删除的地区:


剩下的新长三角地区输出为shp文件,就可以用stata空间计量模块处理了:

比某种需要join数据的快很多倍的时间!
提示:可视化地图不适合中文论文,但是做数据框需要这个地图shp文件。
2.stata空间计量命令转换为stata格式数据:
ESRI文件是最常用的。这种格式,一个地图一般包括三个文件:.shp形状文件,.dbf数据库文件和.shx索引文件。我们只需要.shp和.dbf文件。可用shp2dta命令将这些文件转换成两个关联的.dta文件。安装以下三个命令
ssc install spmap
ssc install shp2dta
ssc install mif2dta
//http://blog.sina.com.cn/s/blog_4b86ce47010005xo.html
cd E:\stata_ref\statamap数据\stata作空间分布图\新长三角地区shp
shp2dta using newyangziriverdelta, database(newyangziriverdelta_db) coordinates(newyangziriverdelta_coord) genid(id)
这样得到两个stata数据文件之一,这个数据可以继续添加数据,随后的非洲数据便这样:

刚导出的非洲地图空间数据:

合并后的新非洲地理空间数据:

stata软件制作非洲地图:
spmap odis using "africa5_coord.dta", id(id) ///
ndfcolor(red) ///
title("中国对非洲国家投资存量", size(*0.8)) ///
subtitle("(in US$ million,by the end of 2012)", size(*0.8)) ///
legstyle(2) legend(region(lcolor(black)))
关于做动态图,可进一步参考下面的网文:Stata:空间计量之用 spmap 绘制地图
3.stata空间计量操作指南,一个经管之家的帖子,对陈强教授的内容进行了扩充
汇总在文末:阅读原文,还包括空间权重矩阵的一个特别处理。
4.stata空间面板模型的操作参考以下推文:
5.空间数据的内生性问题处理有了新命令spivregress,stata内help,或参考:
空间计量百科全书式的使用指南的do file公开, 暨空间计量研究小组成立100天
参考文献:
区域经济研究方法Master Program:根据经纬度Stata制空间权重矩阵代码数据
基于ArcGlobe的全球移动电话,中国、美国、“一带一路”沿线夜间灯光可视化分布
空间计量及Stata、ArcGIS及Matlab、GeoDa软件应用
动态空间面板模型教程(一文读懂动态面板空间spregdpd操作应用)
一文读懂空间计量及stata应用(附lr检验、空间可视化、权重矩阵、检验、模型dofile等)
空间计量及应用SPATIAL AUTOREGRESSIVE MODELS USING STATA(一)
根据经纬度Stata制空间权重矩阵代码数据:附空间计量经济学软件代码与参考文献





