在Stata中可以實現空間面板數據的自相關性檢驗,具體代碼如下:
//*安裝Stata命令:
findit spatreg
net install sg162.pkg
//*空間權重矩陣:
spatwmatusing columbusswm.dta, name(W)
*「name(W)」表示將根據數據集columbusswm.dta生成的空間權重矩陣命令為w
//*空間自相關的指標計算與假設檢驗:
//*被解釋變量crime的全域自相關指標及相應檢驗:
use columbusdata.dta, clear
spatgsacrime, weights(W) moran geary go twotail
*「spatgsa」進行全局空間自相關檢驗,「gsa」表示「global spatial autocorrelation」;「weights(W)」指定空間權重矩陣W;「moran geary go」分別表示計算莫蘭指數I,吉爾裡指數C以及Getic-Ord指數G;「twotail」表示進行雙側檢驗,默認為單側檢驗(即認為只可能存在空間正相關性)
//*被解釋變量crime的局域自相關指標及相應檢驗:
spatlsacrime, weights(W) moran twotail
*「spatlsa」進行局域空間自相關檢驗,「lsa」表示「local spatial autocorrelation」;其餘選擇項與spatgsa類似。
//*根據距離計算空間自相關:
spatcorr crime, bands(0(1)5) xcoord(x) ycoord(y)
*距離從0到5,每一個單位分段,分別計算Moran’s I;spatcorr能夠考察空間自相關指標對距離臨界值d的依賴性
spatcorr crime, bands(0(1)5) xcoord(x) ycoord(y) cumulative
*距離從0到5,每一個單位分段,累計計算Moran’s I
spatcorr crime, bands(0(1)5) xcoord(x) ycoord(y) graph
*作圖