下载KEGG某一条通路的基因list

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=",")
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成，浏览时请结合常识与多方信息审慎甄别。
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。
禁止转载，如需转载请通过简信或评论联系作者。

相关阅读更多精彩内容

友情链接更多精彩内容