reads_demo_02: Forward reads and reverse reads in 'integer' 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_02)

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_02)

## Nucleosome positioning
rjmcmc(startPosForwardReads = reads_demo_02$readsForward,
         startPosReverseReads = reads_demo_02$readsReverse,
         nbrIterations = 150, lambda = 3, kMax = 30,
         minInterval = 144, maxInterval = 290, minReads = 6)

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