bayesLife (version 5.3-1)

get.e0.mcmc: Accessing MCMC Results

Description

The function get.e0.mcmc retrieves results of an MCMC simulation and creates an object of class bayesLife.mcmc.set. Function has.e0.mcmc checks the existence of such results. Function e0.mcmc extracts a single chain, and e0.mcmc.list extracts several or all chains from the simulation results.

Usage

get.e0.mcmc(sim.dir = file.path(getwd(), "bayesLife.output"), 
    chain.ids = NULL, low.memory = TRUE, burnin = 0, verbose = FALSE)
    
has.e0.mcmc(sim.dir)
e0.mcmc(mcmc.set, chain.id = 1)
e0.mcmc.list(mcmc.set, chain.ids = NULL)

Value

get.e0.mcmc returns an object of class bayesLife.mcmc.set. has.e0.mcmc returns a logical value.

e0.mcmc returns an object of class bayesLife.mcmc, and e0.mcmc.list returns a list of bayesLife.mcmc objects.

Arguments

sim.dir

Directory where the simulation results are stored.

chain.ids

Chain identifiers in case only specific chains should be included in the resulting object. By default, all available chains are included.

low.memory

If FALSE full MCMC traces are loaded into memory.

burnin

Burnin used for loading traces. Only relevant, if low.memory=FALSE.

verbose

Logical switching log messages on and off.

chain.id

Chain identifier.

mcmc.set

Object of class bayesLife.mcmc.set.

Author

Hana Sevcikova

See Also

bayesLife.mcmc.set

Examples

sim.dir <- file.path(find.package("bayesLife"), "ex-data", "bayesLife.output")
m <- get.e0.mcmc(sim.dir)
summary(m)

# summary of the world parameters for a single chain
# (the same as above since there is only one chain in this toy example)
summary(e0.mcmc.list(m)[[1]], par.names.cs = NULL)

# the same as
summary(e0.mcmc(m, chain.id = 1), par.names.cs = NULL)