stata|生成空间权重矩阵之邻接矩阵（0-1矩阵）

点击上方蓝色字体，关注我们

“

人穷莫入众，言轻莫劝人

”1
经纬度获取

1.网上搜索各城市经纬度坐标，适用于样本量比较少的情况。

2.使用arcgis生成质心坐标

（1）打开shp地图的属性

（2）按照相同的方法生成Y

(3)属性表转化成Excel表

（4）如下根据需要进行更改

其中：pgdp是2006-2018年的平均GDP，X经度，Y纬度

2
Stata生成空间权重矩阵

1.构造0-1矩阵

关于0-1矩阵：接壤赋值为1，不接壤赋值为0，所以得到的0-1矩阵中的数据全是0和1.考虑到仅有经纬度坐标是没有办法判断两个地区是否接壤，所以这里借助geoda软件，生成0-1矩阵再将其导入到stata中，具体如下：

（1）打开geoda,载入shp地图，创建空间权重

（2）添加唯一ID

（3）创建

（4）保存

（5）将.gal文件导入stata中

spwmatrix import using D:\geoda\zhongyuan01.gal,wname(james) xport(zhongyuan001, dat)putexcel set results23,replaceputexcel B02=matrix(james)//没有的命令需要提前安装//常规安装方法是：ssc install xx命令

wname(james)：给矩阵命名

xport(zhongyuan001, dat)：输出成.dat文件,R或者Matlab可直接打开

putexcel set results23,replace：建立一个名字叫 results23 的excel表格，位置在，zhongyuan01.gal 所在的目录下 

putexcel B02=matrix(james)：把矩阵导入到excel中，其中B02代表，excel的B行2列，从B行2列开始输入

（6）结果如下

（6）得到的excel表格

（7）进一步编辑

之所以在B02开始导入，目的就是为了一会手动将ID（1-30）复制到行和列当中

（8）导入stata中生成.dta文件即可

如果本文对您有所帮助，欢迎分享到朋友圈！

感谢