Sample size allowed us to divide our dataset by substrate porosity and construct linear mixed-effects models using data from a single porosity to explore potential drivers of Zr: we regressed Zr values against a combination of quantitative variables (e.g., blood source, temperature, and class of biomolecule) which served as fixed effects and individual studies were included as a random effect.