GO富集分析柱状图

 1 target_gene_id <- unique(read.delim("miRNA-gene interactions.txt")$EntrezID)
 2 # BiocInstaller::biocLite("clusterProfiler")
 3 # BiocInstaller::biocLite("org.Hs.eg.db")
 4 
 5 display_number = c(15, 10, 15)
 6 ## GO enrichment with clusterProfiler
 7 library(clusterProfiler)
 8 ego_MF <- enrichGO(OrgDb="org.Hs.eg.db",
 9              gene = target_gene_id,
10              pvalueCutoff = 0.05,
11              ont = "MF",
12              readable=TRUE)
13 ego_result_MF <- as.data.frame(ego_MF)[1:display_number[1], ]
14 # ego_result_MF <- ego_result_MF[order(ego_result_MF$Count),]
15 
16 ego_CC <- enrichGO(OrgDb="org.Hs.eg.db