常见的基因功能富集分析方法可以认为分两代。

(1)第一代:基于目标基因集预筛选的功能富集分析方法

基本步骤包括两步:

(a)从背景基因集合,按照一定固定阈值(例如,是否差异显著)筛选目标基因集。这属于一个二分类的策略,例如将基因分为差异表达与无差异表达的。

(b)通过一些统计方法(例如,超几何检验)证明目标基因集(例如差异表达基因)的基因更多集中在某些通路中。

第一代富集分析方法,简单易懂,所以我们常见的结题报告基本上都使用这类方法。但第一代富集分析方法最大的不足就是:筛选目标基因集的标准(例如差异分析中的差异倍数大于2 & FDR <0.05)会直接影响富集分析的结果,尤其当差异基因很少的时候,可能无法富集到任何通路。

(2)第二代:基于所有基因评分或排序的功能富集分析方法

第二代功能富集分析方法则可以很好弥补第一代方法的不足,最有代表性的方法就是GSEA(Gene Set Enrichment Analysis)的方法。以GSEA为例,分析也可以分为两步:

(a)利用某个指标(例如差异倍数)对所有基因进行排序(而不是用固定阈值筛选)。

(b)判断某个通路的基因在所有的排序中是否更为靠前(或靠后),来证明实验处理会该通路产生更显著的影响。

GSEA法因为不需要对数据进行预筛选,所以很好地解决了一代方法的不足,对通路富集的检测能力更加敏感(PS:在基迪奥最新版本的转录组结题报告中也加入了GSEA富集分析的内容)。

想了解更多关于GSEA的内容,请到基迪奥Omicshare class搜索关键词“GSEA”收看教学视频:

GSEA分析可以直接使用单机JAVA版的GSEA软件完成,但该软件的一大不足就是绘制出的GSEA图上不了台面。

JAVA版的GSEA软件图形输出格式是png格式,主要有问题包括:

(1)像素偏低,发文章略显清晰度不够(应该没有到达300dpi);

(2)图有些丑,例如周边是灰色,背景加网格线(当然,也看每个人喜好);

(3)输出格式是位图,不利于后期编辑。

那么如何解决呢?比较简单的解决方案是:GSEA  + R的策略:

(1)用GSEA软件完成分析

(2)基于分析结果,用R语言绘图的

在本期Omicshare小课堂中,我们演示如何用一个R脚本(我们提供现成的脚本)绘制一个更简洁清晰的GSEA plot。代码非常简单,一学就会,如下:

>source("gseaplot_modified.r")

>replotGSEA("my_analysis.GseaPreranked_result","HALLMARK_TNFA_SIGNALING_VIA_NFKB","AML")

然后就可以绘制出如下的GSEA plot:画面简洁、像素足够且方便后期修图。

是不是很简单,那就点击文章底部 阅读原文 登录免费学习吧~

最后,

本节视频课程为基迪奥春季转录组培训班的课程一小部分,如果想要学习GSEA更多的内容,可以参加基迪奥线下生信培训班。

如何用r语言搜集报表_基迪奥免费小课堂——如何用R语言绘制GSEA plot相关推荐

  1. 语言迭代发_从今天开始,每天学点R语言~

    首先,在这里先跟大家说声对不起,技术团队最近一直在做课程开发,本着宁缺毋滥的原则,我们的微信文章一再搁置,在编辑部催了无数遍之后,终于可以把课程放出来给大家了. 为了增加课程的实用性,在前期课程的基础 ...

  2. python语言运行效率高_如何评价说PYTHON是最快的语言?

    我是视频的作者,首先感谢朋友们的关注,并且诞生了这个论题,让我可以看到大家对我视频中的观点的反馈.对于大家的讨论,在日常时间允许的条件下,我会尽量参与.下面就我目前已经看到的朋友们的评论,做一下展开. ...

  3. 用python画圆角矩形_详解微信小程序canvas圆角矩形的绘制的方法

    详解微信小程序canvas圆角矩形的绘制的方法 发布时间:2020-10-04 18:20:31 来源:脚本之家 阅读:103 作者:清夜 微信小程序允许对普通元素通过 border-radius 的 ...

  4. canvas 圆角矩形填充_详解微信小程序canvas圆角矩形的绘制的方法

    微信小程序允许对普通元素通过 border-radius 的设置来进行圆角的绘制,但有时候在使用 canvas 绘图的时候,也需要圆角,例如需要将页面上某块区域导出为图片下载到本地的时候,常用的解决方 ...

  5. python语言编写爬虫_自写Python小爬虫一个 - 『编程语言区』 - 吾爱破解 - LCG - LSG |安卓破解|病毒分析|www.52pojie.cn...

    [Python] 纯文本查看 复制代码import requests import re import os #类型列表 Type_list = {1:'qingchun',2:'xiaohua',3 ...

  6. 04747 java语言程序设计(一)_自考软件工程:2017年4月04747Java语言程序设计(一)历年真题及答案...

    Java语言程序设计(一)历年真题试卷及答案由湖南自考生网www.zikaosw.com收集.提供.更多真题及答案,可在我办网站"自考历年真题"栏目中查看.软件工程历年真题建议 同 ...

  7. java语言就业方向_四大就业方向彻底解决你学习Java语言编程的后顾之忧-生活感悟与随笔 -亦是美网络...

    前面的文章中已经给大家介绍了关于Java编程语言的优势.特点以及就业前景等内容,大家之所以这么热衷于Java编程学习,不仅是因为Java是编程市场上最为广泛流行的语言,大家更关心的是学习了Java编程 ...

  8. angular语言前端开发_详解使用angular-cli发布i18n多国语言Angular应用_唇印_前端开发者...

    在模板 Hello world! 使用ng命令产生xlf格式的message.xlf文件 $ ng xi18n --output-path src/i18n 命令执行后,生成 src/i18n/mes ...

  9. 单片机c语言篮球比分_基于单片机的篮球计时计分器设计(C语言编程、含proteus仿真图)...

    摘要 单片机,亦称单片微电脑或单片微型计算机.它是把中央处理器( CPU ). 随机存取存储器( RAM ).只读存储器( ROM ).输入 / 输出端口( I/0 )等主要计 算机功能部件都集成在一 ...

最新文章

  1. 我用24小时、8块GPU、400美元在云上完成训练BERT!
  2. 真菌其实是长歪了的动物
  3. redis在容器里连接不上_Redis服务器被劫持风波,服务器相关知识共享学习
  4. 2.34模型--简单字符串查找(占位).c
  5. 网络获取json数据并解析
  6. Windows下visual studio code搭建golang开发环境
  7. 重构 pdf_三维温度场的重构方法,更准确地监测储能系统的电池堆内部温度
  8. 没有mmu可以跑Linux吗?
  9. (转)C#开发微信门户及应用(5)--用户分组信息管理
  10. android下拉菜单刷新,android下拉刷新列表实现(pull to refresh)
  11. 博客访问者来自15个国家和地区
  12. python sorted函数
  13. BGP Flow Specification原理
  14. IplImage 应用解读
  15. 局域网助手_IP地址管理和局域网远程助手
  16. 为什么装了个visual studio就多了个dvd驱动器?
  17. 桌面上计算机图标移动变成复制,电脑桌面图标都变成lnk后缀怎么办
  18. 一个APP想法,如何做到从想法到上线的全过程操作
  19. 如何看待人工智能的现状
  20. 英飞凌基础学习笔记PSI5(Peripheral Sensor Interface)

热门文章

  1. 蒙特卡洛能解决啥_用蒙特卡洛方法解决“无法解决”的问题
  2. amazon rds 性能_Amazon S3 —云文件存储可提高性能并节省成本
  3. JDBC原理及应用模板
  4. MVC教程第五篇:MVC整合Ajax
  5. C++ C# 中作用域限定符
  6. 漫步最优化十七——点对点映射
  7. 漫步数学分析十七——连续映射上的运算
  8. bat tomcat程序在后台执行_jenkins 中利用 cmd 命令启动 tomcat,进程被杀
  9. leetcode - 56. 合并区间
  10. 西瓜书——极大似然估计和朴素贝叶斯