Projection of COVID-19 intensive care hospitalizations in the Federal District, Brazil: an analysis of the impact of social distancing measures

Ivan Zimmermann Mauro Sanchez Jonas Brant Domingos Alves About the authors

Abstract

Objective:

To build scenarios and analyze the impact of social distancing policies on the spread of COVID-19 and the need for intensive care unit beds.

Methods:

Three dissemination scenarios were built according to level of adherence to social distancing measures in the context of Brazil's Federal District, based on a dynamic transition compartmental model and Monte Carlo simulations. The model's parameter values were based on official sources, indexed bibliographic databases and public data repositories.

Results:

The favorable scenario, with constant 58% adherence to social distancing, estimated a peak of 189 (interquartile range [IQR]: 57 – 394) ICU hospitalizations on March 3rd2021. Absence of social distancing would result in an unfavorable scenario with a peak of 6,214 (IQR: 4,618 – 8,415) ICU hospitalizations probably as soon as July 14th2020.

Conclusion:

The projections indicate the high impact of social distancing measures and emphasize the applicability of public indicators for COVID-19 monitoring.

Keywords:
Coronavirus Infections; Intensive Care Units; Bed Occupancy; Health Evaluation; Public Policy

Introduction

COVID-19 is a disease resulting from Sars-CoV-2 virus infection. Its clinical manifestations can include severe acute respiratory syndrome (SARS), among others. 11. Mohammadi M, Meskini M, Nascimento Pinto AL. 2019 Novel coronavirus (COVID-19) overview. Z Gesundh Wiss [Internet]. 2020 Apr [cited 2020 Sep 25];1-9. Available from: https://doi.org/10.1007/s10389-020-01258-3
https://doi.org/10.1007/s10389-020-01258...
The first case in Brazil was confirmed on 25/2/2020; by 3/5/2020 the country had officially registered over 100,000 cases and 7,000 deaths. 22. Ministério da Saúde (BR). Painel coronavírus [Internet]. Brasília: Ministério da Saúde; 2020 [citado 2020 set 25]. Disponível em: https://covid.saude.gov.br/
https://covid.saude.gov.br/...
Given its far-reaching nature and ability to spread rapidly, its potential to bring the hospital network to a collapse, in particular insufficient availability of intensive care unit (ICU) beds, COVID-19 is considered to be a very serious health condition for Public Health. 11. Mohammadi M, Meskini M, Nascimento Pinto AL. 2019 Novel coronavirus (COVID-19) overview. Z Gesundh Wiss [Internet]. 2020 Apr [cited 2020 Sep 25];1-9. Available from: https://doi.org/10.1007/s10389-020-01258-3
https://doi.org/10.1007/s10389-020-01258...

Although there are reports of experiments with treatment and initiatives to find a vaccine, at the time of writing this manuscript, there was no robust proof of the efficacy of pharmaceutical interventions in preventing or treating this clinical condition. 33. Rothan HA, Byrareddy SN. The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak. J Autoimmun [Internet]. 2020 May [cited 2020 Sep 25];109:102433. Available from: https://doi.org/10.1016/j.jaut.2020.102433
https://doi.org/10.1016/j.jaut.2020.1024...
Ample testing, personnel protection (health services and essential services), isolation of infected individuals and adoption of non-pharmaceutical interventions, such as public policies on social distancing, appear as the most feasible option for addressing COVID-19. 44. Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health [Internet]. 2020 May [cited 2020 Sep 25];5(5):e261-70. Available from: https://doi.org/10.1016/s2468-2667(20)30073-6
https://doi.org/10.1016/s2468-2667(20)30...
On 11/3/2020, the first local social distancing measures came into force throughout Brazil's Federal District in view of the risk of novel coronavirus infection, including the suspension of events with audiences of more than 100 people and educational activities at all schools, universities and faculties. 55. Governo do Distrito Federal. Decreto n° 40.509, de 11 de março de 2020. Dispõe sobre as medidas para enfrentamento da emergência de saúde pública de importância internacional decorrente do novo coronavírus, e dá outras providências [Internet]. Diário Oficial do Distrito Federal, Brasília (DF), 2020 mar 11 [citado 2020 set 25];Edição extra. Disponível em: http://www.sinj.df.gov.br/sinj/Norma/ad0fae78af5f4e50b46c7357b7ee8597/Decreto_40509_11_03_2020.html
http://www.sinj.df.gov.br/sinj/Norma/ad0...
Moreover, bars and restaurants had to ensure a minimum distance of two meters between tables. On 21/4/2020, in addition to the measures already taken, the Federal District government began testing symptomatic people, increasing sample collection and screening points exclusively for these people in places other than health centers.

With the aim of monitoring and informing public policy decisions on COVID-19, a series of initiatives and networks have joined forces and made available case prediction tools including projections of the need for hospital beds and equipment. 66. Rede CoVida. Os impactos das medidas de distanciamento social e redução de fluxo intermunicipal na Bahia. Bol CoVida [Internet]. 2020 abr [citado 2020 set 25];1(2):19. Disponível em: https://covid19br.org/main-site-covida/wp-content/uploads/2020/04/BoletimCovida_2ed_Final.pdf
https://covid19br.org/main-site-covida/w...

7. COVID-19 Brasil. Monitoramento e análises da situação do Coronavírus no Brasil [Internet]. São Paulo; 2020 [citado 2020 set 25]. Disponível em: https://ciis.fmrp.usp.br/covid19/
https://ciis.fmrp.usp.br/covid19/...
-88. Organização Pan-Americana de Saúde - OPAS. Ministério da Saúde (BR). Universidade de Brasília. Butantã. Pressão hospitalar por COVID-19 [Internet]. [S.l];2020. Disponível em: https://covid-calc.org/
https://covid-calc.org/...
In this respect, analysis of data on public policies, with regard to evaluation of their impact, is a field of possible inferences, for which robust data are essential. In scenarios with temporal limitation of data availability, building models and counterfactuals (possible but not necessarily observed scenarios) is of great value for identifying trends and projections of possible impacts of a given public policy. 99. Presidência da República (BR). Casa Civil. Avaliação de políticas públicas: guia prático de análise ex post [Internet]. Brasília: Presidência da República; 2018 [citado 2020 set 25]. 301 p. Disponível em: https://www.ipea.gov.br/portal/index.php?option=com_content&view=article&id=34504
https://www.ipea.gov.br/portal/index.php...

Based on the context of the occurrence of COVID-19 cases among the Federal District population, the objective of this study was to build scenarios to analyze the impact of social distancing policies on the spread of COVID-19 and the need for ICU beds in the Federal District, with the aid of a mathematical model and Monte Carlo simulations.

Methods

Taking as our reference the bases and designs of studies in the field of public policy evaluation, 99. Presidência da República (BR). Casa Civil. Avaliação de políticas públicas: guia prático de análise ex post [Internet]. Brasília: Presidência da República; 2018 [citado 2020 set 25]. 301 p. Disponível em: https://www.ipea.gov.br/portal/index.php?option=com_content&view=article&id=34504
https://www.ipea.gov.br/portal/index.php...
this is a case study based on simulations. Following data and indicator acquisition, we built scenarios of projections and counterfactuals relevant for analysis of the impact of social distancing on the demand for ICU beds in the context observed.

The case used for this analysis is comprised of all the administrative regions within the Federal District, the total population of which is estimated as being 3,223,048 inhabitants in 2020, according to projections made by the Brazilian Institute of Geography and Statistics (IBGE). 1010. Instituto Brasileiro de Geografia e Estatística - IBGE. Diretoria de Pesquisas. Coordenação de População e Indicadores Sociais. Gerência de Estudos e Análises da Dinâmica Demográfica. Projeção da população do Brasil e Unidades da Federação por sexo e idade para o período 2000-2030 [Internet]. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística; 2020 [citado 2020 set 25]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?ibge/cnv/projpopuf.def
http://tabnet.datasus.gov.br/cgi/tabcgi....
With the aim of obtaining a reference as to existing hospital capacity, we consulted the database of the National Registry of Health Establishments as at June 2020 to find the number of adult ICU beds in the Federal District. In addition, with the aim of contextualizing the adherence of the population to the social distancing measures recommended by the policies in force, we consulted the public data contained in the Federal District ‘social isolation index’ for the period 1/3/2020 to 31/7/2020, made available by the Inloco technology company. 1111. Inloco. Mapa brasileiro da COVID-19 [Internet]. São Paulo: Inloco; 2020 [citado 27 abr 2020]. Disponível em: https://www.inloco.com.br/pt/covid-19
https://www.inloco.com.br/pt/covid-19...
This index has assisted initiatives monitoring the spread of COVID-19, undertaken by public universities and governmental institutions. By using telephone device geolocation data for part of the local population - with prior authorization -, it is possible to estimate on a daily basis the percentage of that population which has adhered to public policy advice on social isolation. 1111. Inloco. Mapa brasileiro da COVID-19 [Internet]. São Paulo: Inloco; 2020 [citado 27 abr 2020]. Disponível em: https://www.inloco.com.br/pt/covid-19
https://www.inloco.com.br/pt/covid-19...

Model structure

A mathematical model was built in order to estimate the need for ICU beds as a result of the spread of COVID-19. Taking the methodological guidelines published by the Health Ministry's Science, Technology and Strategic Supplies Secretariat 1212. Ministério da Saúde (BR). Secretaria de Ciência Tecnologia e Insumos Estratégicos. Departamento de Ciência e Tecnologia. Diretrizes metodológicas: diretriz de avaliação econômica [Internet]. Brasília: Ministério da Saúde; 2014 [citado 2020 set 25]. 132 p. Disponível em: http://bvsms.saude.gov.br/bvs/publicacoes/diretrizes_metodologicas_diretriz_avaliacao_economica.pdf
http://bvsms.saude.gov.br/bvs/publicacoe...
and the practical bases for modeling communicable diseases, 1313. Vynnycky E, White RG, Fine P. An introduction to infectious disease modelling. Oxford: Oxford University Press; 2010. the concept used for the model was a dynamic transition compartmental model. In this model, derived from the classic SIR model (‘susceptible, infected and recovered’) proposed by Kermack & McKendrick, transition between compartments (or states of health) is represented by a system of dependent differential equations, among other parameters, containing the number of infected people ( I ), transmission coefficient ( β ), latency ( α ) and recovery ( γ ).

The transition and compartment structure used by this study was based on models already provided by other authors who also sought to project hospital demand, 1414. Rocha Filho TM, Santos FSG, Gomes VB, Rocha TAH, Croda JHR, Ramalho WM, et al. Expected impact of COVID-19 outbreak in a major metropolitan area in Brazil. medRxiv Prepr [Internet]. 2020 Mar [cited 2020 sep 25]. Available from: https://doi.org/10.1101/2020.03.14.20035873
https://doi.org/10.1101/2020.03.14.20035...

15. Hill A. Modeling COVID-19 spread vs healthcare capacity [Internet]. [S.l.]; 2020 [cited 2020 Sep 25]. Available from: https://alhill.shinyapps.io/COVID19seir/
https://alhill.shinyapps.io/COVID19seir/...
-1616. Penn University. CHIME model: discrete-time SIR modeling of infections/recovery [Internet]. [S.l.]; 2020 [cited 2020 Sep 25]. Available from: https://code-for-philly.gitbook.io/chime/what-is-chime/sir-modeling
https://code-for-philly.gitbook.io/chime...
as well as the impact of non-pharmaceutical interventions - such as use of facemasks - on community transmission of COVID-19. 1717. Eikenberry SE, Mancuso M, Iboi E, Phan T, Eikenberry K, Kuang Y, et al. To mask or not to mask: modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic. Infect Dis Model [Internet]. 2020. In this proposal, the model for the daily correlation between the number of susceptible ( S ), exposed ( E ), infected ( I ) and removed ( R ) individuals - SEIR - is represented in Figure 1 .

Figure 1
Representation of the structure of the compartmental dynamic model used

Within the infected compartment, there are also subdivisions of those who do not need to be hospitalized ( I1 ), those who progress to needing hospitalization ( I2 ) and patients in a critical state ( I3 ) who need to be hospitalized in an ICU.

Main premises of the model

  1. No change in the total size of the population is expected (overall birth and mortality rates), except due to death from COVID-19. Given the short time period analyzed, this restriction was considered to simplify the model with little impact on the results.

  2. Similarly to previous approaches, based on COVID-19 SEIR models, 1717. Eikenberry SE, Mancuso M, Iboi E, Phan T, Eikenberry K, Kuang Y, et al. To mask or not to mask: modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic. Infect Dis Model [Internet]. 2020. it is assumed that hospitalized individuals are not effectively exposed to the general population and, therefore, do not contribute to the calculation of the force of infection in the community.

  3. The possibility of a recovered individual being reinfected is not foreseen, even though there is no conclusive evidence as to the extent to which acquired immunity following infection protects against reinfection. 1818. World Health Organization - WHO. “Immunity passports” in the context of COVID-19 [Internet]. Genebra: World Health Organization; 2020 [cited 2020 Sep 25]. Available from: https://www.who.int/news-room/commentaries/detail/immunity-passports-in-the-context-of-covid-19
    https://www.who.int/news-room/commentari...

  4. As it is not an objective of the analysis, this model does not take into consideration potential increase in mortality due to lack of hospital beds. It is important to note that lethality found by the model may be underestimated in a scenario of lack of hospital beds.

  5. The estimates of effective contact assume random mixing, 1313. Vynnycky E, White RG, Fine P. An introduction to infectious disease modelling. Oxford: Oxford University Press; 2010. whereby each infected individual may have contact with a susceptible individual according to the same average probability.

The complete equations and further details of the model are provided in Supplementary Material 1.

Scenarios

Three scenarios of the spread of COVID-19 were built with the aim of analyzing the potential impact of the disease on the need for hospital care. The scenarios were based on a one-year time period which was considered to be sufficient to take into account both a complete wave of infections and also uncertainty as to availability or otherwise of specific interventions (vaccines and treatment). The start date was 22/3/2020, providing a ten-day margin between the beginning of the social distancing measures and the peak of social isolation levels obtained in the Federal District.

Scenario A - COVID-19 spread profile with no adherence to social distancing measures

Definition

This scenario represents the contrafactual that would occur if no social distancing measures were taken or adhered to by the population.

There is therefore no reduction in social contact over the course of the model, taking an average estimate of social isolation observed in the last ten days before the local measures were brought into force on 11/3/2020.

The model therefore uses a constant transmission coefficient ( β ) value.

Scenario B - COVID-19 spread profile maintaining adherence to social distancing

Definition

This scenario represents the profile of the spread of infection and the profile of need for hospital beds maintaining the same level of average social isolation obtained during the first ten days after the peak in social isolation.

In this scenario, the model also uses a constant transmission coefficient ( β ) value.

Scenario C - COVID-19 spread profile based on the observed dynamics of a falling level of social distancing

Definition

This scenario represents the dynamics of social isolation observed in the Federal District and the possible projection of the need for hospital beds when maintaining the last average levels of isolation observed.

Seeking to simulate an interruption or return of distancing measures, this context also estimates population isolation targets according to demand for ICU beds.

As mentioned in the definitions presented, the scenarios did not include a complete simulation of adherence to distancing above the average observed following the implementation of distancing polices in the Federal District. Despite this, we emphasize that the analysis of the isolation measures in Scenario C was not subject to this restriction.

Parameter values and uncertainties

With the aim of bringing the analysis model closer to local reality and to the size of the population, we used the Federal District age structure to adjust the expected indicators of lethality and hospitalization for the region. Federal District Health Department official bulletins 1919 18. Governo do Distrito Federal. Secretaria de Saúde. Boletins Informativos DIVEP/CIEVES (COE). Boletins informativos sobre coronavirus (COVID-19). Brasília: GDF; 2020. Disponível em: http://www.saude.df.gov.br/boletinsinformativos-divep-cieves/
http://www.saude.df.gov.br/boletinsinfor...
were used to obtain numbers of cases by date of symptom onset, numbers of ICU hospitalizations and COVID-19 deaths in the Federal District. In addition to the estimates of infection and lethality by age group for the Federal District (Supplementary Material 2), the complete series of ICU hospitalizations, deaths and percentage social isolation available for the period analyzed (data from 1/3/2020 to 31/7/2020) is shown in Supplementary Material 3. In addition to the parameters mentioned, the daily evolution of the instantaneous reproduction number (Rt) was estimated for each scenario.

Given the absence of local data for the remaining parameters, we also consulted publications indexed on Pubmed ( http://www.ncbi.nlm.nih.gov/pubmed ), preprints on MedRvix ( www.medrxiv.org ) and public repositories of modeling parameters, such as the initiative of the MIDAS network ( http://midasnetwork.us/COVID-19 ), using an adaptation of the basic search structure with the terms (COVID-19 OR coronavirus) AND (model OR modelling). The choice of the data sources of available parameters prioritized those used in models validated by other institutions considered to be references on the subject, such as the Imperial College London models, seeking comparability of the model proposed here. The retrieved data included point estimations and their parametric uncertainties (standard error, confidence intervals or probability distributions). The values of the main parameters considered in model proposed here, along with their parametric uncertainties, description and sources, are shown in Table 1 .

Table 1
List of values, parametric uncertainties, distributions and sources of the main parameters of the model of the spread of COVID-19 in the Federal District, Brazil

In order to take into consideration the impact of the parametric uncertainties of the nine parameters described in Table 1 on the model results, the analyses were performed using a stochastic (probabilistic) approach, using Monte Carlo simulations, with multiple draws on expected value levels (2ndorder). 2020 19. Briggs A, Claxton K, Sculpher M. Decision modelling for health economic evaluation. Oxford: Oxford University Press; 2011. 256 p. Based on the characteristics of the variables, a priori probability distributions were adjusted for the following parameters: incubation time (days), duration of mild clinical picture (days), fraction of individuals with a severe clinical picture (%), time hospitalized in a ward (days), fraction of individuals who progressed to ICU hospitalization UTI (%), time hospitalized in an ICU (days), lethality (%) and time from onset of symptoms to death (days), as illustrated in Supplementary Material 4. In the absence of intervals in the data sources, we used bootstrap resampling to estimate parametric uncertainty (confidence intervals). The results were summarized into medians and coverage amplitudes of the simulations (interquartile range [IQR] and confidence interval [CI]).

Validation

In order to calibrate the model, we used as a reference the scenario that reflects the spread of COVID-19 during the observed levels of social distancing (Scenario C). We recorded the values of the minimization statistic of the root-mean-square error (RMSE) of the number of deaths predicted by the model, as observed in the time series under analysis. The adequacy of the number of Monte Carlo simulations was assessed by visual inspection of the convergence of the main parameters of the model.

All the analyses were performed with the aid of Microsoft Excel® and the Visual Basic for Applications (VBA) programming language. The adjustment of the calibration parameters was performed using the non-linear least squares method, with the support of the non-linear Generalized Reduced Gradient (GRG) optimization algorithm, using the SOLVER function. The complete model and the history of its versions are available at a public data repository. 2121 20. Zimmermann I. Modelo de projeção da demanda por leitos de UTI por COVID-19. Mendeley Data [Internet]. 2020 set [citado 2020 set 25]. Disponível em: http://dx.doi.org/10.17632/bypzkmwnrv
http://dx.doi.org/10.17632/bypzkmwnrv...

Results

According to the official data consulted, between 1/3/2020 and the final date of data acquisition for the model (31/7/2020), 106,292 COVID-19 cases and 1,469 deaths had been confirmed in the Federal District. 1919 18. Governo do Distrito Federal. Secretaria de Saúde. Boletins Informativos DIVEP/CIEVES (COE). Boletins informativos sobre coronavirus (COVID-19). Brasília: GDF; 2020. Disponível em: http://www.saude.df.gov.br/boletinsinformativos-divep-cieves/
http://www.saude.df.gov.br/boletinsinfor...
Bearing in mind the date of the first local decree on distancing measures, there was important initial adherence to the measures by the population, given the gradual increase in the rate of social isolation over the days following the decree, from an average of 30.0% between 2/3/2020 and 11/3/2020, to an average of 57.4% between 22/3/2020 and 30/4/2020 ( Figure 2 ).

Figure 2
Evolution of percentage social isolation in relation to the implementation date of COVID-19 combat measures in the Federal District, Brazil

Considering the records of cases by date of symptom onset, all the scenarios began simulations on 22/3/2020, with an average estimate of 366 infected (I1) individuals and 732 exposed (E) individuals - given the existence of a minimum instantaneous reproduction number (Rt) of 2 without social isolation. With an RMSE of 83.5 in the predicted death values, the value of parameter τ was estimated as 0.6979 for calibrating the model. According to the National Registry of Health Establishments database, as at June the Federal District had a total of 1534 adult ICU beds, 378 of which were public. Below are the results of each of the scenarios.

Scenario A - COVID-19 spread profile with no adherence to social distancing measures

In the first scenario, the transmission coefficient is calculated assuming the level of social isolation to be 30.0%, this being coherent with the average values observed up until 11/3/2020. As illustrated in Figure 3 , given the need to hospitalize critical cases, if there were no adherence to social distancing by the population, according to the median of the simulations, there would be a peak of 6,214 people needing to be hospitalized in ICUs (IQR: 4,618 – 8,415) probably on 14/7/2020 (IQR: 9/7/2020 – 20/7/2020). However, according to the National Registry of Health Establishments data, the total number of public ICU beds registered would already have been surpassed on 17/5/2020 (IQR: 12/5/2020 – 22/5/2020). This demand, greater than the number of registered public ICU beds, would continue for at least 125 days (IQR: 115 – 135).

Figure 3
Projections of the demand for ICU beds due to COVID-19 and the instantaneous reproduction number (Rt) in Scenario (A) COVID-19 spread profile with no adherence to social distancing measures, Scenario (B) COVID-19 spread profile maintaining adherence to social distancing and Scenario (C) COVID-19 spread profile based on the observed dynamics of a falling level of social distancing, in the Federal District, Brazil

Over the period of a year, there would probably be a peak of 178,615 active infections (IQR: 153,820 – 203,955) on 1/7/2020 (IQR: 26/6/2020 – 6/7/2020). At this rate, by 31/3 of the following year (2021), up to 73.2% (IQR: 70.3 – 75.9%) of the Federal District population would have been infected. Even if it were possible to provide care for all infected people, 16,143 individuals would have died (IQR: 10,948 – 22,326) - expected lethality would be 0.67% (IQR: 0.45 – 0.94%).

Scenario B - COVID-19 spread profile maintaining adherence to social distancing

Considering the hypothesis of maintaining the average level of social isolation at 57.38%, this being the level reached in the initial period of the distancing measures, this scenario projected a peak of 189 ICU hospitalizations (IQR: 57 – 394) only on 7/3/2021 (IQR: 6/1/2021 – 31/1/2021) and, in this context, it would be feasible to meet the demand with the number of beds registered in the Federal District.

Over the one-year simulated period, the total number of infected individuals could reach a peak of 5,269 active infections (IQR: 1,538 – 11,549). By 31/3/2021, the cumulative total of infections would affect only 7.6% of the Federal District population DF (IQR: 3.0 – 16.4%) and, with 0.70% lethality (IQR: 0.47 – 0.94%), there would be up to 1,438 deaths (IQR: 583 – 3,234).

Scenario C - COVID-19 spread profile based on the observed dynamics of a falling level of social distancing

Using the proposed methodology described for this scenario, the transmission coefficient followed the dynamic variation of the level of isolation recorded as at 31/7/2020. With effect from 1/8/2020, the average level of 39.7% isolation was assumed to continue, based on the levels recorded on the last five days available. 1111. Inloco. Mapa brasileiro da COVID-19 [Internet]. São Paulo: Inloco; 2020 [citado 27 abr 2020]. Disponível em: https://www.inloco.com.br/pt/covid-19
https://www.inloco.com.br/pt/covid-19...
As such, the fall observed in the level of distancing was projected to lead to a peak of up to 3,460 ICU hospitalizations (IQR: 2,634 – 4,779) probably on 10/9/2020 (IQR: 31/8/2020 – 22/9/2020) ( Figure 3 ). In this scenario, even raising the level of isolation to an average of 70% (a level close to the maximum level recorded on 22/3/2020) with effect from 1/8/2020, there would still be a peak of 2,097 hospitalizations (IQR: 1.245 – 3.312) on 19/8/2020 (IQR: 16/8/2020 – 21/8/2020). On the other hand, complete interruption of distancing, returning to the 30% level of isolation observed before the distancing policies were implemented, would lead to a peak of 4,812 ICU hospitalizations (IQR: 3,643 – 6,283) on 12/9/2020 (IQR: 04/9/2020 – 21/9/2020).

Discussion

Based on public indicators and data used in this analysis, it is possible to identify adherence to social distancing following the publication of the first decree establishing measures to address COVID-19 in the Federal District. Likewise, linking the temporal evolution of the level of isolation with the consequent progressive reduction in the spread of the disease (Rt), it would be possible to reach reference values for the control of outbreaks and, at the same time, delay and reduce the peak of ICU hospitalizations to a plausible estimate in relation to the capacity of the local healthcare network.

Differently to the scenario in which social distancing was maintained, both the other scenarios with absence of social distancing, whether because of no measures being implemented (Scenario A), or because of low adherence (Scenario C), provide a projection compatible with the collapse of the local healthcare network due to high demand for ICU hospitalizations. Notwithstanding the uncertainty about the data and the magnitude of demand, the alternative projections of Scenario C, which seeks to incorporate the observed percentages of social isolation, pointed to peaks of hospitalization in August and September of this year. These projections are consistent with the ICU bed waiting list situation in the Federal District reported in August. Also coherent with the findings of a recent study of other Brazilian states carried out by researchers from Imperial College London, 2222 21. Mellan TA, Hoeltgebaum HH, Mishra S, Whittaker C, Schnekenberg RP, Gandy A, et al. Report 21: estimating COVID-19 cases and reproduction number in Brazil. Imperial College London [Internet]. 2020 May [cited 2020 Sep 25]. Available from: https://doi.org/10.25561/78872
https://doi.org/10.25561/78872...
notwithstanding significant reductions in the spread of cases, the levels of isolation achieved would perpetuate an Rt >1, indicating the continuation of epidemic growth of COVID-19 in the Federal District.

Despite this analysis being restricted to the context of the Federal District, it shows the magnitude of the potential effect of distancing on the need for hospital beds, it being unlikely that this relationship, even with a different magnitude, cannot be generalized to other contexts. The results are also consistent with other initiatives analyzing social distancing, 44. Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health [Internet]. 2020 May [cited 2020 Sep 25];5(5):e261-70. Available from: https://doi.org/10.1016/s2468-2667(20)30073-6
https://doi.org/10.1016/s2468-2667(20)30...
including national contexts such as the analysis by Rede CoVida of the positive effects of social distancing and the reduction in movement of people between municipalities on the control of COVID-19 in the state of Bahia, 66. Rede CoVida. Os impactos das medidas de distanciamento social e redução de fluxo intermunicipal na Bahia. Bol CoVida [Internet]. 2020 abr [citado 2020 set 25];1(2):19. Disponível em: https://covid19br.org/main-site-covida/wp-content/uploads/2020/04/BoletimCovida_2ed_Final.pdf
https://covid19br.org/main-site-covida/w...
the available data of which strongly suggest that the measures taken managed to flatten the epidemic curve in that state.

The model adopted here follows methodological lines compatible with the population dynamics of infections 1313. Vynnycky E, White RG, Fine P. An introduction to infectious disease modelling. Oxford: Oxford University Press; 2010. and, as much as possible, was built on local parameters, including the adjustment of the infection and lethality rates to the Federal District age structure. Moreover, we sought to make transparent use of available evidence and deal in a robust manner with the parametric uncertainties by using Monte Carlo simulations. 2020 19. Briggs A, Claxton K, Sculpher M. Decision modelling for health economic evaluation. Oxford: Oxford University Press; 2011. 256 p. As demonstrated in this analysis of the Federal District, monitoring dynamic indicators of isolation and spread of infection, such as estimated Rt, along with their uncertainties, is an efficient and pragmatic strategy for monitoring policies intended to address COVID-19. 2323 22. Tariq A, Lee Y, Roosa K, Blumberg S, Yan P, Ma S, et al. Real-time monitoring the transmission potential of COVID-19 in Singapore, March 2020. medRxiv Prepr [Internet]. 2020 Jun [cited 2020 Sep 25]. Available from: https://dx.doi.org/10.1186%2Fs12916-020-01615-9
https://dx.doi.org/10.1186%2Fs12916-020-...

Despite compartmental models at expected levels being able to represent average behavior of epidemics in populations, they are not capable of portraying individual histories. 1212. Ministério da Saúde (BR). Secretaria de Ciência Tecnologia e Insumos Estratégicos. Departamento de Ciência e Tecnologia. Diretrizes metodológicas: diretriz de avaliação econômica [Internet]. Brasília: Ministério da Saúde; 2014 [citado 2020 set 25]. 132 p. Disponível em: http://bvsms.saude.gov.br/bvs/publicacoes/diretrizes_metodologicas_diretriz_avaliacao_economica.pdf
http://bvsms.saude.gov.br/bvs/publicacoe...
Likewise, estimating the dynamic values of the rates of effective contacts is discussed according to varied approaches in the literature, with different potential impacts on the results, depending on the calculation method. 2424 23. Chowell G, Sattenspiel L, Bansal S, Viboud C. Mathematical models to characterize early epidemic growth: a review. Phys Life Rev [Internet]. 2016 Sep [cited 2020 sep 25];18:66-97. Available from: https://doi.org/10.1016/j.plrev.2016.07.005
https://doi.org/10.1016/j.plrev.2016.07....
Although the adjustments relating to local age structure and calibration based on available records of hospitalizations were taken into consideration, other factors capable of interfering in the projection scenarios were not taken into consideration, such as comorbidities related to greater severity of the disease, including cardiovascular, respiratory and endocrine diseases. 2525 24. Wang X, Fang X, Cai Z, Wu X, Gao X, Min J, et al. Comorbid chronic diseases and acute organ injuries are strongly correlated with disease severity and mortality among COVID-19 patients: a systemic review and meta-analysis. Research (Wash D C) [Internet]. 2020 Apr [cited 2020 Sep 25];2020:2402961. Available from: https://doi.org/10.34133/2020/2402961
https://doi.org/10.34133/2020/2402961...
And still in relation to demographic and social characteristics, we emphasize that the model assumes the premise of random contact, and the values of its estimates may diverge from measurements obtained when considering the specific contact patterns of each group (influence of heterogeneity). When interpreting the number of deaths expected by the model, it is important to remember that we did not take into consideration the possibility of increased deaths because of lack of hospital beds in the event of a possible collapse of the health care network.

Finally, we highlight that the data shown here on registered beds should also be interpreted paying attention to the reference of the targets for reducing impact on the hospital network. This is because they do not take into consideration the real proportion of hospital bed occupancy due to other clinical conditions, nor additional beds that are being put in place by health authorities.

Despite its limitations, this analysis reinforces conclusions as to the positive effect of social distancing measures on reducing the impact of the spread of COVID-19 on demands on the hospital network. Given the alarming scenario projected in the event of interruption or low adherence to social distancing, it is essential that local health authorities consider daily monitoring of distancing indicators, spread of the disease and availability of hospital beds, when planning and implementing public policies to combat COVID-19.

References

  • 1
    Mohammadi M, Meskini M, Nascimento Pinto AL. 2019 Novel coronavirus (COVID-19) overview. Z Gesundh Wiss [Internet]. 2020 Apr [cited 2020 Sep 25];1-9. Available from: https://doi.org/10.1007/s10389-020-01258-3
    » https://doi.org/10.1007/s10389-020-01258-3
  • 2
    Ministério da Saúde (BR). Painel coronavírus [Internet]. Brasília: Ministério da Saúde; 2020 [citado 2020 set 25]. Disponível em: https://covid.saude.gov.br/
    » https://covid.saude.gov.br/
  • 3
    Rothan HA, Byrareddy SN. The epidemiology and pathogenesis of coronavirus disease (COVID-19) outbreak. J Autoimmun [Internet]. 2020 May [cited 2020 Sep 25];109:102433. Available from: https://doi.org/10.1016/j.jaut.2020.102433
    » https://doi.org/10.1016/j.jaut.2020.102433
  • 4
    Prem K, Liu Y, Russell TW, Kucharski AJ, Eggo RM, Davies N, et al. The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study. Lancet Public Health [Internet]. 2020 May [cited 2020 Sep 25];5(5):e261-70. Available from: https://doi.org/10.1016/s2468-2667(20)30073-6
    » https://doi.org/10.1016/s2468-2667(20)30073-6
  • 5
    Governo do Distrito Federal. Decreto n° 40.509, de 11 de março de 2020. Dispõe sobre as medidas para enfrentamento da emergência de saúde pública de importância internacional decorrente do novo coronavírus, e dá outras providências [Internet]. Diário Oficial do Distrito Federal, Brasília (DF), 2020 mar 11 [citado 2020 set 25];Edição extra. Disponível em: http://www.sinj.df.gov.br/sinj/Norma/ad0fae78af5f4e50b46c7357b7ee8597/Decreto_40509_11_03_2020.html
    » http://www.sinj.df.gov.br/sinj/Norma/ad0fae78af5f4e50b46c7357b7ee8597/Decreto_40509_11_03_2020.html
  • 6
    Rede CoVida. Os impactos das medidas de distanciamento social e redução de fluxo intermunicipal na Bahia. Bol CoVida [Internet]. 2020 abr [citado 2020 set 25];1(2):19. Disponível em: https://covid19br.org/main-site-covida/wp-content/uploads/2020/04/BoletimCovida_2ed_Final.pdf
    » https://covid19br.org/main-site-covida/wp-content/uploads/2020/04/BoletimCovida_2ed_Final.pdf
  • 7
    COVID-19 Brasil. Monitoramento e análises da situação do Coronavírus no Brasil [Internet]. São Paulo; 2020 [citado 2020 set 25]. Disponível em: https://ciis.fmrp.usp.br/covid19/
    » https://ciis.fmrp.usp.br/covid19/
  • 8
    Organização Pan-Americana de Saúde - OPAS. Ministério da Saúde (BR). Universidade de Brasília. Butantã. Pressão hospitalar por COVID-19 [Internet]. [S.l];2020. Disponível em: https://covid-calc.org/
    » https://covid-calc.org/
  • 9
    Presidência da República (BR). Casa Civil. Avaliação de políticas públicas: guia prático de análise ex post [Internet]. Brasília: Presidência da República; 2018 [citado 2020 set 25]. 301 p. Disponível em: https://www.ipea.gov.br/portal/index.php?option=com_content&view=article&id=34504
    » https://www.ipea.gov.br/portal/index.php?option=com_content&view=article&id=34504
  • 10
    Instituto Brasileiro de Geografia e Estatística - IBGE. Diretoria de Pesquisas. Coordenação de População e Indicadores Sociais. Gerência de Estudos e Análises da Dinâmica Demográfica. Projeção da população do Brasil e Unidades da Federação por sexo e idade para o período 2000-2030 [Internet]. Rio de Janeiro: Instituto Brasileiro de Geografia e Estatística; 2020 [citado 2020 set 25]. Disponível em: http://tabnet.datasus.gov.br/cgi/tabcgi.exe?ibge/cnv/projpopuf.def
    » http://tabnet.datasus.gov.br/cgi/tabcgi.exe?ibge/cnv/projpopuf.def
  • 11
    Inloco. Mapa brasileiro da COVID-19 [Internet]. São Paulo: Inloco; 2020 [citado 27 abr 2020]. Disponível em: https://www.inloco.com.br/pt/covid-19
    » https://www.inloco.com.br/pt/covid-19
  • 12
    Ministério da Saúde (BR). Secretaria de Ciência Tecnologia e Insumos Estratégicos. Departamento de Ciência e Tecnologia. Diretrizes metodológicas: diretriz de avaliação econômica [Internet]. Brasília: Ministério da Saúde; 2014 [citado 2020 set 25]. 132 p. Disponível em: http://bvsms.saude.gov.br/bvs/publicacoes/diretrizes_metodologicas_diretriz_avaliacao_economica.pdf
    » http://bvsms.saude.gov.br/bvs/publicacoes/diretrizes_metodologicas_diretriz_avaliacao_economica.pdf
  • 13
    Vynnycky E, White RG, Fine P. An introduction to infectious disease modelling. Oxford: Oxford University Press; 2010.
  • 14
    Rocha Filho TM, Santos FSG, Gomes VB, Rocha TAH, Croda JHR, Ramalho WM, et al. Expected impact of COVID-19 outbreak in a major metropolitan area in Brazil. medRxiv Prepr [Internet]. 2020 Mar [cited 2020 sep 25]. Available from: https://doi.org/10.1101/2020.03.14.20035873
    » https://doi.org/10.1101/2020.03.14.20035873
  • 15
    Hill A. Modeling COVID-19 spread vs healthcare capacity [Internet]. [S.l.]; 2020 [cited 2020 Sep 25]. Available from: https://alhill.shinyapps.io/COVID19seir/
    » https://alhill.shinyapps.io/COVID19seir/
  • 16
    Penn University. CHIME model: discrete-time SIR modeling of infections/recovery [Internet]. [S.l.]; 2020 [cited 2020 Sep 25]. Available from: https://code-for-philly.gitbook.io/chime/what-is-chime/sir-modeling
    » https://code-for-philly.gitbook.io/chime/what-is-chime/sir-modeling
  • 17
    Eikenberry SE, Mancuso M, Iboi E, Phan T, Eikenberry K, Kuang Y, et al. To mask or not to mask: modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic. Infect Dis Model [Internet]. 2020.
  • 18
    World Health Organization - WHO. “Immunity passports” in the context of COVID-19 [Internet]. Genebra: World Health Organization; 2020 [cited 2020 Sep 25]. Available from: https://www.who.int/news-room/commentaries/detail/immunity-passports-in-the-context-of-covid-19
    » https://www.who.int/news-room/commentaries/detail/immunity-passports-in-the-context-of-covid-19
  • 19
    18. Governo do Distrito Federal. Secretaria de Saúde. Boletins Informativos DIVEP/CIEVES (COE). Boletins informativos sobre coronavirus (COVID-19). Brasília: GDF; 2020. Disponível em: http://www.saude.df.gov.br/boletinsinformativos-divep-cieves/
    » http://www.saude.df.gov.br/boletinsinformativos-divep-cieves/
  • 20
    19. Briggs A, Claxton K, Sculpher M. Decision modelling for health economic evaluation. Oxford: Oxford University Press; 2011. 256 p.
  • 21
    20. Zimmermann I. Modelo de projeção da demanda por leitos de UTI por COVID-19. Mendeley Data [Internet]. 2020 set [citado 2020 set 25]. Disponível em: http://dx.doi.org/10.17632/bypzkmwnrv
    » http://dx.doi.org/10.17632/bypzkmwnrv
  • 22
    21. Mellan TA, Hoeltgebaum HH, Mishra S, Whittaker C, Schnekenberg RP, Gandy A, et al. Report 21: estimating COVID-19 cases and reproduction number in Brazil. Imperial College London [Internet]. 2020 May [cited 2020 Sep 25]. Available from: https://doi.org/10.25561/78872
    » https://doi.org/10.25561/78872
  • 23
    22. Tariq A, Lee Y, Roosa K, Blumberg S, Yan P, Ma S, et al. Real-time monitoring the transmission potential of COVID-19 in Singapore, March 2020. medRxiv Prepr [Internet]. 2020 Jun [cited 2020 Sep 25]. Available from: https://dx.doi.org/10.1186%2Fs12916-020-01615-9
    » https://dx.doi.org/10.1186%2Fs12916-020-01615-9
  • 24
    23. Chowell G, Sattenspiel L, Bansal S, Viboud C. Mathematical models to characterize early epidemic growth: a review. Phys Life Rev [Internet]. 2016 Sep [cited 2020 sep 25];18:66-97. Available from: https://doi.org/10.1016/j.plrev.2016.07.005
    » https://doi.org/10.1016/j.plrev.2016.07.005
  • 25
    24. Wang X, Fang X, Cai Z, Wu X, Gao X, Min J, et al. Comorbid chronic diseases and acute organ injuries are strongly correlated with disease severity and mortality among COVID-19 patients: a systemic review and meta-analysis. Research (Wash D C) [Internet]. 2020 Apr [cited 2020 Sep 25];2020:2402961. Available from: https://doi.org/10.34133/2020/2402961
    » https://doi.org/10.34133/2020/2402961
  • 26
    25. Sanche S, Lin YT, Xu C, Romero-Severson E, Hengartner N, Ke R. The novel Coronavirus, 2019-nCoV, is highly contagious and more infectious than initially estimated. medRxiv Prepr [Internet]. 2020 Feb [cited 2020 Sep 25]. Available from: https://doi.org/10.1101/2020.02.07.20021154
    » https://doi.org/10.1101/2020.02.07.20021154
  • 27
    26. Linton NM, Kobayashi T, Yang Y, Hayashi K, Akhmetzhanov AR, Jung S, et al. Incubation period and other epidemiological characteristics of 2019 novel coronavirus infections with right truncation: a statistical analysis of publicly available case data. J Clin Med [Internet]. 2020 Feb [cited 2020 Sep 25];9(2):538. Available from: https://dx.doi.org/10.3390%2Fjcm9020538
    » https://dx.doi.org/10.3390%2Fjcm9020538
  • 28
    27. Verity R, Okell LC, Dorigatti I, Winskill P, Whittaker C, Imai N, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis [Internet]. 2020 Mar [cited 2020 Sep 25];20(6):669-77. Available from: https://doi.org/10.1016/S1473-3099(20)30243-7
    » https://doi.org/10.1016/S1473-3099(20)30243-7
  • 29
    28. Neil M Ferguson, Laydon D, Nedjati-Gilani G, Imai N, Ainslie K, Baguelin M, et al. Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand. Imperial College London [Internet]. 2020 [cited 2020 Sep 25]. Available from: https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-9-impact-of-npis-on-covid-19/
    » https://www.imperial.ac.uk/mrc-global-infectious-disease-analysis/covid-19/report-9-impact-of-npis-on-covid-19/
  • 30
    29. Grasselli G, Zangrillo A, Zanella A, Antonelli M, Cabrini L, Castelli A, et al. Baseline characteristics and outcomes of 1591 patients infected with SARS-CoV-2 admitted to ICUs of the Lombardy Region, Italy. JAMA [Internet]. 2020 Apr [cited 2020 Sep 25];323(16):1574-81. Available from: https://dx.doi.org/10.1001%2Fjama.2020.5394
    » https://dx.doi.org/10.1001%2Fjama.2020.5394

Supplementary Material 1 - Transitions and equations

Following contact with an infected person, an individual who up to then was susceptible ( S ) becomes exposed ( E ) during a period of latency. This period is related to average disease incubation time. Following the latency period, the individual becomes infectious ( I1 ) and can spread infection before being hospitalized ( I2 or I3 ) or being removed ( R ) from circulation by being placed in quarantine. Hospitalization was foreseen for individuals who, following a mild clinical picture ( I1 ), progressed to a severe clinical picture ( I2 ) or a critical clinical picture ( I3 ). In keeping with the follow-up of patients with COVID-19, in the case of patients who progress to a critical clinical picture, transition to death is foreseen according to specific lethality of the disease ( µ ). In each of the compartments of hospitalized individuals, the possibility is also foreseen of them being removed from the model depending on length of hospitalization ( γ−1 ). The transitions between the compartments were solved by a system of equations, whereby the total population ( N ) is equal to and represents the sum of susceptible ( S ), exposed ( E ), infected ( I ), removed ( R ) and dead ( M ) individuals at each moment in time t of the model. The model's differential equation system is shown in equations 1 to 8:

(1)S(t+1)=S(t)(βI1)S(t)N
(2)E(t+1)=E(t)+(βI1)S(t)NαE(t)
(3)I1(t+1)=I1(t)+αE(t)(γ1+p1)I1(t)
(4)I2(t+1)=I2(t)+p1I1(t)(γ2+p2)I2(t)
(5)I3(t+1)=I2(t)+p2I2(t)(γ3+μ)I3(t)
(6)R(t+1)=R(t)+γ1I1(t)+γ2I2(t)+γ3I3(t)
(7)M(t+1)=M(t)+μI3(t)
(8)N=S(t)+E(t)+I1(t)+I2(t)+I3(t)+R(t)+M(t)

Below we present equations that are complementary to the model's differential equation system, where fhosp refers to the fraction of individuals who progressed to hospitalization, futi refers to the fraction of hospitalized individuals who progressed to a critical clinical picture, and ifr refers to lethality:

(9)α=(incubation time)1
(10)δ=ifr(fhospfuti)
(11)μ=δ1(time unit_deathtime unit isolation_or hospitalizationtime unit isolationor hospitalization)
(12)γ1=fhosp(time until isolation or hospitalization)1
(13)γ2=(1futi)(time in a regular bed)1
(14)γ3=(1δ)(time in an ICU bed)1
(15)p1=(time until isolation or hospitalization)1γ1
(16)p2= futi(time in a regular bed)1γ1
Relationship between social distancing and transmission coefficient

With the aim of simulating the spread of the disease in view of the effects of social distancing, we calculated the transmission coefficient ( β ), based on available epidemic estimates and social isolation estimates. 1717. Eikenberry SE, Mancuso M, Iboi E, Phan T, Eikenberry K, Kuang Y, et al. To mask or not to mask: modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic. Infect Dis Model [Internet]. 2020.,2424 23. Chowell G, Sattenspiel L, Bansal S, Viboud C. Mathematical models to characterize early epidemic growth: a review. Phys Life Rev [Internet]. 2016 Sep [cited 2020 sep 25];18:66-97. Available from: https://doi.org/10.1016/j.plrev.2016.07.005
https://doi.org/10.1016/j.plrev.2016.07....
As such, we take into consideration here the dynamic relationship whereby the effect contact rate ( β ) at each moment in time follows the function of a transmissibility factor ( τ ), this being a more stable factor strongly associated with pathogen virulence, and the function of a exposed individual contact factor ( c ), this being a more dynamic factor and dependent on social isolation, i.e. β = τc . As such, the isolation levels recorded daily in the Federal District time series, are imputed in the model as a complement of factor c . Likewise, similarly to what was proposed by Eikenberry et al., 1717. Eikenberry SE, Mancuso M, Iboi E, Phan T, Eikenberry K, Kuang Y, et al. To mask or not to mask: modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic. Infect Dis Model [Internet]. 2020. the value of τ was estimated by means of calibration using the values of the time series of the number of ICU hospitalizations due to COVID-19 in the Federal District between 22/3/2020 and 21/5/2020.

Once the transmission coefficient β was obtained, we also calculated the instantaneous reproduction number ( Rt ) for each day in the period analyzed, with the aim of illustrating the dynamics of the spread of COVID-19 as per the proposed scenarios. In keeping with previous methodological approaches, the instantaneous reproduction number ( Rt ) considers the specific rates of entering and leaving the infectious compartment by means of equation 17:

(17)Rt=β(t)S(t)(p1+γ1)

Supplementary Material 2 - Infection, hospitalization and lethality rates expected for COVID-19, according to age structure, in the Federal District, Brazil

Supplementary Material 3 - Time series of cumulative confirmed COVID-19 cases, intensive care hospitalizations, number of deaths and percentage social isolation in the Federal District, Brazil

Supplementary Material 4 - Illustration of the probability distribution functions of the model's main parameters

Publication Dates

  • Publication in this collection
    14 Dec 2020
  • Date of issue
    2020

History

  • Received
    29 May 2020
  • Accepted
    24 Aug 2020
Secretaria de Vigilância em Saúde e Ambiente - Ministério da Saúde do Brasil Brasília - Distrito Federal - Brazil
E-mail: ress.svs@gmail.com