gb2_HD {bayespm}

The Highest Density (HD) interval of Generalized Beta of the second kind distribution.

Description

gb2_HD is used to derive the narrowest interval of a predetermined coverage from a Generalized Beta of the second kind distribution, containing the Highest Density (HD) values.

Usage

gb2_HD( cover = NULL, a = NULL, c = NULL, d = NULL, plot = FALSE,
        xlab = "x", ylab = "Density" )

Arguments

cover

scalar (between 0 and 1); The coverage (probability) of the Highest Density (HD).

a

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

c

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

d

scalar (positive); It is the scale parameter of the Generalized Beta of the second kind distribution and needs to be a number.

plot

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

xlab, ylab

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

Details

gb2_HD provides the Highest Density (HD) of a Generalized Beta of the second kind distribution. HD is the narrowest interval with a pretermined coverage. The values included in HD have density 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

gb2_HD(0.95, 4, 6, 6, plot = TRUE)