as.bugs.array: Convert to bugs object

In MultiBUGS/R2MultiBUGS: Running MultiBUGS from R

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)

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 MultiBUGS 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

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.

If DIC=TRUE then DIC can be either already passed to argument DICOutput or calculated from deviance values in sims.array.

Value

A bugs object is returned

Author(s)

Jouni Kerman, [e-mail] with modification by Andrew Gelman,
[e-mail], packaged by Uwe Ligges, [e-mail].

See Also

bugs

MultiBUGS/R2MultiBUGS documentation built on Aug. 14, 2019, 3:15 p.m.