1

我正在使用包中的HSD.test功能进行 TukeyHSD 事后测试。该函数工作正常，但我不确定 p 值隐藏在哪里。中的字母表示重要性，但实际的 p 值在哪里？谢谢agricolaeRgroups

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
4

1 回答 1

6

更改group为FALSE返回比较而不是紧凑型字母显示。

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
于 2017-11-07T22:37:49.063 回答