DBAmmd-class: Class DBAmmd

In dhelekal/mmdiff3: Statistical Testing for ChIP-Seq data sets

Description

The DBAmmd Class defines a container for differential binding analysis using MMDiff2. For this class a number of methods is foreseen, among which accessors for every slot. As MetaData, it needs to contain the path to the data directory and the name of a sampleSheet csv file.

Value

DBAmmd Object

Constructor

DBAmmd()returns an empty DBAmmd Object.
DBAmmd(MetaData) initializes a DBAmmd Object for a new Experiment.
(See below and the package vignette for more details.)

Slots

MetaData:

List containing an ExpData and an AnaData compartment. "ExpData" needs a dataDir and a SampleSheet entry. A genome entry, which should be a valid BSGenome name, is useful to find sequence motifs. (Note the genome version needs to correspond to the one used for the read alignment. Use available.genomes() to find the right name.) The AnaData entry is used to store and access parameters for the MMDiff2 Analysis, like the sigma of the RBF Kernel.

rowRanges:

GRanges object containing Regions of Interests (Peaks)

Reads:

List containing positions of mapped reads, i.e. exact start and end positions of mapped fragments. In the case of single-end reads, the left most postions of fragments mapping to the positive strands and the right most positions of fragments mapping to the negative strands are stored in "Left.p" and "Right.n". Use getPeakReads to fill this slot and estimateFragmentCenters to add the (estimated) positions of fragment centers.

RawTotalCounts:

m x n matrix containing total counts of reads mapping to m peaks in n samples (including input samples)

RawCounts.p:

m x n matrix containing counts of reads mapping to positive (forward) strand

RawCounts.n:

m x n matrix containing counts of reads mapping to negative (reverse) strand

Hists:

List of lists, each of length m (number of Peaks). Compartments could be 'Left.p','Right.n','Left.n','Right.p','Center.n', 'Center.p','Center','Left','Right', defining whether left or right ends or centers of fragments should be considered for positive ('p') or negative ('n') strand, or both strands combined. For a given compartment there is one entry per peak, which is a n x L_i matrix, where n is the number of samples and L_i is the number of bins used to cover the extend of the peak. Note, L_i varies between peaks of different lengths. See compHists() for more details.

DISTs:

List with compartments for different methods to compute distances (e.g. MMD). Each compartment contains a m x N matrix with computed distances for each Peak between N pairs of samples. See compDists() for more details.

mCounts:

(for internal use only)

Contrasts:

List of lists. Each entry contains a contrast i.e. the definition of two groups that should be compared to each other in a differential analysis. A Contrast needs entries "name1", "name2" for group names, as well as group memberships given in "group1" and "group2". Results of a differential test for this contrast are stored in an entry given by the method name, e.g. "MMD.locfit"

Author(s)

Gabriele Schweikert

See Also

DBAmmd-Accessors,getPeakReads

Examples

## Example using a small data set provided in the MMDiffBamSubset package

# setting the Experiment meta data:
ExpData <- list(dataDir=system.file("extdata", package="MMDiffBamSubset"),
           sampleSheet="Cfp1.csv")

MetaData <- list('ExpData' = ExpData)

# Creating a DBAmmd data set:
MMD <- DBAmmd(MetaData)

dhelekal/mmdiff3 documentation built on July 25, 2019, 8:18 p.m.