Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling

Analyses of ecological data should account for the uncertainty in the process(es) that generated the data. However, accounting for these uncertainties is a difficult task, since ecology is known for its complexity. Measurement and/or process errors are often the only sources of uncertainty modeled when addressing complex ecological problems, yet analyses should also account for uncertainty in sampling design, in model specification, in parameters governing the specified model, and in initial and boundary conditions. Only then can we be confident in the scientific inferences and forecasts made from an analysis. Probability and statistics provide a framework that accounts for multiple sources of uncertainty. Given the complexities of ecological studies, the hierarchical statistical model is an invaluable tool. This approach is not new in ecology, and there are many examples (both Bayesian and non-Bayesian) in the literature illustrating the benefits of this approach. In this article, we provide a baseline for concepts, notation, and methods, from which discussion on hierarchical statistical modeling in ecology can proceed. We have also planted some seeds for discussion and tried to show where the practical difficulties lie. Our thesis is that hierarchical statistical modeling is a powerful way of approaching ecological analysis in the presence of inevitable but quantifiable uncertainties, even if practical issues sometimes require pragmatic compromises.

Stability in The Proportion of Harbor Seals Hauled Out Under Locally Ideal Conditions

We monitored the haul-out behavior of 68 radio-tagged harbor seals (Phoca vitulina) during the molt season at two Alaskan haul-out sites (Grand Island, August-September 1994; Nanvak Bay, August-September 2000). For each site, we created a statistical model of the proportion of seals hauled out as a function of date, time of day, tide, and weather covariates. Using these models, we identified the conditions that would result in the greatest proportion of seals hauled out. Although those “ideal conditions” differed between sites, the proportion of seals predicted to be hauled out under those conditions was very similar (81.3% for Grand Island and 85.7% for Nanvak Bay). The similar estimates for both sites suggest that haul-out proportions under locally ideal conditions may be constant between years and geographic regions, at least during the molt season.