Hur läser jag resultatet av Dunn ' s test?
On februari 1, 2021 by adminHur läser jag resultaten från Dunn ”s test ? Specifikt, vad säger värdena i tabellen nedan?
Jag har icke-parametriska data i fyra grupper, och jag gjorde först en Kruskal-Wallis test för att bekräfta att gruppernas distributioner var olika från varandra och den sammanlagda datamängden. Jag använde sedan Dunns test för att se vilka grupper som skilde sig från varandra och vilka inte.
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
Utgången efter Kruskal-Wallis-testet ger alla möjliga parvisa jämförelser (sex i fallet med fyra grupper). Så den på första raden jämför grupp B med grupp A, den första på andra raden jämför grupp C med grupp A, etc.).
Det övre numret för varje jämförelse är Dunns parvisa z teststatistik. Det lägre numret är i det här exemplet den råa p -värdet som är associerat med testet (dvs. du skulle jämföra med $ \ alpha / 2 $, även om denna p -värde ändras beroende på om familjefel eller felaktig upptäckthastighet flera jämförelser justeringsalternativ. För stegvisa flera jämförelsejusteringar (t.ex. Holm, Benjamini-Hochberg, etc.) kommer justerade p -värden att ha en asterisk nästa till dem om du skulle avvisa nollhypoteserna på den angivna signifikansnivån (vilket inte nödvändigtvis indikeras direkt av de justerade p -värdena eftersom avslag beror på beställning … se dokumentationen och citaten däri för mer detaljer.).
Jag är författare till detta paket (mejla mig, vilket uttryckligen anges i dokumentationen, skulle sannolikt vara det bästa sättet att komma i kontakt med mig direkt).
kruskal.test
skiljer sig något från sådana som erhållits med alternativetkw=TRUE
i ovanstående paket. Beror det på några approximationer eller är det två olika saker?kruskal.test
justerar inte för band, medandunn.test
justerar inte för slipsar. Som det står i dokumentationen … precis ovanför där det ger e-postmeddelandet att kontakta för supportfrågor.