"aggregate" function in R give non-unique values
1
0
Entering edit mode
Wendy Qiao ▴ 360
@wendy-qiao-4501
Last seen 10.3 years ago
Hi all, When I analyse my microarray data, I used "aggregate" to calculate the mean of the probes that have the same gene symbols. Following is the code that I use: data.lumi.aggre<-aggregate(data.lumi,by=list(data.lumi$SYMBOL),FUN=mea n) I think that "aggregate" function gives unique gene symbols; however in my results, multiple rows have the same gene symbols. I wonder if anybody encountered the same problem. Is this a bug in R? How could I solve it? Thank you, Wendy [[alternative HTML version deleted]]
Microarray Microarray • 1.3k views
ADD COMMENT
0
Entering edit mode
@martin-morgan-1513
Last seen 4 months ago
United States
On 09/07/2011 02:49 PM, Wendy Qiao wrote: > Hi all, > > When I analyse my microarray data, I used "aggregate" to calculate the mean > of the probes that have the same gene symbols. Following is the code that I > use: > > data.lumi.aggre<-aggregate(data.lumi,by=list(data.lumi$SYMBOL),FUN=m ean) Hi Wendy -- can you provide a simple reproducible example that we can copy from your email and paste into our R session to replicate the problem? Also, can you provide the output of the sessionInfo() command? Martin > > I think that "aggregate" function gives unique gene symbols; however in my > results, multiple rows have the same gene symbols. I wonder if anybody > encountered the same problem. Is this a bug in R? How could I solve it? > > Thank you, > Wendy > > [[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 -- Computational Biology Fred Hutchinson Cancer Research Center 1100 Fairview Ave. N. PO Box 19024 Seattle, WA 98109 Location: M1-B861 Telephone: 206 667-2793
ADD COMMENT
0
Entering edit mode
Hi Martin, Thank you for your reply. I realized that the problem that I am encountering is not because of the "aggregate" function, but the conversion between gene symbols from my array to official symbols: I am trying to convert the gene symbols from my array to official gene symbols using data.lumi.officialSymbol<-alias2SymbolTable(data.lumi.aggre.Symbol,sp ecies="Hs") Then I replaced the gene symbols on the array with the official symbols. I just realized that multiple gene symbols on the array were matched to the same official symbol, which I thought was an error of the aggregate function. I am sorry for confusing you. Thanks, Wendy [[alternative HTML version deleted]]
ADD REPLY

Login before adding your answer.

Traffic: 816 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6