R语言 | GO富集GOCircle弦表图绘制

文摘   2024-11-21 08:02   云南  

一边学习,一边总结,一边分享!

由于微信改版,一直有同学反映。存在长时间接收不到公众号的推文。那么请跟随以下步骤,将小杜的生信筆記设置为星标,不错过每一条推文教程。

欢迎关注《小杜的生信笔记》!!

如何加入社群

小杜的生信笔记仅有微信社群

1. 微信群:付费社群。添加小杜好友,加友请知:加友须知!!,加入社群请查看小杜的生信笔记付费加友入群声明

入群声明

2. 小杜个人微信:若你有好的教程或想法,可添加小杜个人微信。值得注意的是,小杜个人微信并不支持免费咨询长时间咨询,但支持小问题2-3个免费咨询。

小杜微信:

知识星球:

Bioinfo Notes学习交流群② 成立啦

基于长期积累和各位大佬、同学的支持,我们Bioinfo Notes学习交流群①已经满员,群②已成立,详情可看新的征程,学习永远不晚!Bioinfo Notes生信交流群②成立!

本期教程

引言

在以前的教程中,出过多个GO富集或KEGG富集分析及绘图教程,如R语言绘图 | GO富集柱状图,R语言绘图 | 多组GO和KEGG富集结果,分面堆积柱状图,富集桑基气泡图绘制教程,R语言绘图 | GO、KEGG等富集棒棒图,aPEAR包绘制功能富集网络图,绘制差异基因富集柱状图等。但是,对于GOCircle弦表图的绘制教程是空缺的。

读图

左图:灰色背景标识DEG的LogFC值,包含up和down基因。中间表示每个GO term的Z-score值。右图:表示GO term。

Z-score值计算:(up-down)/squrt(count)。

Code

# install.packages("GOplot")
library(GOplot)

使用事例数据

#'@来自与R包的事例数据
data(EC)
#使用 circle_dat 函数生成用于绘图的数据对象
circ<-circle_dat(EC$david, EC$genelist)
# 绘制 circ# 绘制 GOCircle 图
##'@绘制基础图形
p <- GOCircle(circ,
nsub = 5, # 显示前5个GO条目
label.size = 3, # 标签字体大小
table.legend = TRUE, # 显示图例
zsc.col = c("blue", "red"), # LogFC值的颜色梯度
lfc.col = c("blue", "red")) # LogFC值的颜色梯度

ggsave("./20241121.GOplot.jpg",plot = p,width = 8,height = 4)

使用自己的富集结果进行会绘图

我们在使用R绘图时,会存在使用事例数据,可以正常分析和绘图,但是一旦使用自己数据,就无法运行,以及很多教程也只给你事例数据的代码。

加载富集数据集

setwd("D:\\小杜的生信笔记\\2024\\20241121_Goplot绘制富集图")
# 导入对应的数据
GO_data <- read.table("01.GO.Results.txt",header = T, sep = "\t")
GO_data[1:5,1:5]
FC_data <- read.table("02.DE_logFC.txt",header = T, sep = "\t")
FC_data[1:5,1:2]

生成用于绘图的数据对象

小杜的生信筆記
小杜的生信筆記,生信小白,初来乍到请多指教。 主要学习分享,转录组数据分析,基于R语言数据分析和绘制图片等,以及相关文献的分享。
 最新文章