Tarea Calculo Indice de Moran

DATOS

library(clhs)
set.seed(12345)
MO=rnorm(n = 150,mean = 3,sd = 0.5)
xy=expand.grid(x=seq(1,10),y=seq(1,15))

Se crea un plot con los datos de distribucion espacial de Materia organica un data frame y un ggplot para ver la distribucion

plot(xy,col=MO,pch=19, main = "Distribucion espacial de Materia Organica")

dim(xy)
## [1] 150   2
dfmo<-data.frame(MO, xy)
head(dfmo)
##         MO x y
## 1 3.292764 1 1
## 2 3.354733 2 1
## 3 2.945348 3 1
## 4 2.773251 4 1
## 5 3.302944 5 1
## 6 2.091022 6 1
library(ggplot2)
ggplot(dfmo, aes(x = x, y=y, fill = MO))+
  geom_tile( size = 3)

Se genera la matriz con luego el inverso de la misma y como se genera un infinito entonces lo igualamos a cero para posteriormente calcular el Indice de Moran

library(ape)
mdistancias <- as.matrix(dist(cbind(xy$x, xy$y)))#Matriz de distancias
mdistanciasinv <- 1/mdistancias #inverso de la matriz de las distancias
diag(mdistanciasinv) <- 0 #pasar la diagona de infinito a cero
 
mdistanciasinv[1:5, 1:5]
##           1         2   3         4         5
## 1 0.0000000 1.0000000 0.5 0.3333333 0.2500000
## 2 1.0000000 0.0000000 1.0 0.5000000 0.3333333
## 3 0.5000000 1.0000000 0.0 1.0000000 0.5000000
## 4 0.3333333 0.5000000 1.0 0.0000000 1.0000000
## 5 0.2500000 0.3333333 0.5 1.0000000 0.0000000
Moran.I(MO, mdistanciasinv) #indice de moran
## $observed
## [1] -0.009650003
## 
## $expected
## [1] -0.006711409
## 
## $sd
## [1] 0.007694112
## 
## $p.value
## [1] 0.7025151

De acuerdo al pvalor se acepta la hipotesis nula y se dice que no hay dependencia espacial