在使用stata进行空间关系分析时，moran散点图的绘制往往是不可避免的。stata也提供了比较方面的程序包spatlsa来帮助我们进行moran散点图的绘制。但是spatlsa绘制的图形对其进行编辑，这是经常困扰我们的一个问题。因为我们可能需要判断散点图内每个象限具体都包含哪些地区，判断他们的空间关系。这在地区很少的情况下还可以通过肉眼观察进行判断，但是地区一多就无能为力了，这就需要我明白moran散点图的绘制原理，有需要的话进行自我绘制与自定义。
moran散点图横轴表示标准化的被检验变量，纵轴则表示其空间滞后项
spatlsa 绘制
use moran.dta, clear
spatwmat using w01n.dta, name(w01n) standardize
spatlsa wuliu2008, weights(w01n) moran graph(moran) symbol(id) id(names2008)
image
命令运行完成后图形为黑底，直接右击图片复制-黏贴到word自动为白底
自我绘制
- 变量的标准化
egen sd_wuliu2008 = std(wuliu2008)
- 空间滞后项的生
splagvar sd_wuliu2008, wname(w01n) wfrom(Stata)
- 散点图的绘制
scatter wy_sd_wuliu2008 sd_wuliu2008, xline(0, lp(shortdash) lwidth(0.5)) yline(0, lp(shortdash) lwidth(0.5))
image
往期回顾
- R中的颜色设置(Ⅱ)-自定义颜色排序
- R中的颜色设置(Ⅰ)
- 【前文速递】中国能否通过提升全球价值链减少制造业出口的碳排放？
- ArcGis专题地图绘制
- 使用QGIS绘制流向图
- 投入产出表可视化——弦图
- 投入产出表可视化——桑基图
- 使用R进行空间计量模型估计-面板篇
- 空间计量边际效应可视化
- 使用R进行空间计量模型估计-横截面篇
- （反）阈值空间权重矩阵批量计算与导出
image
