TukeyHSD
Tukey's Honest Significant Difference test is een simultane toetsings procedure. Wordt gebruikt om, na afloop van een variantieanalyse met een significant resultaat, te weten te komen welke van de betrokken steekproeven voor dat resultaat hebben gezorgd.
In R gaat die als volgt:
> bron <- "http://www.mzandee.net/~zandee/statistiek/data/haemo.txt"
> ziek <- read.table(bron, header=T)
> attach(ziek)
> names(ziek)
[1] "haemo" "ziekte"
> TukeyHSD(aov(haemo ~ ziekte, ziek))
Tukey multiple comparisons of means
95% family-wise confidence level
Fit: aov(formula = haemo ~ ziekte, data = ziek)
$ziekte
diff lwr upr p adj
B-A -3.25 -5.042354 -1.457646 0.0014748
C-A -6.50 -8.389307 -4.610693 0.0000074
C-B -3.25 -5.042354 -1.457646 0.0014748
In geen enkel lwr-upr interval zit het getal 0 ingesloten. Dat wil zeggen dat alle opgesomde paarsgewijze verschillen tussen de 3 groepen als significant zijn te beschouwen.
We kunnen het resultaat van TukeyHSD toekennen aan een variabele en vervolgens plotten, waarmee we deze paarsgewijze verschillen zichtbaar kunnen maken:
> ziekHSD <- TukeyHSD(aov(haemo ~ ziekte, ziek))
> plot(ziekHSD)
