[BioC] LIMMA - Intensity normalization methods

ninoochka at free.fr ninoochka at free.fr
Thu Jun 5 16:51:38 CEST 2008


Hi,

i want to work on normalized intensity values in my experiment (oligo arrays).
To do this, I first need to normalize data with a "within array normalization",
followed by a "between array normization".
So, I used printtiploess function followed by the Aquantile function. Just for
testing purpose,  I also tried with with "quantile" method and I was surprised
to see that toptables were really different: resulting p-values were better with
this second method. And, if I take the first 100 genes in the toptables, 50% are
different. As quantile function normalizes again ratios, does it mean that
printiploess method did not correctly normalized my ratios ?

Here are my first genes:

PrinTipLoess + Quantile:


        "logFC" "AveExpr" "P.Value"
"8802"         2,44         10,86         4,20E-17
"9786"         1,96         10,77         7,51E-17
"8906"         2,24         10,75         7,21E-16
"11808" 2,14         10,74         8,86E-16
"5276"         1,93         10,71         1,55E-15

PrinTipLoess + AQuantile:


        "logFC" "AveExpr" "P.Value"
"8802"         2,73         10,84         1,69E-15
"7678"         1,51         7,96         4,45E-15
"7657"         1,48         8,59         6,59E-15
"9254"         1,71         8,66         8,18E-15
"9786"         2,41         10,75         2,04E-14


Sincerely,

Lucie



More information about the Bioconductor mailing list