Kevin2015的个人博客分享 http://blog.sciencenet.cn/u/Kevin2015

博文

clusterProfiler 举例

已有 3702 次阅读 2016-7-7 17:39 |个人分类:知识点专题|系统分类:科研笔记

http://www.bioconductor.org/packages/release/bioc/vignettes/clusterProfiler/inst/doc/clusterProfiler.html  
brca <- scan("brca.txt",what="")
kirc<- scan("kirc.txt",what="")
prad <- scan("prad.txt",what="")
thca <- scan("thca.txt",what="")
ucec <- scan("ucec.txt",what="")

library("org.Hs.eg.db")
library(clusterProfiler)
db <- org.Hs.eg.db
entriz_brca <- select(db,brca,"ENTREZID","SYMBOL")
entriz_kirc <- select(db,kirc,"ENTREZID","SYMBOL")
entriz_prad <- select(db,prad,"ENTREZID","SYMBOL")
entriz_thca <- select(db,thca,"ENTREZID","SYMBOL")
entriz_ucec <- select(db,ucec,"ENTREZID","SYMBOL")

enrichGO_brca <- enrichGO(gene = na.omit(entriz_brca[,2]),  organism = "human", ont = "BP", pvalueCutoff = 0.05, readable= TRUE)
enrichGO_kirc <- enrichGO(gene = na.omit(entriz_kirc[,2]), organism = "human", ont = "BP", pvalueCutoff = 0.05, readable= TRUE)
enrichGO_prad <- enrichGO(gene = na.omit(entriz_prad[,2]), organism = "human", ont = "BP", pvalueCutoff = 0.05, readable= TRUE)
enrichGO_thca <- enrichGO(gene = na.omit(entriz_thca[,2]), organism = "human", ont = "BP", pvalueCutoff = 0.05, readable= TRUE)
enrichGO_ucec <- enrichGO(gene = na.omit(entriz_ucec[,2]), organism = "human", ont = "BP", pvalueCutoff = 0.05, readable= TRUE)
 
summary(enrichGO_brca)
 summary(enrichGO_kirc)
 summary(enrichGO_prad)
 summary(enrichGO_thca)
 summary(enrichGO_ucec)




mylist <- list(summary(enrichGO_brca)$ID,summary(enrichGO_kirc)$ID,summary(enrichGO_prad)$ID,summary(enrichGO_thca)$ID,summary(enrichGO_ucec)$ID)
intersect_GO <-  Reduce(intersect,mylist)



https://wap.sciencenet.cn/blog-2609994-989248.html

上一篇:Ven图 5 set
下一篇:正则
收藏 IP: 14.204.63.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-20 20:04

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部