coda.mcmc: Create an 'mcmc.list' object for 'coda' users

In ArchaeoPhases: Post-Processing of the Markov Chain Simulated by 'ChronoModel', 'Oxcal' or 'BCal'

Create an mcmc.list object for coda users

Description

This wrapper function extracts parallel chains from a data frame to create an mcmc.list object for use with coda diagnostic tools

Usage

coda.mcmc(data, numberChains = 1, iterationColumn = NULL)

Arguments

data	Data frame containing the output of the MCMC algorithm.
numberChains	Number of parallel chains, default = 1.
iterationColumn	Column number corresponding to the iteration values, default = NULL.

Value

An mcmc.list object.

Author(s)

Anne Philippe, [e-mail] and

Marie-Anne Vibet, [e-mail]

See Also

mcmc

mcmc.list

Examples

  data(Events)
  mcmcList = coda.mcmc(data = Events, numberChains = 3, iterationColumn = 1)
  plot(mcmcList)
  gelman.diag(mcmcList)
  # The multivariate criterion can not be evaluated when a phase
  # contains only one date. This induces colinearity problems.
  gelman.diag(mcmcList, multivariate = FALSE)

ArchaeoPhases documentation built on June 22, 2022, 1:05 a.m.