Dk: Death Submove Probability

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

Description

Calculation of the death submove probability of a randomly selected nucleosome using a truncated Poisson distribution.

Usage

Dk(k, lambda, kMax = 30)

Arguments

k	a positive integer, the number of nucleosomes.
lambda	a positive numeric, the theorical mean of the Poisson distribution.
kMax	a positive integer, the maximum number of nucleosomes authorized. When k is equal or superior to kMax, the returned value is 0. Default: 30.

Value

a numeric value representing the calculated death submove probability. The value 0 when k is equal or superior to kMax or when k is equal to 1.

Author(s)

Rawane Samb

Examples

## Return the death submove probability
RJMCMC:::Dk(k = 12L, lambda = 2L, kMax = 30L)

## Zero is returned when k = 1
RJMCMC:::Dk(k = 1L, lambda = 3L, kMax = 30L)

## Zerio is returned when k is superior to kMax
RJMCMC:::Dk(k = 31L, lambda = 2L, kMax = 30L)

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