Caterpillar Plots of opm Model Parameters

Share:

Description

Creates side-by-side plots of equal-tailed credible intervals of opm model parameters. The intervals are displayed as horizontal lines, with 90% interval using a thicker line width and 95% interval a thinner one. The posterior median is indicated with a dot.

Usage

1
2
caterplot(x, parm, main = paste("Caterpillar plot of", xname),
  xlab = "Range of parameter samples", labels = colnames(ranges))

Arguments

x

an instance of class opm

parm

a specification of which parameters are to be plotted, either a vector of names ("rho", "sig2" and "beta" are the only legal values) or a vector of positional indices. If missing, all parameters are considered.

main, xlab

useful defaults for the plot title and X-axis label

labels

labels for each parameter's interval: see axis

Value

A matrix of 2.5%, 5%, 50%, 95%, and 97.5% quantiles for each of the desired parameters, with parameters arranged in columns.

Examples

1
2
3
4
## Not run: 
caterplot(o, main = NULL, labels = expression(alpha, beta, sigma^2))

## End(Not run)