Articulo de referencia

Moran I

Los cuadrados blancos y negros están perfectamente distribuidos, de modo que el índice de Moran sería -1 si se utilizara la definición de vecinos de torre . Si los cuadrados bla...

Los cuadrados blancos y negros están perfectamente distribuidos, de modo que el índice de Moran sería -1 si se utilizara la definición de vecinos de torre . Si los cuadrados blancos se colocaran en una mitad del tablero y los negros en la otra, el índice de Moran se aproximaría a +1 a medida que N aumenta. Una disposición aleatoria de los colores de los cuadrados haría que el índice de Moran se aproximara a 0. 

En estadística, el índice I de Moran es una medida de autocorrelación espacial desarrollada por Patrick Alfred Pierce Moran . [ 1 ] [ 2 ] La autocorrelación espacial se caracteriza por una correlación en una señal entre ubicaciones cercanas en el espacio. La autocorrelación espacial es más compleja que la autocorrelación unidimensional porque la correlación espacial es multidimensional (es decir, de 2 o 3 dimensiones espaciales) y multidireccional.

La Moran global I

El índice I de Moran global es una medida de la agrupación general de los datos espaciales. Se define como

I=norteS0i=1nortej=1nortewij(incógnitaiincógnita¯)(incógnitajincógnita¯)i=1norte(incógnitaiincógnita¯)2{\displaystyle I={\frac {N}{S_{0}}}{\frac {\sum _{i=1}^{N}\sum _{j=1}^{N}w_{ij}(x_{i}-{\bar {x}})(x_{j}-{\bar {x}})}{\sum _{i=1}^{N}(x_{i}-{\bar {x}})^{2}}}}

dónde

  • norte{\displaystyle N}es el número de unidades espaciales indexadas pori{\displaystyle i}yj{\displaystyle j};
  • incógnita{\displaystyle x}es la variable de interés;
  • incógnita¯{\displaystyle {\bar {x}}}es la media deincógnita{\displaystyle x};
  • wij{\displaystyle w_{ij}}son los elementos de una matriz de pesos espaciales con ceros en la diagonal (es decir,wii=0{\displaystyle w_{ii}=0});
  • yS0{\displaystyle S_{0}}es la suma de todoswij{\displaystyle w_{ij}},S0=i=1nortej=1nortewij{\displaystyle S_{0}=\sum _{i=1}^{N}\sum _{j=1}^{N}{w_{ij}}}. [ 3 ] [ 4 ]
Estadística I de Moran calculada para diferentes patrones espaciales. Usando vecinos ' torre ' para cada celda de la cuadrícula, estableciendowij=1{\displaystyle w_{ij}=1}para los vecinos j{\displaystyle j}dei{\displaystyle i}y luego normalizando por filas la matriz de pesos. Arriba a la izquierda se muestra una anticorrelación que da un I negativo . Arriba a la derecha se muestra un gradiente espacial que da un I positivo grande . Abajo a la izquierda se muestran datos aleatorios que dan un valor de I cercano a 0 (o1/(norte1)0,04{\displaystyle -1/(N-1)\simeq -0.04}). En la parte inferior derecha se muestra una "mancha de tinta" o patrón de dispersión con autocorrelación positiva.

Definición de la matriz de ponderaciones espaciales

El valor deI{\displaystyle I}puede depender bastante de las suposiciones incorporadas en la matriz de ponderaciones espacialeswij{\displaystyle w_{ij}}La matriz es necesaria porque, para abordar la autocorrelación espacial y modelar la interacción espacial, debemos imponer una estructura que restrinja el número de vecinos a considerar. Esto se relaciona con la primera ley de geografía de Tobler , que establece que todo depende de todo lo demás, pero las cosas más cercanas aún más; en otras palabras, la ley implica una función de decaimiento de la distancia espacial , de modo que, si bien todas las observaciones influyen en todas las demás, después de cierto umbral de distancia, esa influencia puede despreciarse.

La idea es construir una matriz que refleje con precisión sus suposiciones sobre el fenómeno espacial particular en cuestión. Un enfoque común es dar un peso de 1 si dos zonas son vecinas, y 0 en caso contrario, aunque la definición de "vecinos" puede variar. Otro enfoque común podría ser dar un peso de 1 ak{\displaystyle k}vecinos más cercanos, 0 en caso contrario. Una alternativa es utilizar una función de decaimiento de distancia para asignar pesos. A veces, la longitud de una arista compartida se utiliza para asignar diferentes pesos a los vecinos. La selección de la matriz de pesos espaciales debe guiarse por la teoría sobre el fenómeno en cuestión. El valor deI{\displaystyle I}Es bastante sensible a los pesos y puede influir en las conclusiones que se extraen sobre un fenómeno, especialmente al utilizar distancias.

Valor esperado

El valor esperado del índice I de Moran bajo la hipótesis nula de ausencia de autocorrelación espacial es

mi(I)=1norte1{\displaystyle E(I)={\frac {-1}{N-1}}}

La distribución nula utilizada para esta expectativa es queincógnita{\displaystyle x}La entrada se permuta mediante una permutación.π{\displaystyle \pi }elegido uniformemente al azar (y la esperanza es sobre elegir la permutación).

Con tamaños de muestra grandes (es decir, cuando N tiende a infinito), el valor esperado tiende a cero.

Su varianza es igual a

Var(I)=norteS4S3S5(norte1)(norte2)(norte3)S02(mi(I))2{\displaystyle \operatorname {Var} (I)={\frac {NS_{4}-S_{3}S_{5}}{(N-1)(N-2)(N-3)S_{0}^{2}}}-(E(I))^{2}}

dónde

S0=ijwij{\displaystyle S_{0}=\sum _{i}\sum _{j}w_{ij}}[ 3 ]
S1=12ij(wij+wji)2{\displaystyle S_{1}={\frac {1}{2}}\sum _{i}\sum _{j}(w_{ij}+w_{ji})^{2}}[ 3 ]
S2=i(jwij+jwji)2{\displaystyle S_{2}=\sum _{i}\left(\sum _{j}w_{ij}+\sum _{j}w_{ji}\right)^{2}}[ 3 ]
S3=norte1i(incógnitaiincógnita¯)4(norte1i(incógnitaiincógnita¯)2)2{\displaystyle S_{3}={\frac {N^{-1}\sum _{i}(x_{i}-{\bar {x}})^{4}}{(N^{-1}\sum _{i}(x_{i}-{\bar {x}})^{2})^{2}}}}
S4=(norte23norte+3)S1norteS2+3S02{\displaystyle S_{4}=(N^{2}-3N+3)S_{1}-NS_{2}+3S_{0}^{2}}
S5=(norte2norte)S12norteS2+6S02{\displaystyle S_{5}=(N^{2}-N)S_{1}-2NS_{2}+6S_{0}^{2}}[ 5 ]

Los valores significativamente inferiores a -1/(N-1) indican autocorrelación espacial negativa y los valores significativamente superiores a -1/(N-1) indican autocorrelación espacial positiva. Para la comprobación de hipótesis estadísticas, los valores I de Moran se pueden transformar en puntuaciones z .

Valores deI{\displaystyle I}rango entrenorteS0wmetroinorte{\displaystyle {\frac {N}{S}}_{0}w_{min}}ynorteS0wmetroaincógnita{\displaystyle {\frac {N}{S}}_{0}w_{max}}[ 6 ] dondewmetroinorte{\displaystyle w_{min}}ywmetroaincógnita{\displaystyle w_{max}}son los valores propios mínimo y máximo correspondientes de la matriz de pesos. Para una matriz normalizada por filasnorteS0=1{\displaystyle {\frac {N}{S}}_{0}=1}.

El índice I de Moran está inversamente relacionado con el índice C de Geary , pero no son idénticos. El índice I de Moran mide la autocorrelación espacial global, mientras que el índice C de Geary es más sensible a la autocorrelación espacial local.

Moran local I

Agrupaciones del porcentaje estimado de personas en situación de pobreza por condado en los Estados Unidos continentales en 2020, calculado utilizando el índice I de Moran local de Anselin.

El análisis de autocorrelación espacial global proporciona una única estadística que resume toda el área de estudio. En otras palabras, el análisis global presupone homogeneidad. Si esta suposición no se cumple, entonces tener una sola estadística carece de sentido, ya que la estadística debería variar espacialmente.

Además, incluso si no hay autocorrelación global ni agrupamiento, aún podemos encontrar agrupaciones a nivel local utilizando el análisis de autocorrelación espacial local. El hecho de que el I de Moran sea una suma de productos cruzados individuales es aprovechado por los "indicadores locales de asociación espacial" (LISA) para evaluar el agrupamiento en esas unidades individuales calculando el I de Moran local para cada unidad espacial y evaluando la significancia estadística para cada I i . A partir de la ecuación del I de Moran global , podemos obtener:

Ii=incógnitaiincógnita¯metro2j=1nortewij(incógnitajincógnita¯){\displaystyle I_{i}={\frac {x_{i}-{\bar {x}}}{m_{2}}}\sum _{j=1}^{N}w_{ij}(x_{j}-{\bar {x}})}

dónde:

metro2=i=1norte(incógnitaiincógnita¯)2norte{\displaystyle m_{2}={\frac {\sum _{i=1}^{N}(x_{i}-{\bar {x}})^{2}}{N}}}

entonces,

I=i=1norteIiS0{\displaystyle I=\sum _{i=1}^{N}{\frac {I_{i}}{S_{0}}}}

Yo soy el Moran GlobalI{\displaystyle I}midiendo la autocorrelación global, I i es local y S_0 es nuevamente la suma de los pesos.

Los LISA se pueden calcular en GeoDa y ArcGIS Pro , que utiliza el I de Moran local , [ 7 ] [ 8 ] propuesto por Luc Anselin en 1995. [ 9 ]

Usos

Informe de autocorrelación espacial generado por ArcGIS Pro para la población de los condados de EE. UU. en 2022.

El índice de Moran se utiliza ampliamente en los campos de la geografía y la ciencia de la información geográfica . Algunos ejemplos incluyen:

  • El análisis de las diferencias geográficas en las variables de salud. [ 10 ]
  • Caracterización del impacto de las concentraciones de litio en el agua potable sobre la salud mental. [ 11 ]
  • En dialectología , para medir la importancia de la variación lingüística regional. [ 12 ]
  • Definición de una función objetivo para la segmentación significativa del terreno para estudios geomorfológicos [ 13 ]

Véase también

Referencias

  1. Moran, PAP (1950). " Notas sobre fenómenos estocásticos continuos". Biometrika . 37 (1): 17– 23. doi : 10.2307/2332142 . JSTOR 2332142. PMID 15420245 .  
  2. Li, Hongfei; Calder, Catherine A. ; Cressie, Noel (2007). "Más allá del I de Moran : prueba de dependencia espacial basada en el modelo autorregresivo espacial". Análisis geográfico . 39 (4): 357– 375. doi : 10.1111/j.1538-4632.2007.00708.x .
  3. 1 2 3 4 Cliff, AD; Ord, JK "Medidas de autocorrelación en el plano". Procesos espaciales: modelos y aplicaciones . Londres : Pion Limited. 1.5.1 Procedimientos generales (pág. 19). ISBN 0-85086-081-4.
  4. "Calcular el I de Moran" . Consultado el 31 de marzo de 2026 .
  5. Cliff y Ord (1981), Procesos espaciales, Londres
  6. ^ de Jong, P.; Sprenger, C.; van Veen, F. (1984). "Sobre valores extremos de I de Moran y c de Geary". Análisis Geográfico . 16 (1): 17– 24. doi : 10.1111/j.1538-4632.1984.tb00797.x .
  7. Anselin, Luc (2005). "Explorando datos espaciales con GeoDa: un libro de trabajo" (PDF) . Laboratorio de análisis espacial. pág. 138. 
  8. "Análisis de clústeres y valores atípicos (I de Moran local de Anselin) (Estadística espacial)" . ESRI . Consultado el 28 de mayo de 2024 .
  9. Anselin, Luc (1995). "Indicadores locales de asociación espacial—LISA" . Análisis geográfico . 27 (2): 93– 115. doi : 10.1111/j.1538-4632.1995.tb00338.x .
  10. Getis, Arthur (3 de septiembre de 2010). "Análisis de la asociación espacial mediante el uso de estadísticas de distancia" . Análisis geográfico . 24 (3): 189– 206. doi : 10.1111/j.1538-4632.1992.tb00261.x .
  11. Helbich, M; Leitner, M; Kapusta, ND (2012). "Examen geoespacial del litio en el agua potable y la mortalidad por suicidio" . Int J Health Geogr . 11 (1): 19. doi : 10.1186/1476-072X-11-19 . PMC 3441892. PMID 22695110 .  
  12. Grieve, Jack (2011). "A regional analysis of contraction rate in written Standard American English". International Journal of Corpus Linguistics. 16 (4): 514–546. doi:10.1075/ijcl.16.4.04gri.
  13. Alvioli, M.; Marchesini, I.; Reichenbach, P.; Rossi, M.; Ardizzone, F.; Fiorucci, F.; Guzzetti, F. (2016). "Automatic delineation of geomorphological slope units with r.slopeunits v1.0 and their optimization for landslide susceptibility modeling". Geoscientific Model Development. 9: 3975–3991. doi:10.5194/gmd-9-3975-2016.