Wie lese ich die Ergebnisse des Dunn ' -Tests?
On Februar 1, 2021 by adminWie lese ich die Ergebnisse von Dunns test ? Was sagen mir die Werte in der folgenden Tabelle?
Ich habe nicht parametrische Daten in 4 Gruppen und habe zuerst einen Kruskal-Wallis durchgeführt Test, um zu bestätigen, dass die Gruppenverteilungen voneinander und vom aggregierten Datensatz verschieden waren. Ich habe dann Dunns Test verwendet, um festzustellen, welche Gruppen sich voneinander unterschieden und welche nicht.
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
Antwort
Die Ausgabe nach dem Kruskal-Wallis-Test bietet alle möglichen paarweisen Vergleiche (sechs bei vier Gruppen). Die in der ersten Zeile vergleicht also Gruppe B mit Gruppe A, die erste In der zweiten Zeile wird Gruppe C mit Gruppe A usw. verglichen.
Die obere Zahl für jeden Vergleich ist Dunns paarweise z -Teststatistik. Die niedrigere Zahl ist in diesem Beispiel der rohe p -Wert, der dem Test zugeordnet ist (dh Sie würden ihn mit $ \ alpha / 2 $ vergleichen, obwohl sich dieser p -Wert abhängig ändert Bei der Anpassungsoption für die familienbezogene Fehlerrate oder die Falschentdeckungsrate für Mehrfachvergleiche. Für schrittweise Mehrfachvergleichsanpassungen (z. B. Holm, Benjamini-Hochberg usw.) werden die angepassten p -Werte als nächstes mit einem Sternchen versehen an sie, wenn Sie die Nullhypothesen mit dem angegebenen Signifikanzniveau ablehnen würden (was nicht unbedingt direkt durch die angepassten p -Werte angezeigt wird, da die Ablehnung von der Bestellung abhängt … weitere Informationen finden Sie in der Dokumentation und den Zitaten darin Details.).
Ich bin der Autor dieses Pakets (eine E-Mail an mich, wie in der Dokumentation ausdrücklich angegeben, ist wahrscheinlich der beste Weg, um direkt mit mir in Kontakt zu treten).
Kommentare
- Entschuldigung, wenn dies die Frage in eine andere Richtung lenkt, aber wann würden Sie den Test anpassen? Welche Kriterien Sind Ihre Daten erforderlich, damit Sie eine Anpassung auswählen können?
- @pocketlizard Siehe familienbezogene Fehlerrate und Falschentdeckungsrate .
- Ich habe gesehen, dass die mit
kruskal.test
erhaltenen Kruskal-Wallis-Werte geringfügig von der abweichen diejenigen, die mit der Optionkw=TRUE
im obigen Paket erhalten wurden. Liegt dies an einigen Annäherungen oder handelt es sich um zwei verschiedene Dinge? - @GennaroTedesco
kruskal.test
passt sich nicht an Bindungen an, währenddunn.test
passt sich den Bindungen an. Wie in der Dokumentation angegeben … direkt darüber, wo die E-Mail-Adresse für Support-Fragen angegeben wird.
Schreibe einen Kommentar