2019-08-26 biomaRt 进行基因转换

library("biomaRt")

human = useMart("ensembl", dataset = "hsapiens_gene_ensembl")

mouse = useMart("ensembl", dataset = "mmusculus_gene_ensembl")

gene<- read.table("geneinrepeat1000kb.txt")

  gene<-gene[,1]

  aa<-getBM(attributes = c("refseq_mrna","external_gene_name"),filters = "refseq_mrna",values = gene,mart = mouse)

  write.table(aa,file="geneinreat1000kb-genename.txt")

©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。

推荐阅读更多精彩内容