rbacon (version 2.4.2)

scissors: Remove the first n iterations.

Description

Removes iterations of the MCMC time series, and then updates the output file.

Usage

scissors(burnin, set = get("info"))

Arguments

burnin

Number of iterations to remove of the iterative time series. If this value is higher than the amount of remaining iterations, a warning is given and the iterations are not removed. If the provided number is negative, the iterations will be removed from the end of the run, not from the start. If a range is given, this range of iterations is removed.

set

Detailed information of the current run, stored within this session's memory as variable info.

Value

NA

Details

Bacon will perform millions of MCMC iterations for each age-model run by default, although only a fraction of these will be stored. In most cases the remaining MCMC iterations will be well mixed (the upper left panel of the fit of the iterations shows no undesirable features such as trends or sudden systematic drops or rises). If the run has a visible remaining burn-in, scissors can be used to remove them. To remove, e.g., the first 300 iterations, type scissors(300). To remove the last 300 iterations, type scissors(-300). To remove iterations 300 to 600, type scissors(300:600).

References

Blaauw, M. and Christen, J.A., Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6 (2011), no. 3, 457--474. https://projecteuclid.org/euclid.ba/1339616472

See Also

http://www.qub.ac.uk/chrono/blaauw/manualBacon_2.3.pdf

Examples

# NOT RUN {
  Bacon(ask=FALSE, coredir=tempfile())
  scissors(100)
  agedepth()
# }
# NOT RUN {
# }