bonferroni') 8另一种多重比较方法 TukeyHSD(result.aov) plot(TukeyHSD(result.aov)) 9非参数检验 kruskal.test(y$V1,a1) #两两比较 pairwise.wilcox.test print(kruskal.test(y$V1,a1)) print(pairwise.wilcox.test(y$V1,a1,p.adjust.method = 'bonferroni')) sink