R语言作方差分析

#单因素方差分析
rice_yields <- read.csv(file = './data/rice_yields.csv',header = TRUE)
head(rice_yields)

R语言中的回归诊断-- car包

#通过QQ图来检验正态性
install.packages('car',repos = 'https://mirrors.ustc.edu.cn/CRAN/')

library(car)
with(rice_yields,qqPlot(lm(产量~品种)))

=#检查方差齐性
bartlett.test(rice_yields$产量~rice_yields$品种)

data:  rice_yields$产量 by rice_yields$品种
Bartlett's K-squared = 2.0554, df = 7, p-value = 0.9567

#进行方差分析
result_rice<-aov(rice_yields$产量~rice_yields$品种)
summary(result_rice)

install.packages('gplots',repos = 'https://mirrors.ustc.edu.cn/CRAN/')

#绘制均值图
library(gplots)
plotmeans(产量~品种,data = rice_yields)

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！