mlda {MLDA} | R Documentation |
A function to transform the log ratios of signal intensities into log likelihood ratios using three linear models on DMH array
mlda(dataset, trim = 0.05, norm = TRUE)
dataset |
data object from mlda.data() |
trim |
the trim degree when estimate the variance of the residuals to the linear model |
norm |
normalisation data using unmethylated model as reference by default. The intercept of the unmethylation linear model is adjusted around 0 |
Three linear models were constructed in this function. One is the unmethylation model constructed using mitochondrial sequences or other unmethylated controls. One is the intermediate model used to select the loci that can be used as the methylated reference. The other one is the methylation model constructed using the loci identified by the intermediate model. We estimated the distance of each locus to the methylated and unmethylated lines using standardised residuals. The likelihood of the locus being methylated was estimated by the difference of the squared standardised residual from the methylated line to that from the unmethylated line. Consequently, we had a measure of methylation status for each locus based on the log likelihood ratio of methylation rather than log ratio of cy3 (undigested) and cy5(digested) intensities.
a list containing:
ratio1 |
A matrix of log likelihood ratios |
ratio2 |
A matrix of log likelihood ratios on dye-swapped/duplicated array |
Wei Dai w.dai@imperial.ac.uk
data(hcgi12k) mlda.result<-mlda(hcgi12k,norm=FALSE,trim=0.05)