Predicting the number of total COVID-19 cases and deaths in Brazil by the Gompertz model.

In this work, we estimate the total number of infected and deaths by COVID-19 in Brazil and two Brazilian States (Rio de Janeiro and Sao Paulo). To obtain the unknown data, we use an iterative method in the Gompertz model, whose formulation is well known in the field of biology. Based on data collected from the Ministry of Health from February 26, 2020, to July 2, 2020, we predict, from July 3 to 9 and at the end of the epidemic, the number of infected and killed for the whole country and for the Brazilian states of Sao Paulo and Rio de Janeiro. We estimate, until July 9, 2020, a total of 1,709,755 cases and 65,384 deaths in Brazil, 331,718 cases and 15,621 deaths in Sao Paulo, 134,454 cases and 11,574 deaths in Rio de Janeiro. We also estimate the basic reproduction number R 0 for Brazil and its two states. The estimated values ( R 0 ) were 1.3, 1.3, and 1.4 for Brazil, Sao Paulo, and Rio de Janeiro, respectively. The results show a good fit between the observed data and those obtained by the Gompertz. The proposed methodology can also be applied to other countries and Brazilian states, and we provide an executable as well as the source code for a straightforward application of the method on such data. © Springer Nature B.V. 2020.

Introduction

The World Health Organization declared the coronavirus disease 2019 a pandemic on March 11, pointing to the over 118,000 cases in over 110 countries and territories around the world at that time [18]. The COVID-19 pandemic in Brazil began on February 26, 2020, when a man from Sao Paulo who returned from Italy tested positive for the virus. The first case of COVID-19 in Rio de Janeiro was confirmed on March 5, 2020. Nowadays, Brazil is considered the epidemic center of Latin America, occupying the second place in the total number of cases and, more recently, in the total number of deaths. Currently (July 9, 2020), Sao Paulo and Rio de Janeiro are the states with the highest number of deaths by the new coronavirus in Brazil, according to data from the Brazilian Health Ministry.

In the past few months, a considerable number of studies related to the evolution of COVID-19 in the world have been submitted and published. In the following, we describe some of these works.

Ahmadi et al. [1] developed mathematical models to predict the number of COVID-19 cases in Iran from April 3, 2020, to May 13, 2020. The unknown parameters in these models were estimated by running the fminsearch, a MATLAB function, which is a least-squares algorithm. Torrealba et al. [22] analyzed the modeling and prediction of COVID-19 in Mexico, from an initial approximation, and using the Gauss–Newton algorithm, the authors estimated parameters in the Gompertz and Logistic models. Articles [12, 15, 16] estimated the number of total COVID-19 cases and deaths in world using the Gompertz model. In [11], they analyzed the number of deaths by COVID-19 with social distancing in Brazil, Sao Paulo. Using the SEIR model, the authors recommend temporary lockdowns, however, great the economic and social costs. In [2], the authors studied the impact on the evolution of cases in Rio de Janeiro, Brazil, considering a Susceptible-Infectious-Quarantined-Recovered (SIQR) model with containment.

In this work, we model the near-future trajectory of the cumulative number of infections and deaths by COVID-19 for Brazil and two Brazilian states, given the cumulative number of infected cases and deaths from February 26, 2020 to July 2, 2020. Data were obtained from the Brazilian Ministry of Health until July 2 (https://covid.saude.gov.br), and we considered the cumulative reports from the date on which the first case was notified in Brazil and in each analyzed state. We also estimate the basic reproduction number , which represents the average number of secondary infections generated by each infected person. To estimate the unknown parameters, we use the Gompertz model and propose a gradient type iterative method (Minimal error method).

The Gompertz model is one of the particular cases of the Richards model, as well as Brody, negative exponential, logistic and von Bertalanffy models, see [20, 21].

The growth functions can be grouped in three main categories: those without inflection point (Brody and negative models), those with sigmoidal shape and a fixed inflection point (Gompertz, logistic, and von Bertalanffy models), and those with a flexible inflection point (Richards model). The logistic, Gompertz, and von Bertalanffy models exhibit inflection points at about 50, 37, and 30% of the upper asymptote, respectively, which means that the Gompertz and von Bertalanffy processes are asymmetric, whereas the logistic is a symmetric process, see [19].

The Gompertz and the logistic models are the most frequently used sigmoid functions, and the literature on these models is extensive [3, 21]. In general, the cumulative number of deaths and cases by COVID-19 presents an asymmetrical sigmoidal growth curve. Therefore, using an inappropriate growth curve can have a substantial impact on forecasting [3].

The Gompertz model [5] was proposed by Benjamin Gompertz in 1825. Since then, this exponential model has been used to describe growth in plants, animals, bacteria, and cancer cells [22]. The Gompertz differential equation has the following form,where t represents time; N(t) is the cumulative population size at time t; r the intrinsic growth rate of model; K is the maximum value of model (N) when t goes to infinity; represents the initial time; and is the initial population or condition.

The analytic solution from (1) iswhere .

The turning point is the time at which the rate of accumulation changes from increasing to decreasing or vice versa and can be easily located by finding the inflection point of the epidemic curve, that is, the moment at which the trajectory begins to decline. Clearly, this quantity is of epidemiological importance, indicating either the beginning (i.e., the moment of acceleration after deceleration) or end (i.e., the moment of deceleration after acceleration) of a phase [6]. From Eq. (2), it is trivial to show that the inflection point is given byAt the inflection point, the number of infected or killed is given by .

In this paper, we consider the discrete function, from Eq. (2) and for ,where i represents time in days and is the cumulative number of infected cases (or deaths) at day i.

In this paper, we consider that equality (3) is the inflection point of Eq. (4).

We denote . The inverse problem is to estimate , from (4), given .

Method: inverse problem to the Gompertz model

Here, we consider an application of the minimal error method (MEM) to the inverse problem at the Gompertz model. The MEM is a variant of the conjugate gradient method [4, 10, 17], and this approach was used to estimate unknown parameters in a computational neuroscience models [24, 25]. Also, [14] solves the Cauchy problem in linear elasticity with MEM.

Knowing the cumulative number of infections or deaths from COVID-19 , we want to determine assuming that Eq. (4) holds.

The transpose of a vector is denoted by . Let the nonlinear operator defined bywhere are multivariable functions with real value and solves (4).

To obtain an approximation for , given and , we used the minimal error iterationwhere is the directional derivative of computed at , and is its adjoint. We also definewhere solves (4), replacing by . Note that .

It is possible to show that, under certain conditions (we assume that is the case), converges to a solution of ; see [9, Theorem 3.21].

From Eq. (5), we obtain an approximation for , but the adjoint is not known. In this paper, we obtain the unknown operator (see Appendix) together with the following iterationwhereThe regression coefficient is used to evaluate the fitting ability of various methods and can be obtained by the following equation,whereis the average of the cumulative confirmed COVID-19 cases (or deaths).

Fig. 1

Plots for Brazil. Subplot a shows extrapolated infection trajectory from February 26 to July 9. Subplot b presents the number of COVID-19 deaths from March 17 to July 9. The blue curve shows the exact data. The red line represents the estimate of the blue curve with the Gompertz model. The green asterisk is the inflection point

Remark 1

From differential Eq. (1) and since is known, we consider the following initial approximation for

The numerical scheme for our method poposed is in Algorithm 1.

Remark 2

The basic reproduction number is the most important parameter to analyze any epidemic model for any disease. In a practical sense, is the average number of secondary infections caused by one infected individual during his/her entire infectious period. An epidemic occurs only if a single individual can spread its illness to more than one individual . If , on average an infectious individual infects less than one person and the contagion is expected to stop spreading.

Our result can be used to compute the basic reproduction numbersee [6–8, 13], where r denotes the intrinsic growth rate in the Gompertz model and is the generation time of disease transmission. In this paper, we consider .

Results

In this section, we will provide the results obtained for the whole country (Brazil) and for the Brazilian states of Sao Paulo and Rio de Janeiro. Table 1 describes the estimated values for the Gompertz model and also presents the regression coefficient, the inflection point, and the basic reproduction number.

Table 1

Estimates of parameters unknown

	K	r	R^2	Peak day(t_i.p)	R_0
Brazil
Cases	4,618,893	0.0204	0.999	July 9	1.3
Deaths	91,440	0.0309	0.999	June 4	–
Sao Paulo
Cases	820,202	0.0202	0.997	July 4	1.3
Deaths	20,802	0.0311	0.994	May 30	–
Rio de Janeiro
Cases	255,854	0.0235	0.998	June 20	1.4
Deaths	17,585	0.0275	0.998	June 7	–

Second and third columns are estimated values for the Gompertz model (4). Note that K is the number of infected (or killed) at the end of the pandemic. The fourth column describes the regression coefficient according to Eq. (8). The fifth column presents the date when the number of cases or deaths reaches its peak (i.e., the moment of deceleration after acceleration), see Eq. (3). Finally, the sixth column is the basic reproduction number, see Remark (2)

The estimated daily values for the cumulative number of confirmed cases and deaths, from July 3 to 9, are illustrated in Table 2. Figures 1, 2, 3 illustrate the relationship between observed values and those predicted by the Gompertz model. The tables and figures show the efficiency of the proposed iterative method.

Table 2

Estimated number of cases and deaths by days in Brazil and its two states

Date	Total cases		Total deaths
	Prediction	Real	Prediction	Real
Brazil
July 3	1,501,970	1,539,081	61,059	63,174
July 4	1,536,474	1,577,004	61,815	64,265
July 5	1,571,050	1,603,055	62,557	64,867
July 6	1,605,681	1,623,284	63,285	65,487
July 7	1,640,353	1,668,589	63,998	66,741
July 8	1,675,049	1,713,160	64,698	67,964
July 9	1,709,755	1,755,779	65,384	69,184
Sao Paulo
July 3	295,141	310,702	14,729	15,694
July 4	301,245	312,530	14,886	15,996
July 5	307,349	320,179	15,039	16,078
July 6	313,451	323,070	15,190	16,134
July 7	319,548	332,708	15,337	16,475
July 8	325,637	341,365	15,480	16,788
July 9	331,718	349,715	15,621	17,118
Rio de Janeiro
July 3	121,965	118,956	10,738	10,500
July 4	124,083	120,440	10,882	10,624
July 5	126,188	121,292	11,025	10,667
July 6	128,278	121,879	11,165	10,698
July 7	130,353	124,086	11,304	10,881
July 8	132,412	126,329	11,440	10,970
July 9	134,454	128,324	11,574	11,115

The cases and deaths are cumulative. The mean absolute percentage error (MAPE) values for the cumulative number of cases and deaths in Brazil are 2.0% and 3.9%, respectively. Considering the data from Sao Paulo, the MAPEs for the cumulative number of cases and deaths are 4.3% and 6.9%, respectively. Finally, the MAPEs for the cumulative number of cases and deaths in Rio de Janeiro are 4.4% and 3.4%, respectively

Fig. 2

Results for Sao Paulo. Subplots a and b are the cumulative number of infected individuals and deaths from February 26 to July 17 and March 17 to July 17, respectively. See Figure 1 for a description of each variable

Fig. 3

Plots for Rio de Janeiro. Subplots a and b show the cumulative number of infected individuals and deaths from March 5 to July 17 and March 19 to July 17, respectively. See Figure 1 for a description of each variable

Conclusion

In this article, we propose an iterative method to estimate the unknown parameters in the Gompertz model. The results of the modeling show a good fit between the estimated and the observed data. We obtained an estimate of for the cumulative number of infected and killed in Brazil. For Sao Paulo, we obtained an approximation of and for the cumulative number of infected and deaths, respectively. For Rio de Janeiro, the estimate was for the cumulative number of infected and deaths.

The basic reproduction number ( or the average number of infections caused by one typically infectious individual ) depends on the growth rate (r), which was estimated in this work, and the duration of infectiousness . Considering , the basic reproduction number was 1.3, 1.3, and 1.4 for Brazil, Sao Paulo, and Rio de Janeiro, respectively. Applying our proposed methodology for China, we obtained , given the number of infected cases accumulated from February 2 to July 2, based on data collected from the world health organization. We did the same for Italy and obtained . Many works obtained these values for China and Italy. Therefore, our estimates for the basic reproduction number are acceptable.

The inflection point of the curve provides vital information about changing trends in the epidemic and may possibly indicate changes in intervention and control. In this paper, the inflection points for the cumulative number of infected cases were July 9, July 4, and June 21 for Brazil, Sao Paulo, and Rio de Janeiro, respectively.

Regarding the methodology used, there are several methods for estimating or determining parameters in a mathematical model. We can divide these methods into two groups, iterative and non-iterative methods. In this paper, we used a gradient type iterative method. This method can recover unknown parameters with a non-uniform distribution (non-constant functions) given noisy data. On the other hand, like all iterative methods, it depends on an initial estimate. If this initial approximation is far from the solution, the iteration will diverge. Although it is difficult to determine an initial estimate for some problems, in this work we found an appropriate initial guess, see Remark 1.