boxplot.true: Compare Boxplots to True Values

Description

Plots side-by-side boxplots of the columns of the matrix x. Each boxplot can have its own reference line (truth) and standard error lines se.truth, if desired. This function was originally written to display MCMC output, where the reference lines were true values used to test an MCMC simulation.

Usage

BoxplotTrue(x, truth = NULL, vnames = NULL, center = FALSE,
            se.truth = NULL, color = "white", truth.color = "black",
            ylim = NULL, ...)

Value

called for its side effect

Arguments

x: The matrix whose columns are to be plotted.
truth: (optional) A vector of reference values with length equal to ncol(x).
vnames: (optional) character vector giving the column names of x.
center: (optional) logical. If truth is supplied then center=TRUE will center each column of x around truth to show the variation around the reference line.
se.truth: (optional) numeric vector of length ncol(x). If truth is supplied then additional reference lines will be drawn at truth +/- 2*se.truth.
color: (optional) vector of colors for each boxplot.
truth.color: A color (or vector of colors) to use for the segments representing true values.
ylim: Limits for the vertical axis. If NULL then the axis will be scaled to fit x, truth, and truth + c(2, -2) * se.truth
...: additional arguments to boxplot.

Author

Steven L. Scott

Examples

Run this code

x <- t(matrix(rnorm(5000, 1:5, 1:5), nrow=5))
BoxplotTrue(x, truth=1:5, se.truth=1:5, col=rainbow(5), vnames =
  c("EJ", "TK", "JT", "OtherEJ", "TJ") )