Acessibilidade / Reportar erro

PREDICTIVE MODEL FOR COVID-19 INCIDENCE IN A MEDIUM-SIZED MUNICIPALITY IN BRAZIL (PONTA GROSSA, PARANÁ)

MODELO PREDICTIVO DE LA OCURRENCIA DE COVID-19 EN UN MUNICIPIO DE TAMAÑO MEDIO EN BRASIL (PONTA GROSSA-PARANÁ)

ABSTRACT

Objective:

to produce a predictive model for the incidence of COVID-19 cases, severity and deaths in Ponta Grossa, state of Paraná.

Methods:

this is an ecological study with data from confirmed cases of COVID-19 reported between March 21, 2020 and May 3, 2020 in Ponta Grossa and proportion of severity, hospitalization and lethality in the literature. A susceptible-infected-recovered (SIR) epidemic model was developed, and reproduction rate (R0), duration of epidemic, peak period, number of cases, hospitalized patients and deaths were estimated. Deaths were calculated by age group and in three scenarios: at day 24, at day 34, and at day 44 of the epidemic.

Results:

in the three scenarios assessed in this study, the variation in the number of cases was explained by an exponential curve (r2=0.74, 0.79 and 0.89, respectively, p<0.0001 in all scenarios). The SIR model estimated that, in the best scenario, the peak period will be around 120 days after the first case (between July 11, 2020 and July 25, 2020), estimated R0 will be 1.07 and will infect 0.23% of the population. In the worst scenario, peak period will involve 4,375 (95% CI; 4156-4594) cases and 825 (95% CI; 700-950) cases in the best scenario. Most cases and hospital admissions will involve patients aged 20 to 39 years, the number of deaths will be higher among the elderly and more pronounced among patients aged ≥80 years.

Conclusion:

this is the first study that provides COVID-19 projections for a municipality that is not a large capital. It shows a peak period at a later moment; therefore, the municipality will have more time to prepare and adopt protective measures to reduce the number of simultaneous cases.

DESCRIPTORS
Estimate; Epidemics; Basic reproduction number; Coronavirus; SARS-CoV-2; Mortality; Hospitalization

RESUMEN

Objetivo:

obtener un modelo predictivo para la ocurrencia de casos, severidad y muertes por COVID-19 en Ponta Grossa-Paraná.

Métodos:

estudio ecológico con datos de casos confirmados de COVID-19 notificados del 21/03/2020 al 3/3/2020 en Ponta Grossa y proporción de severidad, hospitalización y letalidad en la literatura. Se construyó un modelo epidemiológico (SIR) infectado-recuperado susceptible y tasa de reproducción estimada (R0), duración de la epidemia, fecha pico, número de casos, hospitalizaciones y muertes. Este último por grupo de edad y en tres escenarios: a los 24 días, a los 34 días y a los 44 días de epidemia.

Resultados:

en los tres escenarios evaluados, la variación en el número de casos se explicó por una curva exponencial (r2 = 0.74, 0.79 y 0.89, respectivamente y p <0.0001 en total). El modelo SIR estimó que, en el mejor escenario, el pico ocurrirá alrededor de 120 días después del primer caso (entre el 7/11/2020 y el 25/7/2020), el R0 estimado será de 1.07 y alcanzará 0.23 % de habitantes infectados. En el peor de los casos, el pico estimado será de 4375 (IC del 95%: 4156-4594) y 825 (IC del 95%: 700-950) en el mejor de los casos. El mayor número estimado de casos y hospitalizaciones estará en el rango entre 20 y 39 años, el número de muertes será mayor entre los ancianos y más pronunciado entre ≥ 80 años.

Conclusión:

este es el primer estudio con proyecciones para COVID-19 en un municipio fuera de las grandes capitales y demostró que el pico llegará tarde, por lo tanto, el municipio tendrá más tiempo de preparación y que las medidas de protección pueden reducir el número simultáneo de casos.

DESCRIPTORES
Estimación; Epidemias Número reproductivo basal; Coronavirus; SARS-CoV-2; Letalidad; Hospitalización

RESUMO

Objetivo:

obter um modelo preditivo da ocorrência de casos, gravidade e óbitos por COVID-19 em Ponta Grossa-Paraná.

Métodos:

estudo ecológico com dados de casos confirmados de COVID-19 notificados de 21/03/2020 a 03/05/2020 em Ponta Grossa e proporção de gravidade, hospitalização e letalidade da literatura. Um modelo epidemiológico suscetível-infectado-recuperado (SIR) foi construído e estimadas taxa de reprodução (R0), duração da epidemia, data do pico, número de casos, hospitalizações e óbitos. Estas últimas por faixa etária e em três cenários: aos 24 dias, aos 34 dias e aos 44 dias de epidemia.

Resultados:

nos três cenários avaliados, a variação no número de casos foi explicada por uma curva exponencial (r2=0,74, 0,79 e 0,89, respectivamente e p<0,0001 em todos). O modelo SIR estimou que, no melhor cenário, o pico ocorrerá em torno de 120 dias após o primeiro caso (entre 11/07/2020 e 25/07/2020), o R0 estimado será 1,07 e chegará a 0,23% dos habitantes infectados. No pior cenário, o pico estimado será de 4375 (IC 95% 4156-4594) casos e 825 (IC 95% 700-950) no melhor cenário. O maior número estimado de casos e hospitalizações será na faixa entre 20 e 39 anos, o número de óbitos será maior entre idosos e mais acentuado entre ≥ 80 anos.

Conclusão:

este é o primeiro estudo com projeções para a COVID-19 em um município fora das grandes capitais e mostrou que o pico será tardio, portanto, o município terá mais tempo de preparo e que medidas protetivas podem reduzir o número simultâneo de casos.

DESCRITORES
Previsão; Epidemia; Número reprodutivo basal; Coronavirus; SARS-CoV-2; Letalidade; Hospitalização

INTRODUCTION

In 2020, a respiratory syndrome was described, presenting different characteristics from previously diagnosed syndromes (SARS-CoV, MERS-CoV), caused by a new microorganism initially called 2019-nCoV. It causes a more serious respiratory disease in populations at risk, such as adults over 60 years and people with comorbidities.11. Zhu N, Zhang D, Wang W, Li X, Yang B, et al. A novel coronavírus from patients with pneumonia in China, 2019. N Engl J Med [Internet]. 2020 [cited 2020 Apr 28];382:727-33. Available from: https://dx.doi.org/10.1056/nejmoa2001017
https://dx.doi.org/10.1056/nejmoa2001017...
-22. Perlman S. Another decade, another coronavirus. N Engl J Med [Internet]. 2020 [cited 2020 Apr 28];382:8. Available from: https://dx.doi.org/10.1056/nejme1917479
https://dx.doi.org/10.1056/nejme1917479...
The virus was first identified in Wuhan, in the province of Hubei, China, in individuals exposed to a market of live animals and seafood. China’s sanitary authorities reported a fast increase in the number of confirmed cases and many deaths due to this new disease, called COVID-19. Since the disease appeared in China, it soon started to be identified in other countries.33. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) - China, 2020. CCDC Weekly [Internet]. 2020 [cited 2020 Apr 03];2(x):1-10. Available from: https://newsus.cgtn.com/news/files/COVID-19.pdf
https://newsus.cgtn.com/news/files/COVID...

On March 11, 2020, a pandemic was declared by the World Health Organization (WHO) due to the new coronavirus.44. Organização das Nações Unidas News [Internet]. Organização Mundial da Saúde declara novo coronavírus uma pandemia;2020 [cited 2020 Apr 24]. Available from: https://news.un.org/pt/story/2020/03/1706881
https://news.un.org/pt/story/2020/03/170...
The Johns Hopkins University reported 2,744,511 global cases and 192,982 deaths by COVID-19 worldwide55. Johns Hopkins University & Medicine[Internet]. Coronavirus resource center; 2020 [cited 2020 Apr 24]. Available from: https://coronavirus.jhu.edu/map.html
https://coronavirus.jhu.edu/map.html...
until April 24, 2020, and the estimated lethality rate varied from 1.18% in Germany up to 12.07% in Italy.66. European Centre for Disease Prevention and Control [Internet]. Coronavirus disease 2019 (COVID-19) pandemic: increased transmission in the EU/EEA and the UK - seventh update; 2020 [cited 2020 Apr 03]. Available from: https://www.ecdc.europa.eu/en/publications-data/rapid-risk-assessment-coronavirus-disease-2019-covid-19-pandemic#no-link
https://www.ecdc.europa.eu/en/publicatio...

In Brazil, until April 23, 2020, 49,492 cases had been reported, and a lethality rate of 6.7% for any age group. All regions of the country have confirmed cases and deaths by COVID-19, but distinct distribution, behavior and mortality rates are observed in the different Brazilian regions, with a greater number of confirmed cases and deaths in the Southeast and Northeast regions, respectively. On the same day, the South region presented the second lowest number of cases of COVID-19, behind the Central West region.77. Ministério da Saúde [Internet]. Coronavirus Brasil: Painel Coronavirus (Atualizado 23/04/2020); 2020 [cited 2020 Apr 24]. Available from: https://covid.saude.gov.br/
https://covid.saude.gov.br/...
Paraná had 1,082 confirmed cases and a lethality rate of 5.5% on the same day.88. Secretaria da Saúde do Paraná [Internet]. Informe epidemiológico COVID-19 (Atualizado 23/04/2020); 2020 [cited 2020 Apr 24]. Available from: http://www.saude.pr.gov.br/modules/conteudo/conteudo.php?conteudo=3507
http://www.saude.pr.gov.br/modules/conte...

To control the pandemic, gradual social distancing was adopted in the country, but not simultaneously in all states, in addition to closed borders, investments in public health and health infrastructure, and hiring of new health professionals, among other actions. However, such measures were adopted in general at the national level, without taking into account local specificities and structures of states and municipalities.99. Ministério da Saúde [Internet]. Portal Coronavirus (COVID-19);2020 [cited 2020 Apr 24]. Available from: https://coronavirus.saude.gov.br/
https://coronavirus.saude.gov.br/...

Considering Brazil in general and large centers such as São Paulo and Rio de Janeiro, studies have been conducted to find predictive models and help health services develop coping actions.55. Johns Hopkins University & Medicine[Internet]. Coronavirus resource center; 2020 [cited 2020 Apr 24]. Available from: https://coronavirus.jhu.edu/map.html
https://coronavirus.jhu.edu/map.html...
,1010. Cidade de São Paulo: Saúde [Internet]. Vigilância em saúde: boletins COVID-19; 2020 [cited 2020 Apr 24]. Available from: https://www.prefeitura.sp.gov.br/cidade/secretarias/saude/vigilancia_em_saude/index.php?p=295572
https://www.prefeitura.sp.gov.br/cidade/...
-1111. Ribeiro RSM. Previsões a partir do modelo epidemiológico SIR para os casos de infecção pelo COVID-19: Uma aplicação para os estados brasileiros. Belo Horizonte (MG): Grupo de Políticas Públicas e Desenvolvimento, Centro de Desenvolvimento e Planejamento Regional da UFMG; Nota técnica de fevereiro de 2020 [Internet]; [cited 2020 Apr 16]. Available from: https://pesquisas.face.ufmg.br/gppd/wp-content/uploads/sites/24/2020/03/30.03-previsões-de-infeccção-pelo-covid-19-aplicação-para-os-estados-brasileiros.pdf
https://pesquisas.face.ufmg.br/gppd/wp-c...
However, no studies have been developed showing estimated time and intensity of epidemic incidence for small and medium-sized municipalities. Given the need to implement measures for the prevention of virus dissemination and acquisition of patient care infrastructure, this study aimed to produce a predictive model for the disease incidence, severity and deaths by COVID-19 in 2020 in Ponta Grossa, a medium-sized municipality in the state of Paraná, Brazil.

METHOD

Study design

This is an ecological study that used data available on government digital platforms (as described below) without analysis on humans or animals; therefore, not requiring an assessment by a research ethics committee. An epidemic curve with an epidemiological model was built and the concepts, methods and scenarios evaluated are described as follows.

Data sources and study variables

Daily numbers of confirmed cases of COVID-19 in the municipality of Ponta Grossa were used from March 21, 2020 (when the first case was reported) to May 3, 2020, obtained from official public records of the Health Division of the State of Paraná.88. Secretaria da Saúde do Paraná [Internet]. Informe epidemiológico COVID-19 (Atualizado 23/04/2020); 2020 [cited 2020 Apr 24]. Available from: http://www.saude.pr.gov.br/modules/conteudo/conteudo.php?conteudo=3507
http://www.saude.pr.gov.br/modules/conte...
A ‘confirmed case’ means a positive result in RT-PCR testing for COVID-19, or by clinical-epidemiological criteria, close contact at home with a laboratory-confirmed case, with fever and/or any respiratory symptom within 14 days after the last contact with the patient and for which specific laboratory investigation was not possible.1212. Secretaria da Saúde do Paraná [Internet]. Como é definido um caso suspeito do coronavírus?; 2020 [cited 2020 Apr 19]. Available from: http://www.saude.pr.gov.br/modules/conteudo/conteudo.php?conteudo=3514
http://www.saude.pr.gov.br/modules/conte...
Therefore, the first 24 days of the epidemic were used as the primary input for the SIR (susceptible-infected-recovered) model described below. The population of Ponta Grossa considered in this study was 352,749 inhabitants.1313. Instituto Paranaense de Desenvolvimento Econômico e Social [Internet]. Projeção da população dos municípios do Paraná - 2018 - 2040; 2020 [cited 2020 Apr 16]. Available from:http://www.ipardes.pr.gov.br/pagina/projecao-da-populacao-dos-municipios-do-parana-2018-2040
http://www.ipardes.pr.gov.br/pagina/proj...

The following data were considered: 1.2% asymptomatic cases, 80.3% mild and moderate cases, 13.8% severe cases, and 4.7% critical cases, according to a descriptive study on COVID-19 cases in China published by the country’s disease control center (CDC). Given the fact that it is the territory where the epidemic started and where the first epidemiological studies on the disease were published, and that the local CDC was used as the source, these were considered reliable data.33. Novel Coronavirus Pneumonia Emergency Response Epidemiology Team. The epidemiological characteristics of an outbreak of 2019 novel coronavirus diseases (COVID-19) - China, 2020. CCDC Weekly [Internet]. 2020 [cited 2020 Apr 03];2(x):1-10. Available from: https://newsus.cgtn.com/news/files/COVID-19.pdf
https://newsus.cgtn.com/news/files/COVID...
For the percentage of infected individuals, hospital admissions and deaths in each age group, an intense search was conducted in the literature and several studies were combined to obtain the parameters. These and other references are presented in detail in the Supplementary Material 1 .

Statistical methods: SIR (susceptible-infected-recovered) epidemiological model

The model used in this study was developed by Kermack and McKendrick in 1927 and consists of explaining the behavior of an epidemic by dividing a population into 3 groups: susceptible (S), infected (I) and recovered (R); then, the model calculates the system of differential equations:1414. Bacaër N. McKendrick and Kermack on epidemic modelling (1926-1927). In: A short history of mathematical population dynamics. London: Springer [Internet]; 2011. p. 89-96 [cited 2020 Apr 16]. Available from: https://dx.doi.org/10.1007/978-0-85729-115-8_16
https://dx.doi.org/10.1007/978-0-85729-1...

d S d t = - β I N S d I d t = β I N S - γ I d R d t = γ I

Where:

N= population of Ponta Grossa 20201313. Instituto Paranaense de Desenvolvimento Econômico e Social [Internet]. Projeção da população dos municípios do Paraná - 2018 - 2040; 2020 [cited 2020 Apr 16]. Available from:http://www.ipardes.pr.gov.br/pagina/projecao-da-populacao-dos-municipios-do-parana-2018-2040
http://www.ipardes.pr.gov.br/pagina/proj...

S= number of susceptible individuals (all population was assumed as susceptible in the beginning of the epidemic)

I= number of infected individuals

R= number of recovered individuals

β= contagion coefficient

γ= recovery coefficient

t= time in days

Contagion and recovery coefficients were obtained by adjusting data (input) and by minimizing the total sum of squares (TSS) of the difference between the absolute number of infected individuals and the number of infected individuals adjusted from the system of differential equations above. The final values of β and γ were obtained by an iterative optimization process using the method developed by Byrd et al.1616. Byrd RH, Lu P, Nocedal J, and Zhu C. A limited memory algorithm for bound constrained optimization. SIAM Journal on Scientific Computing [Internet]. 1995 [cited 2020 Apr 16];16: 1190-1208. Available from: https://dx.doi.org/10.1137/0916069
https://dx.doi.org/10.1137/0916069...
that allows restrictions on the parametric domain of data. This analysis is more robust for non-linear optimization problems than Monte Carlo simulations, for instance. In the end, the reproduction rate (R0) was obtained, which is an estimated number of individuals infected by one infected person. The following values were then generated: SIR curves from the beginning to the end of the epidemic, estimated duration of the epidemic, estimated peak period, curve to estimate the number of cases over time, and percentage of infected individuals at the peak period.

Three scenarios were simulated to estimate the curves of expected cases:

Scenario 1: estimated curve based on confirmed data from the first 24 days (March 21, 2020 to April 13, 2020).

Scenario 2: estimated curve based on confirmed data from the first 34 days (March 21, 2020 to April 23, 2020).

Scenario 3: estimated curve based on confirmed data from the first 44 days (March 21, 2020 to May 3, 2020).

All statistical analyses were performed by the study authors using R software.1515. R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available from: https://www.R-project.org/
https://www.R-project.org/...

RESULTS

There were 5 confirmed cases in the first 24 days since the first case reported on March 21, 2020 in the municipality, characterizing a cumulative incidence in the period of 1.42 cases per 100,000 inhabitants. Dispersion of cases over time and the estimated curve during the first 24, 34 and 44 days are illustrated in Figures 1A, 1B and 1C, respectively. The curves that best explained data variation in the three scenarios were exponential curves of r2=0.74 (p<0.0001), r2=0.79 (p<0.0001) and r2=0.89 (p<0.0001), respectively. Variation in the scenarios showed a growing curve and decline of the growth angle over the period. The coefficient started at 0.069 and ended at 0.054; that is, the speed during the first days was more accentuated than the speed in the whole period.

Figure 1 -
Curves of confirmed cases during the first 24, 34 and 44 days of the COVID-19 epidemic (day 0 = March 21, 2020) in the municipality of Ponta Grossa, Paraná, Brazil, 2020.

A SIR model was built with data from the first 24 days of the disease and the city population (Figure 2A). The graph in Figure 2A starts on day 1 (March 21, 2020) with the susceptible population on the day the first case was reported (N - 1), number of infected individuals = 1 and number of recovered individuals = 0. Over time, a fluctuation occurred among susceptible, infected and recovered individuals, reaching its peak around 120 days after the first case, that is, between July 10, 2020 and July 24, 2020, and a balance after about 170 days, with zero infected individuals and N = recovered + susceptible individuals. At the end of the period, the reproduction rate (R0) of the municipality would be 1.18 and, at the peak, it would infect around 1.24% of the population. Figures 2B and 2C show the SIR models generated for the first 34 and 44 days. The reproduction rate would be 1.14 and 1.07, respectively, and would infect 0.88% and 0.23% of the population at the peak. Given that the real exponential growth of cases at day 44 was slower than in the first 24 days, it also impacted the estimates of the SIR model.

Figure 2 -
Estimates generated by the SIR (S = susceptible, I = infected and R = recovered) model for 200 days of the epidemic based on confirmed cases during the first 24, 34 and 44 days of the COVID-19 epidemic (day 0 = March 21, 2020 ) in the municipality of Ponta Grossa, Paraná, Brazil, 2020.

Figure 3 shows the estimates for the simulation of three scenarios regarding the total number of cases - mild and moderate, and severe and critical cases. The scenario with data from the first 24 days was the worst scenario, with an estimated peak of infected individuals of 4,375 (5% lower, corresponding to 4,156, and 5% higher, corresponding to 4,594), of these, 3,513 (±5%, 3338 -3689) were mild and moderate cases, and 809 (±5%, 769-850) were severe and critical cases. The scenario with data from the first 34 days was an intermediate scenario, with an estimated peak of infected individuals of 3,117 (±5%, 2961-3273), of these, 2,503 (±5%, 2378-2628) were mild and moderate cases, and 577 (±5%, 548-605) were severe and critical cases. The scenario with data from the first 44 days was the best scenario, with an estimated peak of infected individuals of 825 (±5% 783-866), of these, 662 (±5%, 629-695) were mild and moderate cases, and 153 (±5%, 145-160) were severe and critical cases. Along the first 44 days of the epidemic, it slowed down when compared to estimates at the beginning of the epidemic.

Figure 3 -
Graphs of the SIR model showing estimated numbers of COVID-19 cases, with the total number of cases (red curve), number of mild and moderate cases (green curve) and severe and critical cases (blue curve) in three scenarios: confirmed cases during the first 24, 34 and 44 days of the COVID-19 epidemic (day 0 = March 21, 2020) in the municipality of Ponta Grossa, Paraná, Brazil, 2020.

In all scenarios estimated by age group (Table 1), the highest number of cases and hospital admissions due to COVID-19 were calculated for the 20-29 and 30-39 age groups. However, the estimated absolute number of deaths was higher among the elderly, and the number of deaths among the oldest group (≥80 years), in all scenarios, was higher than the number of deaths in the age group of 40 years, and at least twice the number of individuals aged 70-79 years.

Table 1 -
Estimated peak of COVID-19 cases in each age group and 95% confidence interval (95% CI) in the three scenarios: confirmed cases in the first 24, 34 and 44 days of the COVID-19 epidemic (day 0 = March 21, 2020) in the municipality of Ponta Grossa, Paraná, Brazil, 2020.

DISCUSSION

This is the first study to characterize COVID-19 projections for a municipality that is not a large capital city of Brazil. The characteristics of this municipality are similar to the reality of many others; therefore, the methodology used in this study can be adapted and extrapolated to other locations. The municipality population is 352,749 (fourth in the state of Paraná), covering an area of 2,055 km2, with a demographic density of 171.6 inhabitants/km2.1313. Instituto Paranaense de Desenvolvimento Econômico e Social [Internet]. Projeção da população dos municípios do Paraná - 2018 - 2040; 2020 [cited 2020 Apr 16]. Available from:http://www.ipardes.pr.gov.br/pagina/projecao-da-populacao-dos-municipios-do-parana-2018-2040
http://www.ipardes.pr.gov.br/pagina/proj...
According to the Brazilian Institute of Geography and Statistics (IBGE), medium-sized municipalities are those with a population between 100,000 and 500,000; therefore, this classification includes 276 municipalities in Brazil, along with Ponta Grossa, the 22nd in Paraná.1717. IBGE [Internet]. Indicadores sociais municipais [cited 2020 Apr 28]. Available from: https://biblioteca.ibge.gov.br/visualizacao/livros/liv54598.pdf
https://biblioteca.ibge.gov.br/visualiza...

In terms of number of cases and cumulative incidence, Ponta Grossa was below Curitiba and São Paulo, where at day 24, the incidence was 5.89 and 8.89 cases/100,000 inhabitants, respectively, and similar to the state of Paraná in general, where the incidence was 2.52 cases/100,000 inhabitants in the same period.88. Secretaria da Saúde do Paraná [Internet]. Informe epidemiológico COVID-19 (Atualizado 23/04/2020); 2020 [cited 2020 Apr 24]. Available from: http://www.saude.pr.gov.br/modules/conteudo/conteudo.php?conteudo=3507
http://www.saude.pr.gov.br/modules/conte...
,1010. Cidade de São Paulo: Saúde [Internet]. Vigilância em saúde: boletins COVID-19; 2020 [cited 2020 Apr 24]. Available from: https://www.prefeitura.sp.gov.br/cidade/secretarias/saude/vigilancia_em_saude/index.php?p=295572
https://www.prefeitura.sp.gov.br/cidade/...
This result points to the need for differentiated and proportional estimates that take into account the characteristics of municipalities other than large centers.

Although the small number of cases and low cumulative incidence may indicate a small magnitude of the disease in the municipality, the proportional approximation to the exponential curve in all scenarios shows that, despite the specificities, the variation of cases will be similar to what has been reported in other realities. The behavior of the epidemic in the first stage with exponential growth was also found in several other locations in Brazil and worldwide.55. Johns Hopkins University & Medicine[Internet]. Coronavirus resource center; 2020 [cited 2020 Apr 24]. Available from: https://coronavirus.jhu.edu/map.html
https://coronavirus.jhu.edu/map.html...
,1515. R Core Team (2019). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available from: https://www.R-project.org/
https://www.R-project.org/...
,1818. Governo do Estado do Rio de Janeiro [Internet]. Boletim diário: coronavírus; 2020 [cited 2020 Apr 24]. Available from: https://coronavirus.rj.gov.br/boletins/
https://coronavirus.rj.gov.br/boletins/...
-2020. Remuzzi A, Remuzzi G. COVID-19 and Italy: what next?. Lancet [Internet]. 2020 [cited 2020 Apr 04];395:1225-1228. Available from: https://dx.doi.org/10.1016/S0140-6736(20)30627-9
https://dx.doi.org/10.1016/S0140-6736(20...

The estimated peak of disease cases after around 120 days in all scenarios may indicate the municipality will have more time to adopt measures to handle the disease than other municipalities and that the measures adopted in the first 44 days may have delayed its incidence. In China, the peak period of the disease was observed around 90 days after the first case was reported.2121. Wang H, Wang Z, Dong Y, Chang R, Xu C, Zhang S, et al. Phase-adjusted estimation of the number of coronavirus disease 2019 cases in Wuhan, China. Cell Discov [Internet]. 2020 [cited 2020 Apr 16];6:10. Available from: https://dx.doi.org/10.1038/s41421-020-0148-0
https://dx.doi.org/10.1038/s41421-020-01...
In Brazil, a study using the same SIR model estimated that in the state of São Paulo, the peak period would be within around 70 days, and in Minas Gerais Gerais and Rio de Janeiro, within about 60 days.1111. Ribeiro RSM. Previsões a partir do modelo epidemiológico SIR para os casos de infecção pelo COVID-19: Uma aplicação para os estados brasileiros. Belo Horizonte (MG): Grupo de Políticas Públicas e Desenvolvimento, Centro de Desenvolvimento e Planejamento Regional da UFMG; Nota técnica de fevereiro de 2020 [Internet]; [cited 2020 Apr 16]. Available from: https://pesquisas.face.ufmg.br/gppd/wp-content/uploads/sites/24/2020/03/30.03-previsões-de-infeccção-pelo-covid-19-aplicação-para-os-estados-brasileiros.pdf
https://pesquisas.face.ufmg.br/gppd/wp-c...
Therefore, large cities may present a decline in the epidemic curve, when medium-sized municipalities, such as the one analyzed in this study, will have an increase in cases, admissions and deaths.

While these municipalities will have an opportunity to prepare for the disease dissemination and complications, they will face an increasing incidence of COVID-19 concomitant with the respiratory infections that cyclically appear in winter season. However, regional services that offer health care to more serious cases, including COVID-19 and other respiratory complications caused by infections, are generally located in medium-sized municipalities, so if the predicted scenarios become a reality, they will take place with situations observed every year and which overload the health services in these locations. In this case, surveillance and prevention and planning efforts are highly important, not only for COVID-19, but to ensure a proper response to what has been anticipated for the next months.

Reproduction rates (R0) of 1.18, 1.14 and 1.07 were below those estimated in other national and international studies. In a study published by the response committee of Imperial College London with an estimated incidence of COVID-19 for several countries, including Brazil, R0 ranged from 2.5 in Denmark to 4.8 in Spain.2222. Flaxman S, Mishra S, Gandy A, Unwin H, Coupland H, Mellan T, et al. Report 13: estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 european countries. Imperial College London [Internet]. 2020 [cited 2020 Apr 16]. Available from: https://dx.doi.org/10.25561/77731
https://dx.doi.org/10.25561/77731...
For Brazil, estimates ranged from 1.4 in Rio de Janeiro to 2.6 in Ceará. The number obtained for Ponta Grossa was also reflected in the final proportion of infected individuals of 1.24%, which for other states, ranged from 5.23% in Rio Grande do Sul to 20.5% in Ceará.1111. Ribeiro RSM. Previsões a partir do modelo epidemiológico SIR para os casos de infecção pelo COVID-19: Uma aplicação para os estados brasileiros. Belo Horizonte (MG): Grupo de Políticas Públicas e Desenvolvimento, Centro de Desenvolvimento e Planejamento Regional da UFMG; Nota técnica de fevereiro de 2020 [Internet]; [cited 2020 Apr 16]. Available from: https://pesquisas.face.ufmg.br/gppd/wp-content/uploads/sites/24/2020/03/30.03-previsões-de-infeccção-pelo-covid-19-aplicação-para-os-estados-brasileiros.pdf
https://pesquisas.face.ufmg.br/gppd/wp-c...

In addition, flattening the curve over time may be related to an increase in the number of infected, recovered or immune individuals, a trend already observed in other viral infections. That is, as the epidemic increases and more people are ill and recover, R0 tends to decrease.

When analyzing the behavior of estimated data and simulations with a reduced peak, the focus is on predicting urgent hospital needs for the period and reducing the number of deaths. However, the initial protocol of the Ministry of Health of Brazil must be considered, according to which only hospitalized patients were tested.2323. Brasil, 2020. Boletim epidemiológico 05, Centro de Operações de Emergências em Saúde Pública, COVID-19, Ministério da Saúde, 13/03/2020 [cited 2020 Apr 27]. Available from: https://www.saude.gov.br/images/pdf/2020/marco/24/03--errata---boletim-epidemiologico-05.pdf
https://www.saude.gov.br/images/pdf/2020...
In view of the above, the number of confirmed cases would represent only 30% of the actual infected individuals, as demonstrated in another study.2424. Wynants L, Calster BV, Bonten MMJ, Collins GS, Debray TPA, Vos MD, et al. Prediction models for diagnosis and prognosis of covid-19 infection: systematic review and critical appraisal. BMJ [Internet]. 2020 [cited 2020 Apr 16]; 369:m1328. Available from: https://dx.doi.org/10.1136/bmj.m1328
https://dx.doi.org/10.1136/bmj.m1328...

The limitations of this study refer to estimates based on a simplified traditional epidemiological model (SIR), assuming the behavior of cases during the first days and the estimated population of the municipality. In addition, the scenarios were developed using the parameters of proportion of mild, moderate, severe, and critical cases, hospital admissions and deaths present in the literature until the moment of publication. A systematic review with a critical analysis of the predictive models for COVID-19 published between January and March 2020 analyzed the results of 10 predictive models for progression and mortality. The main factors for bias in estimated numbers were: absent demonstration of equations, incomplete description of assumptions and inadequate exclusion of portions of the population.2424. Wynants L, Calster BV, Bonten MMJ, Collins GS, Debray TPA, Vos MD, et al. Prediction models for diagnosis and prognosis of covid-19 infection: systematic review and critical appraisal. BMJ [Internet]. 2020 [cited 2020 Apr 16]; 369:m1328. Available from: https://dx.doi.org/10.1136/bmj.m1328
https://dx.doi.org/10.1136/bmj.m1328...
In this study, the equations were summarized in the section of material and methods and presented in detail in the references; all parameters of the simulations are presented in the supplementary material 1 , and no portion of the population was excluded for the calculations. However, despite considering the limitations and probability of bias in the estimates presented in this study, this study is relevant due to its (not exact) approximation to the reality that may be experienced in other municipalities in Brazil.

To elucidate the assumptions that could change the scenarios, Figure 4 contains the facts used in the calculations of this study and an illustrative summary of the assumptions causing changes in the scenarios. The effects of protective measures and their impact on the current pandemic in Brazil were also discussed in a recently published documentary study showing the evolution of curves in Brazil and other countries with and without protective measures.2525. Oliveira AC, Lucas TC, Iquiapaza RA. O que a pandemia da COVID-19 tem nos ensinado sobre adoção de medidas de precaução?. Texto Contexto Enferm [Internet]. 2020 [cited 2020 May 04]; 29:e20200106. Available from: https://dx.doi.org/10.1590/1980-265x-tce-2020-0106
https://dx.doi.org/10.1590/1980-265x-tce...

Figure 4 -
Summary of facts used in the calculations of the predictive models for COVID-19 incidence and factors that could influence and change the scenarios anticipated for Ponta Grossa, Paraná, Brazil.

CONCLUSIONS

The total incidence of COVID-19 in Ponta Grossa in the first 44 days of the disease was lower than in large centers and, despite this, the increase in cases can be explained by an exponential curve, just like other locations. According to this study, the peak period will occur around 120 days after the first case and almost four months after the beginning of the epidemic in Brazil. At the end of the period, R0 of the municipality will be 1.07, also below other locations in the world, and at the peak period it will infect about 0.23% of the population.

A reduced number of cases due to protective measures can reduce the simultaneous burden of care for critically ill and severe patients; however, the model takes into account many facts and, when a change occurs, there may be a variation in the estimated scenarios. Anyway, this study showed that the municipality will have more time to prepare and implement actions, and it can be extended to other municipalities of the same size.

REFERENCES

NOTES

  • ACKNOWLEDGEMENTS

    The authors would like to thank the State University of Ponta Grossa, Fundação Municipal de Saúde - Prefeitura Municipal de Ponta Grossa, and AAC&T Consultoria em Pesquisa LTDA.

Publication Dates

  • Publication in this collection
    31 July 2020
  • Date of issue
    2020

History

  • Received
    06 May 2020
  • Accepted
    15 June 2020
Universidade Federal de Santa Catarina, Programa de Pós Graduação em Enfermagem Campus Universitário Trindade, 88040-970 Florianópolis - Santa Catarina - Brasil, Tel.: (55 48) 3721-4915 / (55 48) 3721-9043 - Florianópolis - SC - Brazil
E-mail: textoecontexto@contato.ufsc.br