A probability region is also commonly called a credible region. For
more information on probability regions, see `p.interval`

.

Joint probability regions are plotted only for two variables, and the
regions are estimated with functions modified from the `car`

package. The internal ellipse functions assume bivariate normality.

This function is often used to plot posterior distributions of
samples, such as from the `LaplacesDemon`

function.