发帖
雷达卡
引言： 权重矩阵，是做空间面板回归、空间界面回归必须的文件
从陈强第二版的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
京公网安备 11010802022788号
