- Allow flat structure by setting K_factor > K_fine
- Update stan code to reflect new array syntax
- Hierarchical structure of sections changed so that child sections are offset from parents. This will eventually break backwards compatibility, for now if the old argument "K" is used, "K_fine" is calculated from "K". The change improves the age models by reducing the influence of the positions of the breaks in the low resolution layers.
- 1 sigma quantiles now included in summaries
- calibration function can handle offset uncertainties (e.g. reservoir age uncertainty)
- Hamstr logo added to Readme
- patchwork used instead of ggpubr to reduce dependencies
- some little used arguments moved to hamstr_control and stan_sampler_args to simplify main function
- default K structure changed to powers of 2
- ability to model hiatuses reimplemented
- Smoothing of accumulation rates in plots and output
- Ability added to model age-heterogeneity from bioturbation and its effect on age models.
- Default plot changed to put diagnostic plots at the bottom.
-
Minor change to the way acc_shape is adjusted to account for the number of hierarchical levels. acc_shape is now adjusted to control the shape of the gamma distributed alpha parameters, not the total variance. The difference is minor.
-
New methods to extract and plot accumulation rates
New default hierarchical structure K based on default resolution of 1 cm (up to 1000 cm) and number of levels approximately equal to number of new child sections per parent per level.
- Significant update to behaviour when using hierarchical structure. The acc_shape parameter is now adjusted to control the total variance in the alpha parameters when there are multiple levels
- Default mem_mean and mem_strength values updated to 0.5 and 10 to match updated defaults from bacon 2.5.1 onwards
- Methods written for generic functions plot and summary
- Package and main function renamed to hamstr (hierarchical accumulation models with Stan and R)
- dev branch with hierarchical sections merged to master
- Hiatuses can now be modelled
- Added more information to the summary figure