[BioC] getSYMBOL( ) error
Rao,Xiayu
XRao at mdanderson.org
Tue Jul 9 18:31:39 CEST 2013
Oh, I see. Thank you very much for your kind help, James!
-----Original Message-----
From: James W. MacDonald [mailto:jmacdon at uw.edu]
Sent: Tuesday, July 09, 2013 11:19 AM
To: Rao,Xiayu
Cc: 'bioconductor at r-project.org'
Subject: Re: [BioC] getSYMBOL( ) error
Hi Xaiyu,
On 7/8/2013 6:13 PM, Rao,Xiayu wrote:
> Hello,
>
> I tried to find gene names that match affyIDs, and I tried two ways as below, but neither worked. I got all NAs for the Symbol column. I have no idea what happened. Thank you very much for helping me out!
>
> (1) library(affy)
> mydata<- ReadAffy()
> eset<- rma(mydata)
> library(annotate)
> library(hgu95av2.db)
> ID<- featureNames(eset)
> Symbol<- getSYMBOL(ID,"hgu95av2.db")
> fData(eset)<- data.frame(ID=ID,Symbol=Symbol)
>
> All NAs for Symbol in the topTable
>
> (2) fit$genes$Symbol<- getSYMBOL(fit$genes$ID, "hgu95av2.db")
>
>> fit
> $genes
> ID Symbol
> 1 7892501 NA
> 2 7892502 NA
> 3 7892503 NA
> 4 7892504 NA
> 5 7892505 NA
> 32316 more rows ...
Those aren't hgu95av2 probe IDs. They appear to be hugene 1.0 transcript cluster IDs. So you need to use the hugene10sttranscriptcluster.db package.
Best,
Jim
>
> Thanks,
> Xiayu
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> 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