¿Cómo leer los resultados de la prueba de Dunn ' s?
On febrero 1, 2021 by admin¿Cómo leo los resultados de Dunn «s prueba ? Específicamente, ¿qué me dicen los valores de la siguiente tabla?
Tengo datos no paramétricos en 4 grupos, y primero hice un Kruskal-Wallis prueba para confirmar que las distribuciones de los grupos eran diferentes entre sí y del conjunto de datos agregado. Luego utilicé la prueba de Dunn para ver qué grupos eran diferentes entre sí y cuáles no.
library(dunn.test) dunn.test(data, g=area, kw=TRUE) Kruskal-Wallis rank sum test data: x and area Kruskal-Wallis chi-squared = 1730.4401, df = 3, p-value = 0 Comparison of x by area (No adjustment) Row Mean-| Col Mean | A B C ---------+--------------------------------- B | 20.62310 | 0.0000 | C | 26.66519 -0.087499 | 0.0000 0.4651 | D | 39.09084 7.401256 9.469204 | 0.0000 0.0000 0.0000
Respuesta
El resultado de la prueba de Kruskal-Wallis proporciona todas las posibles comparaciones por pares (seis en el caso de cuatro grupos). Por lo tanto, el de la primera fila compara el grupo B con el grupo A, el primero en la segunda fila compara el grupo C con el grupo A, etc.).
El número superior para cada comparación es la estadística de prueba z por pares de Dunn. El número más bajo es en este ejemplo el valor p sin procesar asociado con la prueba (es decir, se compararía con $ \ alpha / 2 $, aunque este valor p cambia según en la opción de ajuste de comparaciones múltiples de tasa de error familiar o tasa de descubrimiento falso. Para ajustes de comparación múltiple escalonada (por ejemplo, Holm, Benjamini-Hochberg, etc.), los valores ajustados de p tendrán un asterisco a continuación para ellos si rechaza las hipótesis nulas en el nivel de significancia especificado (que no está necesariamente indicado directamente por los valores ajustados de p , ya que el rechazo depende del orden … consulte la documentación y las citas que contiene para obtener más información detalles.).
Soy el autor de este paquete (enviarme un correo electrónico, como se indica explícitamente en la documentación, probablemente sería la mejor manera de ponerse en contacto conmigo directamente).
kruskal.test
difieren ligeramente de los los obtenidos con la opciónkw=TRUE
en el paquete anterior. ¿Esto se debe a algunas aproximaciones o son dos cosas diferentes?kruskal.test
no ajusta los empates, mientras quedunn.test
ajusta los empates. Como dice la documentación … justo encima de donde da el correo electrónico para contactar si tiene preguntas de soporte.