A Novel Tool for Real-time Estimation of Epidemiological Parameters of Communicable Diseases Using Contact-Tracing Data: Development and Deployment

BACKGROUND: The Surveillance Outbreak Response Management and Analysis System (SORMAS) contains a management module to support countries in their epidemic response. It consists of the documentation, linkage, and follow-up of cases, contacts, and events. To allow SORMAS users to visualize data, compute essential surveillance indicators, and estimate epidemiological parameters from such network data in real-time, we developed the SORMAS Statistics (SORMAS-Stats) application. OBJECTIVE: This study aims to describe the essential visualizations, surveillance indicators, and epidemiological parameters implemented in the SORMAS-Stats application and illustrate the application of SORMAS-Stats in response to the COVID-19 outbreak. METHODS: Based on findings from a rapid review and SORMAS user requests, we included the following visualization and estimation of parameters in SORMAS-Stats: transmission network diagram, serial interval (SI), time-varying reproduction number R(t), dispersion parameter k, and additional surveillance indicators presented in graphs and tables. We estimated SI by fitting lognormal, gamma, and Weibull distributions to the observed distribution of the number of days between symptom onset dates of infector-infectee pairs. We estimated k by fitting a negative binomial distribution to the observed number of infectees per infector. Furthermore, we applied the Markov Chain Monte Carlo approach and estimated R(t) using the incidence data and the observed SI computed from the transmission network data. RESULTS: Using COVID-19 contact-tracing data of confirmed cases reported between July 31 and October 29, 2021, in the Bourgogne-Franche-Comté region of France, we constructed a network diagram containing 63,570 nodes. The network comprises 1.75% (1115/63,570) events, 19.59% (12,452/63,570) case persons, and 78.66% (50,003/63,570) exposed persons, including 1238 infector-infectee pairs and 3860 transmission chains with 24.69% (953/3860) having events as the index infector. The distribution with the best fit to the observed SI data was a lognormal distribution with a mean of 4.30 (95% CI 4.09-4.51) days. We estimated a dispersion parameter k of 21.11 (95% CI 7.57-34.66) and an effective reproduction number R of 0.9 (95% CI 0.58-0.60). The weekly estimated R(t) values ranged from 0.80 to 1.61. CONCLUSIONS: We provide an application for real-time estimation of epidemiological parameters, which is essential for informing outbreak response strategies. The estimates are commensurate with findings from previous studies. The SORMAS-Stats application could greatly assist public health authorities in the regions using SORMAS or similar tools by providing extensive visualizations and computation of surveillance indicators.