birthMove: Birth move when the current number of nucleosomes is superior...

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

Description

Attempt to add a component for a new nucleosome. The function is specific to the case when when the current number of nucleosomes is superior to 1.

Usage

birthMove(paramValues, kValue, muValue, sigmafValue, sigmarValue, deltaValue,
  wValue, dfValue, aValue, dimValue)

Arguments

paramValues	a list containing: startPSF a vector of positive integer, the start position of all the forward reads. startPSR a vector of positive integer, the start position of all the reverse reads. kmax a integer the maximum number of nucleosomes allowed. lambda a positive numeric, the theorical mean of the Poisson distribution. minReads a integer, the minimum number of reads in a potential canditate region. y a vector of numeric, the position of all reads ( including forward and reverse reads). nr a integer, the number of reverse reads. nf a integer, the number of forward reads. nbrReads a integer, the total number of reads (including forward and reverse reads). zeta integer, the length of a nucleosome. deltamin a integer, the minimum distance between the maxima of the forward and reverse reads position densities for each nucleosome. deltamax a integer, the maximum distance between the maxima of the forward and reverse reads position densities for each nucleosome. minReadPos a numeric, the minimum position of the reads. maxReadPos a numeric, the maximum position of the reads.
kValue	a integer, the number of nucleosomes.
muValue	a vector of numeric of length kValue, the positions of the nucleosomes.
sigmafValue	a vector of numeric of length kValue, the variance of the forward reads for each nucleosome.
sigmarValue	a vector of numeric of length kValue, the variance of the reverse reads for each nucleosome.
deltaValue	a vector of numeric of length kValue, the distance between the maxima of the forward and reverse reads position densities for each nucleosome.
wValue	a vector of positive numerical of length kValue, the weight for each nucleosome.
dfValue	a vector of positive numerical of length kValue, the number of degrees of freedom.
aValue	a vector of positive numerical of length kValue + 1, the positions, on the chromosome, delimiting the regions of the reads associated with each nucleosome.
dimValue	a vector of positive numerical of length kValue, the number of reads associated to each nucleosome.

Value

a list containing:

• k a integer, the updated number of nucleosomes.

• mu a vector of numeric of length k, the updated positions of the nucleosomes.

• sigmaf a vector of numeric of length k, the updated variance of the forward reads for each nucleosome.

• sigmar a vector of numeric of length k, the updated variance of the reverse reads for each nucleosome.

• delta a vector of numeric of length kValue, the updated distance between the maxima of the forward and reverse reads position densities for each nucleosome.

• w a vector of positive numerical of length k, the updated weight for each nucleosome.

• df a vector of positive numerical of length k, the updated degrees of freedom.

• a a vector of positive numerical of length k + 1, the updated positions, on the chromosome, delimiting the regions of the reads associated with each nucleosome.

• dim a vector of positive numerical of length k, the updated number of reads associated to each nucleosome.

• rho a numeric, the acceptance probability.

Author(s)

Rawane Samb, Pascal Belleau, Astrid Deschenes

Examples

## Load dataset
data(reads_demo)

## Create a list containing the mandatory parameters
paramList <- list(kmax = 5L, nf = length(reads_demo$readsForward),
nr = length(reads_demo$readsReverse),
nbrReads = length(reads_demo$readsForward) + length(reads_demo$readsReverse),
y = sort(c(reads_demo$readsForward, reads_demo$readsReverse)),
startPSF = reads_demo$readsForward,
startPSR = reads_demo$readsReverse, lambda = 2,
zeta = 147, deltamin = 142, deltamax = 152,
minReadPos = min(c(reads_demo$readsReverse, reads_demo$readsForward)),
maxReadPos = max(c(reads_demo$readsReverse, reads_demo$readsForward)))

## Initial position of the nucleosomes
muPosition <- c(72800, 73000)
muPosition

## Birth move when more than one nucleosome
result <- RJMCMC:::birthMove(paramValues = paramList, kValue = 2L,
muValue = muPosition, sigmafValue = c(100, 120), sigmarValue = c(120, 100),
deltaValue = c(200, 210), wValue = c(0.5, 0.5), dfValue = c(3, 4),
aValue = c(min(c(reads_demo$readsReverse, reads_demo$readsForward)), 73150,
max(c(reads_demo$readsReverse, reads_demo$readsForward))),
dimValue = c(length(reads_demo$readsForward),
length(reads_demo$readsReverse)))

## Final position of the nucleosomes
result$mu[1:3]

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