R: clusterProfiler/enrichplot 富集分析与可视化神器

5 篇文章 2 订阅

使用R包clusterProfiler 进行富集分析与可视化

clusterProfiler官方网址:

http://www.bioconductor.org/packages/release/bioc/html/clusterProfiler.html

ID 转换功能

函数:bitr;bitr_kegg
功能:实现GeneSymbol, GeneID, Ensembl,UniprotID, ncbi-proteinid等之间的相互准换。支持约20物种。

在这里插入图片描述

图片说明:上图中x是一个list,里面放的是基因名,也就是SYMBOL这一列,也可以传入其他,只要在fromType一列写上对应的类别即可。

GO/KEGG 富集分析

一条命令即可完成分析
1. GO分析

在这里插入图片描述

2. KEGG分析

在这里插入图片描述

3. 输出结果

在这里插入图片描述

图片说明:第一张图为GO分析,第二张图为KEGG分析,传入的gene 是一个基因列表,该列表可以是geneID GeneSymbol 甚至是ENSEMBL, 只要在keyType 一列说明具体是哪一种类型就可以了。第三张图为输出结果,p.adjust和qvalue都是对p值的校正,p.adjust 在分析时可以选用方法,比如图中选择的是BH, qvalue是FDR校正。

GO/KEGG 富集分析可视化

直接使用上述分析结果ego 和kk,一条命令即可完成可视化

柱状图

barplot(ego, showCategory=20)

在这里插入图片描述

气泡图

dotplot(kk)

在这里插入图片描述

通路-基因网络图

cnetplot(ego, categorySize=“pvalue”, foldChange=geneList)

在这里插入图片描述

cnetplot(ego, foldChange=geneList, circular = TRUE, colorEdge = TRUE)

在这里插入图片描述

热图

heatplot(ego, foldChange=geneList)

在这里插入图片描述

pathway-pathway Network

emapplot(ego)

在这里插入图片描述

GO通路网络图

goplot(ego)

在这里插入图片描述

KEGG通路注释

图片说明:在这里,geneList 是一个耗油基因名及对应的flodchange 的列表,图中绿色即为下调,红色为上调。
其他可视化如GSEA分析等,可参照官方说明。
  • 16
    点赞
  • 121
    收藏
    觉得还不错? 一键收藏
  • 7
    评论

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

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 7
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值