controlCorrection {nucleR} | R Documentation

## Correct experimental profiles with control sample

### Description

This function allows the correction of experimental coverage profiles (usually MNase digested nucleosomal DNAs in this library) with control samples (usually naked DNA sample digested with MNase). This is useful to correct MNase biase.

### Usage

    controlCorrection(exp, ctr, ...)

    ## S4 method for signature 'SimpleRleList'
    controlCorrection(exp, ctr, mc.cores = 1)

    ## S4 method for signature 'Rle'
    controlCorrection(exp, ctr)

    ## S4 method for signature 'list'
    controlCorrection(exp, ctr, mc.cores = 1)

    ## S4 method for signature 'numeric'
    controlCorrection(exp, ctr)

### Arguments

`exp, ctr` | Comparable experimental and control samples (this means same format and equivalent preprocessment)
`...` | Further arguments to be passed to or from other methods.
`mc.cores` | Number of cores available for parallel list processing

### Details

This substracts the enrichment in the control sample respect it's mean from the experimental profile.

This is useful for examinating the effect of the MNase digestion in nucleosome experiments using a nucleosomal DNA and a genomic (naked) DNA sample. Notice that genomic DNA samples cannot be strand-corrected using single end data, so only paired end controls are useful for this proupose, despite they can be compared against extended nucleosomal DNA single end reads. Furthermore, both datasets must be converted to reads per milion.

This process dificults the nucleosome positioning due the lower sharpness of the peaks, but allows a complementary study of the MNase digestion effect.

### Value

Corrected experimental profile

### Author(s)

Oscar Flores oflores@mmb.pcb.ub.es

### Examples

    map = syntheticNucMap(as.ratio=TRUE)
    exp = coverage.rpm(map$syn.reads)
    ctr = coverage.rpm(map$ctr.reads)
    corrected = controlCorrection(exp, ctr)

[Package nucleR version 2.18.4 Index]