[BioC] Different Normalizations

Rafael A. Irizarry ririzarr at jhsph.edu
Sat Feb 7 18:40:48 MET 2004


On Sat, 7 Feb 2004, Anja von Heydebreck wrote:

> Hi Julia,
> 
> >Hi,
> >I am comparing different normalizations of Affymetrix microarrays:
> >1. expresso with rma, vsn, pmonly and medianpolish
> >2. calling rma and vsn with: eset=rma(Data); Set=vsn(eset);
> >
> rma gives you expression levels on the log scale, whereas
> vsn expects input values on the raw scale. So it doesn't make sense
> to call vsn after rma like in 2. (apart from the question why you
> should normalize your data twice - I guess in 1. you used rma only for
> background correction; but calling rma like in 2. involves also
> quantile normalization and medianpolish).

of cource you could do:

exprs(eset) <- 2^exprs(eset)
eset = vsn(eset)
exprs(eset) <- log2(exprs(eset))

this should make things look more alike.

> 
> >I know that medianpolish does log2 transformation and vsn does log 
> >transformation (base e), but when I compare the results, the difference 
> >between the two datasets is huge- not a matter of log-base.
> >Shouldn't the results be very similar?
> >
> >Any help would be much appreciated.
> >
> >Thanks,
> >Julia
> >  
> >
> Best,
> Anja
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://www.stat.math.ethz.ch/mailman/listinfo/bioconductor
>



More information about the Bioconductor mailing list