Ich verwende die HSD.test
-Funktion aus dem agricolae
-Paket in R
, um einen TukeyHSD-Post-hoc-Test durchzuführen. Die Funktion funktioniert gut, außer ich bin mir nicht sicher, wo die p-Werte versteckt sind. Die Buchstaben in groups
bedeuten Signifikanz, aber wo sind die tatsächlichen p-Werte? DankeP-Werte von HSD.test
library(agricolae)
data(sweetpotato)
model<-aov(yield~virus, data=sweetpotato)
out <- HSD.test(model,"virus", group=TRUE,console=TRUE,
main="Yield of sweetpotato\nDealt with different virus")
Study: Yield of sweetpotato
Dealt with different virus
HSD Test for yield
Mean Square Error: 22.48917
virus, means
yield std r Min Max
cc 24.40000 3.609709 3 21.7 28.5
fc 12.86667 2.159475 3 10.6 14.9
ff 36.33333 7.333030 3 28.0 41.8
oo 36.90000 4.300000 3 32.1 40.4
Alpha: 0.05 ; DF Error: 8
Critical Value of Studentized Range: 4.52881
Minimun Significant Difference: 12.39967
Treatments with the same letter are not significantly different.
yield groups
oo 36.90000 a
ff 36.33333 ab
cc 24.40000 bc
fc 12.86667 c
@PoGibas Sicher muss es oder wie würden die Buchstaben bestimmt werden? –
Gute Frage. Ich bin mir nicht sicher, wie man es mit diesem Befehl macht. Versuchen Sie 'TukeyHSD (model)', das ein Basis-R-Befehl ist. Dies gibt Ihnen p-Werte (angepasst) für paarweise Vergleiche. – AntoniosK