葫芦娃大王 发表于 2014-7-5 21:49
非常好,感谢楼主!
从编程的角度看你确定是正确的吗?
我看了相关的mran值期望方差公式,感觉,moran值应该没问题,但是怎么感觉方差有点问题呢?
楼主代码里面的方差计算公式
第一, w1=w1+(w(i,j)+w(j,i))^2;好像应该还要除个2;
第二,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;
好像应该是
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+1)*(n-3))-ei^2;
你觉得呢?欢迎大家讨论啊,严谨对待!