Hello, when processing an affybatch by performing all steps "manually" (i.e. not using expresso or justRMA), when do I call pmcorret? does this order make sense? 1. bgcorrect (rma) 2. normalize (quantiles) 3. pmcorrect (pmonly) ... no summary! In this case (since I'm using rma) do I have to take care that I remove the MM intensities from the affybatch before passing the bgcorrected affybatch to the normalization (because on the PMs were corrected)? Basically I'd just like to skip the last step the "summary" to finally have a matrix with the probe insensities (and not the summarized *probeset* intensities). justMA cannot leave out the summary step, can it? kind regards + thanks for your comments, Arne -- Arne Muller, Ph.D. Toxicogenomics, Aventis Pharma arne dot muller domain=aventis com

In general you'd apply the background correction function to an affybatch, to get a bg corrected affybatch. Then send this to a normaliation routine and so on .... For the case that you are asking for no "pmcorrection" is needed . The recommended way of getting what you want is mentioned on: http://www.stat.berkeley.edu/~bolstad/ComputeRMAFAQ/ComputeRMAFAQ.html "If you are interested in applying the RMA preprocessing to just the PM probes of an affybatch Data. You might find the following code slightly more efficient than applying the standard operations bg.correct(),normalize() to an affybatch. my.PM <- apply(pm(Data), 2, bg.adjust) my.PM <- normalize.quantilesmy.PM)" On Thu, 2004-04-22 at 10:05, Arne.Muller@aventis.com wrote: > Hello, > > when processing an affybatch by performing all steps "manually" (i.e. not using expresso or justRMA), when do I call pmcorret? > > does this order make sense? > 1. bgcorrect (rma) > 2. normalize (quantiles) > 3. pmcorrect (pmonly) > ... no summary! > > In this case (since I'm using rma) do I have to take care that I remove the MM intensities from the affybatch before passing the bgcorrected affybatch to the normalization (because on the PMs were corrected)? > > Basically I'd just like to skip the last step the "summary" to finally have a matrix with the probe insensities (and not the summarized *probeset* intensities). justMA cannot leave out the summary step, can it? > > kind regards + thanks for your comments, > > Arne > > -- > Arne Muller, Ph.D. > Toxicogenomics, Aventis Pharma > arne dot muller domain=aventis com > > _______________________________________________ > Bioconductor mailing list > Bioconductor@stat.math.ethz.ch > https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor -- Ben Bolstad <bolstad@stat.berkeley.edu> http://www.stat.berkeley.edu/~bolstad