Computes the shortest credible interval within which an unobserved parameter value falls with a particular probability.

## Value

A three-columns `numeric`

`matrix`

giving the lower and upper boundaries
of the credible interval and associated probability.

## See also

Other summary statistics:
`confidence_binomial()`

,
`confidence_mean()`

,
`confidence_multinomial()`

,
`interval_hdr()`

## Examples

```
## HDR of the Old Faithful eruption times
interval_hdr(faithful$eruptions)
#> start end p
#> [1,] 1.307160 2.837942 0.33
#> [2,] 3.150567 5.295819 0.62
```