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")
生成的基因列表直接用clusterProfiler和org.Hs.eg.db富集分析
©著作权归作者所有,转载或内容合作请联系作者
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 仅针对人类而言。虽然一直对GO和KEGG不感冒,但流程化的分析还是要做的。主要包括两部分: 将基因编号转为ENTR...
- https://www.jianshu.com/p/0089b1b3cbc4?tdsourcetag=s_pcti...
- 错误: package or namespace load failed for ‘org.Hs.eg.db’:l...
- 1、出现的问题: err:package or namespace load failed for ‘org.Hs...
- 为什么pathway富集分析结果没有我感兴趣的通路? 图1 GSEA原理(图片来自plob.org) GSEA分析...