Instantaneous R calculation for COVID-19 epidemic in Brazil

Introduction

On December 31, 2019, 27 cases of viral pneumonia were reported in the city of Wuhan, China. A new coronavirus, related to SARS-Cov and MERS-Cov, was isolated from the patients’ airways, being initially named 2019-nCov (Zhu et al., 2020). The rapid spread of this new coronavirus (now called SARS-Cov-2) gave rise to the current pandemic, (Bedford et al., 2020). If there is a significant number of asymptomatic patients transmitting the virus isolation of only the high-risk population will be ineffective (Hellewell et al., 2020). Widespread social restrictions, therefore, are necessary for the effective control of the pandemic (Anastas-sopoulou et al., 2020; Hou et al., 2020).

Outbreaks can be quantitativelly described by the reproduction number (R) (Heesterbeek and Dietz, 1996; Delamater et al., 2019) a common measure of pathogen spread. R values higher than 1 are associated with increasing spread and epidemic phase, whereas decline of spread is seen with values of less than 1. An estimate of COVID-19 R₀ value (basic reproduction number on the initial outbreak) was made using Wuhan data, calculating a number between 2.2 and 3.6 (Li et al., 2020). The first brazillian case was registered on February 26, 2020. As of today (April 16), the total number of confirmed cases is 30,425. It is still unclear what the epidemic dynamics are locally and how it has responded to social restrictions.

Methods

Results

The evaluation of R_i at the first time point (7 days after the first case was registered in each state) showed a variation in the range of 1.55 to 2.81. Median value was 1.96 and median 0.95 quartile was 2.75. This indicated that the epidemic was unabated and probably reflected the basic reproduction number of SARS-Cov-2 (table 1). The value of R_i at 21 days after the first case was registered in each state ranged from 1.09 to 1.89 (median 1.32 and 0.95 quartile median 1.5), indicating that the proliferation rate of the disease was diminishing, and the virus spread was slower (table 2). The EpiEstim package output includes an epidemic curve and a graph of the estimated R_i value, with 95% confidence intervals. The inspection of the R_i graphs could give some insight into the regional differences in COVID-19 infection spread (figures 1-3).

• Download figure
• Open in new tab

Figure 1:

Epidemic curves and R_i estimate curves for some brazillian states: AC - Acre, AL - Alagoas, AM - Amazonas, AP - Amapá, BA - Bahia, CE - Ceará, DF - Distrito Federal, ES - Espírito Santo, GO - Goiás.

• Download figure
• Open in new tab

Figure 2:

Epidemic curves and R_i estimate curves for some brazillian states: MA - Maranhão, MG - Minas Gerais, MS - Mato Grosso do Sul, MT - Mato Grosso, PA - Pará, PB - Paraíba, PE - Pernambuco, PI - Piaúi, PR - Paraná.

• Download figure
• Open in new tab

Figure 3:

Epidemic curves and R_i estimate curves for some brazillian states: RJ - Rio de Janeiro, RN - Rio Grande do Norte, RO - Rondônia, RR - Roraima, RS - Rio Grande do Sul, SC - Santa Catarina, SE - Sergipe, SP - São Paulo, TO - Tocantins.

Discussion

Our results showed an initial R_i compatible with the rapid epidemic growth rate in the beginning of the pandemic spread in Brazil, with most values > 2. However, 2 weeks later, the R_i showed a very impressive decline, reaching < 1.5 for the majority of brazillian states. This probably is a direct consequence of control measures instituted by local governments (social distancing, quarantine, lockdown). Furthermore, the highest R_i initial values correlated with the states that experienced high infection rates, like AM (Amazonas), RJ (Rio de Janeiro), and CE (Ceará).

In this study, we used incidence data to derive R_i and looked at brazillian states specific COVID-19 infection dynamics. In a model incorporating multiple variables, instantaneous R sensitively described real-time shifts of COVID-19 incidence, and varied accordingly with epidemic phase. Superspreading events were associated with R_i with high values, tipically much higher than 2, as well as rapid epidemic growth phases. In contrast, epidemic decline stage was characterized by R_i < 2 (Bandoy and Weimer, 2020).

The data on estimated R_i points to specific differences between different states. Some graphs exhibit spikes, like BA (Bahia) or ES (Espírito Santo) that may correlate with non compliance with control measures. Alternatively, local dynamics could play a role in determining R_i variation between different brazillian regions. Even in those regions that depicted variable behavior and R_i spikes, however, the trend was towards decreasing values and, therefore, less viral spread.

The basic reproduction number R₀ is the expected number of infections caused by an individual in the absence of widespread immunity. Once widespread immunity is achieved, the effective reproduction number R will become lower than R₀ and once R is less than 1, the population is said to have developed herd immunity and the epidemic declines. Immunity can only be obtained with certainty by vaccination, and a vaccine or effective treatment is not expected to be available soon. The best (and only) strategy right now is to rely on social distancing measures until a sustained epidemic suppression (R < 1) could be attained (Ferretti et al., 2020).

Conclusion

Instantaneous R estimating showed to be a convenient way to investigate epidemic dynamics. The COVID-19 pandemic figures in Brazil indicate a trend towards the amellioration of the epidemic. However, this is heavily dependent upon regional governments decisions and politics. There is a high risk of new viral spread if the control measures are precociously eased. Authorities must take this into account and decide wisely upon social restrictions in the next few weeks, because they may be critical in this scenario.