betabinom_HM {bayespm}

The Highest Mass (HM) interval of Beta-Binomial distribution.

Description

betabinom_HM is used to derive the narrowest interval of a predetermined coverage from a Beta-Binomial distribution, containing the Highest Mass (HM) values.

Usage

betabinom_HM( cover = NULL, n = NULL, a = NULL, b = NULL, plot = FALSE,
              xlab = "x", ylab = "Probability" )

Arguments

cover

scalar (between 0 and 1); The coverage (probability) of the Highest Mass (HM).

n

scalar (integer); It is the number of trials of the Binomial distribution and needs to be a number.

a

scalar (positive); It is the shape1 parameter of the Beta distribution and needs to be a number.

b

scalar (positive); It is the shape2 parameter of the Beta distribution and needs to be a number.

plot

logical; if TRUE, the plot of the HM region is displayed. It is FALSE by default.

xlab, ylab

The titles of the x-axis and y-axis for the HM plot. The default values are "x" and "Probability".

Details

betabinom_HM provides the Highest Mass (HM) of a Beta-Binomial distribution. HM is the narrowest region, with the minimum absolute difference from a pretermined coverage. The values included in HM have probability larger than every point outside the region.

Value

Returns a data frame with the upper and the lower limits of the Highest Predictive Denisty or Mass (HPrD/M) region, the lower and the achieved coverage. If plot=T, then the HPrD/M is visualized.

Examples

betabinom_HM(0.95, 10, 20, 180, plot = TRUE)