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)總結
- 上一篇: 奉化哪里有笋罐头空瓶批发?
- 下一篇: 卡冻结了怎么办