均值和标准误计算

AB <- read.table(file = "~/Desktop/R/data1.txt",header = TRUE) AB

head(AB);str(AB)

#计算均值和标准误

library(stringr)

mean_se <- vector()

#⾃自定义函数：计算标准误

func.se <- function(data){

n.t <- length(data)

n.na <- length(which(is.na(data)))

n <- n.t - n.na

sd <- sd(data,na.rm = TRUE)

se <- sd/sqrt(n)

return(se)

}

#利利⽤用循环计算均值和标准误并使⽤用±符号拼接

for (i in1:nrow(AB)) {

data.i <- as.numeric(AB[i,2:4])

mean.i <- round(mean(data.i,na.rm = TRUE),3)

se.i <- round(func.se(data.i),3)

mean_se[i] <- str_c(mean.i,se.i,sep = "±")

}

AB$mean_se <- mean_se

#导出数据

write.csv(AB,file = "data_treated.csv")