R2-27-第七次作业-火山图

whk 2017-12-20 00:25:17 阅读: 992

library(ggplot2)
library(reshape2)
library(ggthemes)
mtcars
w=ggplot(mtcars,aes(factor(cyl)))+geom_bar(aes(fill=factor(gear)))+labs(x="cyl",y="count",title="R2-27")+scale_fill_discrete(name="gear")
mytheme=w+theme(panel.grid.major =element_blank(), panel.grid.minor = element_blank(),panel.background = element_blank(),axis.line = element_line(colour = "black"))

library(latex2exp)

data=read.csv("D:\\nature23643-s4.csv")      

head=(data)

data$significant = ifelse(abs(log2(data$FC)) > 1 & data$q_value < 0.05,"yes","no" )

data$mylabels =ifelse(abs(log2(data$FC)) > 1 & data$q_value < 0.05,as.character(data$锘縂ene.ID),NA)

mysig = subset(data,significant == "yes")

mysig

ydrop = -0.6

textdrop = 0.3

myplot = ggplot(data,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 = 18,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-27") + labs(x = expression(log[2])) +labs(x = expression(paste(log[2],"(fold change)")),y = expression(paste(-log[2],"(Q value)")))

myplot

mytheme

Rplot16.pngRplot15.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-2024 科研好助手,京ICP备20005780号-1 建议意见

服务热线

178 0020 3020

微信服务号