数据库代码公开：KEGG富集分析

#引用R包
library(R.utils)
R.utils::setOption("clusterProfiler.download.method","auto")
library("clusterProfiler")
library("org.Hs.eg.db")
library("enrichplot")
library("ggplot2")
#富集的基因列表
geneList = c("PCNA PSMD8 PSMD7 SET SNRPA1 RAN SRSF2 G3BP1 STARD7 NPM1 BUB3 EIF3D XPO1 FBL EIF4A1 CANX NAP1L1 CBX3 CCT3 C1QBP U2AF1 UBE2L3 SSBP1 SRSF1 TCP1 MCM2 EIF3B PSMD14 SNRPA PWP1 APEX1 TXNL4A HNRNPR PSMB2 HPRT1 MCM6 NME1 SNRPD1 EEF1B2 HSPD1 CAD RPL18 PGK1 DDX18 RPS2 LDHA RUVBL2 RNPS1 EIF2S1 RANBP1 MCM5 IARS1 UBE2E1 MYC")
geneList = unlist(strsplit(geneList," "))
#转换基因ID
trans <- bitr(geneList, 
              fromType="SYMBOL",
              toType="ENTREZID", 
              OrgDb="org.Hs.eg.db")
trans = trans[is.na(trans[,"ENTREZID"])==F,]  
genes = trans$ENTREZID 
KEGG = enrichKEGG(gene = genes, organism = "hsa", pvalueCutoff =1, qvalueCutoff =1)
KEGG = as.data.frame(KEGG)
KEGG = KEGG[(KEGG$pvalue < 0.05 & KEGG$qvalue < 0.05),]