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")

你可能感兴趣的:(2019-08-26 biomaRt 进行基因转换)