面板数据（N=31,T=3）的莫兰指数计算MATLAB模拟

   应读者要求，小范编写了面板数据(N=31,T=3，可依据需要调整)的全局、局部莫兰指数MATLAB模拟的代码，请各位同行参考。
   代码基于软件MATLAB R2018a编制。

   全局莫兰指数和局部莫兰指数公式：

   代码运行效果图：

   代码自如下部分开始：
clc;clear all;
% 本代码由重庆科技学院范巧副教授于2018年07月20日17:26编制；% 本代码版权保护，不得用于除学习和研究之外的其他用途；% 31个省份3个年份数据全局莫兰指数和局部莫兰指数Matlab模拟代码（可以依据研究对象进行相应调整）；% 参考文献1：刘贤赵,高长春,张勇,张东水,谢金宁,宋焱,王志强.中国省域碳强度空间依赖格局及其影响因素的空间异质性研究[J].地理科学,2018,38(05):681-690.% 参考文献2：生延超,周玉姣.适宜性人力资本与区域经济协调发展[J].地理研究,2018,37(04):797-813.% 本代码基于MATLAB R2018a软件编写。

% 定义空间权重矩阵的基本形式（根据实际分析进行定义）；N = 31;W0 = randi(2,N)-1;W1 = diag(W0).*eye(N);W2 = W0-W1;

% 空间权重矩阵的标准化；W3 = sum(W2,2);W4 = repmat(W3,1,N);W = W2./W4;

% 检查空间权重矩阵是否已经被标准化；W5 = sum(W,2);for i = 1:N; W6(i) = round(W5(i));end
if round(W6)==ones(N,1); ('Your spatial weight matrix has been row standarded.')else ('Your spatial weight matrix needs to be row standarded.')end

% 录入待分析的数据序列（截面数据或者面板数据）T = 3;Y0 = 100*randn(T,N);Y = abs(Y0);

% 计算T=1年全局莫兰指数；Y1 = Y(1,:);I10 = sum(sum(W));I11 = mean(Y1);I12 = Y1-I11;I13 = sum(I12.^2);I14 = repmat(I12,N,1);I15 = I14';I16 = sum(sum(W.*I14.*I15));Global_Moran_I1 = N*I16/(I10*I13);
if Global_Moran_I1 == 0; ('there is no spatial correlation in the year T=1 !')elseif Global_Moran_I1=Z11)('spatial agglomeration exsits aound region A1 in the year T=1')

% 局部莫兰指数作图；subplot(2,2,2),plot(Z1,'.');axis([1,31,-1,1]);xlabel('T=1');ylabel('Local Moran I');

% 修改第33-75行数据，可以获得各年度的全局和莫兰指数；

% 确定T=2年的全局莫兰指数；Y2 = Y(2,:);I20 = sum(sum(W));I21 = mean(Y2);I22 = Y2-I21;I23 = sum(I22.^2);I24 = repmat(I22,N,1);I25 = I24';I26 = sum(sum(W.*I24.*I25));Global_Moran_I2 = N*I26/(I20*I23);
if Global_Moran_I2 == 0; ('there is no spatial correlation in the year T=2!')elseif Global_Moran_I2=Z11)('spatial agglomeration exsits aound region A2 in the year T=2')

% 局部莫兰指数作图；subplot(2,2,3),plot(Z2,'.');axis([1,31,-1,1]);xlabel('T=2');ylabel('Local Moran I');

% 确定T=3年的全局莫兰指数；Y3 = Y(3,:);I30 = sum(sum(W));I31 = mean(Y3);I32 = Y3-I31;I33 = sum(I32.^2);I34 = repmat(I32,N,1);I35 = I34';I36 = sum(sum(W.*I34.*I35));Global_Moran_I3 = N*I36/(I30*I33);
if Global_Moran_I3 == 0; ('there is no spatial correlation in the year T=3 !')elseif Global_Moran_I3=Z11)
('spatial agglomeration exsits aound region A3 in the year T=3')

% 局部莫兰指数作图；subplot(2,2,4),plot(Z3,'.');axis([1,31,-1,1]);xlabel('T=3');ylabel('Local Moran I');

% 全局莫兰指数作图；Global_Moran_I = [Global_Moran_I1,Global_Moran_I2,Global_Moran_I3];subplot(2,2,1),plot(Global_Moran_I,'-*');xlabel('T');ylabel('Global Moran I');

资源仅供学术交流使用，不得商用！

如有侵权，请联系小编ermao__1128
推文期数：2019121
责任编辑：陈惟刚李崇巍
推文审核：骆丹云王善辉
转载来源：小范空间计量工作室（2018-07-20）
总审核：学术无界顾问团队
下期预告：进展：城镇化与生态环境耦合器（UEC）
往期回顾：MATLAB 实用技术分享（三）
灰色预测GM(1,1)的MATLAB模拟
读研三年，这些时间轴很重要
就业招聘 | 硕博招聘信息三十七

Lucky3141 · 2022-5-6 15:19:01

你好，我觉得代码可能有点错误，可以探讨一下吗？

面板数据（N=31,T=3）的莫兰指数计算MATLAB模拟

1 个回复