跟着Cell学作图 | R语言绘制多组差异火山图

数据和代码获取：请查看主页个人信息！！！

大家好！欢迎来到R语言数据分析视界。相信大家对火山图的绘制方法已经并不陌生，我们前面的推文也有介绍过。

火山图适合展示两个分组之间的差异基因/代谢物等信息；然而，我们在进行实验设计时，往往会纳入多个分组。如果我们想看所有分组之间差异情况，就需要绘制多个火山图。今天我将向大家介绍如何使用ggplot2绘制多组火山图，用于同时展示多组之间差异情况。

相关性热图样式灵感来源于Cell杂志的一篇文章；绘图代码参考：https://zhuanlan.zhihu.com/p/516955474

   

接下来我们来进行分析和可视化展示，首先载入本次绘图数据：

Step1：数据载入

rm(list=ls())pacman::p_load(tidyverse, ggrepel, reshape2, rio)# 载入数据DEG_limma_voom <- read.table("DEG_limma_voom.txt", header=T, sep="\t")

接下来，我们对数据进行处理，剔除不显著的数据。我们使用mutate函数创建一个新的列"change"，根据给定的条件给基因分类为"Up"（上调）或"Down"（下调），并过滤掉"No change"的基因