R/plot.bugs.R

In MultiBUGS/R2MultiBUGS: Running MultiBUGS from R

#' Plotting a bugs object
#'
#' Plotting a \code{bugs} object
#'
#' @param x an object of class `bugs', see \code{\link{bugs}} for details
#' @param display.parallel display parallel intervals in both halves of the
#' summary plots; this is a convergence-monitoring tool and is not necessary
#' once you have approximate convergence (default is \code{FALSE})
#' @param ... further arguments to \code{\link{plot}}
#' @seealso \code{\link{bugs}}
#' @keywords hplot
#' @export
plot.bugs <- function(x, display.parallel = FALSE, ...){
  mar.old <- par("mar")
  pty.old <- par(pty = "m")
  mfrow.old <- par("mfrow")
  layout(matrix(c(1, 2), 1, 2))

  bugs.plot.summary(x, ...)
  bugs.plot.inferences(x, display.parallel, ...)
  header <- ""
  if (!is.null(x$model.file)){
    header <- paste(header, "Bugs model at \"", x$model.file, "\", ", sep = "")
  }
  if (!is.null(x$program)){
    header <- paste(header, "fit using ", x$program, ", ", sep = "")
  }
  header <- paste(header,
                  x$n.chains,
                  " chains, each with ",
                  x$n.iter,
                  " iterations (first ",
                  x$n.burnin,
                  " discarded)",
                  sep = "")
  mtext(header, outer = TRUE, line = -1, cex = 0.7)
  par(pty = pty.old[[1]], mar = mar.old, mfrow = mfrow.old)
}