r语言：空间面板Moran 指数计算_r语言分析不同试验点之间的空间自相关性,计算moran's 指数

Moran指数

衡量空间相关性的一种指标。Moran指数越接近1，空间自相关越明显。判定一定范围内的空间实体相互之间是否存在相关关系，比如：一座座居民楼它们是聚集在一块还是离散分布在各处。

莫兰指数数值分布在[-1，1]，[0，1]说明各地理实体之间存在正相关的关系，[-1，0]之间说明存在负相关的关系，而0值则无相关关系。

r中计算莫兰指数

r语言中spdep包提供了局部和全域莫兰指数计算函数。但是需要注意的是，该函数需要的参数格式为listw，一般来说只能计算截面数据的莫兰指数，无法处理面板数据，在空间计量中受限很大。所以我通过一些修改，实现面板数据的处理，还请大家不吝赐教。

具体代码及过程

我们的数据为27个国家，19年的面板数据，共有16个变量。基础空间权重矩阵为2727的方阵，故不能直接带入r的函数中计算。先定义一个1919的单位矩阵，然后对单位矩阵和基础空间权重矩阵求克罗内克积，得到一个分块对角矩阵。即面板数据莫兰指数可以接受的空间权重矩阵。

// 构造基础空间权重矩阵
p=shapefile("E://地图/Export_Output_2.shp",encoding="UTF-8")
countryname=c("ARE","BGD","BLR","DEU","EGY","GRC","IDN","IND","IRN","IRQ","JPN","KAZ","KHM","KOR","LAO","MMR","MNG","MYS","NLD","PAK","RUS","SAU","SGP","THA"
1
2