identifyGammaPars {bootComb}

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

Description

Finds the best-fit gamma distribution parameters for a given confidence interval and returns the shape, rate parameters.

Usage

identifyGammaPars(qLow, qUpp, alpha = 0.05, initPars = c(1, 1), 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(1,1).

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 shape and rate for use with rgamma/dgamma/pgamma/qgamma.

See Also

ssGammaPars, optim, dgamma