índices de Morán

Nicolle Salamanca

Ejercicio conductividad eléctrica

library(ape)
library(clhs)
library(dplyr)

library(readxl)
datosMoran <- BD_MORAN_1_
datosfin<-sample_n(datosMoran,size = 2800)
as.data.frame(datosfin)

plot(datosfin$X_WGS84,datosfin$Y_WGS84,col=datosfin$CEa_075,xlab = 'Longitud',ylab = 'Latitud',main = 'Conductividad eléctrica aparente a 75cm de profundidad',cex=0.5,pch=19)

ce.matrix<-as.matrix(dist(cbind(datosfin$X_WGS84,datosfin$Y_WGS84)))
dim(ce.matrix)

[1] 2800 2800

ce.matrix.inv<- 1/ce.matrix
diag(ce.matrix.inv)<-0
ce.matrix.inv[1:10, 1:10]

           1         2         3        4           5           6        7         8
1     0.0000  217.6834 1887.8043 541.7625    173.1476    172.9213 443.7094  285.8398
2   217.6834    0.0000  196.6004 291.9138    473.4547    472.6594 408.2945  906.3043
3  1887.8043  196.6004    0.0000 468.8671    161.8803    161.6795 360.0861  250.1293
4   541.7625  291.9138  468.8671   0.0000    243.0179    242.5567 527.1478  401.4828
5   173.1476  473.4547  161.8803 243.0179      0.0000 126840.4495 246.3956  353.8248
6   172.9213  472.6594  161.6795 242.5567 126840.4495      0.0000 245.9987  353.1094
7   443.7094  408.2945  360.0861 527.1478    246.3956    245.9987   0.0000  737.3714
8   285.8398  906.3043  250.1293 401.4828    353.8248    353.1094 737.3714    0.0000
9   271.5035 1022.5411  238.6220 364.1010    348.8009    348.1693 673.8657 3749.5618
10 2815.5312  205.6280 4396.1601 520.7393    168.0069    167.7908 389.1677  264.7842
           9        10
1   271.5035 2815.5312
2  1022.5411  205.6280
3   238.6220 4396.1601
4   364.1010  520.7393
5   348.8009  168.0069
6   348.1693  167.7908
7   673.8657  389.1677
8  3749.5618  264.7842
9     0.0000  251.7553
10  251.7553    0.0000

Moran.I(datosfin$CEa_075,ce.matrix.inv)

$observed
[1] 0.3855839

$expected
[1] -0.0003572705

$sd
[1] 0.002985373

$p.value
[1] 0

Debido a que el índice de Moran fue de 0,38, a 75 cm de profundidad se puede observar que la conductividad eléctrica aparente presenta una dependencia espacial positiva.

plot(datosfin$X_WGS84,datosfin$Y_WGS84,col=datosfin$CEa_150,xlab = 'Longitud',ylab = 'Latitud',main = 'Conductividad eléctrica aparente a 150cm de profundidad',cex=0.5,pch=19)

Moran.I(datosfin$CEa_150,ce.matrix.inv)

$observed
[1] 0.2619223

$expected
[1] -0.0003572705

$sd
[1] 0.00298496

$p.value
[1] 0

Debido a que el índice de Moran fue de 0,26 a 150 cm de profundidad se puede observar que la conductividad eléctrica aparente presenta una dependencia espacial positiva.