搜索
查看: 95|回复: 0

[mRNA-seq] [新手]关于学习y叔的clusterprofiler的gsego分析,有一些疑问

[复制链接]

2

主题

4

帖子

57

积分

注册会员

Rank: 2

积分
57
发表于 2017-10-12 14:10:03 | 显示全部楼层 |阅读模式
学习y叔的clusterprofiler,心说示例给出的CC注释,我就试试MF注释把,结果:
> gsecc <- gseGO(geneList=geneList, ont="MF", OrgDb=org.Hs.eg.db, verbose=F)
#no term enriched under specific pvalueCutoff...
试了很多次都是这样的结果
接下来试试BP:
> gsecc <- gseGO(geneList=geneList, ont="BP",OrgDb=org.Hs.eg.db, verbose=F)
是可以的,对GO:0007009做图(附件)

用网上公开的一批拟南芥数据(GSE44062)进行测试,GO和KEGG富集分析都是可以进行的,到了gseGO分析也是这样的结果:
> gsemf<-gseGO(geneList=genelist,ont="BP",keyType="ENTREZID",OrgDb=org.hs,verbose=F)
#no term enriched under specific pvalueCutoff...
> gsemf<-gseGO(geneList=genelist,ont="CC",keyType="ENTREZID",OrgDb=org.hs,verbose=F)
#no term enriched under specific pvalueCutoff...
> gsemf<-gseGO(geneList=genelist,ont="MF",keyType="ENTREZID",OrgDb=org.hs,verbose=F)
#no term enriched under specific pvalueCutoff...
这里面的genelist和y叔给出的示例geneList格式是一样的,连genename我都特意用bitr函数改成了ENTREZID格式

这是为什么呢..?

本帖子中包含更多资源

您需要 登录 才可以下载或查看,没有帐号?立即注册

x



上一篇:宏基因组分析教程-Analysis of Metagenomic Data
下一篇:参考基因组序列中每条染色体两端的N数量是怎么确定的?
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|手机版|小黑屋|生信技能树    

GMT+8, 2017-10-20 04:02 , Processed in 0.105702 second(s), 27 queries .

Powered by Discuz! X3.2

© 2001-2013 Comsenz Inc.