人大经济论坛 › 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › MATLAB等数学软件专版 › 请问如何matlab计算空间全局Moran's I 及 Z P 值

CDA数据分析研究院

商业数据分析与大数据领航教育品牌



经管云课堂

经管/金融/财会/社科/名师公开课



学术培训

Stata 空间计量 SSCI Python

贵宾：通行论坛特权+数据库权限
+案例库+下载特权 VIP：论坛特权+更多下载次数
+ccerdata数据库+更高阅读权限+……

12 下一页

发帖

楼主: wangx3

10238 13

[实际应用] 请问如何matlab计算空间全局Moran's I 及 Z P 值 [推广有奖]

0关注
0粉丝

初中生

66%

还不是VIP/贵宾

威望: 0 级
论坛币: 9 个
通用积分: 0
学术水平: 0 点
热心指数: 0 点
信用等级: 0 点
经验: 229 点
帖子: 8
精华: 0
在线时间: 15 小时
注册时间: 2015-11-3
最后登录: 2022-5-2

楼主

wangx3 发表于 2016-1-19 14:05:36 |只看作者 |坛友微信交流群|倒序 |AI写论文

是否 +2 论坛币

k人参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

请问谁有计算Moran's I 及 Z、 P 的公式或者程序呢
我除了MATLAB其它都没用过。STATA中有现成的程序可以计算么？还是要自己编程呢。
需要用Geoda么，我已经有空间权重矩阵了，就想校验Moran's I 及 Z、 P 值。
请多指教！

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏9 回帖

关键词：MATLAB moran matla atlab Mat matlab 空间如何

相关帖子

使用道具举报

沙发

zhaozhao16 发表于 2016-1-19 16:33:13 |只看作者 |坛友微信交流群

请问你的空间权重矩阵是如何生成的啊，最近初学MATLAB，还望指教

使用道具举报

藤椅

weifenfen4625 发表于 2016-1-24 19:44:15 |只看作者 |坛友微信交流群

以31个省域的R&D费用内部支出为例，计算其Moran和Z，Ｐ值可以通过函数1--NORMSDIST（Z）得到。下面是计算Moran和Z的代码。
x=[11850469 4280921 2818551 1549799 1171877 4459322 1196882 1647838 7767847 14874466 8172675 3520833 3140589 1354972 11758027 3553246 4462043 3270253 14434527 1076790 148357 1764911 3999702 471850 798371 23033 3427454 669194 137541 209042 454598];
n=31;
s=var(x,1);
m=mean(x);
w=[0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
1 1 0 0 1 1 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0
0 0 1 1 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 1 0
0 0 1 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 1 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 0 1 1 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 1 1 0 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 1 1 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 1 1 1 0 0 0 1 1 1 1 0 0 0 0 0 0 0 0 0 0 0
0 0 1 0 0 0 0 0 0 1 0 1 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
0 0 1 1 0 0 0 0 0 0 0 1 0 0 1 0 1 0 0 0 0 0 0 0 0 0 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 1 0 1 0 0 0 1 0 0 0 0 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 1 0 1 0 1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 1 0 1 1 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 1 1 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 1 1 0 0 1 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 1 0 1 1 1 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 1 1 0 1 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 1 0 1
0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0 0 1 1 0 0 0 0 1 0 1 0
0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 1 0 1 1 1
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 0 0 1
0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 1 0 0
];
y=0;
a=0;
for i=1:1:n
for j=1:1:n
      if (i==j)
         continue
      end
      y=y+w(i,j)*(x(i)-m)*(x(j)-m);
      a=a+w(i,j);
end
end
moran=y/(s*a)
k21=0;
k22=0;
w0=0;
w1=0;
w2=0;
for i=1:1:n
k21=k21+(x(i)-m)^4;
k22=k22+(x(i)-m)^2;
   for j=1:1:n
      w0=w0+w(i,j);
      w1=w1+(w(i,j)+w(j,i))^2;
   end
end
for i=1:1:n
w12(i)=0;
w21(i)=0;
  for j=1:1:n
   w12(i)=w12(i)+w(i,j);
   w21(i)=w21(i)+w(j,i);
  end
w2=w2+(w12(i)+w21(i))^2;
end
k2=(n*k21)/(k22^2);
w1=w1/2;
ei=-1/(n-1);
vari=(n*((n^2+3*n+3)*w1-n*w2+3*w0^2)-k2*((n^2-n)*w1-2*n*w2+6*w0^2))/(w0^2*(n-1)*(n-2)*(n-3))-ei^2;
z=(moran-ei)/vari^(1/2)

使用道具举报