ANOVA-in-R

1
2
3
4
5

#安装一个田间实验设计方面的包，为我们提供LSD.test()等好用的函数
install.packages("agricolae")
library(agricolae)
##读入数据
data(sweetpotato)

1

oneway.test()

1

model <- aov(yield~virus, data = sweetpotato)

1
2
3
4
5
6

#多重比较，不矫正P值
out <- LSD.test(model, "virus", p.adj = "none")
##标记字母法显示比较结果
out$group
#可视化
plot(out)