Estoy trabajando con un análisis espacial exploratorio en R usando el paquete spdep.
Encontré una opción para ajustar los valores p de los indicadores locales de asociación espacial (LISA) calculados usando la localmoran
función. Según los documentos, está dirigido a:
... ajuste del valor de probabilidad para múltiples pruebas.
Además en los documentos de p.adjustSP
He leído que las opciones disponibles son:
Los métodos de ajuste incluyen la corrección de Bonferroni ('"bonferroni"') en la que los valores p se multiplican por el número de comparaciones. Holm (1979) ('"holm"'), Hochberg (1988) ('"hochberg"'), Hommel (1988) ('"hommel"') y Benjamini & Hochberg (1995) también incluyen cuatro correcciones menos conservadoras. ('"fdr"'), respectivamente. También se incluye una opción de paso ('"none"').
Los primeros cuatro métodos están diseñados para dar un fuerte control de la tasa de error familiar. Parece que no hay razón para usar la corrección de Bonferroni no modificada porque está dominada por el método de Holm, que también es válido bajo supuestos arbitrarios.
Los métodos de Hochberg y Hommel son válidos cuando las pruebas de hipótesis son independientes o cuando están asociadas de manera no negativa (Sarkar, 1998; Sarkar y Chang, 1997). El método de Hommel es más poderoso que el de Hochberg, pero la diferencia suele ser pequeña y los valores p de Hochberg son más rápidos de calcular.
El método "BH" (también conocido como "fdr") y "BY" de Benjamini, Hochberg y Yekutieli controlan la tasa de descubrimiento falso, la proporción esperada de descubrimientos falsos entre las hipótesis rechazadas. La tasa de descubrimiento falso es una condición menos estricta que la tasa de error familiar, por lo que estos métodos son más potentes que los demás.
Un par de preguntas que aparecieron:
- En palabras simples, ¿cuál es el propósito de este ajuste?
- ¿Es necesario usar tales correcciones?
- En caso afirmativo, ¿cómo elegir entre las opciones disponibles?