首先同质性检验,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) #分组结果 其实,对比发现,即便是在方差不齐性时,两个方法,结果一致。(由于不是很懂方法原理,所以不知道检验结果一致是不是偶然现象) |
|
来自: 昵称27943849 > 《R》