一、前置环境

1.1 R 语言

下载对应系统的R软件

R: The R Project for Statistical Computing (r-project.org)

以win11为演示

https://mirrors.tuna.tsinghua.edu.cn/CRAN/bin/windows/

1.2 R Studio

安装R Studio

https://download1.rstudio.org/electron/windows/RStudio-2022.12.0-353.exe

安装后启动并关联已安装的R语言

设置工作目录Session ----> Set Working Directory ---> Choose Directory

注意最好是全英文路径，避免奇奇怪怪的Bug

1.3 安装ggplot2 R包

##安装包

install.packages("ggplot2")

##加载包

library(ggplot2)

二、DE基因火山图绘制

2.1 准备文件

主要需要三列信息：Gene Log2FoldChange padj

可自行命名，其中

Log2FoldChange列的数值为Log2(FoldChange)

padj 是校正P值或者可以用原始P value

2.2 读取文件

##文件输入

data <- read.table("V_GV_GV.anno.xls_ggplot2",header = T)

##去除NA行

data_na <-na.omit(data)

##或将NA替换为其他值,如1，但出来的图会在底部多出许多异常点，不建议此操作

data[is.na(data)] <-  "1"

2.3 标记上下调基因

##设置阈值，其中padj< 0.05； Log2FoldChange) > 1；需相应更改

data_na$changed <- factor(ifelse(data_na$padj < 0.05 & abs(data_na$Log2FoldChange) > 1,ifelse(data_na$Log2FoldChange > 1,'Up','Down'),'Stable'))

##查看DE基因数

summary(data_na$changed)

2.4 绘图并保存

##绘图

ggplot(data_na,aes(Log2FoldChange,-log10(padj),
color=factor(changed),
)) +
geom_point() +
labs(x=expression("log2"*"(FoldChange)"), ##x轴Title
y=expression(-Log[10]*" (p value)")) + ##y轴Title
theme_grey(base_size = 10)+
scale_color_manual(values = c('green','grey','red'))+ ##对应分类点的颜色
scale_size_manual(values = c(1,1,1))+ ##对应分类点的Size
theme(legend.title = element_blank(),
legend.position = c(0.3,0.9),
legend.background = element_rect(fill='transparent'))

##保存

ggsave("V_GV_GV_mRNA.png",width=10,height=10,dpi=300)

参考文章：

R语言 ggplot2火山图 - 知乎 (zhihu.com)

(22条消息) R语言如何删除数据集中包含缺失数据NA的行？_ofoliao的博客-CSDN博客_r语言删除缺失值所在一行

ggplot2之二：图像的保存 - 知乎 (zhihu.com)

声明：

代码源于参考文章，进行相应更改，如有侵权，请联系删除，感谢各位的教程！

ggplot2绘制差异表达基因火山图相关推荐

BIC无代码绘制差异基因火山图
无代码绘制差异基因火山图 Volcano plot | 别再问我这为什么是火山图一文解释了火山图如何解读.不太难看懂,而一旦看懂了,图也就知道怎么绘制了. 假设我们已经有了一个差异基因鉴定后的表格文 ...
ImageGP/BIC无代码绘制差异基因火山图
无代码绘制差异基因火山图 Volcano plot | 别再问我这为什么是火山图一文解释了火山图如何解读.不太难看懂,而一旦看懂了,图也就知道怎么绘制了. 假设我们已经有了一个差异基因鉴定后的表格文 ...
差异表达基因火山图（ggplot函数）
1. 读入数据差异表达基因来自limma分析结果. # read the file data <- read.csv("diff_expr_genes.csv",row.n ...
差异表达基因-火山图和聚类图解释
想研究某现象的分子机制,老板豪气的来一句,先测个转录组吧,看下差异表达基因. 是否在心里窃喜,制个样就完事了,太easy有木有.等大堆数据回来的时候,是不是傻眼了? 从何下手挑选差异表达基因呢? 今天 ...
差异表达基因热图怎么看_【生信】憨憨版分析转录组基因表达差异（二）
[简单热图] # 调用下载的pheatmap包.library(pheatmap)#设置当前工作目录,setwd("目标路径")即保存待测文件的目录,使R能够调用.setwd(&q ...
R语言R原生及可视化包ggplot2绘制并排的箱图实战（Side-by-Side Boxplots）
R语言R原生及可视化包ggplot2绘制并排的箱图实战(Side-by-Side Boxplots) 目录 R语言R原生及可视化包ggplot2绘制并排的箱图实战(Side-by-Side Boxpl ...
差异表达基因热图怎么看_差异基因热图绘制：heatmap.2
在RNA-seq数据分析中,差异表达基因分析是一项基本的技能,其中热图又是一种特别常见的用来展示差异表达基因分析结果的方式,今天分享一个非常好用的绘制热图的R函数:heatmap.2.该函数来自gpl ...
一行代码绘制高分SCI火山图
一.概述在近半年中,我读了很多的高分SCI文章,很多文章中都有多种不同的火山图,包括「普通的火山图.渐变火山图.以及包含GO通路信息的火山图」! 经过一段时间的文献阅读和资料查询,终于找到了一个好用 ...
差异表达基因热图怎么看_R绘图雷达图-单基因泛癌差异表达的另类展现形式
往期回顾: R语言学习系列之"多变的热图" 蚂蚁金服在线可视化引擎 G2 R绘图:无与伦比的华丽风行(桑基图) R绘图:相关性分析与作图(单基因相关性) R绘图:相关性分析与作图R ...

ggplot2绘制差异表达基因火山图