hd.beta: HD interval for the beta distribution

In bssbinom: Bayesian Sample Size for Binomial Proportions

HD interval for the beta distribution

Description

Computes the highest density interval for the beta distribution by optimization or simulation.

Usage

hd.beta(c, d, rho = NULL, len = NULL, N = 1000)

Arguments

c	First parameter of the beta distribution.
d	Second parameter of the beta distribution.
rho	A number in (0, 1) representing the fixed probability for the HD interval.
len	A positive real number representing the fixed length for the HD interval.
N	Number of replicates used in the simulation when the length is fixed. Default is 1000.

Value

Lower and upper bounds of the HD interval.

Note

For fixed probability (rho) of the HD interval the function uses betaHPD function from pscl R package. For fixed length of the HD interval the function uses an algorithm proposed in M'Lan et al. (2006). For the case that uses Monte Carlo simulations, the provided interval may vary from one call to the next. The difference is expected to decrease as the number of replicates (N) increases.

References

M’Lan, C.E., Joseph, L., Wolfson, D.B. (2006). Bayesian sample size determination for case-control studies. Journal of the American Statistical Association, 101, 760–772.

Examples

hd.beta(c = 2, d = 3, rho = 0.95)

hd.beta(c = 2, d = 3, len = 0.3)

bssbinom documentation built on June 8, 2025, 10:38 a.m.