莫兰指数stata命令_基于Stata的ESDA（探索性空间数据分析）代码解析与实操验证...

一、引言

探索性空间数据分析(Exploratory Spatial Data Analysis，ESDA)是空间计量经济学的重要技术。

网上调研发现，鲜有基于Herrera Gomez(2016)的方法，使用中国本土数据，比较全面地实现ESDA的文章。本文可能是中国第一篇基于中国本土数据，实现Herrera Gomez(2016)技术，并对其进行深度解读的文章。

本文由五大部分组成：第一部分为引言。第二部分为空间权重矩阵的生成。第三部分为Univariate Spatial Tests，包括：使用不同软件和命令计算单变量全域莫兰指数，绘制单变量莫兰散点图；热点分析(Hot Spot Analysis)的实现。第四部分为Choropleth Map，包括：绘制Quantile Map、Equal Intervals Map、Box Map、Deviation Map及Combine Map。

二、生成空间权重矩阵

ASEF已经推送过，很多关于如何生成各种空间权重矩阵的文章或视频。感兴趣的朋友，可以在ASEF公众号的“文章分类”中进行查找。在这里，不再演示如何生成空间权重矩阵。

需要说明的是：本文所有演示采用的矩阵均为距离倒数空间权重矩阵(inverse of the squared distance weight matrix)。

三、Univariate Spatial Tests

(一)计算Moran's I

1.计算全域Moran's I

. use winv.dta,clear

. spatwmat using winv.dta,name(winv)standardize

. use ti&ts_34variables.dta,clear

.spatgsati_2002 ti_2003 ti_2004 ti_2005 ti_2006 ti_2007 ti_2008 ti_2009 ti_2010 ti_2011 ti_2012ti_2013 ti_2014 ti_2015 ti_2016 ti_2017 ti_2018 ,weights(winv) moran，见图1