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.
DBAmmd()returns an empty DBAmmd Object.
DBAmmd(MetaData) initializes a DBAmmd Object for a new
(See below and the package vignette for more details.)
List containing an
AnaDatacompartment. "ExpData" needs a
genomeentry, which should be a valid
BSGenomename, 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
AnaDataentry is used to store and access parameters for the MMDiff2 Analysis, like the sigma of the RBF Kernel.
GRanges object containing Regions of Interests (Peaks)
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
getPeakReadsto fill this slot and
estimateFragmentCentersto add the (estimated) positions of fragment centers.
m x n matrix containing total counts of reads mapping to m peaks in n samples (including input samples)
m x n matrix containing counts of reads mapping to positive (forward) strand
m x n matrix containing counts of reads mapping to negative (reverse) strand
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.
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.
(for internal use only)
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"
1 2 3 4 5 6 7 8 9 10
## 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)