Reproduction number estimation using smoothing methods
Abstract
The reproduction number is estimated to understand the potential spread of an infectious disease within a population. Current estimation methods can be broadly categorized into mathematical and statistical methodological frameworks. Within the latter, renewal process, exponential growth and smoothing methods infer reproduction numbers from generation time distributions, growth rates, or smoothed epidemic trajectories. Our evaluation was based on three epidemic scenarios with distinct transmissibility levels (R0= 11.52, 2.88, and 1.06), simulated using stochastic SEIR models, and on the Hagelloch measles outbreak data. We compared five statistical methodsthe exponential growth method (EG), the renewal process method (RP), the Bayesian P-spline smoothing method (BPS), the filtering-based Bayesian inference method (FBI), and a smoothing method that integrates spline smoothing, derivative estimation, and exponential growth analysis within a generalized additive model (DSS). The performance criteria comprised the estimated bias, mean squared error (MSE), coverage rate (CR), significance rate (SR), and false negative rate (FNR). Across the three scenarios, no method was uniformly superior; the ranking depended on the criterion prioritized. RP and DSS produced the most accurate and stable point estimates (bias between − 0.45 and 1.69; MSE ≤ 7.31) together with the best-calibrated and narrowest 95% uncertainty intervals (coverage up to 97%). BPS was particularly unstable: its inflated mean drove MSE, bias, and coverage to 13.51%. FBI reached the highest significance rate at R0 = 1.06 (65.49%), although this sensitivity partly reflected their upward bias; EG showed the only substantial FNR (15.11%). The Hagelloch measles data corroborated these patterns. Under the significance constraint (SR), RP maximized risk control (CR/FNR) and DSS maximized accuracy (Bias/MSE).
Related
- Structures of the Basic Reproduction Number R0 Across Compartmental Disease Models
- Empirical Estimation of Ambient Contamination in Combinatorial Single-Cell Methods Using Multi-Reference Mapping
- Competing methods of divergence time estimation recover diverse evolutionary timescales.
- Estimation of Genomic Breeding Values for Economic Traits in Alxa-Type Cashmere Goats Using Bayesian Methods.
- A systematic review of biomass and carbon stock estimation approaches: Methods, uncertainties, and emerging opportunities
BibTeX
@article{barro2026reproduction,
title = {{Reproduction number estimation using smoothing methods}},
author = {Mahamadou Barro and Serge Manituo Aymar Somda and Cédric Stéphane Bationo and Bernard Eric Agodio Daboné and Betty Kazanga and Fabrice Anyirékun Somé and Benjamin Roche and Jean Gaudart},
journal = {Discover Public Health},
year = {2026},
month = oct,
doi = {10.1186/s12982-026-03048-w},
url = {https://inserm.hal.science/inserm-05841760v1},
note = {CC BY-NC-ND 4.0},
}
RIS
TY - JOUR
TI - Reproduction number estimation using smoothing methods
AU - Mahamadou Barro
AU - Serge Manituo Aymar Somda
AU - Cédric Stéphane Bationo
AU - Bernard Eric Agodio Daboné
AU - Betty Kazanga
AU - Fabrice Anyirékun Somé
AU - Benjamin Roche
AU - Jean Gaudart
PY - 2026
DA - 2026/10/04
JO - Discover Public Health
DO - 10.1186/s12982-026-03048-w
UR - https://inserm.hal.science/inserm-05841760v1
AB - The reproduction number is estimated to understand the potential spread of an infectious disease within a population. Current estimation methods can be broadly categorized into mathematical and statistical methodological frameworks. Within the latter, renewal process, exponential growth and smoothing methods infer reproduction numbers from generation time distributions, growth rates, or smoothed epidemic trajectories. Our evaluation was based on three epidemic scenarios with distinct transmissibility levels (R0= 11.52, 2.88, and 1.06), simulated using stochastic SEIR models, and on the Hagelloch measles outbreak data. We compared five statistical methodsthe exponential growth method (EG), the renewal process method (RP), the Bayesian P-spline smoothing method (BPS), the filtering-based Bayesian inference method (FBI), and a smoothing method that integrates spline smoothing, derivative estimation, and exponential growth analysis within a generalized additive model (DSS). The performance criteria comprised the estimated bias, mean squared error (MSE), coverage rate (CR), significance rate (SR), and false negative rate (FNR). Across the three scenarios, no method was uniformly superior; the ranking depended on the criterion prioritized. RP and DSS produced the most accurate and stable point estimates (bias between − 0.45 and 1.69; MSE ≤ 7.31) together with the best-calibrated and narrowest 95% uncertainty intervals (coverage up to 97%). BPS was particularly unstable: its inflated mean drove MSE, bias, and coverage to 13.51%. FBI reached the highest significance rate at R0 = 1.06 (65.49%), although this sensitivity partly reflected their upward bias; EG showed the only substantial FNR (15.11%). The Hagelloch measles data corroborated these patterns. Under the significance constraint (SR), RP maximized risk control (CR/FNR) and DSS maximized accuracy (Bias/MSE).
SN - 3005-0774
ER -