merge_classes: MERGE ILLUMINA PROBE CLASSES

Description

Illumina assigns quite specific functional classes to the probes. 11 classes are defined (Bibikova et al. 2009 & 2011):

As DMRforPairs looks at probes in close proximity to each other within each class, very specific annotations might result in too few probes per region per class and a subsequent drop in the number of identified regions and/or statistical power. This function therefore allows grouping and / or selection of classes of interest.

Usage

1
merge_classes(refgene_class,island_class,recode=1,sep=";")

Arguments

refgene_class

see description and classes_gene parameter in DMRforPairs

island_class

see description and classes_island parameter in DMRforPairs

recode

recoding scheme to use for the functional classes. Can be custom scheme (data frame) or build in scheme (0, 1 or 2). See details.

sep

Separator used in the second column of the recode parameter. Use ";" or do not specify if using the build in schemes.

Details

The recode parameter can be set to use one of the build in recoding schemes:

The recode parameter can also be set to a custom recoding scheme (data.frame). For example:

data.frame(c("gene","tss"),c("Body;5\'UTR;3\'UTR;1stExon","TSS1500;TSS200")).

In this scheme the classes are merged into two categories: TSS or other gene region. Probes solely associated with CpG island-related classes are discarded. Probes not annotated to any of the 11 classes are always discarded by DMRforPairs, except when option 2 is used which collects all probes into one class (i.e. ignores classes).

If classes are unknown an m x 1 character vector with "unknown.gene" and "unknown.island" for all m rows (probes) can be used for the refgene_class and island_class parameters respectively (set "recode" to 2!).

Value

List of objects containing:

$pclass

original classes per probe (gene and island classes are merged) (m x 1 data frame)

$pclass_merged

classes after recoding (m x 1 data frame)

$no.pclass

row indexes of probes with no annotation to any of the classes specified in the recoding scheme.

$u_pclass

unique list of the classes of interest after recoding (i.e. the first column of the recode data fram)

Author(s)

Martin Rijlaarsdam, m.a.rijlaarsdam@gmail.com

References

Examples

1
2
3
#merge_classes() is an integrated part of the DMRforPairs() wrapper and is 
#not usually called by the user directly. Please see DMRforPairs() for 
#an example.

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.