查看原文
其他

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

张家星 星8氪 2021-09-21


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



人穷莫入众,言轻莫劝人



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文件即可



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



: . Video Mini Program Like ,轻点两下取消赞 Wow ,轻点两下取消在看

您可能也对以下帖子感兴趣

文章有问题?点此查看未经处理的缓存