经过借鉴前人成果,结合自己实践,终于完成了,废话不多说:
首先由Geoda生成距离权重矩阵——保存——生成Wd.gwt
打开stata
*安装软件包 ssc install sppack,包含-spmat-, -spreg-, and -spivreg-等命令
*使用-spmat import- 读取.GWT文件
spmat import M using Wd.gwt, geoda
*如果还是想生成一个 Stata 格式的权重矩阵文件,需要从spmat对象释出权重矩阵并将其保存为一般的Mata matrix
*假设其名称为W
*该命令会创建一个名叫W的spmat对象,存储从Geoda导入的权重矩阵信息.
spmat getmatrix M W
*进一步保存为 Stata 格式,数据集中的变量会被命名为x1, x2, x3, ……
getmata (x*)=W
之后会在stata数据编辑器中看到空间距离权重矩阵,但是很多小伙伴发现x1, x2, x3, ……和自己定义的编码是不对应的(换句话说:是矩阵乱的)
对应方法:
Arcgis打开生成空间权重矩阵的.shp面打开属性表,属性表中的FID列的数字是和stata生成的权重矩阵中的x1, x2, x3, ……相对应。之后在同一替换,就会得到你想要的矩阵!


雷达卡







京公网安备 11010802022788号







