caterplot function

Caterpillar Plots of opm Model Parameters

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.

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

Returns

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

Examples

## Not run:

caterplot(o, main = NULL, labels = expression(alpha, beta, sigma^2))
## End(Not run)