Kruskal-Wallis-test mit details über paarweise Vergleiche

Die standard-stats::kruskal.test-Modul ermöglicht eine Berechnung der kruskal-wallis-test auf einem Datensatz:

>>> data(diamonds)
>>> kruskal.test(price~carat, data=diamonds)

Kruskal-Wallis rank sum test

data:  price by carat by color 
Kruskal-Wallis chi-squared = 50570.15, df = 272, p-value < 2.2e-16

Dies ist richtig, es ist mir eine Wahrscheinlichkeit, dass alle Gruppen in den Daten haben den gleichen Mittelwert.

Allerdings würde ich gerne die details für jedes paar ein Vergleich, wie wenn Diamanten der Farben D und E haben den gleichen Mittelwert Preis, wie einige andere softwares tun (SPSS), wenn Sie Fragen, für einen test Kruskal.

Habe ich gefunden kruskalmc aus dem Paket pgirmess die mir erlaubt, zu tun, was ich tun will:

> kruskalmc(diamonds$price, diamonds$color)
Multiple comparison test after Kruskal-Wallis 
p.value: 0.05 
Comparisons
      obs.dif critical.dif difference
D-E  571.7459     747.4962      FALSE
D-F 2237.4309     751.5684       TRUE
D-G 2643.1778     726.9854       TRUE
D-H 4539.4392     774.4809       TRUE
D-I 6002.6286     862.0150       TRUE
D-J 8077.2871    1061.7451       TRUE
E-F 2809.1767     680.4144       TRUE
E-G 3214.9237     653.1587       TRUE
E-H 5111.1851     705.6410       TRUE
E-I 6574.3744     800.7362       TRUE
E-J 8649.0330    1012.6260       TRUE
F-G  405.7470     657.8152      FALSE
F-H 2302.0083     709.9533       TRUE
F-I 3765.1977     804.5390       TRUE
F-J 5839.8562    1015.6357       TRUE
G-H 1896.2614     683.8760       TRUE
G-I 3359.4507     781.6237       TRUE
G-J 5434.1093     997.5813       TRUE
H-I 1463.1894     825.9834       TRUE
H-J 3537.8479    1032.7058       TRUE
I-J 2074.6585    1099.8776       TRUE

Jedoch, dieses Paket ermöglicht nur eine kategorische variable (z.B. ich kann nicht, studieren Sie die Preise gruppierten, die durch Farbe und von carat, wie ich es machen kann mit dem kruskal.test), und ich weiß nichts über die pgirmess Paket, ob es eingehalten wird oder nicht, oder wenn es getestet wird.

Können Sie mir empfehlen ein Paket zum ausführen der Kruskal-Wallis-test gibt details, die für jeden Vergleich? Wie würdest du das problem umgehen?

Es ist kruskal Funktion in agricolae Paket. Man könnte prüfen, ob Ihren Bedürfnissen passt.
Offensichtlich meinst du kruskal.test (ohne den zweiten .- test). Ich denke, Sie sind mit Diamanten dataset von ggplot2-Paket. Gut, ich kann nicht herausfinden, warum, aber wenn ich versuche es zu laden bekomme ich eine Fehlermeldung internal error -3 in R_decompress1 In addition: Warning message: restarting interrupted promise evaluation. (Ich weiß, das sieht aus wie eine Frage, die ich posten sollte, aber hat jemand vor dem gleichen problem?)
Marek: ich danke Ihnen sehr. Ich sah agricolae, aber es scheint das gleiche problem wie pgirmess 🙁 @gd047: ich habe keine Ahnung, was Ihre Fehler sein mögen, wirklich. Vielleicht sollten Sie versuchen, neu zu installieren ggplot2. Sie haben den gleichen Fehler mit anderen Datenbeständen oder Pakete?
Ich Frage mich nur, was macht Sie glauben, dass kruskal.test umgehen kann mehr als eine Gruppierung der Variablen. Es ist ein one-way-test, nachdem alle. Wenn Sie wollen, schauen alle Kombinationen von zwei Faktoren ab, erstellen Sie einfach eine neue Kombination mit interaction.

InformationsquelleAutor dalloliogm | 2010-03-19

Schreibe einen Kommentar