差异表达基因火山图（ggplot函数）

1. 读入数据

差异表达基因来自limma分析结果。

# read the file
data <- read.csv("diff_expr_genes.csv",row.names=1)
# sorting
data = diff_genes[order(diff_genes$adj.P.Val),]
head(data)
class(data)colnames(data)
dim(data)# 如果之前保存的是R对象
save(nrDEG_limma_voom, file = 'diff_expr_genes')
load('diff_expr_genes')
data <- nrDEG_limma_voom

2. ggplot作火山图

# 颜色设定
data$color <- ifelse(data$adj.P.Val<0.05 & abs(data$logFC)>= 1,ifelse(data$logFC > 1,'red','blue'),'gray')#install.packages("ggplot2")
#install.packages("ggrepel")# 导入包
library(ggplot2)
library(ggrepel) # labels 不重叠#tiff(filename = "volcano.tif",width = 720, height = 720)
rownames(data) <- data$symbol# 设定要标出的基因名
data$sign <- ifelse(data$adj.P.Val < 1e-10 & abs(data$logFC) >4,rownames(data),NA)p <- ggplot(data, aes(logFC, -log10(adj.P.Val), col = color)) +geom_point() +theme_bw() +scale_color_manual(values = color) +labs(x="log2 (fold change)",y="-log10 (q-value)") +geom_hline(yintercept = -log10(0.05), lty=4,col="grey",lwd=0.6) +geom_vline(xintercept = c(-1, 1), lty=4,col="grey",lwd=0.6) +theme(legend.position = "none",panel.grid=element_blank(),axis.title = element_text(size = 18),axis.text = element_text(size = 14)) +geom_text_repel(aes(label = sign), size = 3)
print(p)#dev.off()
# ggsave保存
ggsave("ggsave_volcano.tif",dpi=100,plot=p,device = 'tiff')#device:    "eps", "ps", "tex" (pictex), "pdf", "jpeg",
#"tiff", "png", "bmp","svg" or "wmf" (windows only)

差异表达基因火山图（ggplot函数）相关推荐

ggplot2绘制差异表达基因火山图
一.前置环境 1.1 R 语言下载对应系统的R软件 R: The R Project for Statistical Computing (r-project.org) 以win11为演示 http ...
差异表达基因-火山图和聚类图解释
想研究某现象的分子机制,老板豪气的来一句,先测个转录组吧,看下差异表达基因. 是否在心里窃喜,制个样就完事了,太easy有木有.等大堆数据回来的时候,是不是傻眼了? 从何下手挑选差异表达基因呢? 今天 ...
差异表达基因热图怎么看_【生信】憨憨版分析转录组基因表达差异（二）
[简单热图] # 调用下载的pheatmap包.library(pheatmap)#设置当前工作目录,setwd("目标路径")即保存待测文件的目录,使R能够调用.setwd(&q ...
BIC无代码绘制差异基因火山图
无代码绘制差异基因火山图 Volcano plot | 别再问我这为什么是火山图一文解释了火山图如何解读.不太难看懂,而一旦看懂了,图也就知道怎么绘制了. 假设我们已经有了一个差异基因鉴定后的表格文 ...
ImageGP/BIC无代码绘制差异基因火山图
无代码绘制差异基因火山图 Volcano plot | 别再问我这为什么是火山图一文解释了火山图如何解读.不太难看懂,而一旦看懂了,图也就知道怎么绘制了. 假设我们已经有了一个差异基因鉴定后的表格文 ...
差异表达基因热图怎么看_差异基因热图绘制：heatmap.2
在RNA-seq数据分析中,差异表达基因分析是一项基本的技能,其中热图又是一种特别常见的用来展示差异表达基因分析结果的方式,今天分享一个非常好用的绘制热图的R函数:heatmap.2.该函数来自gpl ...
差异表达基因热图怎么看_R绘图雷达图-单基因泛癌差异表达的另类展现形式
往期回顾: R语言学习系列之"多变的热图" 蚂蚁金服在线可视化引擎 G2 R绘图:无与伦比的华丽风行(桑基图) R绘图:相关性分析与作图(单基因相关性) R绘图:相关性分析与作图R ...
差异表达基因热图怎么看_多变的热图1（新手专用）
热图(heatmap)用不同的颜色和颜色的深浅来直观的展示数据之间的差异.在测序类的文章里,几乎必有一幅热图用来展示差异表达基因.很多工具都可以完成热图的制作,今天这篇文章主要介绍利用R语言的 phe ...
差异表达基因热图怎么看_画图专题 | 热图函数pheatmap()
画图专题 | 热图函数pheatamp() 百迈客医学 2018-10-31 先看一眼这个函数的参数,这么多,而且最后还有省略号.那么我们应该怎么合理使用这些参数让你的热图看起来更加高大上呢? phe ...

差异表达基因火山图（ggplot函数）

1. 读入数据

2. ggplot作火山图

差异表达基因火山图（ggplot函数）相关推荐

最新文章

热门文章