2019-08-14 biomaRt 转换

2019-08-20  本文已影响0人  njmujjc

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)

上一篇 下一篇

猜你喜欢

热点阅读