reads_demo: Forward reads and reverse reads in 'numeric' format (for demo...

In ArnaudDroitLab/RJMCMC: Bayesian hierarchical model for genome-wide profiling of nucleosome positions based on high-throughput short-read data (MNase-Seq data)

Description

A group of froward and reverse reads that can be used to test the rjmcmc function.

Usage

data(reads_demo)

Format

A list containing:

• readsFoward a vector of non-negative numeric, the start positions of the forward reads.

• readsReverse a vector of non-negative numeric, the start positions of the reverse reads. Beware that the start position of a reverse read is always higher that the end positition.

See Also

• rjmcmc for profiling of nucleosome positions

Examples

## Loading dataset
data(reads_demo)

## Nucleosome positioning
rjmcmc(startPosForwardReads = reads_demo$readsForward,
         startPosReverseReads = reads_demo$readsReverse,
         nbrIterations = 100, lambda = 3, kMax = 30,
         minInterval = 146, maxInterval = 292, minReads = 5)

ArnaudDroitLab/RJMCMC documentation built on May 5, 2019, 7:06 a.m.