大家好,又见面了,我是你们的朋友全栈君。
之前的博客有介绍过R和Geoda计算莫兰指数的方法,考虑到有时候我们需要自定义空间权重矩阵来计算莫兰指数,那以上两种方法显得有点复杂。所以,今天来分享Stata计算莫兰指数的方法~
本次案例使用的数据为15-19年全国的人均GDP,数据图如下:
Stata中导入数据的方式十分便捷,通常可以分以下两种:
use dta文件地址
(例如,dta文件在D盘,则使用use "D:/data.dta"
即可导入)莫兰指数计算的相关程序包需要预先下载,在命令行输入findit spatwmat
后,点击sg162程序包即可
主要命令如下,注意using后没有引号,standardize表示行标准化
spatwmat using d:/weight.dta,name(W) standardize
权重数据文件需自行构造,示例图为邻接矩阵样例
注意点:
以邻接矩阵W,计算15-19莫兰指数
spatgsa y_2015 y_2016 y_2017 y_2018 y_2019,weights(W) moran
如果想一键将结果生成表格,可使用asdoc+命令
,即可将结果输入word中,如下图所示(asdoc需通过ssc install asdoc
安装)
#熟悉stata编程的朋友,这里可以用循环代替
spatlsa y_2015,weights(W) moran
spatlsa y_2016,weights(W) moran
spatlsa y_2017,weights(W) moran
spatlsa y_2018,weights(W) moran
spatlsa y_2019,weights(W) moran
spatlsa y_2015 ,weight(W) moran id(province) graph (moran) symbol(id)
该图在stata中背景显示是黑底,复制到word中是透明的
findint spatwmat #安装程序包
spatwmat using d:/weight.dta,name(W) standardize #导入权重矩阵
spatgsa y_2015 y_2016 y_2017 y_2018 y_2019,weights(W) moran #计算Global moran'I
spatlsa y_2015,weights(W) moran #计算local moran'I
spatlsa y_2015 ,weight(W) moran id(province) graph (moran) symbol(id) #带汉字的moran'I
发布者:全栈程序员栈长,转载请注明出处:https://javaforall.cn/153076.html原文链接:https://javaforall.cn