全局莫兰指数、局部莫兰指数、莫兰散点图的stata实现
不吃鱼的小小熊
2025年05月13日 21:32

(一)、莫兰指数是什么

莫兰指数(Moran's I) 是衡量空间自相关(Spatial Autocorrelation)的经典统计量,用于检验地理空间中某变量的观测值是否呈现聚集或分散模式。简单来说,它回答的问题是:邻近区域的数值是相似(聚集)还是差异显著(分散)?

1. 核心概念

  • 空间自相关

  • 指地理上邻近的区域在某个变量(如人均收入、犯罪率、气温等)上是否存在相关性。

    • 正相关:高值靠近高值,低值靠近低值(如富裕区聚集、贫困区聚集)。

    • 负相关:高值被低值包围,或反之(如城市中心高房价与周边低房价交替)。

    • 零相关:空间分布完全随机。

  • 莫兰指数的范围

    • 理论值介于 -1 到 1 之间:

      • I > 0:正自相关(相似值聚集)。

      • I ≈ 0:无空间模式(随机分布)。

      • I < 0:负自相关(相异值分散)。

    • 实际计算中,I 可能略微超出 [-1,1] 范围。

2. 莫兰指数的类型

  • 全局莫兰指数(Global Moran's I)

  • 评估整个研究区域内变量的整体空间模式。

    • 示例问题:全国各城市的GDP是否存在空间聚集?

  • 局部莫兰指数(Local Moran's I, LISA)

  • 识别局部区域的聚集或异常点(如热点、冷点)。

    • 示例问题:上海市的某个区是否与周边区形成高-高聚集?

3. 数学公式(全局莫兰指数)算了吧放过自己,会用就行了

4.空间计量矩阵的选择

(1)地理邻接空间权重矩阵

(2)反距离空间权重矩阵

(3)反距离平方空间权重矩阵

(4)经济距离平方倒数矩阵

5. 如何判断显著性?

(1)Moran's I值:介于 -1 到 1 之间。

>0:正自相关(相似值聚集)。

≈0:无空间自相关(随机分布)。

<0:负自相关(相异值聚集)。

(2)Expected I:零假设(无空间自相关)下的期望值,通常接近 -1/(n-1)。

(3)z-score & p-value:

若 p < 0.05,拒绝零假设,表明存在显著空间自相关。

(二)、具体操作方法

1、全局莫兰指数的计算

(1)给出空间权重矩阵的路径

use "C:\Users\11498\OneDrive\硕士研究生\002项目\25-4-22国家重大大健康产业\莫兰指数问题\地理邻接空间权重矩阵.dta"

(2)给出待估计指标的路径

use "C:\Users\11498\OneDrive\硕士研究生\002项目\25-4-22国家重大大健康产业\莫兰指数问题\人口高质量发展水平02年~17年(有地名的).dta"

(3)标准化空间权重矩阵

spatwmat using"C:\Users\11498\OneDrive\硕士研究生\002项目\25-4-22国家重大大健康产业\莫兰指数问题\地理邻接空间权重矩阵.dta",name(W) standardize

(4)计算全局莫兰指数

 spatgsa year2002 year2003 year2004 year2005 year2006 year2007 year2008 year2009 year2010 year2011 year2012 year2013 year2014 year2015 year2016 year2017,weights(W) moran

注意:

(a)spatwmat spatgsa为stata的命令包需要提前安装,b站有教程(推荐up主:吾系奥特曼)

(b) year2002为变量

2、局部莫兰指数的计算

spatlsa year2002,weights(W) moran

3、莫兰散点图的绘制

 spatlsa year2002,weights(W) moran graph (moran) id(provience) symbol(id)

注意:

(a) spatwmat spatgsa为stata的命令包需要提前安装,b站有教程(推荐up主:吾系奥特曼)