Two changes ensure chains don’t get stuck in a loop:
n.walkers argument) to ensure ergodicity
(#11, @scheidan)lower.inits and upper.inits
are deprecated in favour of inits which leave more
flexibility to the user. Please read the detailed blog post for more
background about this change and how to migrate.inits can now be a data.frame or a
matrixd.e.mcmc() and s.m.mcmc() are not exported
any more. Please use the wrapper MCMCEnsemble()
instead.vignette("diagnostic-pkgs", package = "mcmcensemble"))
presenting two different options (coda and bayesplot) to plot and
evaluate the MCMC chains produced by mcmcensemble.max.iter %/% n.walkers == 1)coda = TRUE now correctly prompt the user to use
coda = FALSE if they do not wish to install coda.f. This is useful if you do something
like:p.log.named <- function(x) {
B <- 0.03
return(-x["a"]^2/200 - 1/2*(x["b"]+B*x["a"]^2-100*B)^2)
}MCMCEnsemble()
is now recorded in an additional attribute (accessible via
attr(res, "ensemble.sampler")).lower.inits and upper.inits have the same
namesd.e.mcmc() and
s.m.mcmc() now match those of
MCMCEnsemble()Suggests, instead of being
a hard dependency