这篇文章主要介绍“r语言怎么进行单因素方差分析”的相关知识,小编通过实际案例向大家展示操作过程,操作方法简单快捷,实用性强,希望这篇“r语言怎么进行单因素方差分析”文章能帮助大家解决问题。
library(agricolae)data(sweetpotato)
可以看到, 数据为两列, 第一列为不同病毒的类型, 第二列为产量, 为了研究不同病毒感染对产量的影响.
> head(sweetpotato) virus yield1 cc 28.52 cc 21.73 cc 23.04 fc 14.95 fc 10.66 fc 13.1
因素: virus
变量: yield
model<-aov(yield~virus, data=sweetpotato)summary(model)
可以看到不同病毒达到显著性水平.
> summary(model) Df Sum Sq Mean Sq F value Pr(>F)virus 3 1170.2 390.1 17.34 0.000733 ***Residuals 8 179.9 22.5---Signif. codes:0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
这里多重比较采用LSD方法.
out <- LSD.test(model,"virus", p.adj="bonferroni")out$groups
结果可以看到, oo, ff,cc之间不显著(都有a), oo,ff与fc之间显著(字母没有交集).
> out$groups yield groupsoo 36.90000 aff 36.33333 acc 24.40000 abfc 12.86667 b
plot(out)
关于“r语言怎么进行单因素方差分析”的内容就介绍到这里了,感谢大家的阅读。如果想了解更多行业相关的知识,可以关注亿速云行业资讯频道,小编每天都会为大家更新不同的知识点。
亿速云「云服务器」,即开即用、新一代英特尔至强铂金CPU、三副本存储NVMe SSD云盘,价格低至29元/月。点击查看>>
免责声明:本站发布的内容(图片、视频和文字)以原创、转载和分享为主,文章观点不代表本网站立场,如果涉及侵权请联系站长邮箱:is@yisu.com进行举报,并提供相关证据,一经查实,将立刻删除涉嫌侵权内容。
原文链接:https://my.oschina.net/u/4592498/blog/4473032