2019-08-14 biomaRt 转换

library("biomaRt")

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

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

gene<- read.table("../DEgenelist.txt")

gene<-gene[,1]

##从 gene name 转换成 refseq_mrna,  输入的是gene name, 则  filters ="external_gene_name" ##

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

head(aa)

你可能感兴趣的:(2019-08-14 biomaRt 转换)