莫兰指数stata命令_用R计算moran's I(莫兰指数)【笨办法学空间计量之一】

假设我们将要研究11个省市区从2010年到2016年的某项经济发展指标，原始数据经过数据预处理和特征加权后(这些步骤假定读者已经熟悉，不再赘述)，我们得到如图1所示的数据。

图1.png

接下来的任务是分年度计算莫兰指数，从而检查该经济发展指标的空间相关性。

ARCGIS和Geoda当然是做空间计量的首选，其次stata和R也提供了一些空间计量的包，算个moran's I自然不在话下，不过拿着鼠标点点点怎么能体现Learning the hard way的初衷呢，““你们要进窄门。因为引到灭亡，那门是宽的，路是大的，进去的人也多；引到永生，那门是窄的，路是小的，找着的人也少。”所以我们干脆实实在在地把莫兰指数算出来好啦。

图2.png

回顾一下莫兰指数的定义，其实它是融入了空间拓扑关系的相关系数。要计算莫兰指数，我们只需要经济发展数据 Y 和空间权重矩阵 W 这两个数据，而如果使用Rook权重矩阵的话，W就退化成我们熟知的邻接矩阵。

因此我们先用笨办法，不借助任何其他软件，把空间权重矩阵求出来。假定按照上表的顺序，从1开始对省市区进行编码，则用邻接列表表示的空间邻接关系如下(默认海南与广东、广西接壤)：

1 2

2 3 8 1

3 2

4 5 6

5 8 6 4

6 4 5 7

7 6 9

8 2 5

9 7 10 11

10 9 11

11 9 10

其实用邻接列表已经可以代入公式计算莫兰指数了，但清晰起见，我们还是用两个丑陋的for循环，将其转换为邻接矩阵wmx。

wmx：

[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11]

[1,] 0 1 0 0 0 0 0 0 0 0 0

[2,] 1 0 1 0 0 0 0 1 0 0 0

[3,] 0 1 0 0 0 0 0 0 0 0 0

[4,] 0 0 0 0 1 1 0 0 0 0 0

[5,] 0 0 0 1 0 1 0 1 0 0 0

[6,] 0 0 0 1 1 0 1 0 0 0 0

[7,] 0 0 0 0 0 1 0 0 1 0 0

[8,] 0 1 0 0 1 0 0 0 0 0 0

[9,] 0 0 0 0 0 0 1 0 0 1 1

[10,] 0 0 0 0 0 0 0 0 1 0 1

[11,] 0 0 0 0 0 0 0 0 1 1 0

好了，万事具备，以下这段代码分年度计算莫兰指数

S0=sum(wmx)

mi=rep(0,ncol(sd))

for (t in 1:ncol(sd)) {

y=sd[,t]

yhat=mean(y)

SSR=(sd(y)^2)*10

sum=0

for (i in 1:11) {

for (j in 1:11) {

sum=sum+wmx[i,j]*(y[i]-yhat)*(y[j]-yhat)

}

mi[t]=(sum/SSR)*(11/S0)

}

得到的莫兰指数如下，所有年份均为负数，说明该项经济指标与空间聚集度呈现负相关关系。

图4.png

怎样计算莫兰指数的Z值和P值，未完待续，敬请期待。Q461157910

莫兰指数stata命令_用R计算moran's I(莫兰指数)【笨办法学空间计量之一】相关推荐

莫兰指数stata命令_【第六期】Regional Study 群日报
本期看点: 一.用MATLAB估计SLX模型二.各省能源消费结构的计算问题三.stata中DID共同趋势相关代码讨论四.公式请教五.Stata画莫兰指数散点图的命令六.分析铁路影响的数据查找 ...
笨办法学R编程（1）
在倚天屠龙记中,有一人唤作火工头陀.此人练功不靠心法,只靠模仿他人招式,由外而内,自成一家.练习编程也有如此的法门,不看文字描述,只观察和模仿别人的代码.这样也可以由外而内学会编程.<笨办法学p ...
《笨方法学python》_《笨办法学Python》第46课手记
<笨办法学Python> 第46课手记这节课制作了一个Python的项目骨架,花了我一个晚上和一个早上的时间,原因是我下载的pdf里面只有OX S的命令行,而没有win下的.我为此在知道 ...
倾向得分匹配的stata命令_计量方法的适用条件汇总（二）：倾向得分匹配
独家揭秘: 计量经济学的魅力与激情:陈强老师的高级计量现场班侧记(2019.5.1) 如何学好高级计量:探秘陈强老师的高级计量及Stata现场班(2019.10.1) 接上期推文,本期探讨倾向得分匹配 ...
rstudio查询命令_【R教程】RStudio使用指南
来源:数量经济学回顾:[R教程]R与RStudio简介及安装 1.Rstudio界面简介界面一般分为四个窗口, 包括代码编辑窗口.变量浏览窗口.命令执行窗口(控制台(Console)).文件管理窗 ...
倾向得分匹配的stata命令_培训对工资是否影响显著：倾向得分匹配法(PSM)及stata实现...
第一部分模型背景 1.研究目的 2.基本思想第二部分数据介绍以及语法简介 1.数据介绍 2.语法格式第三部分案例讲解以及stata实现 1.变量介绍以及数据描述性统计 2.倾向匹 ...
linux 命令详解大于号_大于号与小于号_笨办法学Bash Shell编程-基础篇视频课程_Linux视频-51CTO学院...
聪明人下笨功夫.本课程所倡导"笨办法"的核心是: ● 手动输入所有代码,不要复制粘贴! ● 正确地输入所有代码,也包括注释 ● 运行代码并保证产生相同的输出 ● 如果出现了bug, ...
笨办法学python 3 48题_附录练习 8-10 笨办法学Python3
附录练习 8 来回移动 (pushd, popd) 在这个练习中,你将学习如何用 pushd 命令保存你当前的位置然后去到一个新的位置,以及如何用 popd 命令返回之前保存的位置. 55.9.1 跟 ...
笨办法学R编程（2）
本例将介绍R语言中的while循环和if条件.最终用它来解决Project Euler的第二个问题.除了练习之外你还需要了解一些斐波纳契数列的知识.废话不多说了,打开R控制台,跟着输入下面的代码,自行 ...

莫兰指数stata命令_用R计算moran's I(莫兰指数)【笨办法学空间计量之一】

莫兰指数stata命令_用R计算moran's I(莫兰指数)【笨办法学空间计量之一】相关推荐

最新文章

热门文章