空间面板回归——建立权重矩阵的方法-距离倒数的权重矩阵-经管之家官网！

免费学术公开课,扫码加入

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

从陈强第二版的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的根目录下

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image005.png

# 此处 第一个目标，形成地理坐标的shp文件生成；

2、用stata

2、1

spshape2dta Export_Output_1,replace

# 把 Export_Output_1整体转换成stata的文件dta（ 要能创建 CX,CYcoordinate，这样才能成功，这样才能dta 带有stata自身的地理坐标格式）

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image006.png

2.2

spset _ID,modify replace

# 目的是，设置成空间数据，_ID 是 北京—xx城市的序号名词（ LONG 数据类型）；

#错误解决1： 已经set了； 那就 spset,clear 清除掉，再set

# 错误解决2： 有可能spset _ID 就能set 了。

这是我们要出现的结果： coordinate 那里要有 _CX ,CY （此时默认的话，stata这两列没有数据，——2.3步我们自己补上）

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image007.png

2.3 统一说一下原始数据，stata的数据

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image009.jpg

我们要把X Y 数据 自己复制到，CXCY里面。

3、生成矩阵

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image011.jpg

结果如下。

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image013.jpg

4、反思路径

4.1 第三步，为什么不能用邻居矩阵？

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image014.png

显示没有发现neighbor，我觉得用GEODA软件的那个邻居可能有。

一开始我直接用xls 导入的 GEODA,

http://file:///C:/Users/xi/AppData/Local/Temp/msohtmlclip1/01/clip_image016.jpg

在这个 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

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP