数据处理过程中，在计算完各个指标各组的均值和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 |
