关于画火山图（by ggplot2）的一些总结和经验

愿武艺晴小朋友一定得每天都开心！ 

---

文献中常用经典的火山图，是展示差异表达基因的利器。每次测完转录组，做实验组和对照组的比较后，都会用到。

我自己也画了不算太多也不算太少的次数。然后最近画的时候忽然间意识到这个可视化方法我经常用，却没系统的整理过，一些tips散落在各个角落里。值此契机，在此时此地此时间，我要把它们整理一遍，内容如下：

1）加载要用的R包

rm(list=ls())
options(stringsAsFactors = F)
library(readr)
library(gplots)
library(ggplot2)
library(org.Mm.eg.db)
library(enrichplot)
library(readxl)
library(patchwork)
library(dplyr)
library(ggrepel)

2)  完成火山图

DEGs <- **** #加载自己的Degs

# 2.1 定义水平
DEGs[which(DEGs$padj %in% NA),'change'] <- 'Not'
DEGs[which(DEGs$log2FC >= 1 & DEGs$padj < 0.05),'change'] <- 'Up&avg_express'
DEGs[which(DEGs$log2FC <= -1 & DEGs$padj < 0.05),'change'] <- 'Down&avg_express'
DEGs[which(abs(DEGs$log2FC) < 1