¿Cómo leer los resultados de la prueba de Dunn?

12

¿Cómo leo los resultados de la prueba de Dunn ? Específicamente, ¿qué me dicen los valores en la tabla a continuación?

Tengo datos no paramétricos en 4 grupos, y primero hice una prueba de Kruskal-Wallis para confirmar que las distribuciones de los grupos eran diferentes entre sí y con el conjunto de datos agregado. Luego usé 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
bubbalouie
fuente

Respuestas:

17

La salida que sigue a la prueba de Kruskal-Wallis proporciona todas las comparaciones por pares posibles (seis en el caso de cuatro grupos). Entonces, 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 bruto asociado con la prueba (es decir, se compararía con , aunque este valor p cambia según 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 paso a paso (por ejemplo, Holm, Benjamini-Hochberg, etc.), los valores p ajustados tendrán un asterisco junto a ellos si rechaza las hipótesis nulas en el nivel de significancia especificado (que no necesariamente se indica directamente por el p ajustadoα/2-valores ya que el rechazo depende del pedido ... consulte la documentación y las citas para obtener más 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 contactarme directamente).

Alexis
fuente
Lo siento si esto lleva la pregunta en una dirección diferente, pero ¿cuándo ajustarías la prueba? ¿Qué criterios son necesarios en sus datos para que desee elegir un ajuste?
pocketlizard
He visto que los valores de Kruskal-Wallis obtenidos con kruskal.testligeramente difieren de los obtenidos con la opción kw=TRUEen el paquete anterior. ¿Esto se debe a algunas aproximaciones o son dos cosas diferentes?
Gented
@GennaroTedesco kruskal.testno se ajusta a los lazos, mientras que dunn.testsí se ajusta a los lazos. Como dice en la documentación ... justo arriba de donde se le da el correo electrónico para contactar con preguntas de soporte.
Alexis