Interval estimates for diatom-inferred lake pH histories

We investigate ordinary least-squares and Bayesian methods for constructing interval estimates for historical lake pH's inferred from diatom sediments. The Bayesian method explicitly models several forms of variability, including the sampling and classification variability of the diatom records, estimation variability, and measurement error in observed pH's. The two methods produce similar interval estimates, but the Bayesian model allows design recommendations to be made.