adaptMCMC (version 1.5)

convert.to.coda: Converts chain(s) into coda objects.

Description

Converts chain(s) produced by MCMC or MCMC.parallel into coda objects.

Usage

convert.to.coda(sample)

Value

An object of the class mcmc or mcmc.list.

Arguments

sample

output of MCMC or MCMC.parallel.

Author

Andreas Scheidegger, or

Details

Converts chain(s) produced by MCMC or MCMC.parallel so that they can be used with functions of the coda package.

See Also

MCMC, mcmc, mcmc.list

Examples

## ----------------------
## Banana shaped distribution

## log-pdf to sample from
p.log <- function(x) {
  B <- 0.03                              # controls 'bananacity'
  -x[1]^2/200 - 1/2*(x[2]+B*x[1]^2-100*B)^2
}


## ----------------------
## generate 200  samples

samp <- MCMC(p.log, n=200, init=c(0, 1), scale=c(1, 0.1),
               adapt=TRUE, acc.rate=0.234)


## ----------------------
## convert in object of class 'mcmc'
samp.coda <- convert.to.coda(samp)

class(samp.coda)

## ----------------------
## use functions of package 'coda'

require(coda)

plot(samp.coda)
cumuplot(samp.coda)