原创 TS的美梦 KS科研分享与服务 2022-04-28 10:40

这次真的和火山图杠上了，我就不信了还复现不完奇奇怪怪的火山图？在一篇文章中看到了一种呈现方式，其实大同小异，但是还是练习下有些画法。此外，小伙伴们又看到新奇的火山图，可截图联系我们，争取将火山图打下来。

来源：[1]S García-Silva, A Benito-Martín, S Sánchez-Redondo, et al. Correction: Use of extracellular vesicles from lymphatic drainage as surrogate markers of melanoma progression and BRAF V600E mutation[J]. Journal of Experimental Medicine, 2019, 216(5):jem.2018152204162019c.

加载数据和ggplot2

setwd("E:/生物信息学/ggplot火山图")

df <- read.csv("DEG.csv",header = T)

library(ggplot2)

接下来作图即可。需要完成的是画出散点图，加上分割线，还有对不同的区域填充颜色，这个做法我们之前的文章提到过（不好好做图的NSC系列(八):ggplot2重现Nature文章多组柱状图+散点，这是一场硬仗），最后在合适位置加上文字标注即可。

#画图

library(ggplot2)

ggplot(df, aes(x=avg_log2FC, y=-log10(p_val_adj))) +

geom_hline(aes(yintercept=1.3), color = "#999999", linetype="dashed", size=1)+

geom_vline(aes(xintercept= -0.25), color = "#999999", linetype="dashed", size=1)+

geom_vline(aes(xintercept= 0.25), color = "#999999", linetype="dashed", size=1)+

geom_rect(aes(xmin=-Inf,

xmax=-0.25,

ymin=0,

ymax=Inf),

fill='blanchedalmond',color='grey90')+

geom_rect(aes(xmin=Inf,

xmax=0.25,

ymin=0,

ymax=Inf),

fill='blanchedalmond',color='grey90')+

geom_point(size=3, shape=21, color="black") +

labs(x = "Log2 fold change",y = "-Log10(P-value)", title = "") +

theme_bw() +

theme(panel.grid.major=element_blank(),

panel.grid.minor=element_blank(),

panel.border = element_blank(),

axis.line = element_line(colour = "black",size = 0.5),

axis.title =element_text(size = 14),

axis.text =element_text(size = 12, color = "black"))+

annotate("text", label = "Down-regulated DEGs\nin treatment",

x = -2.2, y = 250, size=5)+

annotate("text", label = "Up-regulated DEGs\nin treatment",

x = 2.2, y = 250, size= 5)

今天的分享就到这了，欢迎小伙伴发图，一起把他啃下来！

如有侵权，请联系本站删除！