【懒人必备】批量完成均值和SD计算,并以mean±SD的格式批量
2019-08-07 本文已影响0人
Cdudu
数据处理过程中,在计算完各个指标各组的均值和SD(或者SE)之后,就需要以mean±SD的形式输入到表格中,当一个试验有n张表格,n个指标的时候,一个个输入就会特别麻烦。因此用R写了一小段代码,用于批量完成这个工作。
数据准备
先将要处理的所有数据按列进行整理,将应变量放在最前面
group | BW | SGR | CAT | SOD |
---|---|---|---|---|
1 | 4 | 1.2 | 21 | 2.2 |
1 | 5 | 1.4 | 22 | 2.3 |
1 | 6 | 1.5 | 23 | 2.5 |
2 | 14 | 2.2 | 33 | 3.5 |
2 | 15 | 2.4 | 34 | 3.6 |
2 | 16 | 2.5 | 35 | 3.7 |
3 | 24 | 3.2 | 43 | 4.4 |
3 | 25 | 3.4 | 33 | 4.6 |
3 | 26 | 3.5 | 45 | 4.8 |
批量操作
#############################################
##
## 批量求均值和SD
##
##-----------------------------------------
## Author : Cdudu
## Data : 2019 8/7
############################################
library(readxl)
# 读入需要统计的数据表
dat<-read_excel('Data.xls')
summary(dat)
# 创建新表用于存放Mean和SD,此处as.character()中向量长度由组数决定,
# 例如上表共三组,就输入as.character(1:3)
dat2<-data.frame(t1=as.character(1:3))
for ( i in c(2:5)){ #此处向量根据应变量的列数决定,上表应变量在2-5列,因此输入c(2:5)
means<-tapply(dat[[i]],dat$group,mean)
means<-sprintf('%.2f',round(means,2)) #改变数字可以改变Mean的小数点位数
SD<-tapply(dat[[i]],dat$group,sd)
SD<-sprintf('%.2f',round(SD,2)) #改变数字可以改变SD的小数点位数
a<-paste(means,'±',SD)
dat2[i-1]<-a
names(dat2)[i-1]<-names(dat[,i])
}
dat3<-t(dat2) #行列转置
write.csv(dat3,'Mean&SD.csv')
以上就是全部代码,实际操作中只需要两步操作
1.输入处理组数量
dat2<-data.frame(t1=as.character(1:3))
2.输入应变量数量
for ( i in c(2:5)){
结果展示
如果最后是将dat3输出到EXCEL,那么表格如下
NULL | V1 | V2 | V3 |
---|---|---|---|
BW | 3±0.2 | 4±0.3 | 5±0.4 |
SGR | 3±0.2 | 3±0.2 | 3±0.2 |
CAT | 3±0.2 | 3±0.2 | 3±0.2 |
SOD | 3±0.2 | 3±0.2 | 3±0.2 |
V1;V2;V3即是3个试验组,里面的数值是我瞎编的。
当然也可以将dat2输出到EXCEL,那么表格如下
NULL | BW | SGR | CAT | SOD |
---|---|---|---|---|
V1 | 3±0.2 | 3±0.2 | 3±0.2 | 3±0.2 |
V2 | 3±0.2 | 3±0.2 | 3±0.2 | 3±0.2 |
V3 | 3±0.2 | 3±0.2 | 3±0.2 | 3±0.2 |