R/bgeqq2.R

In wateRmelon: Illumina 450 methylation array normalization and metrics

daten2 <-
function(mn, un, onetwo, fudge=100, ret2=FALSE, ...){

   mnc <- dfsfit(mn, onetwo, ...)
   unc <- dfsfit(un, onetwo, ...)
   normalizeQuantiles(mnc) -> mncn
   normalizeQuantiles(unc) -> uncn
   beta <- mncn/(mncn + uncn + fudge)
   if (ret2) return (list(methylated=mncn, unmethylated=uncn, beta=beta))
   beta

}