DPhil research, University of Oxford

Uncertainty in Epidemiology

Uncertainty is often measured but rarely defined. What does it mean to be uncertain about a made-up quantity? How can we tie it to the real-world decisions it informs?

Defining uncertainty

We developed a first-principles decision-theoretic framework that directly relates uncertainty in epidemic models to real-world practical concerns [1]. The framework is presented as four principles:

Uncertainty arises from making decisions with incomplete information

We treat the modeller’s role as making an estimate aa of some unknown quantity zz, which induces a loss ℓ(a,z)\ell(a,z). We show that the choice of ℓ(a,z)\ell(a,z) defines the measure of uncertainty, while the model defines the missing data.

When planning ICU capacity, each unmet bed costs more than each unused one. With a loss that reflects this, uncertainty about peak demand becomes the expected cost of acting on the best estimate.

The missing data determines how much uncertainty can be reduced

In machine learning, irreducible uncertainty is usually what would remain with unlimited data. An epidemic only happens once, so we instead define the missing data as data that could practically be collected. Uncertainty that these data would remove is reducible, and the rest is irreducible.

When estimating how many people in a town are infected from a sample, the missing data are the untested residents. With a perfect test, all of the uncertainty is reducible. With an imperfect test, some remains even after testing everyone.

Reasoning about the data quantifies reducible uncertainty

The missing data are typically unknown in advance. The chosen model makes implicit assumptions about these, allowing us to quantify the expected uncertainty reduction, and thus estimate the value of collecting more data.

Wastewater samples in New Zealand covered only part of the population. By modelling what samples from the rest might show, we estimated how much more uncertainty full coverage would remove. This allows decision-makers to determine whether there is value in expanding the program.

More data isn’t guaranteed to reduce uncertainty

New data can contradict what has already been observed, thus increasing uncertainty. Ideally, new data would be expected to reduce uncertainty, and we show that coherence (the assumptions in the predictive model match the assumptions in the fitting model) is the required condition.

Existing wastewater sites may have been chosen because they give reliable measurements, so new sites could be noisier. For more data to be expected to reduce uncertainty, this assumption must be built into the model itself, not only into the predictions of new samples.

Case study

We applied this framework to a published model of COVID-19 in New Zealand [2], fitted to reported cases and wastewater data. Fitting the model with and without the wastewater data measured how much uncertainty those data removed. Modelling the samples that were not collected estimated how much more uncertainty sampling the whole population would remove.

We found that wastewater sampling substantially reduced uncertainty in the model on most days, although on a few days (where the wastewater sampling contradicted case data) the uncertainty increased. We also found that the value of expanding the wastewater sampling program would be limited [1].

Reported cases onlyWith observed wastewaterWith full wastewater coverage (expected)

Fixing overconfident and delayed estimates

EpiEstim [3] and EpiFilter [4] are two widely used tools for tracking how fast an epidemic is spreading. Each depends on a smoothing parameter that users usually leave at its default. We developed efficient methods for fitting these from the data, showing marked improvements in calibration and timeliness [5].

We used the fact that the sum of the log one-step-ahead predictive probabilities is equivalent to the model likelihood. Adapting the models to produce one-step-ahead predictions resulted in computationally efficient ways of calculating likelihoods, thus allowing us to find the posterior distribution of the smoothing parameters.

This work improves the robustness and timeliness of existing popular methods.

Coverage of 95% intervals
EpiEstim43%84%
EpiFilter89%93%
0%50%95%
Default smoothingEstimated smoothing
Days behind in detecting changes
EpiEstim9.7 days0.3 days
EpiFilter2.7 days0 days
0 days5 days10 days
Default smoothingEstimated smoothing

See more

Resources

  • Code · GitHubDecisionTheoreticEpi

    Julia and R code for the New Zealand wastewater case study.

  • Code · GitHubRobustRtEstimators

    Julia implementations of EpiEstim and EpiFilter with estimated smoothing, with worked examples.

  • Tutorials · JuliaRobust Rt estimation

    Worked examples of EpiEstim and EpiFilter, including likelihoods and estimated smoothing.

References

  1. Steyn, N., Mills, C., Shirvaikar, V., Bickford Smith, F., Donnelly, C. A., & Parag, K. V. (2026). A decision-theoretic framework for uncertainty quantification in epidemiological modelling. American Journal of Epidemiology. https://doi.org/10.1093/aje/kwag154
  2. Watson, L. M., Plank, M. J., Armstrong, B. A., et al., & Steyn, N. (2024). Jointly estimating epidemiological dynamics of Covid-19 from case and wastewater data in Aotearoa New Zealand. Communications Medicine. https://doi.org/10.1038/s43856-024-00570-3
  3. Cori, A., Ferguson, N. M., Fraser, C., & Cauchemez, S. (2013). A new framework and software to estimate time-varying reproduction numbers during epidemics. American Journal of Epidemiology. https://doi.org/10.1093/aje/kwt133
  4. Parag, K. V. (2021). Improved estimation of time-varying reproduction numbers at low case incidence and between epidemic waves. PLOS Computational Biology. https://doi.org/10.1371/journal.pcbi.1009347
  5. Steyn, N., & Parag, K. V. (2025). Robust uncertainty quantification in popular estimators of the instantaneous reproduction number. American Journal of Epidemiology. https://doi.org/10.1093/aje/kwaf165