go kegg_KEGG分析及可视化

上一篇推文中我们解释了GO富集分析及可视化（GO富集分析及可视化），除了GO富集分析，我们经常在paper中看到KEGG分析，KEGG是什么呢，Kyoto Encyclopedia of Genes and Genomes，京都基因和基因组百科全书，是系统分析基因功能，联系基因组信息和功能信息的知识库，其中包含有大量的通路图。

我们今天来介绍如何绘制KEGG的barplot和dotplot。

library(DOSE)
library(org.Hs.eg.db)
library(clusterProfiler)
geneNames <- names(geneList)[1:100] ##数据来自DOSE这个包
enrich <- enrichKEGG(gene = geneNames,
                       organism = "hsa",
                       pvalueCutoff = .1,
                       qvalueCutoff = .1)

我们先来看看结果

> head(enrich)
               ID                                                   Description GeneRatio  BgRatio       pvalue     p.adjust       qvalue       geneID Count
hsa04110 hsa04110                                                    Cell cycle 8/48 124/8081 5.