Mathematical model for the mitigation of the economic effects of the Covid-19 in the Democratic Republic of the Congo ===================================================================================================================== * Zirhumanana Balike Dieudonné ## Abstract A Mathematical model for the spread of Covid-19 in Democratic Republic of Congo taking into account the vulnerability of congolese economy is proposed. The reproduction number for the Covid-19 is calculated and numerical simulations are performed using Python software. A clear advice for policymakers is deduced from the forecasting of the model. Key words * Covid-19 * Infectious * Quarantine * Model * Virus ## Introduction In December 2019 a novel Coronavirus appeared in Wuhan, China [9]. The international committee for taxonomy of virus has attributed the name SRAS-Cov-2 to that disease [10], [4]. On January 30, 2020; the World Health Organization (WHO) declared it to be an epidemic of international concern [19]. Four months later, the virus was spread out worldwide, only less than 10 countries were not yet touched by the disease. This is why WHO declared it to be a pandemic since March 11, 2020 [19]. Countries affected by this outbreak have envisaged several measures to mitigate its negatives effects on their health care systems. Those measures include national lockdown and cancellation of travels to and from outside their borders [14]. Many low-income countries, including African countries, have also adopted the above measures without taking into account of the vulnerability of their economies which rely mainly on informal system [23]. On 10 March 2020, the first infected individual has been detected in Democratic Republic of Congo (DRC, [20]). On March 18, 2020; the DRC President has announced draconian measures to mitigate the circulation of the pandemic. These measures include cancellation of flights out and from Kinshasa (the capital city) to other provinces, cessation of schools and churches activities, suspension of barrooms and restaurants activities, prohibition of gatherings of more than twenty persons, etc. These measures have come in force since March 19, 2020 [11]. However, these measures have drastically impacted the national economy. A study in [21] shows that hunger could be deadlier than Covid-19 in low-income countries. This is exactly the case for the DRC. Individuals whose daily income rely on the informal system are not able to strictly respect them. While they were set to control the epidemic, this economic aspect makes them unappropriated to Congolese context. In the proposed model, there is no need of shutting down all the country, only suspicious and infected people are quarantined. The meaning of ‘suspicious’ will be precised in the following section. Since the apparition of this virus, several models have been proposed [1], [2], [3], [5], [6], [7], [8] and [13]. At my knowledge, there is no model which has been performed to help understanding the spreading of this pandemic in DRC and taking into account the economic impact of the virus on the populations. To fill in this gap and help addressing this challenge, an SEIR model (see [18], [15] and [17] for details of this model) with additional compartments is proposed in this paper. In this model, full lockdown or national quarantine is not envisaged as it may jeopardize the fragile economic system of the country. The following compartments have been added to the traditional SEIR model : *Quarantined* and *Hospitalized*. ## 1. Description of the model The population is divided into the following subgroups: * **Susceptibles (S):** This compartment comprises individuals who are not yet infected and not immunized against the disease. They are recruited at rate *θ* and transferred into the Exposed group at the rate *β*. In DRC, only one laboratory can declare positive individuals amid a Covid-19 testing. It is the National Institute of Biomedical Research (NIBR) located in Kinshasa, the capital city. In many cases, after testing, individuals have to wait for many days before they can get the result as the NIBR is far away from many provinces (for example, Kinshasa is 2000 km from Bukavu) and the transportation system is bankrupted. To prevent potential infected individuals to spread the disease, authorities should quarantine all suspicious individuals from susceptible population at rate *E*. Suspicious individuals are those who have been in contact with an infected (confirmed by testing) person or who have recently sojourned in a high risky area. Individuals quarantined are therefore not necessarily infected but they present some like Covid-19 symptoms or they have been (or are suspected to have been) in contact with a suspicious or confirmed case. * **Exposed (E):** This group contains individuals who are infected but not yet infectious. They move into the Infectious (I) group at rate *α*. * **Infectious (I)**: In this group we have individuals who are infectious, i.e those who are able to spread the disease. In general such individuals are at the onset of the illness. Thus, they will be transferred into the hospital (H) at rate *ν* and are removed at rate *η*. This removal is due either to death or recovery without being transferred at the hospital. This will particularly happen for young people who are resilient to the disease and some unbelievers who never accept that the virus is deadly and are currently propagating into the population. * **Hospitalized (H)**. Two categories of individuals are in this group: those who come from infectious group (I) and those who come from quarantine (Q). They are respectively recruited at rates *ν* and *σ*. Hospitalized individuals are likewise removed (recovery or death) at rate *µ* amid the virus infection. * **Quarantined (Q)**: Individuals who are deemed suspicious are placed in quarantine from Susceptible group (S) at rate *E* as mentioned above. Then, individuals are transferred either in hospital at rate *σ* or returned into Susceptible group at rate *γ* if they are declared negative after testing. * **Removed (R)**: Removed subjects are either recovered or deceased due to Covid-19. In addition, in all compartments, individuals are removed at the same rate *λ* due to natural death. The Fig 1 depicts the structure of this model. ![Formula][1] ![Fig 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/16/2020.12.14.20248182/F1.medium.gif) [Fig 1.](http://medrxiv.org/content/early/2020/12/16/2020.12.14.20248182/F1) Fig 1. Structure of the proposed model. In this section, an analysis of the model (system 1) will be performed. It is clear that this system under study satisfies the Lipschitz (see [12]) and theorem 1 (see [22], page 343) conditions. This guarantees the existence and the uniqueness of the solution. ### 1.1 Feasible solution and positivity of the solution The model is performed for understanding the spread of the SRAS-Cov-2 into a human population. Thus, the parameters are expected to be positive. Let Ω be the set of all feasible solution to the system 1. Each solution to a single equation lie in ℝ+. We hence have ![Formula][2] where *N* is the total size of the population under study. In this region, I want ![Formula][3] Proposition 1: The set Ω is invariant and attract all solutions in ![Graphic][4]. **Proof**. From the assumption that *N* = *S* + *E* + *I* + *H* + *Q* + *R*, I deduce that ![Formula][5] Substituting each derivative by the corresponding value from system 1 and simplifying, I get : ![Formula][6] This equation can be rewritten in the following way: ![Formula][7] Equation 4 can be solved using an integrating factor. Multiplying both sides of that equation by e *∫λdt*, we get: ![Formula][8] We finally get ![Formula][9] where C is an arbitrary constant. For *t* = 0, we have ![Graphic][10] which we substitute in equation 5 to find ![Formula][11] Thus, we have ![Graphic][12] as *t →* +*∞*, provided equation 2. This establishes that Ω is positively invariant and attracts all solutions in ![Graphic][13]. In addition, it can be proven that the solution of system 1 has a positive solution provided the initial data set (*S*, *E*, *I*, *H*, *Q*, *R*) ≥0 *∈* Ω. For example, we have ![Formula][14] That is ![Formula][15] The right side of this inequality is always positive, that is the function S is also positive and consequently its initial value. A similar computation can show that ![Formula][16] where *K**i*, (*i* ∈ {1, …, 5}) are arbitrary constants. ### 1.2 The disease free equilibrium (DFE) From the equation 1, we set *f**i* = 0, *i* = 1, …, 6 to get the new algebraic system below. ![Formula][17] From equations (a) to (f), one successively get ![Formula][18] If *E* = 0 (that is *β* = 0) and *E* = 0, then *I* = *H* = *Q* = *R* = 0 and accordingly, the DFE is ![Graphic][19]. In this case, there is no disease as no infectious individuals are within the population. If *ϵ* ≠ 0, the situation is such that suspicious are quarantined but they are neither infectious nor confirmed infected. This is consistent with the logic of the model: individuals with like Covid-19 symptoms are quarantined waiting for the testing results from the NIBR. The corresponding DFE is ![Graphic][20], 0 as *λ* never vanishes. ### 1.3 The basic reproduction number The infection components in this model are E, I, H and Q. In accordance with the notations in [16], [13] and [2], the new infection matrix *F* and the transition matrix *V* are given by ![Formula][21] The basic reproduction number of model 1 is then defined as the spectral radius of the next generation matrix *FV* −1, i.e ![Formula][22] If the inequality ℛ 1 is satisfied, then no epidemic outbreak is possible, otherwise an epidemic occurs. If *ϵ* = 0, then the fifth equation of the model 1 is removed from the infection components (that is, the model comprises only three infection components: E, I and H). Summing equations (b) and (c) of the system 6 we get the following: ![Formula][23] After substituting *I**∗* from equation 7 and simplifying, we get: ![Formula][24] as *E**∗*≠ 0. This finally yields ![Formula][25] Provided 2 and its result, I have ![Formula][26] that is ![Formula][27] Hence, I get ![Formula][28] I set ![Formula][29] the new reproduction number. This is exactly the spectral radius of the new next generation matrix obtained after ignoring the fifth equation of the system 1 in the infection compartments. Thus, *ϵ* stands for a control on the spread of the disease within the susceptible population. ### 1.4 Local stability I examine the local stability of the DFE by the mean of the jacobian matrix of the functions *f**i*, *i* = 1, …, 6 (where *f**i* are the functions in the right side of the system 1) at that point. I therefore have the jacobian below ![Formula][30] where *J*1 = *βI**∗* + *λ* + *E, J*2 = *ν* + *η* + *λ, J*3 = *σ* + *γ* + *λ*. This jacobian evaluated at the DFE yields ![Formula][31] has all its eigenvalues negative whether *E* is null or not. Thereby, the DFE is locally asymptotically stable. ## 2. Results and discussion This section is devoted to numerical simulations and discussion of the results. One hindrance here is the lack of data related to the Covid-19 in DRC. I fitted data to a curve estimation using Python. The Fig 2 depicts how my model fits the data collected for the 100 first days of Covid-19 crisis in DRC 25. ![Fig 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/16/2020.12.14.20248182/F2.medium.gif) [Fig 2.](http://medrxiv.org/content/early/2020/12/16/2020.12.14.20248182/F2) Fig 2. Fitting data to the model for parameters estimation Due to lack of available data in DRC, some parameters have been reasonably assumed (without either overlooking or exaggerating them). Provided some parameters after fitting the model with data, some other parameters have been calculated and others were borrowed to [2]. To determine the value of the calculated parameters, I proceed in the following way: At the beginning of the epidemic, *E* = *I* and second equation of system 1 becomes ![Formula][32] whose solution is ![Formula][33] where *k* is an arbitrary constant. Using the value got from data fitting, I get ![Formula][34] Since the value of *λ* is known from the data on life expectancy in DRC (![Graphic][35], see [15] and [24] for details) and provided there some others parameters borrowed from [2], I get the values in Table 1. View this table: [Table 1.](http://medrxiv.org/content/early/2020/12/16/2020.12.14.20248182/T1) Table 1. Parameters used in the model Using the values of parameters like provided in Table 1, a simulation of the model has been performed in Python. The model is designed such that Susceptible stratum collapses as early as possible. This makes the Exposed and infectious groups to explode quickly. The impact of *E* can be evaluated on Fig 3 and Fig 4. When the number of individuals in group Q is high, the number of positive cases drops but the number of individuals hospitalized increases. This is a proof that organized quarantine allows a rapid detection of positive cases and enables the timely management of patients who need to be hospitalized. The fact that the infectious curve reaches its peak without ever affecting the maximum of the population under study is proof that many individuals are freely circulating as long as they have not been deemed suspicious for quarantine. This allows the economy to run normally in the midst of a virus crisis. ![Fig 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/16/2020.12.14.20248182/F3.medium.gif) [Fig 3.](http://medrxiv.org/content/early/2020/12/16/2020.12.14.20248182/F3) Fig 3. A simulation result for the outbreak using data from the 100 first days data from DRC (*ϵ* = 1.5). ![Fig 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/12/16/2020.12.14.20248182/F4.medium.gif) [Fig 4.](http://medrxiv.org/content/early/2020/12/16/2020.12.14.20248182/F4) Fig 4. A simulation result for the outbreak using data from the 100 first days data from DRC(*ϵ* = 3.1). ## 3 Conclusion The mathematical model proposed in this paper has been designed taking into account the vulnerability of the economy of the Democratic Republic of the Congo and other low-income countries. Since the country does not have enough means to control the epidemic from Kinshasa as is the case to date, I propose that the management of this crisis be decentralized. Thus, each provincial governor will be able to set up a team responsible for tracking down suspect people and quarantine them while waiting for the results of the NIBR tests. ## Data Availability Public data, access guranteed [http://www.covid19drc.com](http://www.covid19drc.com) * Received December 14, 2020. * Revision received December 14, 2020. * Accepted December 16, 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/) ## References 1. 1.Mohammed A. A. A., Ahmed A.E et al., Optimization method for forecasting confirmed cases of COVID-19 in China. Journal of clinical medicine, J.Clin.Med.2020,9,674; doi:10.33930/jcm9030674, 9 February 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.33930/jcm9030674&link_type=DOI) 2. 2. Chanyu Yang and Jin Wang A Mathematical model for the novel Coronavirus epidemic in Wuhan, China Mathematical Biosciences and Engineering, MBE, DOI:10.3934/mbe.2020148, 17(3), 2708–2724, 11 March 2020. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3934/mbe.2020148&link_type=DOI) 3. 3.Caramelo F., N. Ferreira and Boliveiros Estimation of risk factors for COVID-19 mortality preliminary results medRxiv preprint, doi:[https://doi.org/10.1101/2020.01.24.20027268](https://doi.org/10.1101/2020.01.24.20027268). 4. 4. Natalie M. Linton et al. Epidemiological characteristics of novel Coronavirus infection: A statistical analysis of publicly available case data medRxiv preprint, doi:[https://doi.org/10.1101/2020.01.26.20018754](https://doi.org/10.1101/2020.01.26.20018754). 5. 5. Hui Wan, Jing-an Cui and Guo-Jing Yang, Risk estimation and prediction by modeling the transmission of the novel Coronavirus (COVID-19) in mainland China excluding Hubei province medRxiv preprint, doi:[https://doi.org/10.1101/2020.03.01.20029629](https://doi.org/10.1101/2020.03.01.20029629). 6. 6.Liu Z., P. Magal, O. Seydi and G. Webb Predicting the cumulative number of cases for the COVID-19 epidemic in China from early data arxiv:2002.12298v1[q-bio.PE]27 February 2020. 7. 7. Liangrong Peng et al., Epidemic analysis of COVID-19 in China by dynamical modeling arxiv:2002.06563v1[q-bio.PE], 16 February 2020. 8. 8. Isaac I. Bogoch et al. Potential for global spread of a novel coronavirus from China International Society of Travel Medicine, Journal of Travel Medicine, 2020, 1–3, doi:10.1093/jtm/taaa011. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jtm/taaa011&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31985790&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F16%2F2020.12.14.20248182.atom) 9. 9. Shi Zhao et al. Estimating the Unreported Number of Novel Coronavirus (2019-nCoV) Cases in China in the First Half of January 2020: A Data-Driven Modelling Analysis of the Early Outbreak Journal of Clinical Medicine, J.Clin.Med 2020,9,388, doi:10.3390/jcm9020388. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/jcm9020388&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32024089&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F16%2F2020.12.14.20248182.atom) 10. 10. Qun Li et al. Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus–Infected Pneumonia The new england journal of medicine, January 29, 2020. 11. 11.Tshisekdi T. President’s speech on coronavirus crisis in DRC Presidency of the DRC, Kinshasa, March 18, 2020 12. 12. Samuel O. Sowole et al., Mathematical model for measles disease with control on the susceptible and exposed compartments Open Journal of Mathematical Analysis, Open J. Math. Anal. 2020, 4(1), 60–75; doi:10.30538/psrp-oma2020.0053. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.30538/psrp-oma2020.0053&link_type=DOI) 13. 13. Tian-Mu Chen et al. A mathematical model for simulating the phase-based transmissibility of a novel coronavirus Infectious Disease of poverty, BMC, [https://doi.org/10.1186/s40249-020-00640-3](https://doi.org/10.1186/s40249-020-00640-3). 14. 14.B. bib19 et al. Mathematical modeling of the spread of the coronavirus disease 2019 (COVID-19) taking into account the undetected infections. The case of China Institudo de Mathemática Interdisciplinar, MOMAT, 1 April 2020 [http://www.doi.org/10.13140/RG.2.2.21543.29604](http://www.doi.org/10.13140/RG.2.2.21543.29604). 15. 15. Michael Y. Li, An Introduction to Mathematical Modeling of Infectious Diseases Mathematics of Planet Earth, Springer, 2010. 16. 16. P. van den Driessche & James Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission Mathematical Biosciences 180 (2002) 29–48, Elsevier Science, 2002. 17. 17. Laura Pullum et al., Disease Prediction Models and Operational Readiness PLoS ONE · March 2014, DOI: 10.1371/journal.pone.0091989. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0091989&link_type=DOI) 18. 18.Amit Huppert& Guy Katriel, Mathematical modelling and prediction in infectious disease epidemiology Clinical Microbiology and Infection · November 2013, DOI: 10.1111/1469-0691.12308. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/1469-0691.12308&link_type=DOI) 19. 19.[https://www.who.int/fr/news-room/detail/29-06-2020-covidtimeline](https://www.who.int/fr/news-room/detail/29-06-2020-covidtimeline) 20. 20.[https://pulitzercenetr.org/reporting/controversial-communication-around-covid-19-drc](https://pulitzercenetr.org/reporting/controversial-communication-around-covid-19-drc) 21. 21.[www.washingtonpost.com](http://www.washingtonpost.com), Hungercouldbedeadlyinpoorcountriesthancoronavirus, mMay 15, 2020. 22. 22. William E. Boyce & Richard C. Diprima Elementary differential and boundary value problems John Wiley & Sons, Inc., 2001. 23. 23.Yolaine S. EPOMA The informal sector in Africa DAUDAPA,24 June 2020, [http://guide.dadupa.com/le-secteur-informel-en-afrique/](http://guide.dadupa.com/le-secteur-informel-en-afrique/). 24. 24.[https://www.who.int/countries/cod/fr/](https://www.who.int/countries/cod/fr/), 12 June 2020 25. 25.[https://www.covid19drc.com](https://www.covid19drc.com), 12 June 2020 [1]: /embed/graphic-1.gif [2]: /embed/graphic-3.gif [3]: /embed/graphic-4.gif [4]: /embed/inline-graphic-1.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/inline-graphic-2.gif [11]: /embed/graphic-10.gif [12]: /embed/inline-graphic-3.gif [13]: /embed/inline-graphic-4.gif [14]: /embed/graphic-11.gif [15]: /embed/graphic-12.gif [16]: /embed/graphic-13.gif [17]: /embed/graphic-14.gif [18]: /embed/graphic-15.gif [19]: /embed/inline-graphic-5.gif [20]: /embed/inline-graphic-6.gif [21]: /embed/graphic-16.gif [22]: /embed/graphic-17.gif [23]: /embed/graphic-18.gif [24]: /embed/graphic-19.gif [25]: /embed/graphic-20.gif [26]: /embed/graphic-21.gif [27]: /embed/graphic-22.gif [28]: /embed/graphic-23.gif [29]: /embed/graphic-24.gif [30]: /embed/graphic-25.gif [31]: /embed/graphic-26.gif [32]: /embed/graphic-28.gif [33]: /embed/graphic-29.gif [34]: /embed/graphic-30.gif [35]: /embed/inline-graphic-7.gif