这篇文章将为大家详细讲解有关如何利用R语言的ggplot包绘制火山图，小编觉得挺实用的，因此分享给大家做个参考，希望大家阅读完这篇文章后可以有所收获。
成都创新互联公司成立十多年来，这条路我们正越走越好，积累了技术与客户资源，形成了良好的口碑。为客户提供成都网站设计、网站建设、网站策划、网页设计、主机域名、网络营销、VI设计、网站改版、漏洞修补等服务。网站是否美观、功能强大、用户体验好、性价比高、打开快等等，这些对于网站建设都非常重要，成都创新互联公司通过对建站技术性的掌握、对创意设计的研究为客户提供一站式互联网解决方案，携手广大客户，共同发展进步。
火山（Volcano Plot）图在一张图中显示了两个重要的指标（Fold change/pvalue），可以非常直观且合理地筛选出在两样本间发生差异表达的基因。检验分析出两样本间显著差异表达的基因后，以log2（fold change）为横坐标，以T检验显著性检验P值的负对数-log10（pvalue）为纵坐标，即可得火山图（Volcano Plot）。
> install.packages("ggplot2")
> library(ggplot2)
> R0_vs_R3_isoforms_fitter <- read_excel("~/Downloads/R0-vs-R3.isoforms.fitter.xlsx")
> data <- R0_vs_R3_isoforms_fitter
> r03 <-ggplot(data,aes(log2FC,-1*log10(FDR)))
> r03 + geom_point()
运行过程数据及图片：
改变点的颜色：
> r03 + geom_point(color="red")
> r03 + geom_point(aes(color="red"))
> r03 + geom_point(aes(color=significant))
设置标题，坐标轴范围
>r03+geom_point(aes(color=significant))+xlim(-4,4)+ylim(0,30)
>r03xy=r03+geom_point(aes(color=significant))+xlim(-4,4)+ylim(0,30)
> r03xy + labs(title="Volcano plot",x="log2(FC)")
> r03xy + labs(title="Volcanoplot",x=expression(log[2](FC),y=expression(-log[10](FDR))))
自定义颜色
> r03xy + labs(title="Volcanoplot",x=expression(log[2](FC),y=expression(-log[10](FDR))))
> r03xyc = r03xy+labs(title="Volcanoplot",x=expression(log[2](FC)), y=expression(-log[10](FDR)))
> r03xyc + scale_color_manual(values =c("green","black", "red"))
>r03xyc + scale_color_manual(values =c("green","blue", "red"))
> volcano = r03xyc +scale_color_manual(values = c("green","blue","red"))
>volcano+geom_hline(yintercept=1.3)+geom_vline(xintercept=c(-1,1))
>volcano+geom_hline(yintercept=1.3,linetype=4)+geom_vline(xintercept=c(-1,1),linetype=4)
保存图片:
>ggsave("volcano.png")
>ggsave("volcano8.png",volcano,width=8,height=8)
关于“如何利用R语言的ggplot包绘制火山图”这篇文章就分享到这里了，希望以上内容可以对大家有一定的帮助，使各位可以学到更多知识，如果觉得文章不错，请把它分享出去让更多的人看到。
重庆分公司
重庆分公司
