Hvordan leser du resultatene av Dunn ' s test?
On februar 1, 2021 by adminHvordan leser jeg resultatene fra Dunn «s test ? Hva forteller verdiene i tabellen nedenfor meg spesifikt?
Jeg har ikke-parametriske data i 4 grupper, og jeg gjorde først en Kruskal-Wallis test for å bekrefte at gruppene «distribusjoner var forskjellige fra hverandre og det samlede datasettet. Jeg brukte deretter Dunns test for å se hvilke grupper som var forskjellige fra hverandre, og hvilke som ikke var det.
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
Svar
Resultatet etter Kruskal-Wallis-testen gir alle mulige parvise sammenligninger (seks når det gjelder fire grupper). Så den på første rad sammenligner gruppe B med gruppe A, den første på den andre raden sammenligner gruppe C med gruppe A osv.).
Det øvre tallet for hver sammenligning er Dunns parvise z teststatistikk. Det lavere tallet er i dette eksemplet den rå p -verdien som er knyttet til testen (dvs. du vil sammenligne med $ \ alpha / 2 $, selv om denne p -verdien endres avhengig av på den familiemessige feilraten eller falsk oppdagelsesfrekvens flere sammenligningsjusteringsalternativer. For trinnvise flere sammenligningsjusteringer (f.eks. Holm, Benjamini-Hochberg, etc.), vil de justerte p -verdiene få en stjerne neste til dem hvis du vil avvise nullhypotesene på det angitte signifikansnivået (som ikke nødvendigvis direkte er angitt av de justerte p -verdiene, siden avvisning avhenger av bestilling … se dokumentasjonen og sitatene der for mer detaljer.).
Jeg er forfatter av denne pakken (e-post til meg, som eksplisitt angitt i dokumentasjonen, vil sannsynligvis være den beste måten å komme i kontakt med meg direkte).
kruskal.test
er litt forskjellige fra de som er oppnådd med alternativetkw=TRUE
i pakken ovenfor. Skyldes dette noen tilnærminger, eller er det to forskjellige ting?kruskal.test
justerer ikke for bånd, mensdunn.test
justerer for bånd. Som det står i dokumentasjonen … rett over der den gir e-post til kontakt for spørsmål om støtte.