Description Usage Format Source References See Also Examples
The object contains a small subset of the CIP-Seq data sets generated to assess the link between the histone modification states of H3K4me3 with respect to the mediator proteins Cfp1 in mouse [1]. The data is available as part of ArrayExpress Experiment E-ERAD-79. For more details see the MMDiffBamSubset data package.
1 |
A DBA object containing the additional component MD. A DBA object is an S3 object (class "DBA"), which is introduced with the DiffBind package and more information on this class can be found in the Diffbind vignette.
The added field MD is a list containing the following elements:
- RawTotalCounts: matrix of total counts per peak and sample (nPeaks x nSamples)
- PeakRawHists: list of length nPeaks, containing for each Peak a Matrix of histograms (nSamples x nbins). Note, as Peaks can vary in length, nbins may be different for each Peak.
See the MMDiffBamSubset data package for more details.
[1] Clouaire T et al (2012). Cfp1 integrates both CpG content and gene activity for accurate H3K4me3 deposition in embryonic stem cells. Genes Dev. August 1, 2012 26: 1714–1728
dba, summary.DBA
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | ## Not run:
#The object has been generated using the following commands:
##### STEP1: create DBA object:
library('MMDiffBamSubset')
oldwd <- setwd(system.file("extdata", package="MMDiffBamSubset"))
Cfp1 <- dba(sampleSheet="Cfp1.csv", minOverlap=3,
config = data.frame(RunParallel=FALSE))
##### STEP2: compute histograms from bam files
bin.length <- 50
Peaks <- dba.peakset(Cfp1,bRetrieve=TRUE)
Peaks <- Peaks[1:1000]
Cfp1Profiles <- getPeakProfiles(Cfp1,Peaks,bin.length=bin.length,
save.files=FALSE,run.parallel=FALSE)
setwd(oldwd)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.