生成矩阵_如何使用Stata生成空间权重矩阵、计算Moran's I及画莫兰散点图

一、生成空间权重矩阵

1.Create a row-standardized binary spatial weights matrix assuming spherical coordinates and a distance cut-off of 10 miles(生成行标准化的二元空间权重矩阵)

. use 31_provinces_longitude&latitude,clear

.spwmatrix gecon latitude longitude, wn(wbin) wtype(bin) db(0 550) r(3958.761) row

.clear

.svmat wbin

.save wbin

2.Generate an inverse distance squared spatial weights matrix using projected latitudes and longitudes(生成距离倒数空间权重矩阵)

. use 31_provinces_longitude&latitude,clear

. spwmatrix gecon latitude longitude, wname(winvsq) wtype(inv) alpha(1) dband(0 100) cart

. clear

. svmat winvsq

. save winvsq

二、计算Moran's I

(一)计算全域Moran's I

1.使用距离倒数空间权重矩阵

. use winv.dta,clear

. spatwmat us