segmentation {RJMCMCNucleosomes} | R Documentation |
GRanges
containing reads in a list of smaller
segments for the rjmcmc
function.Split a GRanges
of reads (as example, the reads from
a chromosome) in a list
of smaller GRanges
sot that the
rjmcmc
function can be run on each segments.
segmentation(reads, zeta = 147, delta, maxLength)
reads |
a |
zeta |
a positive |
delta |
a positive |
maxLength |
a positive |
a GRangesList
containing all the segments.
Pascal Belleau, Astrid Deschenes
## Load synthetic dataset of reads data(syntheticNucleosomeReads) ## Use dataset of reads to create GRanges object sampleGRanges <- GRanges(seqnames = syntheticNucleosomeReads$dataIP$chr, ranges = IRanges(start = syntheticNucleosomeReads$dataIP$start, end = syntheticNucleosomeReads$dataIP$end), strand = syntheticNucleosomeReads$dataIP$strand) # Segmentation of the reads segmentation(reads = sampleGRanges, zeta = 147, delta = 50, maxLength = 1000)