R2-37-第七期-火山图

洪学志 2017-12-20 00:10:58 阅读: 1295

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

8-1-3.png

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

8-2.png

 

 
邀请讨论

附件

{{f.title}} 大小 {{f.file_size}} 下载 {{f.count_download}} 金币 {{f.count_gold}}
{{item.nick_name}} 受邀请回答 {{item.create_time}}
{{item.refer_comment.nick_name}} {{item.refer_comment.create_time}}

附件

{{f.title}} 大小 {{f.file_size}} 下载 {{f.count_download}} 金币 {{f.count_gold}}
切换到完整回复 发送回复
赞({{item.count_zan}}) 踩({{item.count_cai}}) 删除 回复 关闭
科研狗©2015-2025 科研好助手,京ICP备20005780号-1 建议意见

服务热线

178 0020 3020

微信服务号