stata根据经纬度生成空间权重矩阵

实用

clear all
cd C:\desk\2019全国行政区划
 use 284个城市经纬度.dta, clear
   //（1）计算01矩阵，此方法不靠谱，建议gis加geoda
 spatwmat, name(W01) xcoord(lon1) ycoord(lat1) band(0 12) binary 
 //加standardize为标准化
 matlist W01
 mat2txt,matrix(W01) saving(284)replace  //矩阵保存为txt文件
dataout using 284.txt,word excel  //矩阵转换为word和excel文件，省级矩阵
 //（2）反距离矩阵和地理矩阵平方的好方法
  clear all
cd C:\desk\2019全国行政区划
    use 284个城市经纬度.dta, clear
  spwmatrix gecon lat1 lon1, wname(w02) wtype(inv)
          matlist w02
mat2txt,matrix(w02) saving(284距离)replace  //矩阵保存为txt文件
 dataout using 284距离.txt,word excel  //矩阵转换为word和excel文件
 mat w1=hadamard(w02,w02)  //距离倒数的平方
 matlist w1
mat2txt,matrix(w1) saving(284平方)replace  //矩阵保存为txt文件
dataout using 284平方.txt,word excel  //矩阵转换为word和excel文件

实践过得代码

cd C:\Users\Admin\Desktop\2022-1115
use"经纬度数据.dta" ,replace //你的dta文件，注意要带经纬度
spatwmat, name(W01) xcoord(x) ycoord(y) band(0 12) binary standardize 
matrix list W01 //列出矩阵
matrix A = W01
clear
svmat A, names(w) 

mat2txt,matrix(W01) saving(100.txt) replace //矩阵保存为txt文件。mat2txt需要安装
dataout using 100.txt,word excel //矩阵转换为word和excel文件，保存在默认的存储文件夹下

网络资源

clear all
cd D:\statasp\stata制空间权重矩阵代码数据
   set matsize 11000
   use county_jingweidu.dta, clear
   //（1）计算01矩阵
 spatwmat, name(W01) xcoord(x) ycoord(y) band(0 12) binary standardize
 matlist W01
 mat2txt,matrix(W01) saving(3180)replace  //矩阵保存为txt文件
dataout using 3180.txt,word excel  //矩阵转换为word和excel文件
 //（2）反距离矩阵和地理矩阵平方的好方法
  clear all
 // ssc install spwmatrix
cd D:\statasp\stata制空间权重矩阵代码数据
   use county_jingweidu.dta, clear
  spwmatrix gecon y x, wname(winvsq) wtype(inv)
          matlist winvsq
mat2txt,matrix(winvsq) saving(3180d)replace  //矩阵保存为txt文件,县级太多，跑不出来，考虑减少样本
 dataout using 3180d.txt,word excel  //矩阵转换为word和excel文件
 mat w1=hadamard(winvsq,winvsq)  //距离倒数的平方
 matlist w1
mat2txt,matrix(w1) saving(3180d)replace  //矩阵保存为txt文件
dataout using 3180.txt,word excel  //矩阵转换为word和excel文件



clear all
cd D:\statasp
   use city_jingweidu.dta, clear
   //（1）计算01矩阵
 spatwmat, name(W01) xcoord(x) ycoord(y) band(0 12) binary standardize
 matlist W01
 mat2txt,matrix(W01) saving(310)replace  //矩阵保存为txt文件
dataout using 310.txt,word excel  //矩阵转换为word和excel文件，地级城市
 //（2）反距离矩阵和地理矩阵平方的好方法
  clear all
 // ssc install spwmatrix
cd D:\statasp
   use city_jingweidu.dta, clear
  spwmatrix gecon y x, wname(winvsq) wtype(inv)
          matlist winvsq
mat2txt,matrix(winvsq) saving(310d)replace  //矩阵保存为txt文件
 dataout using 310d.txt,word excel  //矩阵转换为word和excel文件
 mat w1=hadamard(winvsq,winvsq)  //距离倒数的平方
 matlist w1
mat2txt,matrix(w1) saving(310d)replace  //矩阵保存为txt文件
dataout using 310.txt,word excel  //矩阵转换为word和excel文件


clear all
cd D:\statasp
   use jingweidu.dta, clear
   （1）计算01矩阵
 spatwmat, name(W01) xcoord(x) ycoord(y) band(0 12) binary standardize
 matlist W01
 mat2txt,matrix(W01) saving(100)replace  //矩阵保存为txt文件
dataout using 100.txt,word excel  //矩阵转换为word和excel文件，省级矩阵
 （2）反距离矩阵和地理矩阵平方的好方法
  clear all
  ssc install spwmatrix
cd D:\statasp
   use jingweidu.dta, clear
  spwmatrix gecon y x, wname(winvsq) wtype(inv)
          matlist winvsq
mat2txt,matrix(winvsq) saving(111)replace  //矩阵保存为txt文件
 dataout using 111.txt,word excel  //矩阵转换为word和excel文件
 mat w1=hadamard(winvsq,winvsq)  //距离倒数的平方
 matlist w1
mat2txt,matrix(w1) saving(130)replace  //矩阵保存为txt文件
dataout using 130.txt,word excel  //矩阵转换为word和excel文件