0

我有以下探针列表(来自鼠标),

1460644_at
1460645_at
1460646_at
1460647_a_at
1460648_at
1460649_at
1460650_at
1460651_at
1460652_at
1460653_at
1460654_at
1460655_a_at
1460656_a_at
1460657_at
1460658_at
1460659_at
1460660_x_at
1460661_at

更长的列表可以在这里找到:http: //dpaste.com/1371949/plain/

我想要做的是使用 R 将该名称与基因符号转换。有什么方法可以做到这一点?

我试过这个但失败了

#source("http://bioconductor.org/biocLite.R")
# biocLite("hgu95av2.db")
# biocLite("affy")
library(hgu95av2.db)
library(affy)
dat<-read.table("http://dpaste.com/1371949/plain/")
probes<-as.vector(as.matrix(dat))
Symbols = unlist(mget(probes, hgu95av2SYMBOL, ifnotfound=NA))
# all returns NA

更新:我试过 BioMart 但无济于事:

library(biomaRt)
dat<-read.table("http://dpaste.com/1371949/plain/")
probes<-as.vector(as.matrix(dat))
mouse = useMart("ensembl", dataset = "mmusculus_gene_ensembl")
g = getGene( id = probes, type = "affy_mg_u74av2", mart = mouse)
show(g)

它改为打印 Ensembl 基因 id (egENSMUSG00000057666)

4

0 回答 0