as.bugs.array: Convert to bugs object

Description

Function converting results from Markov chain simulations, that might not be from BUGS, to bugs object. Used mainly to display results with plot.bugs.

Usage

as.bugs.array(sims.array, model.file=NULL, program=NULL,
    DIC=FALSE, DICOutput=NULL, n.iter=NULL, n.burnin=0, n.thin=1)

Value

A bugs object is returned

Arguments

sims.array: 3-way array of simulation output, with dimensions n.keep, n.chains, and length of combined parameter vector.
model.file: file containing the model written in WinBUGS code
program: the program used
DIC: logical; whether DIC should be calculated, see also argument DICOutput and details
DICOutput: DIC value
n.iter: number of total iterations per chain used for generating sims.array
n.burnin: length of burn in, i.e. number of iterations to discarded at the beginning for generating sims.array
n.thin: thinning rate, a positive integer, used for generating sims.array

Author

Jouni Kerman, kerman@stat.columbia.edu with modification by Andrew Gelman, gelman@stat.columbia.edu, packaged by Uwe Ligges, ligges@statistik.tu-dortmund.de.

Details

This function takes a 3-way array of simulations and makes it into a bugs object that can be conveniently displayed using print and plot and accessed using attach.bugs. If the third dimension of sims() has names, the resulting bugs object will respect that naming convention. For example, if the parameter names are “alpha[1]”, “alpha[2]”, ..., “alpha[8]”, “mu”, “tau”, then as.bugs.array will know that alpha is a vector of length 8, and mu and tau are scalar parameters. These will all be plotted appropriately by plot and attached appropriately by attach.bugs.