Finds the best-fit negative binomial distribution parameters for a given confidence interval and returns the size, prob parameters.
identifyNegBinPars(
qLow,
qUpp,
alpha = 0.05,
initPars = c(10, 0.5),
maxiter = 1000
)The observed lower quantile.
The observed upper quantile.
The confidence level; i.e. the desired coverage is 1-alpha. Defaults to 0.05.
A vector of length 2 giving the initial parameter values to start the optimisation; defaults to c(10,0.5).
Maximum number of iterations for optim. Defaults to 1e3. Set to higher values if convergence problems are reported.
A vector of length 2 giving the 2 parameters size and prob for use with rnbinom/dnbinom/pnbinom/qnbinom.