En anteriores tutoriales hemos aprendido a realizar el Análisis de tendencias espaciales con ArcGIS y a utilizar la herramienta de Estadísticos Zonales con ArcGIS. En esta ocasión veremos cómo realizar un análisis de patrones espaciales con ArcGIS para conocer el tipo de distribución que presentan los datos que queremos analizar, es decir, si se reparten de manera dispersa o si por el contrario hay valores que se concentran en determinadas zonas.
Independientemente del tipo de análisis de patrones espaciales con ArcGIS que vayamos a realizar, en todos ellos obtendremos una serie de parámetros que serán los mismos. Siempre partiremos de una hipótesis nula que establece que las entidades o los valores de las entidades se distribuyen de manera aleatoria.
Al realizar el análisis, obtendremos unos parámetros «z» y «p» que nos indicarán si esta hipótesis se cumple o no, es decir, si debemos rechazar la hipótesis nula.
- Valor p: es una probabilidad; si este valor es muy pequeño, significa que la probabilidad de que se cumpla la hipótesis nula es muy pequeña por lo que podremos rechazarla
- Valor z: representa las desviaciones estándar; en general, si éste valor es muy pequeño (próximo a 0) no tendremos suficiente evidencia estadística como para rechazar la hipótesis nula. Por este motivo, z toma valores muy altos (positivos o negativos) cuando el valor de p es pequeño.
Análisis de patrones espaciales con ArcGIS: Relación espacial (Índice de Morans I)
Este tipo de análisis de patrones espaciales con ArcGIS analiza el patrón de distribución de los valores de la entidad . Nosotros tenemos una serie de estaciones pluviométricas con valores de precipitación media anual por lo que este análisis nos va a permitir conocer cómo se distribuyen los valores de precipitación en nuestra región de estudio.
Esta herramienta calcula el valor medio y la varianza del atributo que estamos evaluando (precipitación) y deduce si el patrón de distribución que presentan los datos es agrupado, disperso o aleatorio midiendo la autocorrelación espacial basada en las ubicaciones y los valores de las entidades simultáneamente.
Si los valores tienden a agruparse espacialmente (es decir, los valores altos se agrupan cerca de otros valores altos o los valores bajos se agrupan cerca de otros valores bajos), el Índice de Morans será positivo; por el contrario, si los valores se dispersan (valores altos tienden a estar cerca de valores bajos), el índice será negativo.
Esto implica que, siempre que los valores z y p obtenidos revelen que podemos rechazar la hipótesis nula, un índice de Morans mayor que 0 indicará una tendencia de los valores hacia la agregación (clustering); por el contrario, un índice menor que 0 indicará una tendencia de los valores hacia la dispersión.
Vamos a realizar nuestro análisis ArcToolbox < Spatial Statistics Tools < Analyzing Patterns < Spatial autocorrelation (Morans I). Seleccionamos nuestra capa de estaciones pluviométricas (‘Estaciones’) y el dato que queremos analizar, es decir, la precipitación (Input Field: «P_mm»). Activamos la casilla para que el programa genere el informe de resultados correspondiente (‘Generate report’) y dejamos el resto de parámetros que aparecen por defecto.
Los análisis de patrones no generan capas si no informes que se pueden visualizar mediante un navegador de internet, seleccionando el archivo HTML de la ventana «Resultados».
En esta ventana se muestran los valores obtenidos; sin embargo, la mejor manera de visualizar los resultados lo conseguimos abriendo el informe completo. Para ello, seleccionamos el archivo (se habrá guardado en la carpeta de destino que hayamos indicado), y lo abrimos con el navegador que utilizamos habitualmente.
Análisis de patrones espaciales con ArcGIS: Clustering Alto/Bajo (G general de Getis Ord)
Según hemos visto, mediante el Índice de Morans I podemos averiguar si los valores analizados presentan un patrón de distribución disperso, aleatorio o agregado. Con el análisis de patrones espaciales con ArcGIS que vamos a realizar ahora mediremos la concentración de valores altos o bajos para un mismo área de estudio. Esta herramienta calcula un índice denominado «G general» observado y lo comprara con el esperado.
En este caso, un valor de z mayor que cero indicará que existe una concentración de valores altos; si por el contrario el valor z es negativo, la concentración que se produce es de valores bajos. Lógicamente, si el valor de z se encuentra muy próximo a cero, no existe agrupación evidente de valores de la entidad.
Lanzamos nuestro análisis Clustering ArcToolbox < Spatial Statistics Tools < Analyzing Patterns < High-Low Clustering (Getis-Ord General G). De nuevo seleccionamos nuestra capa de estaciones pluviométricas (‘Estaciones’) y el dato que queremos analizar (Input Field: «P_mm»). Activamos la casilla para que el programa genere el informe de resultados correspondiente (‘Generate report’) y dejamos el resto de parámetros que aparecen por defecto.
Abrimos nuestro informe de la misma forma que en el caso anterior y lo analizamos:
El análisis de Getis-Ord muestra un valor de p relativamente pequeño (p = 0.012) y un valor de z negativo (z = -2.49); esto implica que existe una evidencia estadística suficientemente significativa como para afirmar que se produce una agrupación de valores bajos de precipitación (Low-Clusters).
Según el mensaje inferior, dado el valor de z obtenido, hay menos de un 5% de probabilidad de que el patrón de agregación de valores bajos de precipitación sea resultado de la casualidad.
