
Hi Sebastien,
 
Thanks for the tip :
 
ls(hgu133aSYMBOL)[mget(ls(hgu133aSYMBOL), hgu133aSYMBOL) %in% gensym]
works fine
 
I would be grateful if you can tell me how to convert Affy probe to gene name? i.e. the opposite way
 
I tried using hgu133aGENENAME but that didn't work as its using the description. 
 
How do I use hgu133aSYMBOL to get gene name from affy probe.
 
Many thanks
 
Ruppert
 
> Date: Tue, 6 May 2008 13:07:22 +1000> From: seb@gerega.net> To: ruppert7@hotmail.com> Subject: Re: [BioC] converting gene name to affy probe> > Ruppert Valentino wrote:> > Hello, I am trying to convert a file of gene names to corresponding affy probe names. I managed to write a script that puts the genes in an array then I use the feat = getFeature(symbol = gensym, type = "affy_hg_u133a", mart = mart) in biomaRt however I seem to hit a snag when there is more than probe for a gene name. Does anyone know of an existing script that can do this? thanks Ruppert> > _________________________________________________________________> > Win Indiana Jones prizes with Live Search> >> > [[alternative HTML version deleted]]> >> > _______________________________________________> > Bioconductor mailing list> > Bioconductor@stat.math.ethz.ch> > https://stat.ethz.ch/mailman/listinfo/bioconductor> > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor> >> > > > I think this should do what you want:> library(hgu133a.db)> ls(hgu133aSYMBOL)[mget(ls(hgu133aSYMBOL), hgu133aSYMBOL) %in% gensym]> > hope that helps,> Sebastien
_________________________________________________________________
Be a Hero and Win with Iron Man

	[[alternative HTML version deleted]]

