anova方差分析

首先同质性检验，Levene

（1）方差齐性时（检验结果sig>0.05）,选用 LSD.test {agricolae} p.adj = "bonferroni"，要进行矫正:

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

out <- LSD.test(model,"virus", p.adj="bonferroni")

out$group 为分组结果

（2）方差不齐性时，选用welch.test {onewaytests}：

out <- welch.test(Sepal.Length ~ Species, data = iris)

paircomp(out) #分组结果

其实，对比发现，即便是在方差不齐性时，两个方法，结果一致。（由于不是很懂方法原理，所以不知道检验结果一致是不是偶然现象）