[BioC] tilingArray normalizeByReference Data Scaling
Wolfgang Huber
whuber at embl.de
Tue May 15 09:53:44 CEST 2012
May/15/12 4:00 AM, Dario Strbenac scripsit::
> Thanks for the clarification. Yes, I was confused which one of the options in
> Figure 5 of the article were implemented in the tilingArray packages.
>
> What I meant by the difference in log bases is :
>
> refSig = rowMeans(log2(exprs(reference)[pm, , drop = FALSE]))
> whereas
> ybg[, j] = tapply(log(exprs(x)[background, j], 2), strata, genefilter::shorth,
> tie.action = "min")
>
> ?log says the default base is e. So is there a mathematical reason to use base
> 2 for reference probes, and base e for background probes ?
Hi Dario
I see. The second argument of 'log' is 'base', which in the ybg<- line
above is set to 2. In fact base=2 is used in all calls to the 'log'
function in the package.
Best wishes
Wolfgang
Wolfgang Huber
EMBL
http://www.embl.de/research/units/genome_biology/huber
More information about the Bioconductor
mailing list