【懒人必备】批量完成均值和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
上一篇下一篇

猜你喜欢

热点阅读