Hvordan læses resultaterne af Dunn ' s test?
On februar 1, 2021 by adminHvordan læser jeg resultaterne fra Dunn “s test ? Specifikt, hvad fortæller værdierne i nedenstående tabel mig?
Jeg har ikke-parametriske data i 4 grupper, og jeg lavede først en Kruskal-Wallis test for at bekræfte, at gruppernes “distributioner var forskellige fra hinanden og det samlede datasæt. Derefter brugte jeg Dunns test for at se, hvilke grupper der var forskellige fra hinanden, og hvilke der ikke var.
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 efter Kruskal-Wallis-testen giver alle mulige parvise sammenligninger (seks i tilfælde af fire grupper). Så den på første række sammenligner gruppe B med gruppe A, den første på anden række sammenligner gruppe C med gruppe A osv.).
Det øverste tal for hver sammenligning er Dunns parvise z teststatistik. Det lavere tal er i dette eksempel den rå p -værdi, der er knyttet til testen (dvs. du vil sammenligne med $ \ alpha / 2 $, selvom denne p -værdien ændrer sig afhængigt af på den familiemæssige fejlfrekvens eller falsk opdagelsesrate flere sammenligningsjusteringsmuligheder. For trinvise flere sammenligningsjusteringer (f.eks. Holm, Benjamini-Hochberg osv.) får de justerede p -værdier en stjerne næste til dem, hvis du ville afvise nulhypoteserne på det angivne signifikansniveau (hvilket ikke nødvendigvis direkte er angivet af de justerede p -værdier, da afvisning afhænger af bestilling … se dokumentationen og citater deri for mere detaljer.).
Jeg er forfatter til denne pakke (e-mail til mig, som eksplicit angivet i dokumentationen, ville sandsynligvis være den bedste måde at komme i kontakt med mig direkte).
kruskal.test
adskiller sig lidt fra dem, der fås med indstillingenkw=TRUE
i ovenstående pakke. Er dette på grund af nogle tilnærmelser, eller er de to forskellige ting?kruskal.test
justerer ikke for bånd, mensdunn.test
justerer for slips. Som det står i dokumentationen … lige over hvor den giver e-mailen til kontakt for supportspørgsmål.