实用

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

Logo

有“AI”的1024 = 2048,欢迎大家加入2048 AI社区

更多推荐