matlab z统计量,请问如何matlab计算空间全局Moran's I 及 Z P 值

以31个省域的R&D费用内部支出为例,计算其Moran和Z,P值可以通过函数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)

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值