服务热线
178 0020 3020
library(ggplot2)
library(ggthemes)
b<-ggplot(mtcars,aes(x=cyl,fill=factor(gear)))+geom_bar()+theme(legend.position="none")+labs(title="Here is title",subtitle="here is subtitle",caption="R2-37")
mytheme1<-theme(plot.title=element_text(face="bold.italic",size="14",color="brown"),axis.title=element_text(face="bold",size=10,color="black"),axis.text=element_text(face="bold",size=9,color="black"),axis.line=element_line(colour="black"),panel.grid.major=element_blank(),panel.border=element_blank(),panel.background=element_rect(fill="white"),panel.grid.minor.y=element_blank(),panel.grid.minor.x=element_blank(),legend.position="none")
b+mytheme1
library(ggplot2)
library(latex2exp)
myvolcano<- read.csv("/Users/Hong/Desktop/data/data8/nature23643-s4.csv",header=T)
myvolcano$significant<- ifelse(abs(log2(myvolcano$FC))>1&myvolcano$q_value<0.05,"yes","no")
myvolcano$mylabels<-ifelse(abs(log2(myvolcano$FC))>1&myvolcano$q_value<0.05,as.character(myvolcano$锘縂ene.ID),NA)
mysig<-subset(myvolcano,significant=="yes")
mysig
ydrop=-0.6
textdrop=0.3
myplot<-ggplot(myvolcano,aes(x=log2(FC),y=-log2(q_value),size=significant))+geom_point(col="grey70",alpha=1)+xlim(-2.5,2.7)+ylim(-0.8,15)+scale_size_manual(values=c(2,0))+geom_point(data=mysig,aes(x=log2(FC),y=-log2(q_value),col=mylabels),size=4,position=position_jitter(height=0.1,width=0))+geom_text(aes(label = mylabels),size =3,hjust = 0,vjust = -1,col = "black",fontface = "italic")+annotate("segment",x=1,xend=2.5,y=ydrop,yend=ydrop,color="black",size=1,arrow=arrow(angle=18,ends="last",type="closed"))+annotate("segment",x=-1,xend=-2.5,y=ydrop,yend=ydrop,color="black",size=1,arrow=arrow(angle=-2.5,ends="last",type="closed"))+annotate("text",x=1.7,y=textdrop,label="Upregulated",size=4.5)+annotate("text",x=-1.7,y=textdrop,label="Downregulated",size=4.5)+labs(title="CMTM6 sgRNA2 vs Control",caption="create by R2-37")+labs(x=expression(log[2]))+labs(x=expression(paste(log[2],"(fold change)")),y=expression(paste(-log[2],"(Q value)")))
请问我这里把labs(x=expression(log[2]))+labs(x=expression(paste(log[2],"(fold change)")),y=expression(paste(-log[2],"(Q value)")))改为labs(x=expression(paste(log[2],"(fold change)")),y=expression(paste(-log[2],"(Q value)")))可以出一样的图,为什么答案要加啊?群主能帮忙解答吗?谢谢
myplot+mytheme1
附件