将GO、Pathway富集结果整合在一张高颜值圆圈图上

富集分析是生物医学论文中非常常见的一类分析,例如GO富集分析,Pathway富集分析等。其结果一般包括以下几个要素:1,名字(GO term或者KEGG description);2,该名字所包含的基因数目;3,该名字所包含的基因与进行富集分析的输入基因的重叠数目;4,富集的P值、FDR值等。富集结果的常见绘图方式包括:气泡图,条形图,弦图等。今天我们利用circlize R包整合富集结果,绘制一个高大上的圈图。

                                    

                                                           图1.GO、Pathway富集圈图

如图1所示,从外到内共4个track,中间的为图例(从上到下依次为:分类,p值,上调或者下调)

1分类track

相同颜色的为同一分类,例如黄色表示BP,蓝色表示MF,最外圈刻度表示总的基因数,10的次方表示。总的基因数即图2中所有黄色圈所包含的基因去重后的基因数,或者推荐使用比所有条目中基因数最大值更大的一个值一般设置为整百或者整千)

2) term里边的基因数

矩形长度表示该term中包含的基因数(例如GO:0031625这个term中包含306个基因),与最外圈的刻度成比例。矩形的颜色表示该term的富集p值(或者FDR值,q值等,经过-log10转化),颜色越深表示P越小。

3)重叠基因数

矩形长度表示该term中包含的基因与进行基因富集分析时所输入的基因的重叠基因数。例如进行富集分析的基因有500个,其中与GO:0031625中306个基因重叠的基因为38个。由于这个数字一般较小,因此矩形的长度与该track中最的基因数成比例。一般我们要么使用up的基因,要么使用down的基因进行富集分析,因此这里仅一种颜色。

4)富集分数bar

Bar的高度即第三个track里的基因数占第二个track里对应基因数的比例(如图2所示)。数值范围为0到1,这里每个小圈表示0.2

                                                      图2. 基因数及比例来源示意图

1,打开绘图页面

首先,使用浏览器(推荐chrome或者edge)打开富集圈图绘制页面。左侧为常见作图与分析导航,中间为数据输入框和可选参数,右侧为描述和结果示例。也可以在主页搜索框中搜索GO,找到绘图页面。

微生信 - 在线绘制GO,pathway富集结果圆形图

3.富集圈图绘图页面

2,示例数据

点击右侧“示例数据”链接下载excel格式的示例数据。

                              

                                                       图4. 输入数据示例

示例数据(仅供参考)包括7列:

第1列:GO或者KEGG条目ID。推荐使用ID,因为名字长短不一,很难统一。

2分类名。可以包含1类,2类,3类,或者4类,例如单独绘制KEGG通路富集结果;同时绘制GO-BPKEGG;同时绘制BP、CC、MF;甚至同时绘制BP、CC、MF、KEGG。同一分类放在一起,不要相互穿插。

3总基因数。总基因数即图2中所有黄色圈所包含的基因去重后的基因数,或者推荐使用比所有条目中基因数最大值更大的一个值。例如第4列中最大的数值是470,这里可以使用1200,用来从整体上控制第二个track的矩形长度,是一个比较松散的值一般设置为整百或者整千)

第4列:数据库中某term包含的基因数。在使用同一个数据库版本的情况下,这个数值是固定不变的,例如GO:0010498里边就包含470个基因。

第5列:某term富集的p值(或者q值,FDR值等)

第6列做富集分析时输入基因与该term中包含基因的重叠基因数(见图2)。

第7列第6列第4列的比值(见图2)。

3,粘贴示例数据

直接复制示例数据中的A-G列数据,然后粘贴到输入框。

                                 

                                                               图5. 必需输入

注意:不是拷贝excel文件,是拷贝excel文件里边的数据。另外粘贴到输入框后,格式乱了没关系,只要在excel中是整齐的就行。并且数据矩阵中不能有空的单元格,中文字符等。

4,修改参数,并提交

我们设置了图片尺寸,字体大小,颜色,字体类型等参数,基本能满足日常绘图使用。如需更高级的定制,请联系我们。

                                

                                                     图6.字体、颜色等可调参数

5,提交出图

粘贴好输入数据,调整好参数(或者全部默认)后,点击提交按钮,约8秒后,会在页面右侧出现预览图。我们提供了4种图片格式供下载使用,两种矢量图(pdf,svg)和两种标量图(600 dpi tiff和300 dpi png)。

                                   

                                                                     图7.预览与下载

 Tips:

1,该图最大限度地利用有限的空间展示更多的数据内容。因此条目不能过多,推荐总共30个以内的GO和KEGG条目,否则太拥挤。

2,第4个track最大值默认为1.0,可以修改成其他值,以整体上降低或者降低柱子的高度。

3,一般我们使用up或者down的基因做富集分析,对于使用up+down富集结果不适合这个图。

4,对于DAVID等富集分析结果,各列的名字可能不一样,需要根据数值代表的意思调整成示例数据即可。

没有预览就是没有出图,这时请参考示例数据,检查自己输入数据的格式。

遇到文字截断,需要修改字体、调整字体大小等,使用 科研作图实操:用inkscape编辑svg矢量图

微生信助力高分文章,用户66000+,引用870+

余丁,微生信
关注 关注
  • 2
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
R语言|GO富集分析
weifanbio的博客
04-19 9322
GO富集分析 —— 差异圈图、聚类圈图 一、数据链接 链接: https://pan.baidu.com/s/1PPUW5YyJHjwxvkjmMi0Vtw?pwd=c9s8 提取码: c9s8 二、安装包介绍 clusterProfiler包:功能也比较强大,主要是做GO和KEGG的功能富集及其可视化。 org.Hs.eg.db包:转换NCBI、ensemble等数据库中基因ID,symbol等之间的转换。 enrichplot包:实现多种可视化方法来解释富集结果。 GOplot:功能富集
R语言GO富集分析
一直在水些技术小文
06-12 4249
GO富集分析 —— 差异圈图、聚类圈图一、数据链接 链接: https://pan.baidu.com/s/1PPUW5YyJHjwxvkjmMi0Vtwpwd=c9s8 提取码: c9s8二、安装包介绍 clusterProfiler包:功能也比较强大,主要是做GO和KEGG的功能富集及其可视化。 org.Hs.eg.db包:转换NCBI、ensemble等数据库中基因ID,symbol等之间的转换。 enrichplot包:实现多种可视化方法来解释富集结果。 GOplot:功能富集。三、绘 1、安
零代码在线绘制GO富集结果GOcircle
微生信
05-14 2533
基因本体论(Gene Ontology,GO,http://www.geneontology.org/)是国际基因功能分类系统,旨在建立一套随着研究动态更新的可控词汇表来限定和描述各个物种的基因和蛋白质功能。该数据库将收录的基因本体按照术语描述的内容不同,将所有基因本体分为三大类:生物过程(Biological Process)、细胞组分(Cellular Component)、分子功能(Molecular Function)。GO的基本单位是term,每个term对应一个功能类别或细胞定位,一个基因可能对
GEO生信数据挖掘(八)富集分析(GO 、KEGG、 GSEA 打包带走)
最新发布
zzh1464501547的博客
10-16 3111
例子中结核类型,包括结核,潜隐进展,对照和潜隐,四个类别。本节对差异基因进行富集分析。
单细胞测序流程(九)单细胞的GO圈图
qq_45478665的博客
07-22 5512
系列文章目录 提示:这里可以添加系列文章的所有文章的目录,目录需要自己手动添加 例如:第一章 Python 机器学习入门之pandas的使用 提示:写完文章后,目录可以自动生成,如何生成可参考右边的帮助文档 文章目录 单细胞测序流程(一)简介与数据下载 单细胞测序流程(二)数据整理 单细胞测序流程(三)质控和数据过滤——Seurat包分析,小提琴和基因离差散点 单细胞测序流程(四)主成分分析——PCA 单细胞测序流程(五)t-sne聚类分析和寻找marker基因 单
五维联合展示GO、Pathway富集分析结果
微生信
08-07 4611
GO和Pathway富集分析是功能分析不可或缺的一部分,GO和Pathway富集分析结果的展示也是五花八门:条形,气泡,山脊等等。由于GO分成三部分:生物过程(Biological process)、细胞组分(Cellular component)和分子功能(Molecular function),因此结果展示中需要进行区分。今天微生信为大家带来一张漂亮的多维GO、Pathway富集分析结果。...
GO和KEGG富集分析详细步骤
tqptr_opqww的博客
05-20 2万+
GO和KEGG富集分析 文章目录GO和KEGG富集分析@[toc]1. 将差异表达结果的基因名称转化为id2. GO富集分析3. GO圈图绘制4. KEGG富集分析5. KEGG圈图绘制 1. 将差异表达结果的基因名称转化为id 因为GO和KEGG分析需要用到id,所以这一步需要将基因名字转换为id。具体步骤如下: 新建空白文件夹,将差异分析得到的diff.xls复制粘贴到文件夹中 因为在这里只需要diff.xls中的基因名称和logFC两列,所以只复制这两列粘贴到新建的文本文件symbol.tx
【单细胞高级绘】10.KEGG富集结果的圆圈
TOP生物信息
09-12 2558
本次教程的figure仍然是读者求助的,算得上是kegg富集的新流派。据我的调查,该应该是基迪奥云平台首创(https://www.omicshare.com/tools/Home/Soft/enrich_circle),之后公众号小白鱼的生统笔记进行了复现(仿一个网,使用circlize包绘制圈图可视化基因集富集分析结果)。 最开始也是跟着上述的帖子学习,之后自己对代码进行了改写,重新安排形的布局,使之(在我看来)更有意义。另一个改动是增加了kegg pathway的注释信息,我在之前的帖子中
富集分析原理和clusterProfiler包进行GO、KEGG富集分析详细说明
热门推荐
Tian問的博客
10-03 2万+
概念: 基因富集分析是指对于给定一组基因根据基因组注释信息(GO、KEGG)对基因进行聚类分析,即给定的基因是不是GO中的一个功能(或KEGG中的一个通路)。 基因的功能富集的目的说明给定的基因集对哪些功能的影响有针对性的,不是随机影响的。 原理: 基因富集分析是通过研究给定的基因集在功能节点上是否过出现来得到关注的基因集显著注释的功能节点。通常利用超几何分布等方法计算给定基因集在某个功能(或通路)上的P值,判断给定的基因集在功能(或通路)的基因数目超过了随机的期望,是一个小概率事件。 变
pathway:有点像一个状态机,用于在 iOS 上浏览不同的视控制器
06-22
有点像一个状态机,用于在 iOS 上浏览不同的视控制器。 有时,您有一组 ViewController,它们通过您的应用程序在各种工作流程中重复使用。 你真的不应该在你的 ViewControllers 中存储所有这些状态,因为这会变得...
pathway-mapper:用于癌症途径的交互式协作形策展工具
04-29
运行本地实例为了部署和运行该工具的本地实例,请按照以下步骤操作: 首先,将PathwayMapper克隆到本地计算机,然后导航到本地存储库: 安装git clone https://github.com/iVis-at-Bilkent/pathway-mapper.gi
pathway:用Go编写的基于路径的API服务器
05-18
path-用Go编写的基于路径的API服务器参见文章。构建并运行 go build./pathway测试方法调用 curl -H "Content-Type: application/json" \ -d '{"test": 1}' http://0.0.0.0:4000/api/echo/message增加服务提供了一个...
EnrichedHeatmap:制作丰富的热,以可视化的方式将基因组信号富集到特定的目标区域
05-06
富集是一种特殊类型的热,它可以可视化特定目标区域上基因组信号的富集。 它广泛用于可视化例如组蛋白标记如何丰富到特定位点。 有几种工具可以制作此类热(例如或 )。 在这里,我们通过包实现了丰富的...
微生信平台桑基气泡绘制教程
03-31
微生信平台桑基气泡以5个维度展示基因富集分析结果。视频介绍了该的意义,并使用示例数据展示了如何绘制,及如何进行后期编辑。
☞GO和KEGG富集倍数(Fold Enrichment)如何计算 enrich factor qvalue
生信小博士的博客
07-27 8761
比如说,某一类功能基因的总数是100个,其中差异基因有10个,那么只有实验只是影响了该功能基因的10%,但是如果该实验一共只鉴定到20个差异基因,那么这10个差异基因就占所有差异基因的50%,因此,需要同时考虑差异基因在所有功能分类中的分布,得到的结果才是准确的。其中,a为功能A中差异基因的数量,b为非差异基因的数量,c为所有功能分类中差异基因的总数量,d为非差异基因的总数量,n为识别到的所有差异基因的数目。今天我们来聊聊怎么来展示GO富集分析的结果,下面是一个GO富集分析的结果表。
R语言-Go分析结果的富集
new_bug1234的博客
11-14 863
这篇文章的Fig4,作者在对进行EOSPE和LOSPE进行GO富集分析后,为了总结过多的GO条目,结合了Cytoscape和Jaccard得分,进行了“二次富集”,分为了5+2组。由于文章中没有细述相关步骤,这里就主要推断下计算Jaccard的步骤,仅抛砖引玉,以及提出相应的解决办法,希望大家可以批评指正。接下来就是Cytoscape部分不再详述,参考部分有更好的文章。至于如何挑选的这四组文章中没有细述,按表达是结合了Cytoscape的绘,或许也有先验知识的参与,各位大佬们有什么见解吗?
单细胞测序流程(十)单细胞的kegg富集分析和圈图(终章)
qq_45478665的博客
07-23 1万+
系列文章目录 文章目录 单细胞测序流程(一)简介与数据下载 单细胞测序流程(二)数据整理 单细胞测序流程(三)质控和数据过滤——Seurat包分析,小提琴和基因离差散点 单细胞测序流程(四)主成分分析——PCA 单细胞测序流程(五)t-sne聚类分析和寻找marker基因 单细胞测序流程(六)单细胞的细胞类型的注释 单细胞测序流程(七)单细胞的细胞类型轨迹分析单细胞测序流程(八)单细胞的marker基因转化和​GO富集分析 单细胞测序流程(九)单细胞的GO圈图
使用R语言绘制富集条形,轻松分析基因表达数据
笑不语的博客
06-24 2821
通俗来说,富集分析通过将基因分类到特定的集合中,然后根据基因在集合中的分布和总体分布的比较,来寻找哪些集合与特定的生物过程、疾病或其他功能相关联。随着大数据时代的到来,富集分析与富集条形成为了分析工具的基石之一,如生物医学信息学中的疾病诊断、药物研发、基因功能验证等等,都需要富集条形进行数据展示,以辅助我们进行快速准确的数据分析与解读。富集条形是一种可视化富集分析结果的工具。此外,绘制富集条形时,还可以根据不同的颜色映射,显示富集通路的基因数目和富集通路的显著水平,更直观地为数据提供了参考。
R语言:GO富集和KEGG富集
08-23
R语言中可以使用许多不同的包来进行GO富集和KEGG富集分析。下面是一个使用clusterProfiler包进行富集分析的示例代码: 1. 安装和加载必要的包: ```R install.packages("clusterProfiler") library(clusterProfiler) ``` 2. 对基因列表进行GO富集分析: ```R # 假设有一个基因列表gene_list,包含感兴趣的基因 result_GO <- enrichGO(gene = gene_list, OrgDb = org.Hs.eg.db, # 数据库,例如Human keyType = "ENSEMBL", # 基因ID类型 ont = "BP", # 生物过程 pAdjustMethod = "BH", # 多重检验校正方法 pvalueCutoff = 0.05, # 显著性水平 qvalueCutoff = 0.2) # FDR阈值 # 查看富集结果 result_GO ``` 3. 对基因列表进行KEGG富集分析: ```R result_KEGG <- enrichKEGG(gene = gene_list, organism = "hsa", # 物种,例如Human pvalueCutoff = 0.05, # 显著性水平 qvalueCutoff = 0.2) # FDR阈值 # 查看富集结果 result_KEGG ``` 这些代码会返回富集分析的结果,包括富集的GO term和KEGG pathway,以及统计学上的显著性。你可以根据需要进一步处理和可视化这些结果。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
写文章

热门文章

  • 数据为什么要进行log2转化,倍数变化(fold change)为什么要求个log2FC? 11388
  • 桑吉气泡图 -- KEGG富集气泡图升级版,5维展示富集结果 10513
  • 输入序列就可以预测 - 在线版miRanda靶基因预测来了 8975
  • 左右双侧条形图展示基因本体论(GO)BP、CC、MF三合一富集分析结果C、MF三合一富集分析结果 8186
  • 微生信在线绘制KEGG Pathway通路分类汇总图 7450

分类专栏

  • pathway通路 3篇
  • 富集分析 3篇
  • PATHWAY 3篇
  • 富集结果弦图 1篇
  • KEGG 1篇
  • GO 2篇
  • 论文字体 1篇

最新评论

  • 推荐4个可用的github国内镜像

    m0_65703986: 为什么我都显示无法打开网页啊表情包表情包表情包

  • 根据蛋白质序列,计算其分子量(molecular weight),在线工具,原理和python代码

    普通网友: 引领技术潮流,是不可多得的好文,十分值得借鉴和参考。期待博主未来能够持续分享更多好文【我也写了一些相关领域的文章,希望能够得到博主的指导,共同进步!】

  • 批量将GOID转成GO term名并添加BP,MF,CC分类信息

    CSDN-Ada助手: 恭喜你这篇博客进入【CSDN每天值得看】榜单,全部的排名请看 https://bbs.csdn.net/topics/618632051。

  • 输入序列就可以预测 - 在线版miRanda靶基因预测来了

    lyllyloooo: 请问填入的mRNA3'UTR是不需要反向互补吗?

  • TCGAplot在线版:输入基因,一键绘制TCGA基因表达与MSI相关性雷达图

    CSDN-Ada助手: 推荐 Python入门 技能树:https://edu.csdn.net/skill/python?utm_source=AI_act_python

大家在看

  • 在AMD GPUs上构建解码器Transformer模型 1303
  • 美丽世界(首尔学联组织三校联谊、马拉松、雪岳山、MT) 218
  • 星空探索(原创)

最新文章

  • ToxinPred – 多肽毒性预测、突变设计和理化性质预测
  • piRNA-疾病研究最新综述 – 数据库、软件工具、webserver和计算模型
  • circlize R包在线绘制圆形热图,bug修复版
2024年15篇
2023年6篇
2022年60篇
2021年5篇

目录

目录

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43元 前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值

聚圣源男孩瑞字起名字好听宠物猪起名起名网智能取名大全古代给男孩子起名手机归属给快餐鱼饭店起名龚自珍的诗句山水民宿起名给自己起个小名字公司起名子查询abc3602019魏姓女孩起名新生儿宝宝起名网意外的幸运签舟起名的含义起名字后面什么字好屌丝道士中餐厅起名游戏加速器哪个好起打分100分的的好名字官场新贵下载airpods起名电视家apk重生之第二帝国千百鲁金石良缘剧情介绍伍姓女宝宝起名当代中学生报非正式探险笔记服装起商标名字淀粉肠小王子日销售额涨超10倍罗斯否认插足凯特王妃婚姻让美丽中国“从细节出发”清明节放假3天调休1天男孩疑遭霸凌 家长讨说法被踢出群国产伟哥去年销售近13亿网友建议重庆地铁不准乘客携带菜筐雅江山火三名扑火人员牺牲系谣言代拍被何赛飞拿着魔杖追着打月嫂回应掌掴婴儿是在赶虫子山西高速一大巴发生事故 已致13死高中生被打伤下体休学 邯郸通报李梦为奥运任务婉拒WNBA邀请19岁小伙救下5人后溺亡 多方发声王树国3次鞠躬告别西交大师生单亲妈妈陷入热恋 14岁儿子报警315晚会后胖东来又人满为患了倪萍分享减重40斤方法王楚钦登顶三项第一今日春分两大学生合买彩票中奖一人不认账张家界的山上“长”满了韩国人?周杰伦一审败诉网易房客欠租失踪 房东直发愁男子持台球杆殴打2名女店员被抓男子被猫抓伤后确诊“猫抓病”“重生之我在北大当嫡校长”槽头肉企业被曝光前生意红火男孩8年未见母亲被告知被遗忘恒大被罚41.75亿到底怎么缴网友洛杉矶偶遇贾玲杨倩无缘巴黎奥运张立群任西安交通大学校长黑马情侣提车了西双版纳热带植物园回应蜉蝣大爆发妈妈回应孩子在校撞护栏坠楼考生莫言也上北大硕士复试名单了韩国首次吊销离岗医生执照奥巴马现身唐宁街 黑色着装引猜测沈阳一轿车冲入人行道致3死2伤阿根廷将发行1万与2万面值的纸币外国人感慨凌晨的中国很安全男子被流浪猫绊倒 投喂者赔24万手机成瘾是影响睡眠质量重要因素春分“立蛋”成功率更高?胖东来员工每周单休无小长假“开封王婆”爆火:促成四五十对专家建议不必谈骨泥色变浙江一高校内汽车冲撞行人 多人受伤许家印被限制高消费

聚圣源 XML地图 TXT地图 虚拟主机 SEO 网站制作 网站优化