R语言ape包计算莫兰指数和自己根据公式写的莫兰指数计算生成了不同的结果

有大佬可以帮我检查一下是为什么吗

自己根据公式写的代码计算莫兰指数

1. >w <- matrix(c(0,1,1,1,1,0,1,0,1,1,0,1,1,0,1,0),nrow=4,ncol=4,dimnames=list(c("千","昭","天","名"),c("千","昭","天","名")))
   
2. >a <- c(80,95,67,78)
   
3. >mean(a)
   
4. >b <- a-mean(a)
   
5. >i <- (length(a)/sum(w))*(t(b)%*%w%*%b)/(t(b)%*%b)
   
6. >i

复制代码

           [,1]
[1,] -0.3396985
用R语言ape包 计算莫兰指数

1. >install.packages("ape")
   
2. >library(ape)
   
3. >Moran.I(a,w)

复制代码

$observed
[1] -0.3538526


$expected
[1] -0.3333333


$sd
[1] 0.1521452


$p.value
[1] 0.8927175

顺便问一下，有大佬自己写了moran i 的计算和检定公式吗？谢谢！

请问你莫兰指数不显著最后是怎么解决的啊?谢谢啦，我的qq1421760831@qq.com，可以加个qq问你一下吗？我看你这个p值0.89好大啊，这样是不是不显著啊，谢谢啦

我也遇到这个问题了，请问楼主解决了吗？？

是因为空间权重矩阵的缘故。ape包中在计算莫兰指数时，选取的空间权重矩阵是经过行标准化处理的；而作者您所计算的空间权重矩阵并没有做相关的处理。当您把空间权重矩阵做完行标准化处理之后，莫兰指数的计算结果就与ape包中一致了。

涵紫 发表于 2022-1-27 11:31
是因为空间权重矩阵的缘故。ape包中在计算莫兰指数时，选取的空间权重矩阵是经过行标准化处理的；而作者您所 ...

能看看怎么标准化吗，试了好几个标准化还是不对