Plot Overlaying Densities

Share:

Description

Creates a plot containing overlaying kernel density estimates from different MCMC simulations. This function is used in the denoverplot function to produce plots of overlaying densities for parameters in common from two different MCMC simulations.

Usage

1
2
3
denoverplot1(..., ci = NULL, col = NULL, lty = 1, xlim = NULL,
ylim = NULL, xlab = "", ylab = "Density", main = NULL, style = c("gray",
"plain"), gpar = NULL)

Arguments

...

one or more vectors or a list containing one or more vectors to be plotted.

ci

if non NULL, plots (100*ci)% credible interval limits on the density plots. The default (NULL) is not to plot the intervals.

col

one or more colors for the densities. Default is mcmcplotsPalette(n), where n is the number of elements in the ... argument.

lty

types of lines to plot.

xlim

limits for the x axis.

ylim

limits for the y axis.

xlab

label for the x axis.

ylab

label for the y axis.

main

main title for plot.

style

if "gray", then the plotting region is printed with a gray background, otherwise the default plotting region is used.

gpar

a list of additional graphical parameters to be passed to the plotting function. See help for par.

Value

Creates a plot.

Author(s)

S. McKay Curtis with contributions from Ilya Goldin

See Also

denoverplot, denplot, traplot1

Examples

1
2
3
4
5
denoverplot1(rnorm(1000), rnorm(1000))
denoverplot1(rnorm(1000, 0.0, 1.0), rnorm(1000, 0.1, 1.0),
             style="plain", col=mcmcplotsPalette(2, type="grayscale"),
             ci=0.95)
denoverplot1(list(rgamma(1000, 1, 1), rgamma(1000, 1, 1)))