This repository contains R scripts for Bayesian inference applied to evaluate the best-fit parameters for ice flow laws. It supports the paper:
Sheng Fan et al. (2025). "Flow laws for ice constrained by 70 years of laboratory experiments". Nature Geoscience
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"))