Automatically calculate and apply burnin value

Automatically calculate and apply burnin value

autoburnin(jags_out, return.burnin = FALSE, ...)

Arguments

jags_out

JAGS output

return.burnin

Logical. If TRUE, return burnin value in addition to samples (as list). Default = FALSE.

...

Additional arguments for getBurnin, gelman_diag_mw, and gelman.diag.

Author

Michael Dietze, Alexey Shiklomanov

Examples

     z1 <- coda::mcmc(c(rnorm(2500, 5), rnorm(2500, 0)))
     z2 <- coda::mcmc(c(rnorm(2500, -5), rnorm(2500, 0)))
     z <- coda::mcmc.list(z1, z2)
     z_burned <- autoburnin(z)