1. library(tidyverse)
    2. library(clusterProfiler)
    3. library(org.Hs.eg.db)
    4. setwd("")
    6. gene <- data.frame(gene = c(
    7.   "COA3", 
    8.   "CCDC107",
    9.   "HSPB8",
    10.   "KIAA1671",
    11.   "SCD",
    12.   "KCTD10"
    13. ))
    14. #Kegg
    15. #基因名转换
    16. eg <- bitr(geneID = gene$gene,fromType = "SYMBOL",toType=c("ENTREZID","ENSEMBL"), OrgDb="org.Hs.eg.db")
    18. #kegg分析
    19. kegg <- enrichKEGG(gene = eg$ENTREZID, 
    20.                    organism = "hsa", 
    21.                    keyType = "kegg", 
    22.                    pvalueCutoff = 0.05,
    23.                    pAdjustMethod = "BH",
    24.                    minGSSize = 10, 
    25.                    maxGSSize = 500,
    26.                    qvalueCutoff = 0.2,
    27.                    use_internal_data = FALSE)
    30. dotplot(kegg)
    31. go <- enrichGO(eg$ENTREZID, OrgDb = org.Hs.eg.db, ont='ALL',pAdjustMethod = 'BH',pvalueCutoff = 0.05, 
    32.                qvalueCutoff = 0.2,keyType = 'ENTREZID')
    33. dotplot(go,showCategory = 30)

    image.png