chainsSummary {basicMCMCplots} | R Documentation |
Compare summary statistics from multiple MCMC chains
Description
Plots median and 95
Usage
chainsSummary(
samplesList,
var = NULL,
nrows = NULL,
scale = FALSE,
width = 7,
height = NULL,
legend = !is.null(names(samplesList)),
legend.location = "topright",
jitter,
buffer = NULL,
buffer.right = NULL,
buffer.left = NULL,
cex = 1,
file = NULL
)
Arguments
samplesList |
List of arrays of MCMC samples from different chains |
var |
Parameter names to plot |
nrows |
Number of rows in the resulting plot |
scale |
Logical, whether to normalize each posterior chain |
width |
Width of figure |
height |
Height of figure |
legend |
Logical, whether to include a legend of chain names |
legend.location |
Legend location |
jitter |
Scale factor for spreading out lines from each chain |
buffer |
Buffer margin on both sides. Overrides buffer.right and buffer.left |
buffer.right |
Additional buffer on left side of plot |
buffer.left |
Additional buffer on right side of plot |
cex |
Expansion coefficient for text |
file |
Filename for saving figure to a file |
Examples
samples1 <- cbind(rnorm(1000, 1), rgamma(1000, 1), rpois(1000, 1))
colnames(samples1) <- c('alpha', 'beta', 'gamma')
samples2 <- cbind(rnorm(1000, 2), rgamma(1000, 2), rpois(1000, 2))
colnames(samples2) <- c('alpha', 'beta', 'gamma')
samplesList <- list(chain1 = samples1, chain2 = samples2)
chainsSummary(samplesList, nrow = 1, jitter = .3, buffer.left = .5, buffer.right = .5)
[Package basicMCMCplots version 0.2.7 Index]