Я вычисляю локальное значение Морана I для SpatialPolygonsDataFrame, как показано ниже, но значения не находятся в диапазоне -1 / + 1. Я исследовал это и нашел ниже, что предполагает, что что-то пошло не так с моей стандартизацией (https://stats.stackexchange.com/questions/160459/why-is-morans-i-coming-out-greater-than-1).
Я стремлюсь стандартизировать все Моран I, чтобы они все были в пределах стандартного диапазона.
Любая помощь будет очень признательна!
neighbours.queen.2016<-poly2nb(OA.2016, queen=TRUE)
listw.queen.2016<-nb2listw(neighbours.queen.2016, style="W",zero.policy=TRUE)
LM_Total.Contribution <- localmoran(x = OA.2016@data$Total.Contribution, listw = nb2listw(neighbours.queen.2016, style = "W", zero.policy=TRUE), zero.policy=TRUE)
colnames(LM_Total.Contribution) <- c("LM_Total.Contribution","Standard Error","Variance","Z-Score","P-Value")
LM_Total.Contribution <- LM_Total.Contribution[, c(1,5)]
OA.2016@data<-cbind(OA.2016@data,LM_Total.Contribution)