El conjunto de datos adjunto muestra aproximadamente 6000 árboles jóvenes en aproximadamente 50 huecos forestales de tamaño variable. Estoy interesado en aprender cómo crecen estos retoños dentro de sus respectivos huecos (es decir, agrupados, aleatorios, dispersos). Como saben, un enfoque tradicional sería ejecutar Global Moran's I. Sin embargo, las agregaciones de árboles dentro de agregaciones de huecos parecen ser un uso inapropiado de Moran's I. Realicé algunas estadísticas de prueba con Moran's I usando una distancia de umbral de 50 metros, que produjo resultados sin sentido (es decir, valor p = 0.0000000 ...). Es probable que la interacción entre las agregaciones de brechas produzca estos resultados. He considerado crear una secuencia de comandos para recorrer las brechas de dosel individuales y determinar la agrupación dentro de cada brecha, aunque mostrar estos resultados al público sería problemático.
¿Cuál es el mejor enfoque para cuantificar la agrupación dentro de los grupos?
fuente
Respuestas:
No tiene un campo aleatorio uniforme, por lo que intentar analizar todos sus datos a la vez violará los supuestos de cualquier estadística que elija lanzar al problema. No está claro en su publicación si sus datos son un proceso de puntos marcados (es decir, diámetro o altura asociados con la ubicación de cada árbol). Si estos datos no representan un proceso de puntos marcados, no tengo idea de cómo aplicó un Moran's-I. Si los datos solo representan ubicaciones espaciales, recomendaría usar un Ripley's-K con la transformación Besag-L para estandarizar la expectativa nula en cero. Esto permite una evaluación multiescala del agrupamiento. Si sus datos tienen un valor asociado, entonces su mejor opción es un Moran's-I (LISA) local. Realmente lo miraría con ambas estadísticas. Independientemente de su elección, aún deberá recorrer cada sitio individualmente para producir resultados válidos. Aquí hay un ejemplo de código R para una simulación de Monte Carlo de Ripley's-K / Besag's-L utilizando el conjunto de datos de árbol de redwood incorporado. Debería ser bastante sencillo modificar esto para recorrer sus sitios y generar un gráfico para cada uno.
fuente
Lo que tiene es un patrón de puntos con una ventana que es un número de pequeñas regiones poligonales desconectadas.
Debería poder utilizar cualquiera de las pruebas
package:spatstat
para CSR siempre que la alimente con una ventana correcta. Esto puede ser una serie de conjuntos de pares (x, y) que definen cada compensación o una matriz binaria de valores (0,1) sobre el espacio.Primero definamos algo que se parezca un poco a sus datos:
y pretendamos que nuestros claros son celdas cuadradas que resultan ser esto:
Entonces podemos trazar la función K de esos puntos en esa ventana. Esperamos que esto no sea CSR porque los puntos parecen agrupados dentro de las celdas. Tenga en cuenta que tengo que cambiar el rango de distancias para que sea pequeño, del orden del tamaño de la celda, de lo contrario, la función K se evalúa a distancias del tamaño de todo el patrón.
Si generamos algunos puntos CSR en las mismas celdas, podemos comparar las gráficas de la función K. Este debería ser más como CSR:
Realmente no puede ver los puntos agrupados en las celdas en el primer patrón, pero si lo traza solo en una ventana gráfica, está claro. Los puntos en el segundo patrón son uniformes dentro de las celdas (y no existen en la región negra) y la función K es claramente diferente de
Kpois(r)
la función K de CSR para los datos agrupados y similar para los datos uniformes.fuente
Además de la publicación de Andy:
Lo que desea calcular es una medida de homogeneidad espacial (por ejemplo, la hipótesis: "¿Están agrupados sus puntos?") Como la función L y K de Ripley .
Esta publicación de blog explica bastante bien cómo hacerlo en R. Según el código descrito, primero etiquetaría cada grupo en su conjunto de datos y luego calcularía en un bucle para cada grupo el sobre crítico a través de K de Ripley
fuente