crossrunsymm {crossrun}

Joint Probabilities for Crossings and Runs, Symmetric Case

Description

Joint probability distribution for the number of crossings C and the longest run L in a sequence of n independent Bernoulli observations with success probability p. To enhance precision, results are stored in mpfr arrays and the probabilities are multiplied by m^{n-1} for a multiplier m. This is for the symmetric case with success probability 0.5, in which the multiplied probabilities are integers for the default value 2 of the multiplier.

Usage

crossrunsymm(nmax = 100, mult = 2, prec = 120, printn = FALSE)

Arguments

nmax

; max sequence length.

mult

; multiplier for joint probabilities. Default 2.

prec

; mpft precision.

printn

; logical for including progress output.

Value

pt, list of joint probabilities, multiplied with m^{n-1}. In addition cumulative probabilities qt within each row are also included.

Examples

crs10 <- crossrunsymm(nmax=10,printn=TRUE)