credible.region function

Compute a simultaneous credible region (rectangle) from a sample for a vector valued parameter.