MCMCout: Organize mcmc.list objects from coda/R2jags

In vlahm/MCMCwrap: Organize MCMC output for analysis and plotting

Description

Get posterior draws, means, and modes output by MCMC simulations.

Usage

MCMCout(mcmc.list)

Arguments

x	an mcmc.list object.
...	additional arguments passed to 'density', which determines posterior modes.

Details

Takes mcmc.list objects, which are created by various functions in packages 'R2jags' and 'coda'. Returns posterior draws as arrays with dim(observations, output params, draws). Also returns means and X values of modes (corresponding to highest posterior peaks) as matrices with dim(observations, output params). Works across all parallel MCMC chains (if running several in tandem).

Value

a list containing arrays of 1) all posterior draws from all parallel chains, 2) posterior means, 3) posterior modes.

Author(s)

Mike Vlah, [e-mail]

See Also

MCMCout_mix for organized mixing model output.

Examples

mod1 <- jags.model(model_script, data = list(<data>), 
        inits = function()list(<inits>), n.chains = 3, n.adapt = 2000)
update(mod1, n.iter=5e3)
mod1_out <- coda.samples(mod1, c(<variables out>, n.iter=1e5, thin=100))
MCMCout(mod1_out)

vlahm/MCMCwrap documentation built on May 3, 2019, 6:16 p.m.