Cum să citiți rezultatele testului Dunn '?
On februarie 1, 2021 by adminCum citesc rezultatele din Dunn „s test ? Mai exact, ce îmi spun valorile din tabelul de mai jos?
Am date non-parametrice în 4 grupuri și am făcut mai întâi un Kruskal-Wallis test pentru a confirma că distribuțiile grupurilor au fost diferite între ele și setul de date agregat. Apoi am folosit testul lui Dunn pentru a vedea ce grupuri erau diferite între ele și care nu.
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
Răspuns
Rezultatul care urmează testului Kruskal-Wallis oferă toate comparațiile perechi posibile (șase în cazul a patru grupuri). Deci, cel din primul rând compară grupul B cu grupul A, primul pe al doilea rând compară grupul C cu grupul A etc.).
Numărul superior pentru fiecare comparație este statistica testului Dunn z în perechi. Numărul inferior este în acest exemplu valoarea p brută asociată testului (adică ați compara cu $ \ alpha / 2 $, deși această valoare p se modifică în funcție de pe opțiunea de ajustare a comparației multiple a ratei de eroare în funcție de familie sau a ratei de descoperire falsă. Pentru ajustări multiple de comparație în trepte (de exemplu, Holm, Benjamini-Hochberg etc.), valorile p ajustate vor avea un asterisc pentru ei dacă respingeți ipotezele nule la nivelul de semnificație specificat (ceea ce nu este neapărat indicat în mod direct de valorile p ajustate, deoarece respingerea depinde de ordonare … consultați documentația și citările din aceasta pentru mai multe detalii.).
Sunt autorul acestui pachet (prin e-mail, așa cum este indicat explicit în documentație, ar fi probabil cel mai bun mod de a lua legătura direct cu mine).
kruskal.test
diferă ușor de cele obținute cu opțiuneakw=TRUE
în pachetul de mai sus. Acest lucru se datorează unor aproximări sau sunt două lucruri diferite?kruskal.test
nu se ajustează pentru legături, în timp cedunn.test
se ajustează pentru legături. Așa cum se spune în documentație … chiar mai sus, unde trimite e-mailul de contact pentru întrebări de asistență.