[BioC] RMA + loess normalisation and filtering
naomi at stat.psu.edu
Tue Apr 19 17:51:43 CEST 2005
This topic keeps coming up. Any references supported by data that show
some method is preferable?
Wolfgang, how many replicates do you think are needed for your method to be
(I've gotten burned by reviewers on this one.)
At 12:01 PM 4/19/2005, Wolfgang Huber wrote:
>>>question 1: I have performed *RMA normalisation *of my Affymetrix data.
>>>However, for further analysis I think it is necessary to *filter* the
>>>data (non-expressed genes or below background). However I don't know the
>>>best way to filter the genes that are not expressed or very low
>>>expressed (below the background), based on the RMA normalisation data.
>>My preference is to select genes based on their overall variability,
>>using a criterions such as
>> z = apply(exprs(x), 1, IQR)
>Just remembered - a discussion (with data) on whether to select on mean
>level or variability can be found here:
>There are also additional ideas on pre-filtering, in order to alleviate
>the loss of power from multiple testing.
>European Bioinformatics Institute
>European Molecular Biology Laboratory
>Cambridge CB10 1SD
>Phone: +44 1223 494642
>Fax: +44 1223 494486
>Bioconductor mailing list
>Bioconductor at stat.math.ethz.ch
Naomi S. Altman 814-865-3791 (voice)
Bioinformatics Consulting Center
Dept. of Statistics 814-863-7114 (fax)
Penn State University 814-865-1348 (Statistics)
University Park, PA 16802-2111
More information about the Bioconductor