楼主您好，你文中的局部空间自相关有局部莫兰I的计算程序，想请教下，局部莫兰I的Z值检验程序可以给出下吗？我自己尝试了好久，算出的结果都不显著，非常小，觉得应该是我的程序出问题了，但是不知道怎么调，就是下面这样，菜鸟新手一枚，还望楼主指教，跪求答案：
for k=1:1:n
w0=0;
w1=0;
w12=0;
w21=0;
k21=(X(k)-m)^4;
k22=(X(k)-m)^2;
fori=1:1:n
for j=1:1:n
if (i==j)
continue
end
w0=w0+M(i,j);
w1=w1+(M(i,j)+M(j,i))^2;
w12=w12+M(i,j);
w21=w21+M(j,i);
end
w2=w2+(w12+w21)^2;
end
w1=w1/2;
k2=(n*k21)/(k22^2);
e(k)=-1/(n-1);
var(k)=(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))-e(k)^2;
z(k)=(moran(k)- e(k))/var(k)^(1/2);
end
for k=1:1:n
fprintf('%2.4f\n',z(k))end
扫码加好友，拉您进群
全部版块
我的主页
不错
