A comparative analysis between a SIRD compartmental model and the Richards growth model ======================================================================================= * Antônio M. S. Macêdo * Arthur A. Brum * Gerson C. Duarte-Filho * Francisco A. G. Almeida * Raydonal Ospina * Giovani L. Vasconcelos ## Abstract We propose a compartmental SIRD model with time-dependent parameters that can be used to give epidemiological interpretations to the phenomenological parameters of the Richards growth model. We illustrate the use of the map between these two models by fitting the fatality curves of the COVID-19 epidemic data in Italy, Germany, Sweden, Netherlands, Cuba, and Japan. Keywords * COVID-19 * Fatality curve * SIRD model * Richards growth model * Intervention strategies ## 1. Background The pandemic of the novel coronavirus disease (COVID-19) has created a major worldwide sanitary crisis [1, 2]. Developing a proper understanding of the dynamics of the COVID-19 epidemic curves is an ongoing challenge. In modeling epidemics, in general, compartmental models [3] have been to some extent the tool of choice. However, in the particular case of the COVID-19 epidemic, standard compartmental models, such as SIR, SEIR, and SIRD, have so far failed to produce a good description of the empirical data, despite a great amount of intensive work [4, 5, 6, 7, 8, 9, 10, 11]. In this context, phenomenological growth models have met with some success, particularly in the description of cumulative death curves [12, 13, 14]. The recent discovery, within the context of a generalized growth model known as the beta logistic model [15], of a slow, power-law approach towards the plateau in the final stage of the epidemic curves is another remarkable example of this qualitative success. Growth models, however, have the drawback that their parameters may not be easily interpreted in terms of standard epidemiological concepts [16], as can the parameters of the usual compartmental models. As a concrete example, consider the transmission rate parameter *β* of the SIR model [3]. It can be easily interpreted as the probability that a contact between a susceptible individual and an infective one leads to a transmission of the pathogen, times the number of contacts per day. Although the value of *β* cannot be measured directly in a model independent way, and it is probably not even constant in the COVID-19 epidemic curves, the epidemiological meaning of the parameter is nonetheless easy to grasp conceptually. As a result, models that incorporate such parameters in their basic equations are sometimes regarded as “more epidemiological,” so to speak, than others that do not use similar parameters. This state of affairs creates a somewhat paradoxical scenario, in which we have, on the one hand, the striking empirical success of phenomenological growth models sometimes being downplayed, owing to the lack of a simple epidemiological picture of the underlying mechanism [16], and, on the other hand, the failure of traditional epidemiological compartmental models to produce good quantitative agreement with the empirical COVID-19 data. A glaring instance of the inadequacy of standard compartmental models for the COVID-19 epidemic is their inability to predict the power-law behavior often seen in the early-growth regime as well as in the saturation phase of the accumulated death curves—a feature that is well captured by growth models [15], as already mentioned. It is clear that a kind of compromise is highly desirable, in which we get the benefits of the accuracy of the growth models in describing the epidemic, along with a reasonable epidemiological interpretation of their free parameters. An attempt in this direction was presented by Wang [16], where an approximate map between the Richards growth model [16] and the accumulated number of cases of a SIR model was proposed. The two free parameters of the Richards model were expressed as a function of the epidemiological based parameters of the SIR model. Here we improve on this analysis in two ways: (i) we extend the SIR model to a SIRD model by incorporating the deceased compartment, which is then used as the basis for the map onto the Richards model; (ii) the parameters of the SIRD model are allowed to have a time dependence, which is crucial to gain some efficacy in describing realistic cumulative epidemic curves of COVID-19. ## 2. Data It is in general very hard to estimate the actual number of infected people within a given population, simply because a large proportion of infections go undetected. This happens largely because many carriers of the coronavirus are either asymptomatic or develop only mild symptoms, which in turn makes the number of confirmed cases for COVID-19 a poor proxy for the actual number of infections. This issue is well known in the literature and referred to as the “under-reporting problem” [17, 18]. With this in mind, and in the absence of more reliable estimates for the number of infected cases, we shall here focus our analysis on the fatality curves, defined as the cumulative number of deaths as a function of time. In the present study we considered the mortality data of COVID-19 from the following countries: Italy, Germany, Sweden, Netherlands, Cuba, and Japan. The data used here were obtained from the database made publicly available by the Johns Hopkins University [19], which lists in automated fashion the number of the confirmed cases and deaths attributed to COVID-19 per country. We have used data up to July 30, 2020. ## 3. Methods ### 3.1. The Richards Growth Model The time evolution of the number of cases/deaths in an epidemy can be modelled by means of the Richards model (RM), defined by the following ordinary differential equation [20, 21, 22]: ![Formula][1] where *C*(*t*) is the cumulative number of cases/deaths at time *t, r* is the growth rate at the early stage, *K* is the final epidemic size, and the parameter *α* measures the asymmetry with respect to the s-shaped curve of the standard logistic model, which is recovered for *α* = 1. In the present paper we shall apply the RM to the fatality curves of COVID-19, so that *C*(*t*) will always represent the cumulative numbers of deaths at time *t*, where *t* will be counted in days from the first death. Equation (3.1) must be supplemented with a boundary condition, which can be either the initial time, *t* = 0, or the inflection point, *t* = *t**c*, defined by the condition *C**′′*(*t**c*) = 0, where *C**′′*(*t*) = *d*2*C*(*t*)*/dt*2. A direct integration of (3.1) yields the following explicit formula: ![Formula][2] which will be the basis of our analysis. ### 3.2. SIRD model with constant parameters We start by recalling the standard Susceptible (*S*)-Infected (*I*)-Recovered (*R*)-Deceased (*D*) epidemiological model ![Formula][3] ![Formula][4] ![Formula][5] ![Formula][6] where *S*(*t*), *I*(*t*), *R*(*t*), and *D*(*t*) are the number of individuals at time *t* in the classes of susceptible, infected, recovered, and deceased respectively; whereas *N* is the total number of individuals in the population. i.e., *N* = *S*(*t*)+*I*(*t*)+*R*(*t*)+*D*(*t*). The initial values are chosen to be *S*(0) = *s*, *I*(0) = *i*, with *s* + *i* = *N*, and *R*(0) = 0 = *D*(0). The parameters *γ*1 and *γ*2 are the rates at which infected individual becomes recovered or deceased, respectively. We then consider the following modified SIRD model, where in (3.3) and (3.4) we replace *N* with only the partial population in the *S* and *I* compartments, which takes into account the fact that the recovered (assuming they become immune) and the deceased cannot contribute to the transmission. We thus find ![Formula][7] ![Formula][8] ![Formula][9] ![Formula][10] A fundamental quantity in epidemiology is the basic reproductive ratio, *R*, which is defined as the expected number of secondary infections caused by an infected individual during the period she (or he) is infectious in a population consisting solely of susceptible individuals. In this model, *R* can be calculated using the next generation method [23, 24] and is given by ![Formula][11] Next, we define *y*(*t*) = *S*(*t*) + *I*(*t*) and divide (3.8) by (3.7) to obtain ![Formula][12] Integrating both sides of (3.12), and inserting the result into (3.7), yields a growth equation of the Richards type: ![Formula][13] where *α* = 1 #x2212; 1*/R* and ![Graphic][14]. We now seek to approximate the curve of accumulated death *D*(*t*), obtained from the SIRD model, with the Richards function *C*(*t*), as defined in (3.2). To this end, we first impose the boundary conditions *K* = *D*(*∞*) and *t**c* = *t**i*, where *t**i* is the inflection point of *D*(*t*). By definition ![Graphic][15], which implies from (3.10) that *İ*(*t**i*) = 0 and thus ![Formula][16] Furthermore, we require that at *t* = *t**i* both *C*(*t*) and its derivative *Ċ* (*t*) respectively match *D*(*t*) and ![Graphic][17], thus ![Formula][18] ![Formula][19] Using the condition *İ*(*t**i*) = 0 in the SIRD equations, we find ![Formula][20] ![Formula][21] Using equations (3.15) and (3.16), we finally obtain the connection between the parameters (*r, α*) of the RM and the parameters (*β, γ*1, *γ*2) of the SIRD model ![Formula][22] ![Formula][23] which are the central equations of this paper. We can estimate the precision of the above ‘map’ between the RM and the SIRD model via the relative error function: ![Formula][24] We have verified numerically that ![Formula][25] where *ε* is typically of order 0.1. A typical example of the agreement between the SIRD model, for a given set of parameters (*β, γ*1, *γ*2), and the RM with the parameters obtained from the map described by (3.19) and (3.20), is illustrated in Fig. 1. In Fig. 2 we show the simple monotonic dependence of the Richards parameters (*r, α*) on the parameter *β* of the SIRD model, for the biologically relevant interval 0≤*r, α ≤*1. We also show, for comparison, the behavior of the basic reproduction number *R*. ![Figura 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/06/2020.08.04.20168120/F1.medium.gif) [Figura 1:](http://medrxiv.org/content/early/2020/08/06/2020.08.04.20168120/F1) Figura 1: Illustrative picture (left panel) of the map between the Richards model and the SIRD model with constant parameters. The relative error function between the two models is shown in the right panel. ![Figura 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/06/2020.08.04.20168120/F2.medium.gif) [Figura 2:](http://medrxiv.org/content/early/2020/08/06/2020.08.04.20168120/F2) Figura 2: Behavior of the Richards parameters (*r, α*) and *R* as a function of the parameter *β* of the SIRD model. ### 3.3. SIRD model with time-dependent parameters The SIRD model with constant parameters proved to be insufficient to accommodate properly the human intervention biased dynamics of the COVID-19 epidemics. The simplest solution to this problem is to allow the epidemiological parameter *β* to change in time according to the simple exponential decay function [4] ![Formula][26] where *τ* is the starting time of the intervention and *τ*1 is the average duration of interventions. Here *β* is the initial transmission rate of the pathogen and the product *β**β*1 represents the transmission rate at the end of the epidemic. Remarkably, the central map equations, (3.19) and (3.20), are still valid, although *t**i* is no longer given by (3.14) and should be determined from the maximum of the curve *I*(*t*) obtained from the numerical solution of the SIRD equations, with the parameter *β* replaced by the function *β*(*t*). ## 4. Applications and Discussion In Fig. 3 we demonstrate some applications of the SIRD-RM map by showing the cumulative number of deaths (red circles) attributed to COVID-19 for the following countries: Italy, Germany, Netherlands, Sweden, Japan and Cuba. In all figures shown, the continuous (black) curve is the numerical fit to the empirical data, as produced by the SIRD model with the time-dependent parameter *β*(*t*) given in (3.23), and the dashed (bright green) curve is the corresponding theoretical curve predicted by the RM, with the parameters as obtained from the map (3.19) and (3.20). The statistical fits were performed using the Levenberg-Marquardt algorithm [25], as implemented by the lmfit Python package [26], to solve the corresponding non-linear least square optimization problem. In other words, the *lmfit* package was applied to each empirical dataset to determine the parameters (*β*, *β*1, *γ*1, *γ*2, *τ*, *τ*1) of the SIRD model described in Secs. 3.2 and 3.3. ![Figura 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/08/06/2020.08.04.20168120/F3.medium.gif) [Figura 3:](http://medrxiv.org/content/early/2020/08/06/2020.08.04.20168120/F3) Figura 3: Cumulative number of deaths (red circles) attributed to COVID-19, up to July 30, 2020, for (a) Netherlands, (b) Italy, (c) Germany, (d) Cuba, (e) Sweden, and (f) Japan. The solid black curves are the best fits by the SIRD model with a time-dependent *β*(*t*), where the parameter estimates are given in the inset. The bright green dashes curve is the theoretical curve obtained from the Richards model, with the parameters computed from the numerically determined parameters of the SIRD model via the map (3.19) and (3.20). One can see from Fig. 3 that the agreement between the RM and the SIRD model is very good in all cases considered, which satisfactorily validates the map between these two models. This result thus shows, quite convincingly, that the parameters of the Richards model do bear a direct relationship to epidemiological parameters, as represented, say, in compartmental models of the SIRD type. Although the interpretation of the Richards parameters (*r, α*) are less obvious, in that they involve a nonlinear relation with the probability rates used in compartmental models, these parameters should nonetheless be regarded as bonafide epidemiological parameters. Furthermore, it is important to emphasize the flexibility of the RM: this model, which has only two time-independent parameters, is equivalent (in the sense of the map discussed above) to a SIRD model with *time dependent* parameters. In other words, the two constant parameters of the RM are sufficient to characterize, to a rather good extent, the entire evolution of the COVID-19 epidemic in a given location. It is worth pointing out that the discovery of power-law behaviors in the early-growth regime as well as in the saturation phase of the accumulated death curves, both of which are well described by the beta logistic model (BLM) [15], brings about the challenge to accommodate power laws into a compartmental model. A preliminary analysis [15] shows that substantial modifications in the SIRD equations may be required to achieve power law-behavior in the short-time and long-time regimes of the epidemic curves. The possibility of a map between the BLM and a modified SIRD model with time-dependent parameters is currently under investigation. ## 5. Conclusion The present paper provides a map between a SIRD model with time dependent parameters and the Richards growth model. We illustrated the use of this map by fitting the fatality curves of the COVID-19 epidemics data for Italy, Germany, Sweden, Netherlands, Cuba and Japan. The results presented here are relevant in that they showcase the fact that phenomenological growth models, such as the Richards model, are valid epidemiological models not only because they can successfully describe the empirical data but also because they capture, in an effective way, the underlying dynamics of an infectious disease. In this sense, the free parameters of growth models acquire a biological meaning to the extent that they can be put in correspondence (albeit not a simple one) with parameters of compartmental model, which have a more direct epidemiological interpretation. ## Data Availability The data for countries as well as for states/cities in the US are fetched from the COVID-19 Open Data Repository by CSSE/Johns Hopkins University ([https://systems.jhu.edu](https://systems.jhu.edu) ## Acknowledgements This work was partially supported by the National Council for Scientific and Technological Development (CNPq) in Brazil, through the grants Nos. 303772/2017-4 (GLV), 312612/2019-2 (AMSM), and 305305/2019-0 (RO). AAB also thanks CNPq for its support through a PhD Fellowship (grant No. 167348/2018-3). RO thanks D.A.D.O. and CASTLab. The authors have declared no competing interest. ## Footnotes * antonio.smacedo{at}ufpe.br, arthurbrum{at}df.ufpe.br, gcdf{at}ufs.br, falmeida{at}ufs.br, raydonal{at}de.ufpe.br, giovani.vasconcelos{at}ufpr.br * Received August 4, 2020. * Revision received August 4, 2020. * Accepted August 6, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## Referências 1. [1].WHO, “World Health Organization (WHO). Director-General ‘s opening remarks at the media briefing on COVID-19 - 30 July 2020.” [https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19\---|30-july-2020](https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19\---|30-july-2020), 2020. Accessed: 2020-07-30. 2. [2].WHO, “World Health Organization (WHO). Director-General ‘s opening remarks at the media briefing on COVID-19 - 30 March 2020.” [https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19-30-march-2020](https://www.who.int/dg/speeches/detail/who-director-general-s-opening-remarks-at-the-media-briefing-on-covid-19-30-march-2020), 2020. Accessed: 2020-03-30. 3. [3]. W. O. Kermack and A. G. McKendrick, “A contribution to the mathematical theory of epidemics,” Proceedings of the royal society of london. Series A, Containing papers of a mathematical and physical character, vol. 115, no. 772, pp. 700–721, 1927. 4. [4]. D. Caccavo, “Chinese and Italian COVID-19 outbreaks can be correctly described by a modified SIRD model,” medRxiv medRxiv:2020.03.19.20039388, 2020. 5. [5]. M. Fukui and C. Furukawa, “Power laws in superspreading events: Evidence from coronavirus outbreaks and implications for SIR models,” medRxiv medRxiv:2020.06.11.20128058, 2020. 6. [6]. U. Tirnakli and C. Tsallis, “Epidemiological model with anomalous kinetics - the Covid-19 pandemics,” medRxiv medRxiv:2020.06.24.20139287, 2020. 7. [7]. A. Atkeson, “What will be the economic impact of covid-19 in the us? rough estimates of disease scenarios,” tech. rep., National Bureau of Economic Research, 2020. 8. [8]. E. B. Postnikov, “Estimation of COVID-19 dynamics “on a back-of-envelope”: Does the simplest SIR model provide quantitative parameters and predictions?,” Chaos, Solitons & Fractals, vol. 135, p. 109841, 2020. 9. [9]. Z. Yang, Z. Zeng, K. Wang, S.-S. Wong, W. Liang, M. Zanin, P. Liu, X. Cao, Z. Gao, Z. Mai, et al., “Modified SEIR and AI prediction of the epidemics trend of COVID-19 in China under public health interventions,” Journal of Thoracic Disease, vol. 12, no. 3, p. 165, 2020. 10. [10]. S. Zhao and H. Chen, “Modeling the epidemic dynamics and control of COVID-19 outbreak in China,” Quantitative Biology, pp. 1–9, 2020. 11. [11]. I. Cooper, A. Mondal, and C. G. Antonopoulos, “A SIR model assumption for the spread of COVID-19 in different communities,” Chaos, Solitons & Fractals, p. 110057, 2020. 12. [12]. K. Wu, D. Darcet, Q. Wang, and D. Sornette, “Generalized logistic growth modeling of the COVID-19 outbreak in 29 provinces in China and in the rest of the world.” [https://www.medrxiv.org/content/early/2020/03/16/2020.03.11.20034363.full.pdf](https://www.medrxiv.org/content/early/2020/03/16/2020.03.11.20034363.full.pdf), 2020. 13. [13]. G. L. Vasconcelos, A. M. Macêdo, R. Ospina, F. A. Almeida, G. C. Duarte-Filho, A. A. Brum, and I. C. Souza, “Modelling fatality curves of COVID-19 and the effectiveness of intervention strategies,” PeerJ, vol. 8, p. e9421, June 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7717/peerj.9421&link_type=DOI) 14. [14]. M. Gillman and N. Crokidakis, “Dynamics and future of SARS-CoV-2 in the human host,” medRxiv medRxiv 2020.07.14.20153270, 2020. 15. [15]. G. L. Vasconcelos, A. M. Macêdo, G. C. Duarte-Filho, A. A. Araújo, R. Ospina, and F. A. Almeida, “Complexity signatures in the COVID-19 epidemic: power law behaviour in the saturation regime of fatality curves,” medRxiv medRxiv 2020.07.12.20152140, 2020. 16. [16]. X.-S. Wang, J. Wu, and Y. Yang, “Richards model revisited: Validation by and application to infection dynamics,” Journal of Theoretical Biology, vol. 313, pp. 12–19, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2012.07.024&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22889641&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F06%2F2020.08.04.20168120.atom) 17. [17]. R. Veiga, R. Murta, and R. Vicente, “Age-structured estimation of COVID-19 ICU demand from low quality data,” arXiv preprint arXiv:2006.06530, 2020. 18. [18]. M. F. do Prado, B. Brandão, S. L. Bastos, I. T. Peres, A. d. A. B. da Silva, L. F. Dantas, F. Araújo, and F. A. Bozza, “Analysis of COVID-19 under-reporting in Brazil,” Revista Brasileira de terapia intensiva, vol. 32, no. 2, pp. 224–228, 2020. 19. [19].JHU, “Coronavirus COVID-19 Global Cases by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University (JHU).” [https://coronavirus.jhu.edu/map.html](https://coronavirus.jhu.edu/map.html), 2020. Accessed: 2020-07-30. 20. [20]. F. Richards, “A flexible growth function for empirical use,” Journal of experimental Botany, vol. 10, no. 2, pp. 290–301, 1959. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jxb/10.2.290&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1959XE91800011&link_type=ISI) 21. [21]. X.-S. Wang, J. Wu, and Y. Yang, “Richards model revisited: Validation by and application to infection dynamics,” Journal of Theoretical Biology, vol. 313, pp. 12–19, 2012. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2012.07.024&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22889641&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F06%2F2020.08.04.20168120.atom) 22. [22]. Y.-H. Hsieh, “Richards model: a simple procedure for real-time prediction of outbreak severity,” in Modeling and dynamics of infectious diseases, pp. 216–236, World Scientific, 2009. 23. [23]. O. Diekmann, J. A. P. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations,” Journal of mathematical biology, vol. 28, no. 4, pp. 365–382, 1990. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00178324&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2117040&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F08%2F06%2F2020.08.04.20168120.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990DG35000001&link_type=ISI) 24. [24]. O. Diekmann, J. Heesterbeek, and M. G. Roberts, “The construction of next-generation matrices for compartmental epidemic models,” Journal of the Royal Society Interface, vol. 7, no. 47, pp. 873–885, 2010. 25. [25]. J.J. Moré, “The Levenberg-Marquardt algorithm: implementation and theory,” in Numerical analysis, pp. 105–116, Springer, 1978. 26. [26]. M. Newville, T. Stensitzki, D. Allen, and A. Ingargiola, “Non-linear least-squares minimization and curve-fitting for Python,” Chicago, IL, 2015. [1]: /embed/graphic-1.gif [2]: /embed/graphic-2.gif [3]: /embed/graphic-3.gif [4]: /embed/graphic-4.gif [5]: /embed/graphic-5.gif [6]: /embed/graphic-6.gif [7]: /embed/graphic-7.gif [8]: /embed/graphic-8.gif [9]: /embed/graphic-9.gif [10]: /embed/graphic-10.gif [11]: /embed/graphic-11.gif [12]: /embed/graphic-12.gif [13]: /embed/graphic-13.gif [14]: /embed/inline-graphic-1.gif [15]: /embed/inline-graphic-2.gif [16]: /embed/graphic-14.gif [17]: /embed/inline-graphic-3.gif [18]: /embed/graphic-15.gif [19]: /embed/graphic-16.gif [20]: /embed/graphic-17.gif [21]: /embed/graphic-18.gif [22]: /embed/graphic-19.gif [23]: /embed/graphic-20.gif [24]: /embed/graphic-21.gif [25]: /embed/graphic-22.gif [26]: /embed/graphic-25.gif