README.md

Gaussian Mixture Quantile Normalization (GMQN)

GMQN is a reference based method that removes unwanted technical variations at signal intensity level between samples for 450K and 850K DNA methylation array. It can also easily combined with Subset-quantile Within Array Normalization(SWAN) or Beta-Mixture Quantile (BMIQ) Normalisation to remove probe design bias.

Installation

You can install GMQN R package by following steps. We will submit the R package to Bioconductor soon.

devtools::install_github('MengweiLi-project/gmqn')
library(gmqn)

Dependencies

GMQN depends on the following packages, all available in CRAN.

Demos

when you have raw(.idat) data

# I recommend using minfi to read the raw data and do preprocess.
# But other packages can also be used. 
library(minfi)
RGSet = read.metharray.exp("idat/")
MSet <- preprocessRaw(RGSet) # Other preprocess methods can also be used.
m = data.frame(getMeth(MSet))
um = data.frame(getUnmeth(MSet))

library(gmqn)
# You can skip this line if you want to use default reference.
ref = set_reference(m, um)
beta.GMQN.swan = gmqn_swan_parallel(m, um, ncpu = 45, ref = ref)
beta.GMQN.bmiq = gmqn_bmiq_parallel(m, um, ncpu = 45, ref = ref)



MengweiLi-project/gmqn documentation built on Feb. 26, 2020, 11:35 a.m.