identifyBetaPars {bootComb}

Determine the parameters of the best-fit beta distribution for a given confidence interval for a probability parameter.

Description

Finds the best-fit beta distribution parameters for a given confidence interval for a probability parameter and returns the shape1, shape2 parameters.

Usage

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

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 shape1 and shape1 for use with rbeta/dbeta/pbeta/qbeta.

See Also

ssBetaPars, optim, dbeta