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 of some unknown quantity , which induces a loss . We show that the choice of 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].
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.
See more
Resources
- DecisionTheoreticEpi
Julia and R code for the New Zealand wastewater case study.
- RobustRtEstimators
Julia implementations of EpiEstim and EpiFilter with estimated smoothing, with worked examples.
- Robust Rt estimation
Worked examples of EpiEstim and EpiFilter, including likelihoods and estimated smoothing.
References
- 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
- 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
- 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
- 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
- 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