{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,1,6]],"date-time":"2026-01-06T04:12:07Z","timestamp":1767672727776,"version":"3.48.0"},"reference-count":36,"publisher":"The Royal Society","issue":"233","license":[{"start":{"date-parts":[[2025,12,17]],"date-time":"2025-12-17T00:00:00Z","timestamp":1765929600000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"},{"start":{"date-parts":[[2025,12,17]],"date-time":"2025-12-17T00:00:00Z","timestamp":1765929600000},"content-version":"vor","delay-in-days":0,"URL":"http:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/100025287","name":"Center for Forecasting and Outbreak Analytics","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100025287","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000060","name":"National Institute of Allergy and Infectious Diseases","doi-asserted-by":"publisher","id":[{"id":"10.13039\/100000060","id-type":"DOI","asserted-by":"publisher"}]},{"DOI":"10.13039\/100000057","name":"National Institute of General Medical Sciences","doi-asserted-by":"crossref","id":[{"id":"10.13039\/100000057","id-type":"DOI","asserted-by":"crossref"}]},{"name":"National Science Foundation"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":[],"published-print":{"date-parts":[[2025,12,17]]},"abstract":"<jats:title>Abstract<\/jats:title>\n                  <jats:p>The time-varying basic reproduction number, R0(t), is a key epidemiological metric that quantifies the transmissibility of an infectious pathogen at time t. Accurate estimation and uncertainty quantification of R0(t) are crucial for understanding disease dynamics and informing public health decision-making. In this study, we evaluate six methods for estimating R0(t) using synthetic data generated from a stochastic Susceptible-Infected-Recovered (SIR) model with imposed changes to pathogen transmissibility and empirical COVID-19 case data. The methods include ensemble filter methods and inflation techniques, which are employed to mitigate covariance underestimation and filter divergence. For synthetic data, we compare the ensemble adjustment Kalman filter (EAKF) with no inflation, fixed inflation, and adaptive inflation, and the ensemble square root smoother (EnSRS) with adaptive inflation. For empirical data, we also compare with EpiEstim and EpiFilter. Our results demonstrate that the EAKF and EnSRS methods with adaptive inflation outperform other approaches in accurately estimating R0(t), particularly in scenarios with abrupt changes in transmission rates. The adaptive inflation techniques effectively address covariance underestimation and filter divergence, leading to more robust and reliable estimates of R0(t). These findings highlight the potential of adaptive inflation methods for improving the accuracy of time-varying parameter inference, contributing to more effective public health responses.<\/jats:p>","DOI":"10.1098\/rsif.2025.0131","type":"journal-article","created":{"date-parts":[[2025,12,17]],"date-time":"2025-12-17T09:38:41Z","timestamp":1765964321000},"source":"Crossref","is-referenced-by-count":0,"title":["Data assimilation for estimating time-varying reproduction numbers"],"prefix":"10.1098","volume":"22","author":[{"ORCID":"https:\/\/orcid.org\/0009-0000-6592-2987","authenticated-orcid":false,"given":"Han Yong","family":"Wunrow","sequence":"first","affiliation":[{"name":"Columbia University 1 , , ,","place":["New York, NY, USA"],"department":["Department of Applied Physics and Applied Mathematics"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Sen","family":"Pei","sequence":"additional","affiliation":[{"name":"Columbia University 2 , , ,","place":["New York, NY, USA"],"department":["Department of Environmental Health Sciences"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Jeffrey","family":"Shaman","sequence":"additional","affiliation":[{"name":"Columbia University 2 , , ,","place":["New York, NY, USA"],"department":["Department of Environmental Health Sciences"]},{"name":"Columbia Climate School 3 , , ,","place":["New York, NY, USA"]}],"role":[{"role":"author","vocabulary":"crossref"}]},{"given":"Marc","family":"Spiegelman","sequence":"additional","affiliation":[{"name":"Columbia University 1 , , ,","place":["New York, NY, USA"],"department":["Department of Applied Physics and Applied Mathematics"]},{"name":"Columbia University 4 , , ,","place":["New York, NY, USA"],"department":["Department of Earth and Environmental Sciences"]}],"role":[{"role":"author","vocabulary":"crossref"}]}],"member":"175","published-online":{"date-parts":[[2025,12,17]]},"reference":[{"key":"2026010523072043300_R1","doi-asserted-by":"crossref","DOI":"10.3201\/eid2501.171901","article-title":"Complexity of the basic reproduction number (R0)","volume":"25","author":"Delamater","year":"2019","journal-title":"Emerging Infect. Dis."},{"key":"2026010523072043300_R2","doi-asserted-by":"crossref","first-page":"e32","DOI":"10.2105\/AJPH.2013.301704","article-title":"Unraveling R0: considerations for public health applications","volume":"104","author":"Ridenhour","year":"2014","journal-title":"Am. J. Public Health"},{"key":"2026010523072043300_R3","doi-asserted-by":"crossref","first-page":"1505","DOI":"10.1093\/aje\/kwt133","article-title":"A new framework and software to estimate time-varying reproduction numbers during epidemics","volume":"178","author":"Cori","year":"2013","journal-title":"Am. J. Epidemiol."},{"key":"2026010523072043300_R4","doi-asserted-by":"crossref","first-page":"509","DOI":"10.1093\/aje\/kwh255","article-title":"Different epidemic curves for severe acute respiratory syndrome reveal similar impacts of control measures","volume":"160","author":"Wallinga","year":"2004","journal-title":"Am. J. Epidemiol."},{"key":"2026010523072043300_R5","doi-asserted-by":"crossref","DOI":"10.1126\/science.abb9789","article-title":"Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions","volume":"369","author":"Dehning","year":"2020","journal-title":"Science"},{"key":"2026010523072043300_R6","first-page":"1002","article-title":"Estimating the instantaneous reproduction number (Rt) by using particle filter","volume":"8","author":"Won","year":"2023","journal-title":"Infect. Dis. Model."},{"key":"2026010523072043300_R7","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pone.0002185","article-title":"Real time Bayesian estimation of the epidemic potential of emerging infectious diseases","volume":"3","author":"Bettencourt","year":"2008","journal-title":"PLoS One"},{"key":"2026010523072043300_R8","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1008409","article-title":"Practical considerations for measuring the effective reproductive number, Rt","volume":"16","author":"Gostic","year":"2020","journal-title":"PLoS Comput. Biol."},{"key":"2026010523072043300_R9","doi-asserted-by":"crossref","first-page":"873","DOI":"10.1098\/rsif.2009.0386","article-title":"The construction of next-generation matrices for compartmental epidemic models","volume":"7","author":"Diekmann","year":"2010","journal-title":"J. R. Soc. Interface"},{"key":"2026010523072043300_R10","doi-asserted-by":"crossref","first-page":"599","DOI":"10.1098\/rspb.2006.3754","article-title":"How generation intervals shape the relationship between growth rates and reproductive numbers","volume":"274","author":"Wallinga","year":"2007","journal-title":"Proc. R. Soc. B"},{"key":"2026010523072043300_R11","doi-asserted-by":"crossref","first-page":"2741","DOI":"10.1175\/1520-0493(1999)127<2741:AMCIOT>2.0.CO;2","article-title":"A Monte Carlo implementation of the nonlinear filtering problem to produce ensemble assimilations and forecasts","volume":"127","author":"Anderson","year":"1999","journal-title":"Mon. Weather Rev."},{"key":"2026010523072043300_R12","doi-asserted-by":"crossref","first-page":"72","DOI":"10.1111\/j.1600-0870.2008.00361.x","article-title":"Spatially and temporally varying adaptive covariance inflation for ensemble filters","volume":"61","author":"Anderson","year":"2009","journal-title":"Tellus"},{"key":"2026010523072043300_R13","first-page":"16","article-title":"Data assimilation in the geosciences: an overview of methods, issues, and perspectives","volume":"9","author":"Carrassi","journal-title":"Wiley Interdiscip. Rev."},{"key":"2026010523072043300_R14","doi-asserted-by":"crossref","first-page":"2884","DOI":"10.1175\/1520-0493(2001)129<2884:AEAKFF>2.0.CO;2","article-title":"An ensemble adjustment Kalman filter for data assimilation","volume":"129","author":"Anderson","year":"2001","journal-title":"Mon. Weather Rev."},{"key":"2026010523072043300_R15","unstructured":"Whitaker\n              J\n            , CompoG. 2002An ensemble Kalman smoother for reanalysis. InProc. Symp. on Observations, Data Assimilation and Probabilistic Prediction. Amer. Meteor. Soc, pp. 144\u2013147. Orlando, FL: American Meteorological Society."},{"key":"2026010523072043300_R16","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1009347","article-title":"Improved estimation of time-varying reproduction numbers at low case incidence and between epidemic waves","volume":"17","author":"Parag","year":"2021","journal-title":"PLoS Comput. Biol."},{"key":"2026010523072043300_R17","first-page":"2024","article-title":"Robust uncertainty quantification in popular estimators of the instantaneous reproduction number","author":"Steyn","year":"2024","journal-title":"MedRxiv"},{"key":"2026010523072043300_R18","doi-asserted-by":"crossref","first-page":"616","DOI":"10.1093\/jrsssa\/qnad043","article-title":"A sequential Monte Carlo approach to estimate a time-varying reproduction number in infectious disease models: the Covid-19 case","volume":"186","author":"Storvik","year":"2023","journal-title":"J. R. Stat. Soc. Ser."},{"key":"2026010523072043300_R19","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1008679","article-title":"Estimating the time-varying reproduction number of COVID-19 with a state-space method","volume":"17","author":"Koyama","year":"2021","journal-title":"PLoS Comput. Biol."},{"key":"2026010523072043300_R20","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1009807","article-title":"Bayesian data assimilation for estimating instantaneous reproduction numbers during epidemics: applications to COVID-19","volume":"18","author":"Yang","year":"2022","journal-title":"PLoS Comput. Biol."},{"key":"2026010523072043300_R21","doi-asserted-by":"crossref","first-page":"2745","DOI":"10.1002\/sim.9382","article-title":"Tracking the transmission dynamics of COVID-19 with a time-varying coefficient state-space model","volume":"41","author":"Keller","year":"2022","journal-title":"Stat. Med."},{"key":"2026010523072043300_R22","doi-asserted-by":"crossref","first-page":"115008","DOI":"10.1088\/1361-6420\/ac2cdc","article-title":"Bayesian particle filter algorithm for learning epidemic dynamics","volume":"37","author":"Calvetti","year":"2021","journal-title":"Inverse Probl."},{"key":"2026010523072043300_R23","doi-asserted-by":"crossref","first-page":"105005","DOI":"10.1088\/1361-6420\/aad3e0","article-title":"An approach to periodic, time-varying parameter estimation using nonlinear filtering","volume":"34","author":"Arnold","year":"2018","journal-title":"Inverse Probl."},{"volume-title":"Applied stochastic analysis","year":"2021","author":"Weinan","key":"2026010523072043300_R24"},{"key":"2026010523072043300_R25","doi-asserted-by":"crossref","DOI":"10.1371\/journal.pcbi.1003583","article-title":"Comparison of filtering methods for the modeling and retrospective forecasting of influenza epidemics","volume":"10","author":"Yang","year":"2014","journal-title":"PLoS Comput. Biol."},{"key":"2026010523072043300_R26","doi-asserted-by":"crossref","first-page":"2752","DOI":"10.1073\/pnas.1708856115","article-title":"Forecasting the spatial transmission of influenza in the United States","volume":"115","author":"Pei","year":"2018","journal-title":"Proc. Natl Acad. Sci. USA"},{"key":"2026010523072043300_R27","doi-asserted-by":"crossref","first-page":"489","DOI":"10.1126\/science.abb3221","article-title":"Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2)","volume":"368","author":"Li","year":"2020","journal-title":"Science"},{"key":"2026010523072043300_R28","doi-asserted-by":"crossref","DOI":"10.1016\/j.epidem.2025.100818","article-title":"Modelling COVID-19 in the North American region with a metapopulation network and Kalman filter","volume":"50","author":"Perini","year":"2025","journal-title":"Epidemics"},{"key":"2026010523072043300_R29","doi-asserted-by":"crossref","DOI":"10.1098\/rsos.241275","article-title":"SARS-CoV-2 transmission dynamics in Mozambique and Zimbabwe during the first 3 years of the pandemic","volume":"12","author":"Kaondera-Shava","year":"2025","journal-title":"R. Soc. Open Sci."},{"key":"2026010523072043300_R30","first-page":"1810","article-title":"Inverse problems and data assimilation","author":"Sanz-Alonso","year":"2018","journal-title":"arXiv E-Prints"},{"key":"2026010523072043300_R31","doi-asserted-by":"crossref","first-page":"634","DOI":"10.1175\/1520-0493(2003)131<0634:ALLSFF>2.0.CO;2","article-title":"A local least squares framework for ensemble filtering","volume":"131","author":"Anderson","year":"2003","journal-title":"Mon. Weather Rev."},{"key":"2026010523072043300_R32","doi-asserted-by":"crossref","first-page":"623","DOI":"10.1175\/MWR-D-17-0187.1","article-title":"Enhanced adaptive inflation algorithm for ensemble filters","volume":"146","author":"El Gharamti","year":"2018","journal-title":"Mon. Weather Rev."},{"key":"2026010523072043300_R33","doi-asserted-by":"crossref","first-page":"533","DOI":"10.1016\/S1473-3099(20)30120-1","article-title":"An interactive web-based dashboard to track COVID-19 in real time","volume":"20","author":"Dong","year":"2020","journal-title":"Lancet Infect. Dis."},{"key":"2026010523072043300_R34","unstructured":"Shaman-Lab\n          . 2023shaman-lab\/COVID-19Projection. Seehttps:\/\/github.com\/shaman-lab\/COVID-19Projection."},{"article-title":"epykalman (v1.0.0). Zenodo","year":"2025","author":"Wunrow","key":"2026010523072043300_R35"},{"article-title":"Supplementary material from: Data assimilation for estimating time-varying reproduction numbers. Figshare.","year":"2025","author":"Wunrow","key":"2026010523072043300_R36"}],"container-title":["Journal of the Royal Society Interface"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/royalsocietypublishing.org\/rsif\/article-pdf\/doi\/10.1098\/rsif.2025.0131\/4439579\/rsif.2025.0131.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/royalsocietypublishing.org\/rsif\/article-pdf\/doi\/10.1098\/rsif.2025.0131\/4439579\/rsif.2025.0131.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,1,6]],"date-time":"2026-01-06T04:07:35Z","timestamp":1767672455000},"score":1,"resource":{"primary":{"URL":"https:\/\/royalsocietypublishing.org\/rsif\/article\/22\/233\/20250131\/364698\/Data-assimilation-for-estimating-time-varying"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2025,12,17]]},"references-count":36,"journal-issue":{"issue":"233","published-print":{"date-parts":[[2025,12,17]]}},"URL":"https:\/\/doi.org\/10.1098\/rsif.2025.0131","relation":{},"ISSN":["1742-5662"],"issn-type":[{"type":"electronic","value":"1742-5662"}],"subject":[],"published":{"date-parts":[[2025,12,17]]},"article-number":"20250131"}}