[BioC] Affy: probeset to gene expression with expresso

Martin Preusse martin.preusse at gmail.com
Tue Aug 13 17:46:17 CEST 2013


Thanks! This clarifies everything ;) I misunderstood the 'summary' part of expresso. 



Am Dienstag, 13. August 2013 um 17:41 schrieb James W. MacDonald:

> Hi Martin,
> 
> None of the functions in the affy package know anything about the 
> annotation of the various probesets, so there is no facility there to 
> summarize at anything but the probeset level.
> 
> However, you could use one of the MBNI remapped probesets, which 
> pre-aggregate the probes into probesets based on a few different 
> annotation databases.
> 
> Best,
> 
> Jim
> 
> 
> 
> On 8/13/2013 11:38 AM, Martin Preusse wrote:
> > Hi Jim,
> > 
> > thanks, this is a very clever (and R like) way to average the expressions on the expression matrix. In my example this would work for 'exprs(eset)' when I replace the probe set ids with gene symbols.
> > 
> > Still, it would be great to know if can be achieved with a convenience function (e.g. expresso).
> > 
> > Martin
> > 
> > 
> > Am Dienstag, 13. August 2013 um 17:28 schrieb James W. MacDonald:
> > 
> > > Hi Martin,
> > > 
> > > I just answered a very closely related question. See if this helps:
> > > 
> > > https://stat.ethz.ch/pipermail/bioconductor/2013-August/054353.html
> > > 
> > > Best,
> > > 
> > > Jim
> > > 
> > > 
> > > 
> > > On 8/13/2013 9:47 AM, Martin Preusse wrote:
> > > > I am trying to get the gene level expression values from an Affy micro array, i.e. merge the values for probe sets representing the same gene.
> > > > 
> > > > I tried to use the 'expresso' function from the affy package, but I always end up with an ExpressionSet containing probe sets, not genes.
> > > > 
> > > > What is an easy way to summarize/merge probe sets to (entrez) genes?
> > > > 
> > > > 
> > > > library(affydata)
> > > > library(affy)
> > > > 
> > > > # get the 'Dilution' affy batch
> > > > data(Dilution)
> > > > 
> > > > eset<- expresso(Dilution, bgcorrect.method='rma',
> > > > normalize.method='constant',
> > > > pmcorrect.method='pmonly',
> > > > summary.method='avgdiff')
> > > > 
> > > > 
> > > > write.exprs(eset,'testfile.txt')
> > > > 
> > > > 
> > > > P.S.: I know it might not be the best idea to average probe sets, but I would like to try ;)
> > > > 
> > > > Cheers
> > > > Martin
> > > > 
> > > > _______________________________________________
> > > > Bioconductor mailing list
> > > > Bioconductor at r-project.org (mailto: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
> > 
> 
> 
> 
> -- 
> 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