library("KEGGREST")
listDatabases() #查看KEGG数据库的19个子数据库
gs<-keggGet('mmu00010') #Retrieves all entries from the KEGG database for a set of KEGG identifers
gs[[1]]$GENE #查找所有基因
genes<-unlist(lapply(gs[[1]]$GENE,function(x) strsplit(x,';')))
genelist <- genes[1:length(genes)%%3 ==2]
genelist <- data.frame(genelist)
View(genelist)
#把结果写入表格
write.table(genelist, "m mu00010.csv",
row.names=FALSE,col.names=TRUE,sep=",")
下载KEGG某一条通路的基因list
最后编辑于 :
©著作权归作者所有,转载或内容合作请联系作者
禁止转载,如需转载请通过简信或评论联系作者。
- 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
- 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
- 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
推荐阅读更多精彩内容
- 实例:如何拿到KEGG数据库中多巴胺通路相关的基因集 一、确定目标通路 打开KEGG选择pathway,在搜索框前...
- 有时候我们想看某一特定通路的富集情况,这时就需要获取该通路的基因集。 一、从 GO 中获取 下载脚本 getGoT...
- 差异分析完以后,就有了基因列表和差异倍数,有了这两个东西,就可以用clusterprofiler做GSEA,然后有...
- 1.kegg数据库下载 通常我们如果已知某个基因,想对这个基因做通路很简单,直接用kegg的网页注释就可以;如果想...
- 我们经常要用到KEGG数据库来对基因做功能分析。经常有长得好看的朋友问:如何获得整个通路的基因信息? 其实我们有多...