先上图
image.png
数据形式
1、表格
image.png
2、分组（由于R识别不了“—”，自动换为"."）
image.png
*****代码*****
library(reshape2)
library(Rmisc)
library(ggplot2)
library(reshape2)
qinya <- read.csv("xiao.CSV",header = T)
group <- read.csv("group.CSV",header = T)
##############################
chouqinya <- function(data= data1,group = group,num = 1, aa="a1",bb="b1", cc="c1",dd="d1",ee="e1" ,ff="f1" , x1=1,x2=1,x3=1,x4=1,x5=1,x6=1,x7= 1){
list1 <- list()
namespaste <- c()
names <- c(aa,bb,cc,dd,ee,ff)
date <- c(x1,x2,x3,x4,x5,x6,x7)
b = 0
a = 0
c = 0
for (i in 1:num) {
namespaste[i] <- paste0(names[i],".csv")
if (i == 1) {
list1[[i]] <- data[1:date[i],]
}
if(i >1 ){
a = date[i-1]
b = a + b
c = b+ date[i]
list1[[i]] <- data[(b+1):c,]
}
}
liste <- list()
picture <- c()
for (i in 1:length(list1)) {
row.names(list1[[i]]) <- list1[[i]]$date
eg <- list1[[i]][,-1]
eg <- as.data.frame(t(eg))
eg$sample <- factor(rownames(eg), levels = rev(rownames(eg)))
eg <- melt(eg, id = 'sample')
eg <- merge(eg,group,by="sample")
eg$value <- as.numeric(eg$value)
liste[[i]] <- summarySE(eg, measurevar="value", groupvars=c("group","variable"))
write.csv(liste[[i]],file=namespaste[i],row.names = F)
plota <- liste[[i]]
p1 <- ggplot(plota, aes(x=factor(variable), y=value, colour=group,group=group,shape=group,fill=group)) +
geom_line() +geom_point(size=4)+geom_errorbar(aes(ymin=value-se, ymax=value+se), width=.1)+
theme(panel.grid.major =element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(),axis.line = element_line(colour = "black"))
ggsave(paste0(names[i],".pdf"),p1,device = "pdf",width = 13,height = 9,dpi = 400)
}
}
chouqinya(data = qinya, group = group ,num = 5,aa ="jd" ,bb ="cod" ,cc ="ad" ,dd ="ph" ,ee = "orp",
x1 =3 ,x2 =3 ,x3 =3 ,x4 =22 ,x5 = 23)
#data = qinya 代表数据框
#num =1 代表数据种类
#aa,bb等代表数据名字
#x1,x2等代表每种数据的行数
