aPEAR包绘制功能富集网络图

本期教程

前言

今天学习aPEAR包,绘制KEGG和GO功能富集网络图,用起来还是比较方便的,直接将clusterProfiler富集结果进行绘制,对人类、动物等分析结果非常方便。对于模式植物,使用自己制作的GO或KEGG背景文件进行富集分析,理论上也是可行的,今天尝试一下,没成功,后面继续进行尝试。

教程原文网址aPEAR包绘制功能富集网络图(点击链接)

学习网址

  1. Github (点击即可进入)

    https://github.com/kerseviciute/aPEAR

  1. https://cloud.r-project.org/(点击即可进入)

基础教程

根据的官方教程,使用aPEAR包绘制富集网络图是很方便的,尤其是人类、动物等。

  1. 安装aPEAR包

    ##'@安装aPEAR包
    library(devtools)
    install_github('ievaKer/aPEAR')

  2. 加载所需的包

    library(aPEAR)
    library(clusterProfiler)
    library(org.Hs.eg.db)
    library(DOSE)
    library(ggplot2)

  3. 加载数据

    ##'@加载测试数据
    data(geneList)
    ##'@查看数据
    head(geneList)

  4. 功能富集

    ##'@进行功能富集
    enrich <- gseGO(geneList, OrgDb = org.Hs.eg.db, ont = 'CC')
    ##'@查看富集结果

  5. 绘图

    enrichmentNetwork(enrich@result, drawEllipses = TRUE, fontSize = 2.5)

  1. 根据p值修改颜色

    enrichmentNetwork(enrich@result,
    colorBy = 'pvalue',
    colorType = 'pval',
    pCutoff = -5)

  1. 根据ggplot修改颜色

    enrichmentNetwork(enrich@result,
    colorBy = 'pvalue',
    colorType = 'pval')+
    scale_color_gradientn(colours = c("#B83D3D",'white','#1A5592'),
    name = "pvalue")

  1. Visualize pathway clusters with plotPathClusters()

    set.seed(238923)
    plotPathClusters(
    enrichment = enrich@result,
    sim = clusterssimilarity, clusters = clustersclusters,
    fontSize = 4,
    outerCutoff = 0.01, # Decrease cutoff between clusters and show some connections
    drawEllipses = TRUE
    )

参考:

  1. https://github.com/kerseviciute/aPEAR
  2. https://cloud.r-project.org/web/packages/aPEAR/vignettes/aPEAR-vignette.html
  3. https://mp.weixin.qq.com/s/GVOyfghR0ELLRZLZmTrzew
  4. https://mp.weixin.qq.com/s/wnNhYP5QJ7afs_X28j7nnA

往期文章:

1. 复现SCI文章系列专栏

2. 《生信知识库订阅须知》,同步更新,易于搜索与管理。

3. 最全WGCNA教程(替换数据即可出全部结果与图形)


4. 精美图形绘制教程

5. 转录组分析教程

转录组上游分析教程[零基础]

小杜的生信筆記 ,主要发表或收录生物信息学的教程,以及基于R的分析和可视化(包括数据分析,图形绘制等);分享感兴趣的文献和学习资料!!

相关推荐
一个响当当的名号1 天前
lectrue9 索引并发控制
java·开发语言·数据库
进阶小白猿1 天前
Java技术八股学习Day30
java·开发语言·学习
liu****1 天前
4.Qt窗口开发全解析:菜单栏、工具栏、状态栏及对话框实战
数据库·c++·qt·系统架构
近津薪荼1 天前
优选算法——双指针6(单调性)
c++·学习·算法
三水不滴1 天前
Redis缓存更新策略
数据库·经验分享·redis·笔记·后端·缓存
修修修也1 天前
【无标题】技术欲望是怎样渐进增长的?
学习
企鹅侠客1 天前
第35章—内核解析篇:Redis内存淘汰机制
数据库·redis
西柚小萌新1 天前
【人工智能:Agent】--12.多智能体(二)
数据库
whale fall1 天前
celery -A tool.src.main worker --loglevel=info --queues=worker1_queue & 什么意思
python·学习·apache
快乐非自愿1 天前
【面试题】MySQL 的索引类型有哪些?
数据库·mysql·面试