CpGFilter (version 1.0)

CpGFilterICC: CpG filtering method based on intra-class correlation coefficients.

Description

Filter CpGs based on Intra-class Correlation Coefficients (ICCs). ICCs are calculated by fitting linear mixed effects models to all samples including the un-replicated samples. Including the large number of un-replicated samples improves ICC estimates dramatically. The method accommodates any replicate design.

Usage

CpGFilterICC(dat, rep.design, logit.transform = TRUE, verbose = TRUE)

Arguments

dat

a matrix of CpG beta-values, row - CpG, column - sample

rep.design

a vector indicating the replicate design, it could be factor, character or numeric vectors. Example - c(1, 2, 3, 4, 4, 4, 5, 5) OR c('S1', 'S2', 'S2', 'S2', 'S1')

logit.transform

If TRUE, beta-value will be converted into M-value; Default is TRUE.

verbose

If TRUE, print run information

Value

ICCs for all probes

References

CpGFilter:Model-based CpG probe filtering with replicates for epigenome-wide association studies (2014). Submitted.

Examples

# NOT RUN {
require(CpGFilter)
rep.design <- c(1:900, 901:910, 901:910, 911:915, 911:915, 911:915, 911:915)
rho <- CpGFilterICC(matrix(rnorm(940*10000), 10000, 940), rep.design, logit=FALSE)
# }