A Note on Observation Processes in Epidemic Models

Ground reality versus model-based computation of basic reproductive numbers in epidemics

Computation of basic reproductive numbers is one of the primary goals of epidemic modelers. There are several challenges in such computations, especially when the data from the virus transmission networks are not so easy to collect; this makes model validation almost impossible. We provide a technical comment on the precautions to be taken while computing model-based basic reproductive numbers so that the ground realities of such computation are maintained. Basic reproductive numbers need to be adjusted retrospectively to compensate for reporting errors within the epidemic spread networks. Such an adjustment would lead to revised pandemic preparedness and mitigation plans.

Basic Reproduction Number of the 2019 Novel Coronavirus Disease in the Major Endemic Areas of China: A Latent Profile Analysis

Objective: The aim of this study is to analyze the latent class of basic reproduction number (R₀) trends of the 2019 novel coronavirus disease (COVID-19) in the major endemic areas of China.

Methods: The provinces that reported more than 500 cases of COVID-19 till February 18, 2020 were selected as the major endemic areas. The Verhulst model was used to fit the growth rate of cumulative confirmed cases. The R₀ of COVID-19 was calculated using the parameters of severe acute respiratory syndrome (SARS) and COVID-19. The latent class of R₀ was analyzed using the latent profile analysis (LPA) model.

Results: The median R₀ calculated from the SARS and COVID-19 parameters were 1.84–3.18 and 1.74–2.91, respectively. The R₀ calculated from the SARS parameters was greater than that calculated from the COVID-19 parameters (Z = −4.782 to −4.623, p < 0.01). Both R₀ can be divided into three latent classes. The initial value of R₀ in class 1 (Shandong Province, Sichuan Province, and Chongqing Municipality) was relatively low and decreased slowly. The initial value of R₀ in class 2 (Anhui Province, Hunan Province, Jiangxi Province, Henan Province, Zhejiang Province, Guangdong Province, and Jiangsu Province) was relatively high and decreased rapidly. Moreover, the initial R₀ value of class 3 (Hubei Province) was in the range between that of classes 1 and 2, but the higher R₀ level lasted longer and decreased slowly.

Conclusion: The results indicated that the overall R₀ trend is decreased with the strengthening of comprehensive prevention and control measures of China for COVID-19, however, there are regional differences.

Analyzing the effects of observation function selection in ensemble Kalman filtering for epidemic models

The Ensemble Kalman Filter (EnKF) is a popular sequential data assimilation method that has been increasingly used for parameter estimation and forecast prediction in epidemiological studies. The observation function plays a critical role in the EnKF framework, connecting the unknown system variables with the observed data. Key differences in observed data and modeling assumptions have led to the use of different observation functions in the epidemic modeling literature. In this work, we present a novel computational analysis demonstrating the effects of observation function selection when using the EnKF for state and parameter estimation in this setting. In examining the use of four epidemiologically-inspired observation functions of different forms in connection with the classic Susceptible-Infectious-Recovered (SIR) model, we show how incorrect observation modeling assumptions (i.e., fitting incidence data with a prevalence model, or neglecting under-reporting) can lead to inaccurate filtering estimates and forecast predictions. Results demonstrate the importance of choosing an observation function that well interprets the available data on the corresponding EnKF estimates in several filtering scenarios, including state estimation with known parameters, and combined state and parameter estimation with both constant and time-varying parameters. Numerical experiments further illustrate how modifying the observation noise covariance matrix in the filter can help to account for uncertainty in the observation function in certain cases.

Addressing the COVID-19 transmission in inner Brazil by a mathematical model

In 2020, the world experienced its very first pandemic of the globalized era. A novel coronavirus, SARS-CoV-2, is the causative agent of severe pneumonia and has rapidly spread through many nations, crashing health systems and leading a large number of people to death. In Brazil, the emergence of local epidemics in major metropolitan areas has always been a concern. In a vast and heterogeneous country, with regional disparities and climate diversity, several factors can modulate the dynamics of COVID-19. What should be the scenario for inner Brazil, and what can we do to control infection transmission in each of these locations? Here, a mathematical model is proposed to simulate disease transmission among individuals in several scenarios, differing by abiotic factors, social-economic factors, and effectiveness of mitigation strategies. The disease control relies on keeping all individuals’ social distancing and detecting, followed by isolating, infected ones. The model reinforces social distancing as the most efficient method to control disease transmission. Moreover, it also shows that improving the detection and isolation of infected individuals can loosen this mitigation strategy. Finally, the effectiveness of control may be different across the country, and understanding it can help set up public health strategies.

Alternative graphical displays for the monitoring of epidemic outbreaks, with application to COVID-19 mortality

BACKGROUND

Classic epidemic curves - counts of daily events or cumulative events over time -emphasise temporal changes in the growth or size of epidemic outbreaks. Like any graph, these curves have limitations: they are impractical for comparisons of large and small outbreaks or of asynchronous outbreaks, and they do not display the relative growth rate of the epidemic. Our aim was to propose two additional graphical displays for the monitoring of epidemic outbreaks that overcome these limitations.

METHODS

The first graph shows the growth of the epidemic as a function of its size; specifically, the logarithm of new cases on a given day, N(t), is plotted against the logarithm of cumulative cases C(t). Logarithm transformations facilitate comparisons of outbreaks of different sizes, and the lack of a time scale overcomes the need to establish a starting time for each outbreak. Notably, on this graph, exponential growth corresponds to a straight line with a slope equal to one. The second graph represents the logarithm of the relative rate of growth of the epidemic over time; specifically, log₁₀(N(t)/C(t-1)) is plotted against time (t) since the 25th event. We applied these methods to daily death counts attributed to COVID-19 in selected countries, reported up to June 5, 2020.

RESULTS

In most countries, the log(N) over log(C) plots showed initially a near-linear increase in COVID-19 deaths, followed by a sharp downturn. They enabled comparisons of small and large outbreaks (e.g., Switzerland vs UK), and identified outbreaks that were still growing at near-exponential rates (e.g., Brazil or India). The plots of log₁₀(N(t)/C(t-1)) over time showed a near-linear decrease (on a log scale) of the relative growth rate of most COVID-19 epidemics, and identified countries in which this decrease failed to set in in the early weeks (e.g., USA) or abated late in the outbreak (e.g., Portugal or Russia).

CONCLUSIONS

The plot of log(N) over log(C) displays simultaneously the growth and size of an epidemic, and allows easy identification of exponential growth. The plot of the logarithm of the relative growth rate over time highlights an essential parameter of epidemic outbreaks.