bayesTFR (version 1.2-0)

get.tfr.parameter.traces: Accessing MCMC Parameter Traces

Description

Functions for accessing traces of the MCMC parameters, either country-independent or country-specific.

Usage

get.tfr.parameter.traces(mcmc.list, par.names = tfr.parameter.names(), 
    burnin = 0, thinning.index = NULL, thin = NULL)
    
get.tfr.parameter.traces.cs(mcmc.list, country.obj, 
    par.names=tfr.parameter.names.cs(), 
    burnin=0, thinning.index=NULL, thin=NULL)

Value

Both functions return a matrix with columns being the parameters and rows being the MCMC values, attached to one another in case of multiple chains. get.tfr.parameter.traces contains country-independent parameters, get.tfr.parameter.traces.cs contains country-specific parameters.

Rdversion

1.1

See Also

coda.mcmc.list for another way of retrieving parameter traces.

Examples

sim.dir <- file.path(.find.package("bayesTFR"), "ex-data", 'bayesTFR.output')
m <- get.tfr.mcmc(sim.dir)
tfr.values <- get.tfr.parameter.traces(m$mcmc.list, burnin=10, par.names="sigma0")
hist(tfr.values, main=colnames(tfr.values))
tfr.values.cs <- get.tfr.parameter.traces.cs(m$mcmc.list, 
                    get.country.object('Canada', meta=m$meta),
                    burnin=10, par.names="Triangle_c4")
hist(tfr.values.cs, main=colnames(tfr.values.cs))