Bk: Birth 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 birth submove probability of adding a new nucleosome using a truncated Poisson distribution.

Usage

Bk(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. 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 birth submove probability
RJMCMC:::Bk(k = 14L, lambda = 1L, kMax = 30L)

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

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

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