生成的基因列表直接用clusterProfiler和org.Hs.eg.db富集分析

library(clusterProfiler)
library(org.Hs.eg.db)
gene <- bitr(final1$gene, fromType="SYMBOL", toType=c("ENTREZID","ENSEMBL"), OrgDb="org.Hs.eg.db")
go <- enrichGO(gene$ENTREZID, OrgDb = org.Hs.eg.db, ont='ALL',pAdjustMethod = 'BH',pvalueCutoff = 0.05, 
               qvalueCutoff = 0.2,keyType = 'ENTREZID')
resultgo <- go@result
resultgo <- arrange(resultgo,resultgo$p.adjust)
kegg <- enrichKEGG(gene$ENTREZID, organism = "hsa", pvalueCutoff = 0.05)
resultkegg <- kegg@result
resultkegg_methylation1 <- arrange(resultkegg,resultkegg$p.adjust)
write.csv(resultkegg_methylation1,file="resultkegg_methylation1.csv")
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容