Diagnostic graphics for MCMC output

Description

Display trace plots, density plots or autocorrelation plots for the chains in the MCMC output. Each chain is plotted with a different colour.

Usage

1
2
3
4
5
tracePlot(object, ask=TRUE, ...)

densityPlot(object, ask=TRUE, ...)

acfPlot(object, lag.max=NULL, ask=TRUE, ...)

Arguments

object

An object of class Bwiqid.

ask

If TRUE, the user will be prompted before the next page of output is displayed.

lag.max

Maximum lag at which to calculate the acf; see acf.

...

Additional graphical parameters.

Value

Return NULL invisibly. Used for their plotting side-effects.

Author(s)

Mike Meredith

See Also

plotPost for a histogram and summary statistics.

Examples

1
2
3
4
5
6
7
8
9
# Create a fake Bwiqid object:
fake <- data.frame(mu = rnorm(3000), sigma=rlnorm(3000))
class(fake) <- c("Bwiqid", "data.frame")
attr(fake, "n.chains") <- 3
fake

tracePlot(fake)
densityPlot(fake, col=1:5, lwd=2)
acfPlot(fake, lag.max=10)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.