什么是倾斜45度的火山图?
生信宝典
共 4804字,需浏览 10分钟
· 2021-04-23
今天是转录组培训最后一天,剩下一点时间自由探索。朋友们来自山南海北,还是北京的最多。
一位老师聊起火山图(Volcano plot | 别再问我这为什么是火山图 (在线轻松绘制)),说见过倾斜45度的类似图,可否演示怎么画?想了下,可能是下面这种图,绘起来看看。
检查和安装包
a = rownames(installed.packages())
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager", repos = site)
a = rownames(installed.packages())
install_bioc <-
c(
"egg",
"ggpubr",
"ggplot2"
)
for (i in install_bioc) {
if (!i %in% a)
BiocManager::install(i, update = F)
}
if (!"ImageGP" %in% a){
# devtools::install_github("Tong-Chen/ImageGP")
devtools::install_git("https://gitee.com/ct5869/ImageGP.git")
}
library(ImageGP)
library(ggplot2)
library(ggpubr)
library(egg)
读入数据并标记差异基因
数据是基于DESeq2分析的差异基因数据。
file <- "data/ehbio.simplier.DESeq2.trt._vs_.untrt.results.txt"
diffexpr <- sp_readTable(file, row.names=1)
# 做一个log转换
diffexpr$trt <- log2(diffexpr$trt+1)
diffexpr$untrt <- log2(diffexpr$untrt+1)
head(diffexpr)
## trt untrt baseMean log2FoldChange pvalue padj
## ENSG00000152583 10.88130 6.354646 983.042 4.546 1.219e-91 2.149e-87
## ENSG00000189221 12.09256 8.706410 2391.559 3.387 9.955e-61 8.779e-57
## ENSG00000179094 10.40316 7.343044 757.249 3.065 2.435e-54 1.432e-50
## ENSG00000116584 10.50370 11.567470 2242.427 -1.064 3.957e-49 1.745e-45
## ENSG00000120129 12.55314 9.579082 3386.426 2.975 1.930e-48 6.807e-45
## ENSG00000134686 12.04738 10.587715 2884.847 1.460 1.846e-45 5.427e-42
标记差异基因
diffexpr$level <- ifelse(diffexpr$padj<0.05,
ifelse(diffexpr$log2FoldChange>=1, "trt up",
ifelse(diffexpr$log2FoldChange<=-1, "untrt up", "NoSig")),"NoSig")
head(diffexpr)
## trt untrt baseMean log2FoldChange pvalue padj level
## ENSG00000152583 10.88130 6.354646 983.042 4.546 1.219e-91 2.149e-87 trt up
## ENSG00000189221 12.09256 8.706410 2391.559 3.387 9.955e-61 8.779e-57 trt up
## ENSG00000179094 10.40316 7.343044 757.249 3.065 2.435e-54 1.432e-50 trt up
## ENSG00000116584 10.50370 11.567470 2242.427 -1.064 3.957e-49 1.745e-45 untrt up
## ENSG00000120129 12.55314 9.579082 3386.426 2.975 1.930e-48 6.807e-45 trt up
## ENSG00000134686 12.04738 10.587715 2884.847 1.460 1.846e-45 5.427e-42 trt up
基于ImageGP绘图
p <- sp_scatterplot(diffexpr, xvariable = "trt", yvariable = "untrt",
color_variable = "level",
color_variable_order = c("NoSig","trt up", "untrt up"),
manual_color_vector = c("grey","firebrick","dodgerblue"),
legend.position = c(0.8,0.3)) + coord_fixed(1)
p
到这满足需求了,又有老师说能不能加上数据分布展示?
拼上Marginal plot
xplot <- ggplot(diffexpr, aes(x=trt)) + geom_histogram(fill="firebrick") +
theme_classic() +
theme(axis.line.x=element_blank(),
axis.ticks.x=element_blank(),
axis.text.x = element_blank(),
axis.title.x = element_blank(),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
panel.border = element_blank(),
panel.background = element_blank(),
legend.title = element_blank(),
legend.position = c(0.75,0.85),
#legend.title = element_text(size = 5),
legend.text = element_text(size = 8),
legend.key.size = unit(0.5, "lines"),
legend.spacing = unit(0.3, "cm"),
) + ylab("Density")
xplot
yplot <- ggplot(diffexpr, aes(x=untrt)) + geom_histogram(fill="dodgerblue") +
theme_classic() +
theme(axis.line.y=element_blank(),
axis.ticks.y=element_blank(),
axis.text.y = element_blank(),
axis.title.y=element_blank(),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
panel.border = element_blank(),
panel.background = element_blank(),
legend.title = element_blank(),
legend.position = c(0.7,0.8),
#legend.title = element_text(size = 5),
legend.text = element_text(size = 8),
legend.key.size = unit(0.5, "lines"),
legend.spacing = unit(0.3, "cm"),
) + ylab("Density") +
rotate()
yplot
合并起来
white <- ggplot() + theme_void()
egg::ggarrange(xplot, white, p, yplot,
widths=c(5,2),heights = c(2,5),padding=unit(0,"line"))
也可以用现成的工具
用ggpubr绘制
ggscatterhist(diffexpr, x="trt", y="untrt", color="level",
palette=c("grey","firebrick","dodgerblue"), margin.plot="histogram" )
输出数据用在线工具绘制
见次条
http://www.ehbio.com/Cloud_Platform/front/#/analysis?page=b%27MTI%3D%27
sp_writeTable(diffexpr, file="DE_gene_label.txt")
数据和代码见 https://gitee.com/ct5869/shengxin-baodian/Plot
往期精品(点击图片直达文字对应教程)
后台回复“生信宝典福利第一波”或点击阅读原文获取教程合集
评论
多人同时导出 Excel 干崩服务器!新来的阿里大佬给出的解决方案太优雅了!
点击关注公众号,Java 干货及时推送↓推荐阅读:面试辅导,我们出大成果了!来源:juejin.cn/post/7259249904777838629前言 业务诉求:考虑到数据库数据日渐增多,导出会有全量数据的导出,多人同时导出可以会对服务性能造成影响,导出涉及到mysql查询的io操作,
Java技术栈
1
有意思!一个关于 Spring 历史的在线小游戏
发现 Spring One 的官网上有个好玩的彩蛋,分享给大家!进到Spring One的官网,可以看到右下角有个类似马里奥游戏中的金币图标。点击该金币之后,会打开一个新的页面,进入下面这样一个名为:The History Of Spring 的在线小游戏你可以使用上下左右的方向键来控制Spring
公众号程序猿DD
1
了解加密货币到加密货币的互换
1、什么是加密货币互换?加密货币到加密货币的互换是指以现行市场汇率将一种加密货币直接兑换为另一种加密货币。与需要法定货币存款和较长流程的传统交易所不同,加密货币到加密货币的互换可以无缝地促进交换。掉期在提高加密货币的流动性和效率方面发挥着重要作用。该功能使用户能够将他们的加密货币与钱包中的其他代币进
区块链头条
0
李彦宏:开源大模型不如闭源,后者会持续领先;周鸿祎:“开源不如闭源” 的言论是胡说八道
架构师大咖
架构师大咖,打造有价值的架构师交流平台。分享架构师干货、教程、课程、资讯。架构师大咖,每日推送。
公众号该公众号已被封禁0、李彦宏:开源大模型不如闭源,后者会持续领先当今
源码共读
0
【第129期】程序员的新宠:三款终端工具,让你告别Xshell!
概述 WindTerm:跨平台的SSH利器 首先介绍的是WindTerm,这是一款使用C语言开发的跨平台SSH客户端。它不仅完全免费,而且没有商业使用的限制。WindTerm支持SSH v2、Telnet、Raw Tcp等协议,而且性能出色,甚至超过了FinalShell和Electerm。功能
前端微服务
0
字节员工:35岁以后被裁员的,后来都走了哪条路?现在2-2,要不要利用最后一年拼命上个岸。
架构师大咖
架构师大咖,打造有价值的架构师交流平台。分享架构师干货、教程、课程、资讯。架构师大咖,每日推送。
公众号该公众号已被封禁在当今竞争激烈的职场环境中,年龄并不总是一个决定性
源码共读
0
顶级 Javaer 都在用的 20 个类库,真香!
点击关注公众号,Java 干货及时推送↓推荐阅读:投了 100 多份简历后…优秀且经验丰富的Java开发人员的特征之一是对API的广泛了解,包括JDK和第三方库。我花了很多时间来学习API,尤其是在阅读了Effective Java 3rd Edition之后 ,Joshua Bloch建
Java技术栈
1
985学历真好用,一面答的再差也能过?
JavaGuide官方网站:javaguide.cn相对其他行业如律师、金融等,程序员求职面试对学历的包容度还是相对较高的,尤其是在过去几年。通常情况下,如果候选人在其他方面(如实习经历、项目经历、竞赛经历)表现突出,可以在一定程度上弥补学历的不足。不过,由于现在面试越来越卷,一些大厂、国企和研究所
JavaGuide
1