ego <- enrichGO(gene = eg$ENTREZID,
OrgDb = org.Hs.eg.db,
ont = "BP",
pAdjustMethod = "BH",
pvalueCutoff = 0.01,
qvalueCutoff = 0.05,
readable = TRUE)
ego2=simplify(ego)
ego3=data.frame(ego2)
ego2@result$Description=paste0(toupper(substring(ego2@result$Description, 1, 1)),substring(ego2@result$Description, 2))
p1 <- dotplot(ego2, showCategory=20,color = "pvalue") + ggtitle("dotplot for GO")
p2 <- dotplot(kegg, showCategory=20,color = "pvalue") + ggtitle("dotplot for KEGG")
pdf("combinde.pdf",width = 15,height = 9)
plot_grid(p2, p1, ncol=2, labels = c("A","B"))
dev.off()

image.png
网友评论