在使用stata進(jìn)行空間關(guān)系分析時(shí)，moran散點(diǎn)圖的繪制往往是不可避免的。stata也提供了比較方面的程序包spatlsa來幫助我們進(jìn)行moran散點(diǎn)圖的繪制。但是spatlsa繪制的圖形對其進(jìn)行編輯，這是經(jīng)常困擾我們的一個(gè)問題。因?yàn)槲覀兛赡苄枰袛嗌Ⅻc(diǎn)圖內(nèi)每個(gè)象限具體都包含哪些地區(qū)，判斷他們的空間關(guān)系。這在地區(qū)很少的情況下還可以通過肉眼觀察進(jìn)行判斷，但是地區(qū)一多就無能為力了，這就需要我明白moran散點(diǎn)圖的繪制原理，有需要的話進(jìn)行自我繪制與自定義。

moran散點(diǎn)圖橫軸表示標(biāo)準(zhǔn)化的被檢驗(yàn)變量，縱軸則表示其空間滯后項(xiàng)

spatlsa 繪制

use moran.dta, clear
spatwmat using w01n.dta, name(w01n) standardize
spatlsa wuliu2008, weights(w01n) moran graph(moran) symbol(id) id(names2008)

image

命令運(yùn)行完成后圖形為黑底，直接右擊圖片復(fù)制-黏貼到word自動為白底

自我繪制

1. 變量的標(biāo)準(zhǔn)化

egen sd_wuliu2008 = std(wuliu2008)

2. 空間滯后項(xiàng)的生

splagvar sd_wuliu2008, wname(w01n) wfrom(Stata)

3. 散點(diǎn)圖的繪制

scatter wy_sd_wuliu2008 sd_wuliu2008, xline(0, lp(shortdash) lwidth(0.5)) yline(0, lp(shortdash) lwidth(0.5))

image

往期回顧

• R中的顏色設(shè)置(Ⅱ)-自定義顏色排序
• R中的顏色設(shè)置(Ⅰ)
• 【前文速遞】中國能否通過提升全球價(jià)值鏈減少制造業(yè)出口的碳排放？
• ArcGis專題地圖繪制
• 使用QGIS繪制流向圖
• 投入產(chǎn)出表可視化——弦圖
• 投入產(chǎn)出表可視化——桑基圖
• 使用R進(jìn)行空間計(jì)量模型估計(jì)-面板篇
• 空間計(jì)量邊際效應(yīng)可視化
• 使用R進(jìn)行空間計(jì)量模型估計(jì)-橫截面篇
• （反）閾值空間權(quán)重矩陣批量計(jì)算與導(dǎo)出

image