[BioC] summarizing probe intensites before or after normalization- 1. how to do with RMA 2. Opinions?

k. brand k.brand at erasmusmc.nl
Mon Sep 11 16:52:20 CEST 2006


Ben,

My apologies, i see my error now.

library(affyPLM)
dat <- ReadAffy()
datrma <- rma(dat, normalize=FALSE)
datrma.postqnorm <- normalize(datrma)
boxplot(datrma.postqnorm)

The above worked a treat!

Id still like very much to hear your opinion on whether i might be 
better to normalize before or after summarizing given the variation i 
have to deal with.

thanks again,

Karl

on 9/11/2006 3:51 PM Ben Bolstad said the following:
>> 1. Help to effect Bolstad normalization of the RMA preprocessed and 
>> summarized data. Whilst I succeed in generating unnormalized RMA 
>> preprocessed data with-
>>
>> library(affy)
>> data <- ReadAffy()
>> datarma <- rma(data, normalize=FALSE)
>>
>> As a result of my limited R experience, I failed in finding a method to 
>> effect Bolstad (quantile) normalization on this output.
> 
> library(affyPLM)
> datarma.postqnorm <- normalize(datarma)
> 
> 
> 
> 

-- 
Karl Brand <k.brand at erasmusmc.nl>
Department of Cell Biology and Genetics
Erasmus MC
Dr Molewaterplein 50
3015 GE Rotterdam
lab +31 (0)10 408 7409 fax +31 (0)10 408 9468



More information about the Bioconductor mailing list