cpg.qc {CpGassoc}R Documentation

Performs quality control on Illumina data.

Description

cpg.qc is designed to perform quality control on Illumina data prior to analysis. In addition to the matrix of beta values, this function requires as input matrices of Signal A, Signal B, and detection p-values. It will remove samples that have low intensity (mean signal intensity less than half of the overall median or 2000). It can also set to NA datapoints with detection p-values exceeding a user-specified cutoff, and can remove samples or sites that have a missing rate above a user-specified value. Finally, users can opt to compute beta values as M/(U+M) or M/(U+M+100).

Usage

cpg.qc(beta.orig,siga,sigb,pval,p.cutoff=.001,cpg.miss=NULL,sample.miss=NULL,
constant100=FALSE,sig.return=FALSE)

Arguments

beta.orig

The original beta values matrix obtained from GenomeStudio.

siga

The unmethylated signals matrix obtained from GenomeStudio.

sigb

The methylated signals matrix obtained from GenomeStudio.

pval

A matrix of detection p-values obtained from GenomeStudio. pval should have the same dimension as the beta values and signals: one row for each site and one column for each individual.

p.cutoff

The user-specified cutoff for detection p-values (default=.001).

cpg.miss

Optional cutoff value. If specified, cpg.qc will remove cpg sites where the proportion of missing values exceeds this cutoff.

sample.miss

Optional cutoff value. If specified, cpg.qc will remove samples where the proportion of missing values exceeds this cutoff.

constant100

Logical. If true, the new beta values will be calculated as M/(U+M+100); if false (default) they will be calculated as M/(U+M).

sig.return

Logical. If true, cpg.qc returns the qced signal data as well.

Details

It is important that all the matrices or data frames listed above (pval, siga, sigb, beta.orig) are ordered similarly with respect to samples and CpG sites.

Value

cpg.qc returns a new matrix of beta values that has been subjected to the specified quality control filters. This matrix can be input directly into cpg.assoc.

Author(s)

Barfield, R.; Conneely, K.; Kilaru,V.
Maintainer: R. Barfield: <rbarfield01@fas.harvard.edu>

See Also

cpg.work cpg.perm plot.cpg scatterplot cpg.combine cpg.assoc


[Package CpGassoc version 2.60 Index]