This repository contains R scripts for Bayesian inference applied to evaluate the best-fit parameters for ice flow laws. It supports the paper:
Fan, S., Wang, T., Prior, D. J., Breithaupt, T., Hager, T. F., & Wallis, D. (2025). Flow laws for ice constrained by 70 years of laboratory experiments. Nature Geoscience, 18(4), 296β304.
DOI: 10.1038/s41561-025-01661-z
R/: Core R scripts for Bayesian analysis, data processing, and visualisation.data/input/: Input datasets.notebooks/: RMarkdown files with step-by-step analysis.
Ensure you have the necessary R packages installed:
install.packages(c("R2jags", "ggplot2", "coda", "bayesplot", "lattice", "gridExtra"))