coda.mcmc {ArchaeoPhases}

R Documentation

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, Anne.Philippe@univ-nantes.fr and

Marie-Anne Vibet, Marie-Anne.Vibet@univ-nantes.fr

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)

---