synlik (version 0.1.0)

smcmc-class: smcmc-class

Description

Object representing the results of MCMC estimation on an object of class synlik, from which it inherits.

References

Vihola, M. (2011) Robust adaptive Metropolis algorithm with coerced acceptance rate. Statistics and Computing.

Examples

# Load "synlik" object
data(ricker_sl)

plot(ricker_sl)

# MCMC estimation
set.seed(4235)
ricker_sl <- smcmc(ricker_sl,
                   initPar = c(3.2, -1, 2.6),
                   niter = 50,
                   burn = 3,
                   priorFun = function(input, ...) 1,
                   propCov = diag( c(0.1, 0.1, 0.1) )^2,
                   nsim = 200,
                   multicore = FALSE)

# Continue with additional 50 iterations
ricker_sl <- continue(ricker_sl, niter = 50)

plot(ricker_sl)