clusterprofiler进行富集分析

library(readr)
library(topGO)
library(KEGGREST)
library(clusterProfiler)
library(org.Hs.eg.db)
library(tidyverse)


test1 = bitr(gene, fromType="SYMBOL", toType=c("ENSEMBL", "ENTREZID"), OrgDb="org.Hs.eg.db")


ego_ALL <- enrichGO(gene = test1$ENTREZID, 
                    OrgDb = org.Hs.eg.db, 
                    ont = "ALL", 
                    pAdjustMethod = "BH", 
                    pvalueCutoff = 1,
                    qvalueCutoff = 1,
                    readable = TRUE)

barplot(ego_ALL, showCategory=20,title="EnrichmentGO_ALL")


## kegg
kegg_all <- enrichKEGG(gene = test1$ENTREZID,
                       organism = "hsa",
                       pvalueCutoff = 1,
                       qvalueCutoff = 1)