This is a helper function that compute the sum of squares between two theoretical and observed quantiles of a beta distribution (typically the lower and upper bounds of a confidence interval). This function is for internal use to find the best-fit beta distribution for a given confidence interval.
ssBetaPars(abPars, qLow, qUpp, alpha = 0.05)
The shape1 and shape2 parameters of the theoretical beta distribution.
The observed lower quantile.
The observed upper quantile.
The confidence level; i.e. the desired coverage is 1-alpha. Defaults to 0.05.
A single number, the sum of squares.