引言： 权重矩阵，是做空间面板回归、空间界面回归必须的文件

从陈强第二版的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