identifyNegBinPars {bootComb}

Determine the parameters of the best-fit negative binomial distribution for a given confidence interval.

Description

Finds the best-fit negative binomial distribution parameters for a given confidence interval and returns the size, prob parameters.

Usage

identifyNegBinPars(
  qLow,
  qUpp,
  alpha = 0.05,
  initPars = c(10, 0.5),
  maxiter = 1000
)

Arguments

qLow

The observed lower quantile.

qUpp

The observed upper quantile.

alpha

The confidence level; i.e. the desired coverage is 1-alpha. Defaults to 0.05.

initPars

A vector of length 2 giving the initial parameter values to start the optimisation; defaults to c(10,0.5).

maxiter

Maximum number of iterations for optim. Defaults to 1e3. Set to higher values if convergence problems are reported.

Value

A vector of length 2 giving the 2 parameters size and prob for use with rnbinom/dnbinom/pnbinom/qnbinom.

See Also

ssNegBinPars, optim, dnbinom