我想要每个阶段的输出_生存图。特别是假根生命阶段的'iii',我已经尝试过但只有在excel中分离才能得到。现在很费时间我有三种治疗的所有数据。
我知道我必须从舞台上重新创建另一个列表——这就是我在浪费时间。 然后我基本上想要进行营养处理 (PES) 并确认阶段 - 哪一天 PES 的治疗(光照强度)是第一天,存活率是多少?
附上的图片是我一直在各个阶段使用的,但一直给我不同的情节。在我逐步提取信息之前,我现在已经回到原始数据集。
enter code here
数据集<- summarySE(DataAll,na.rm=TRUE, measurevar='Survival', groupvars=c("Location","Stage","time"))
Datase$Location<-factor(Datase$Location, labels=c("P0","P5","P10","P15","P20","PR"))
Datase$Stage<-factor(Datase$Stage, labels=c("i","ii","iii","iv","v","vi"))
Datase$time<-factor(Datase$time, labels=c("1","2","3","4","5","6","7","8","10","14"))
P0 = ggplot(Datase[Datase$Location=="P0",],aes(colour= factor(Stage), x=time, y=Survival,group=(Stage),shape = factor(Stage))) + theme_bw()+xlab("Time")+ylab(expression(paste("Stage Prop", " ", sep="")))+theme(axis.text.x = element_blank())+ geom_line(position=position_dodge(0),size = 0.5)+ theme(text = element_text(size=14))+facet_grid(.~Location, scales="free_y")+ theme(strip.background = element_blank(), strip .text = element_blank()) +geom_errorbar(position=position_dodge(0),aes(ymin=Survival-se,ymax=Survival+se),width=.1,colour="black") + geom_point(position=position_dodge(0),size = 5)+ theme(text = element_text(size=14))+guides(shape = guide_legend(reverse=F))+ theme(legend.position="none",legend. key.size = unit(0.1, "mm"))+facet_grid(.~Location, scales="free_y") +geom_errorbar(position=position_dodge(0),aes(ymin=Survival-se,ymax=Survival+se) ,width=.1,colour="黑色") + facet_wrap(~Location) + ggtitle("PES - 0 umol") + coord_cartesian(ylim=c(0,100)) + 主题(panel.grid.major = element_blank(),panel.grid.minor = element_blank(),panel.background = element_blank(),axis.line = element_line(颜色=“黑色”))