引言: 权重矩阵,是做空间面板回归、空间界面回归必须的文件
从陈强第二版的ussaww.spmat文件(做空间面板回归的)[url=],可以看到,[/url][x1] 是一个48个城市的行正规化的文件。
我们要做回归,50%的工作是要生成这样一个权重矩阵。
理论上 ,3种的计算模式,有临近、经济距离和距离倒数(我用的这个,其他两个做不出);
我的素材是,选了特定45个城市,希望做空间面板,首先想生成这样一个权重矩阵文件.
需要的数据是45个城市的X\Y坐标(百度 全国各城市经纬度坐标就找到);
需要的软件有:arcgis10.0; stata15; GEODA; EXCEL
视频教程软件: 华东理工大学,空间面板公开课。
http://mooc1.chaoxing.com/course/101569745.html?edit=false&knowledgeId=101614908&module=2#content
我实现的方法如下:
1 生成shp文件;by arcgis—— 必经之路。
adddata,导入经纬度坐标文件。 xx.xls (excel要是2003版本的文件簿文件最好,这样能识别)
显示坐标, x ,y已经自动的导入
添加投影仪坐标 wgs84( 如果不添加,后期添加到stata中,stata不自动识别经纬度坐标)
导出shp文件
这样我们就拿到了shp文件 ——全部放到stata的根目录下
# 此处 第一个目标,形成地理坐标的shp文件生成;
2、用stata
2、1
spshape2dta Export_Output_1,replace
# 把 Export_Output_1整体转换成stata的文件dta( 要能创建 CX,CYcoordinate,这样才能成功,这样才能dta 带有stata自身的地理坐标格式)
2.2
spset _ID,modify replace
# 目的是,设置成空间数据, _ID 是 北京—xx城市的序号名词( LONG 数据类型);
#错误解决1: 已经set了; 那就 spset,clear 清除掉,再set
# 错误解决2: 有可能spset _ID 就能set 了。
这是我们要出现的结果: coordinate 那里要有 _CX ,CY (此时默认的话,stata这两列没有数据,——2.3步我们自己补上)
2.3 统一说一下原始数据,stata的数据
我们要把X Y 数据 自己复制到,CXCY里面。
3、生成矩阵
结果如下。
4、反思路径
4.1 第三步,为什么不能用邻居矩阵?
显示没有发现neighbor,我觉得用GEODA软件的那个邻居可能有。
一开始我直接用xls 导入的 GEODA,
在这个 4‘ 的地方, 有个最近的neighbor,用xls 生成的文件我看不懂。
可能用shp是另外的生成方式,。
并且GEODA 导出的文件——还通过某一个步骤,转到stata 就能直接用
。我想突破了这个转换的步骤,by用shp; xls , 就能生成新的权重矩阵—— 临近的。
#
[x1]查看方法,
先把 usaww.spmat (陈强空间面板中的 权重矩阵)
导入stata
前置的操作是,spmat的命令导入
(help spmat)——找到该命令,然后click here install 再用即可。)
使用文件: spmat use usaww usingusaww.spmat
导出成txt: spmat export usaww usingdaochu.txt