Instantaneous R calculation for COVID-19 epidemic in Brazil =========================================================== * Francisco H. C. Felix * Juvenia Fontenele ## Abstract COVID-19 pandemic represents a major challenge to health systems of all countries. Brazilian regions habe been showing marked differences in onset and number of cases. Health authorities instituted widespread social distancing and lockdown measures but their implementation has also varied. The authors used data on confirmed cases of COVID-19 in Brazil and its states to calculate the value of instantaneous reproduction number at these regions. The results show a reduction of instantaneous reproduction number with time, probably due to social distancing measures put in place in the last weeks by brazillian authorities. It seems logical to maintain restrictions to social contact until the epidemic peak has occurred 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 ### Data source The data were obtained from official sources. The number of cumulative confirmed total cases and deaths in Ceará was obtained from the Ceará State Government information site, IntegraSUS (do Estado do Ceará, 2020). The number of cumulative confirmed total cases and deaths in the state of São Paulo was obtained from the coronavirus information site in the state of São Paulo (do Estado de São Paulo, 2020). The number of cumulative confirmed total cases and deaths in Brazil and its other states was obtained from the COVID-19 brazillian information website (da Saúde do Brasil, 2020). We used the number of cumulative active cases (total minus deaths) to build the epidemic curves. ### Statistical analysis We used EpiEstim, a web application to estimate disease transmissibility during an infectious disease out-break from incidence time series, computing time-dependent reproduction number (R) (Thompson et al., 2019). The Wallinga and Teunis method, as implemented by Ferguson (Wallinga, 2004; Cori et al., 2013) is a likelihood based estimation procedure that captures the temporal pattern of effective reproduction numbers from an observed epidemic curve. Instantaneous R (*R**i*) was estimated within a 5 day time window. Prior mean and standard deviation values for R were set at 3 and 1. Serial interval was estimated using a parametric distribution with uncertainty (offset gamma). We compared the results at two time points (day 7 and day 21 after the first case was registered at each region) from different brazillian states in order to make inferences about the epidemic dynamics. ## 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). View this table: [Table 1:](http://medrxiv.org/content/early/2020/04/29/2020.04.23.20077172/T1) Table 1: Estimated instantaneous R at day 7 after first case of COVID-19 registered on each state of Brazil. The R values correlate well with the regions that would have higher case numbers and faster epidemic growth. View this table: [Table 2:](http://medrxiv.org/content/early/2020/04/29/2020.04.23.20077172/T2) Table 2: Estimated instantaneous R at day 21 after first case of COVID-19 registered on each state of Brazil. It is clear that R values markedly reduced across all country, paralleling the institution of lockdown. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/29/2020.04.23.20077172/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2020/04/29/2020.04.23.20077172/F1) Figure 1: Epidemic curves and Ri 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. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/29/2020.04.23.20077172/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2020/04/29/2020.04.23.20077172/F2) Figure 2: Epidemic curves and Ri 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á. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/04/29/2020.04.23.20077172/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2020/04/29/2020.04.23.20077172/F3) Figure 3: Epidemic curves and Ri 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. ### Suplement Links to data: - Ceará State - [https://indicadores.integrasus.saude.ce.gov.br/api/casos-coronavirus/export-csv](https://indicadores.integrasus.saude.ce.gov.br/api/casos-coronavirus/export-csv) - São Paulo State - [http://www.seade.gov.br/wp-content/uploads/2020/04/Dados-covid-19-estado.csv](http://www.seade.gov.br/wp-content/uploads/2020/04/Dados-covid-19-estado.csv) - Brazil - [https://covid.saude.gov.br/](https://covid.saude.gov.br/) ## 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. ## Data Availability All data used in this research is public. Links were provided. * Received April 23, 2020. * Revision received April 23, 2020. * Accepted April 29, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. C Anastassopoulou, L Russo, A Tsakris, and C Siettos. Data-based analysis, modelling and forecasting of the COVID-19 outbreak. PLoS One, 15:e0230405, 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0230405&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) 2. DJ Darwin R Bandoy and Bart C Weimer. Pandemic dynamics of COVID-19 using epidemic stage instantaneous reproductive number and pathogen genome identity (GENI) score: modeling molecular epidemiology. mar 2020. doi: 10.1101/2020.03.17.20037481. URL [https://doi.org/10.1101%2F2020.03.17.20037481](https://doi.org/10.1101%2F2020.03.17.20037481). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wMy4xNy4yMDAzNzQ4MXYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDQvMjkvMjAyMC4wNC4yMy4yMDA3NzE3Mi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. J Bedford, D Enria, J Giesecke, DL Heymann, C Ihekweazu, G Kobinger, HC Lane, Z Memish, MD Oh., AA Sall, A Schuchat, K Ungchusak, and LH Wieler. COVID-19: towards controlling of a pandemic. Lancet, 395:1015–1018, Mar 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(20)30673-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) 4. Anne Cori, Neil M. Ferguson, Christophe Fraser, and Simon Cauchemez. A New Framework and Software to Estimate Time-Varying Reproduction Numbers During Epidemics. American Journal of Epidemiology, 178 (9):1505–1512, sep 2013. doi: 10.1093/aje/kwt133. URL [https://doi.org/10.1093%2Faje%2Fkwt133](https://doi.org/10.1093%2Faje%2Fkwt133). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwt133&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24043437&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) 5. Ministério da Saúde do Brasil. Painel de casos de doença pelo coronavírus 2019 (COVID-19) no Brasil pelo Ministério da Saúde. Retrieved from: ¡ [https://covid.saude.gov.br/](https://covid.saude.gov.br/)¿ Access on: 16 abr. 2020, 2020. URL [https://covid.saude.gov.br/](https://covid.saude.gov.br/). 6. Paul L. Delamater, Erica J. Street, Timothy F. Leslie, Y. Tony Yang, and Kathryn H. Jacobsen. Complexity of the Basic Reproduction Number (R0). Emerging Infectious Diseases, 25(1):1–4, jan 2019. doi: 10.3201/eid2501.171901. URL [https://doi.org/10.3201%2Feid2501.171901](https://doi.org/10.3201%2Feid2501.171901). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2501.171901&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30560777&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) 7. Governo do Estado de São Paulo. Coronavírus. Casos em São Paulo. Disponível em: ¡[https://www.seade.gov.br/coronavirus/](https://www.seade.gov.br/coronavirus/)¿ Acesso em: 05 abr. 2020, 2020. URL [https://www.seade.gov.br/coronavirus/](https://www.seade.gov.br/coronavirus/). 8. Governo do Estado do Ceará. Boletim epidemiológico novo coronavírus (COVID-19). Disponível em: ¡[https://indicadores.integrasus.saude.ce.gov.br/indicadores/indicadores-coronavirus/coronavirus-ceara](https://indicadores.integrasus.saude.ce.gov.br/indicadores/indicadores-coronavirus/coronavirus-ceara)¿ Acesso em: 05 abr. 2020, 2020. URL [https://indicadores.integrasus.saude.ce.gov.br/indicadores/indicadores-coronavirus/coronavirus-ceara](https://indicadores.integrasus.saude.ce.gov.br/indicadores/indicadores-coronavirus/coronavirus-ceara). 9. Luca Ferretti, Chris Wymant, Michelle Kendall, Lele Zhao, Anel Nurtay, Lucie Abeler-Dörner, Michael Parker, David Bonsall, and Christophe Fraser. Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing. Science, page eabb6936, mar 2020. doi: 10.1126/science.abb6936. URL [https://doi.org/10.1126%2Fscience.abb6936](https://doi.org/10.1126%2Fscience.abb6936). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjgvNjQ5MS9lYWJiNjkzNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzA0LzI5LzIwMjAuMDQuMjMuMjAwNzcxNzIuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 10. J. A. P. Heesterbeek and K. Dietz. The concept of Roin epidemic theory. Statistica Neerlandica, 50(1):89–110, mar 1996. doi: 10.1111/j.1467-9574.1996.tb01482.x. URL [https://doi.org/10.1111%2Fj.1467-9574.1996.tb01482.x](https://doi.org/10.1111%2Fj.1467-9574.1996.tb01482.x). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1467-9574.1996.tb01482.x&link_type=DOI) 11. J Hellewell, S Abbott, A Gimma, NI Bosse, CI Jarvis, TW Russell, JD Munday, AJ Kucharski, WJ Edmunds, S Funk, and RM Eggo. Feasibility of controlling COVID-19 outbreaks by isolation of cases and contacts. Lancet Glob Health, 8:e488–e496, Apr 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2214-109X(20)30074-7&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32119825&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) 12. Can Hou, Jiaxin Chen, Yaqing Zhou, Lei Hua, Jinxia Yuan, Shu He, Yi Guo, Sheng Zhang, Qiaowei Jia, Chenhui Zhao, Jing Zhang, Guangxu Xu, and Enzhi Jia. The effectiveness of the quarantine of Wuhan city against the Corona Virus Disease 2019 (COVID-19): well-mixed SEIR model analysis. Journal of Medical Virology, apr 2020. doi: 10.1002/jmv.25827. URL [https://doi.org/10.1002%2Fjmv.25827](https://doi.org/10.1002%2Fjmv.25827). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmv.25827&link_type=DOI) 13. Qun Li, Xuhua Guan, Peng Wu, Xiaoye Wang, Lei Zhou, Yeqing Tong, Ruiqi Ren, Kathy S.M. Leung, Eric H.Y. Lau, Jessica Y. Wong, Xuesen Xing, Nijuan Xiang, Yang Wu, Chao Li, Qi Chen, Dan Li, Tian Liu, Jing Zhao, Man Liu, Wenxiao Tu, Chuding Chen, Lianmei Jin, Rui Yang, Qi Wang, Suhua Zhou, Rui Wang, Hui Liu, Yinbo Luo, Yuan Liu, Ge Shao, Huan Li, Zhongfa Tao, Yang Yang, Zhiqiang Deng, Boxi Liu, Zhitao Ma, Yanping Zhang, Guoqing Shi, Tommy T.Y. Lam, Joseph T. Wu, George F. Gao, Benjamin J. Cowling, Bo Yang, Gabriel M. Leung, and Zijian Feng. Early Transmission Dynamics in Wuhan China, of Novel Coronavirus–Infected Pneumonia. New England Journal of Medicine, 382(13): 1199–1207, mar 2020. doi: 10.1056/nejmoa2001316. URL [https://doi.org/10.1056%2Fnejmoa2001316](https://doi.org/10.1056%2Fnejmoa2001316). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/nejmoa2001316&link_type=DOI) 14. RN Thompson, JE Stockwin, Gaalen RD van, JA Polonsky, ZN Kamvar, PA Demarsh, E Dahlqwist, S Li, E Miguel, T Jombart, J Lessler, S Cauchemez, and A Cori. Improved inference of time-varying reproduction numbers during infectious disease outbreaks. Epidemics, 29:100356, Dec 2019. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.epidem.2019.100356&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) 15. J. Wallinga. Different Epidemic Curves for Severe Acute Respiratory Syndrome Reveal Similar Impacts of Control Measures. American Journal of Epidemiology, 160(6):509–516, sep 2004. doi: 10.1093/aje/kwh255. URL [https://doi.org/10.1093%2Faje%2Fkwh255](https://doi.org/10.1093%2Faje%2Fkwh255). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/aje/kwh255&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15353409&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223938000001&link_type=ISI) 16. N Zhu, D Zhang, W Wang, X Li, B Yang, J Song, X Zhao, B Huang, W Shi, R Lu, P Niu, F Zhan, X Ma, D Wang, W Xu, G Wu, GF Gao, and W Tan. A Novel Coronavirus from Patients with Pneumonia in China, 2019. N Engl J Med, 382:727–733, Feb 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa2001017&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F04%2F29%2F2020.04.23.20077172.atom)