[BioC] Converting Affymetrix (*_at) Mouse Probe ID to HUGO gene symbol
James W. MacDonald
jmacdon at uw.edu
Mon Sep 9 17:17:35 CEST 2013
You can do either
library(mgu74av2.db)
thedata <- select(mgu74av2.db, probes, "SYMBOL")
or
g <- getBM( c("affy_mg_u74av2", "mgi_symbol"), "affy_mg_u74av2",
probes, mouse)
But you should note that the two probes you list below are not from the
chip you are using in your query. The probes you supply are from some
version of the Mouse Genome 430 series. If I assume the 430 version 2:
getBM( c("affy_mouse430_2", "mgi_symbol"), "affy_mouse430_2",
probes,mart)
affy_mouse430_2 mgi_symbol
1 1460645_at Chordc1
2 1460644_at Bckdk
> select(mouse4302.db, probes, "SYMBOL")
PROBEID SYMBOL
1 1460644_at Bckdk
2 1460645_at Chordc1
Best,
Jim
On Monday, September 09, 2013 10:31:35 AM, Gundala Viswanath wrote:
> I have the following list of probes (from Mouse),
>
> 1460644_at
> 1460645_at
>
> The longer list can be found here: http://dpaste.com/1371949/plain/
> What I want to do is to convert that name with HUGO gene symbol using
> R. What's the way to do it?
>
> I tried this but failed
>
> __BEGIN__
>
> 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)
>
> __END__
>
> It prints Ensembl gene id instead (e.g.ENSMUSG00000057666)
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at r-project.org
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor
--
James W. MacDonald, M.S.
Biostatistician
University of Washington
Environmental and Occupational Health Sciences
4225 Roosevelt Way NE, # 100
Seattle WA 98105-6099
More information about the Bioconductor
mailing list