用R语言进行ANOVA分析

打开.rda（rdata的简写）文件，可以用R.studio打开

#R
#对于比较不同组之间是否存在差别，用单因素方差分析法（ANOVA）

# ll
# 2019-03-14
#对于安装报错的包，可以修改镜像为china(lanzhou)
install.packages("car", dependencies=TRUE, INSTALL_opts = c('--no-lock'))
install.packages("multcomp", dependencies=TRUE, INSTALL_opts = c('--no-lock'))

library(readxl)
library(gplots)
library(car)
library(multcomp)

# read excel 打开第二个工作薄 数据格式为第一列是类别名，第二列是相对应的值 
data<-read_xlsx('C:\\Users\\guosheng\\Desktop\\name.xlsx',sheet=2)

# 统计各个类别的个数 （行标题为equ，value）
# 利用attach()来加载数据，不然在用table函数时会提示找不到对象'equ'
attach(data)
table(equ)

#各组均值
result_mean <-aggregate(val,by = list(equ),FUN= mean)

#各组标准差
result_sd <-  aggregate(val,by = list(equ),FUN= sd)

#检验组间差异（ANOVA）
fit<-aov(val~equ)
summary(fit)

# 绘制各组均值及其置信区间的图形
plotmeans(val ~ equ,xlab = 'nam