“站长,火山图咋画?”
火山 站长
2023-06-13 09:15:26 时间
火山图 就是比热图更热的图~~~
少废话,直接上代码~
基于R.3.4.4
#载入相关包
library(ggplot2)
library(ggrepel)
library(ggsci)
library(dplyr)
#提取数据,resG是注释过基因名字的DEseq2结果
Vol.TG<-na.omit(resG)
#给结果加上change,其中padj和log2FoldChange后面的值根据自己的需要修改
Vol.TG$change <- as.factor(ifelse(Vol.TG$padj < 0.01 & abs(Vol.TG$log2FoldChange) > 1,ifelse(Vol.TG$log2FoldChange > 1,'UP','DOWN'),'NOT'))
#加个配色
mypal_V = pal_gsea("default", n = 3, alpha = 0.9)(3)
#开始设置ggplot
#其中geom_point中alpha是点的透明度,size是点大小都可以修改。
#其中xlim是x轴范围,自己根据数据调整
#theme_bw中base_size是图中文字字号
#theme中那些事调整格子的,想知道效果删掉一个看一下。
#scale_color_manual中values后面值是配色,可以用c("","","")自己填颜色代号调整。limits代表change的顺序,这个要配合颜色顺序调整。
p <- ggplot(data = Vol.TG, aes(x = log2FoldChange, y = -log10(padj), color = change))
+geom_point(alpha=0.8, size = 1)
+xlim(-9,9)+theme_bw(base_size = 18)
+theme(panel.grid.minor = element_blank(),panel.grid.major = element_blank())
+scale_color_manual(name = "", values = mypal_V, limits = c("DOWN", "NOT", "UP"))
#输出p
p
#加上基因名,其中padj和log2FoldChange后面的值根据自己的需要修改
p+geom_text_repel(data=filter(Vol.TG,padj < 0.01,abs(log2FoldChange)>6), aes(label=genename))
效果图
为了让大更好的掌握这个用法,Chris生命科学小站qq群里有resG的例表给大家练习用。如果想知道这个resG是怎么得到的,请看下面星球介绍。
相关文章
- 从抖音到火山引擎——看流媒体技术演进和机会
- 火山图标记基因_火山地形图
- R语言绘制火山图_r语言画曲线图
- 小站工具网站功能总结:查询基因在肿瘤中的表达情况+火山图+GSEA+Survival+ROC+ENSG注释,
- 云化生长,火山引擎的“云原生”在讲些什么?
- 火山引擎爆发!字节跳动推出800G RoCE网络和自研DPU卡
- 对话火山引擎王悦:多媒体继续向更高清,更交互,更沉浸演进
- 干货 | 如何打造企业专属A/B平台?火山引擎DataTester开放平台技术揭秘
- 基拉韦厄火山正在重塑夏威夷!熔岩大爆发有多恐怖
- 古亚洲人74000年前被超级火山喷发“团灭”?新研究显示:并没有
- 字节跳动火山引擎加入 Linux 云原生计算基金会(CNCF)
- 火山安卓系统攻克MSSQL连接瓶颈(火山安卓连接mssql)