2017-11-07 4 views
0

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 
+0

@PoGibas Sicher muss es oder wie würden die Buchstaben bestimmt werden? –

+0

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

Antwort

3

Ändern group-FALSE kehrt die Vergleiche statt der Kompaktbrief-Display.

out = HSD.test(model, "virus", group = FALSE) 
print(out$comparison) 

     difference pvalue signif.   LCL   UCL 
cc - fc 11.5333333 0.0686  . -0.8663365 23.9330031 
cc - ff -11.9333333 0.0592  . -24.3330031 0.4663365 
cc - oo -12.5000000 0.0482  * -24.8996698 -0.1003302 
fc - ff -23.4666667 0.0014  ** -35.8663365 -11.0669969 
fc - oo -24.0333333 0.0012  ** -36.4330031 -11.6336635 
ff - oo -0.5666667 0.9988   -12.9663365 11.8330031 
+0

Erstaunlich! Ich danke dir sehr. Das ist genau das, was ich brauche. –