Revision as of 03:18, 18 November 2023

Abstract

The goal of this work is a dynamic and numerical study of a compartmental mathematical model in order to know the evolution of state variables of the pandemic caused by COVID-19 on red spotlights in México. The model include 7 compartments and the results were compared with the classic SIR model, resulting that our predictions adjusted better the official data. A parameter sensitivity analysis is included.

keywords: COVID-19 pandemic; local stability; Equilibrium point; reproductive basic number.

2 Introduction

In December 31, 2019, Wuhan Municipal Health Commission in Hubei Province, China reported a cluster of 27 pneumonia cases with unknown etiology. The outbreak is associated with common exposures in one wholesale seafood and living animals market in Wuhan City including seven serious cases [1]. The starting point for the first symptomatic case was notified on December 2019. By the 7th of January 2020, the Chinese authorities identified as the causative agent of this outbreak a new type of coronavirus from the family coronaviridae that was called later SARS-CoV-2 which genetic sequence was shared by the Chinese authorities on January the 12th [2]. The via of transmission between human beings, it is similar to the one described by other coronaviruses, by secretions of infected people, mainly by having direct contact with breathing droplets up to 5 microns (allowing them to be transmitted over distances closer than 2 meters) and the contaminated hands or fomites with such secretions, followed by the direct contact of mucosa emissions from the mouth, nose or eyes [3]. The virus SARS-CoV-2 has been founded into the nasopharynx excretions too, including the saliva [4]. On March 11 2020 OMS declared COVID-19 as a world pandemic, since then it is considered the main health problem not only by the mortality rate but also by the secondary effects presented on people such as: Psychological illness, economic loss and the negative impact on their daily activities [5]. From March 20th to May 31st The Secretariat of Health Mexico implemented the social distancing, best known as Healthy Distance (quarantine) because of the COVID-19 pandemic, then to avoid a catastrophic impact on the economy's country, on June 1st, a “The New Normality” started and the goverment reported weekly health alert per region symilar to the traffic light sytem, so with colors there was indicated the type of authorized activities in all economic, labor, educational and social spheres. In this context, the mathematical models of the dynamic transmission on infected diseases have taken an unprecedented relevance. They were applied to help the government`s actions against COVID-19. These models play an important role so as to quantify some possible strategies for mitigation and control of infected diseases [6,7]. There are several models for the infected diseases, since the Clasical SIR Model to more complex proposals [8]. The applied mathematical models, such as Gompertz and Logistic [9,10], have been used with some success to predict the number of infected people with COVID-19. Unfortunately these ones predict on right acurately in short time periods and a huge quantity of data is requested for a good sickness estimation, as it is shown on a recent work done by Torrealba Rodriguez et al. [11], for a COVID-19 estimation in México. On his contributions, Denaiuru et al. [12], makes the proposal for a compartmental mathematical model to the spread of COVID-19 disease, by giving special emphasis on the super-spreaders individuals transmissibility. The problem by working with extensive compartmental models, it is the lack of information for most of the parameters. In a recent work, Neves Armando G.M et. al [13] propose a slightly generalized version of the A-SIR model and propose a scheme for fitting the parameter for the model to real data using the time series only of the deceased individuals. The scheme is applied to the concrete cases of Lombardy, Italy and So Paulo state, Brazil, showing different aspects of the epidemic. In both cases they see strong evidence that the adoption of social distancing measures contributed to a lower increase in the number of deceased individuals when compared to the baseline of no reduction in the infection rate.

In this work, we consider a adaptation of the model in [12] by discarding the super-infectious class but considering the asymtomatic class as an infectious vector. Also we have adjusted the parameters of the resulting model to the Mexican country situations. So we present a consistent compartmental model involving to susceptible class ${\textstyle S}$, exposed class ${\textstyle E}$, infectious with symptoms class ${\textstyle I}$, asymptomatic infectious class ${\textstyle A}$, Hospitalized class ${\textstyle H}$, recovery class ${\textstyle R}$, and death class ${\textstyle F}$. Then a qualitative analysis of the model is done. We show the basic reproduction number and the local stability of the free disease equilibrium. The model is solved numerically to predict the dynamic's pandemic in two federal entities which were considered “red spotlights” for the COVID-19 : the CDMX and Tabasco State. As in the whole Country, in both the National Health Distance and the “New Normality” were implemented. Finally we present the estimations with the Classic SIR Model in order to be compared with the proposed model predictions.

3 Model formulation

We propose a compartment model that takes into account the hospitalized individuals and died individuals due to the virus infection. Considering that the size of the population is constant, ${\textstyle N}$, we subdivide this population into seven classes: suceptibles class ${\textstyle S}$, exposed class ${\textstyle E}$, infectious and symptomatic class ${\textstyle I}$, infectious and asymptomatic class ${\textstyle A}$, sick hospitalized class ${\textstyle H}$, recovery class ${\textstyle R}$ and fatality class ${\textstyle F}$, We consider that the rate at which susceptible individuals are infected by having contact with the symptomatic class is ${\textstyle \beta _{s}}$ and the rate at which succeptible individuals are infected by having contact with the asymptomatic class is ${\textstyle \beta _{a}}$. On the other hand, the incubation period of the virus is considered to be ${\textstyle {\frac {1}{k}}}$ and that a fraction ${\textstyle \sigma }$ of exposed individuals do not develop symptoms, while a fraction ${\textstyle \phi _{h}}$ of individuals who developed symptoms, require hospitalization. The average recovery period of symptomatic individuals, individuals asymptomatic and hospitalized individuals is represented by ${\textstyle {\frac {1}{\gamma _{s}}}}$, ${\textstyle {\frac {1}{\mu }}}$ and ${\textstyle {\frac {1}{\gamma _{h}}}}$, respectively. Furthermore, we assume that symptomatic individuals and hospitalized individuals die at a rate ${\textstyle \delta _{s}}$ and ${\textstyle \delta _{h}}$ respectively. These relationships between the different classes are represented in the diagram of the Figure 1.

Figure 1: Transference diagram for the proposed COVID-19 model.

According to the diagram in Figure 1, the dynamics of the epidemic can be represented using the following model of ordinary differential equations

The region of interest for the analysis of the model is

4 Stability analysis of disease-free equilibrium.

The model (1) has the disease-free equilibrium point ${\textstyle E_{0}=(N,0,0,0,0,0,0)}$. The next generation operator method of Van den Driessche is used to calculate ${\textstyle R_{0}}$ [14]. According to the method, we consider the group of disease variables to ${\textstyle x=(E,I,A,H)}$ and the disease-free variable group to ${\textstyle y=(S,R,F)}$. Thus,

where

and

represent new infections and transitions of individuals between groups, respectively. Therefore

and

where ${\textstyle h_{1}=\phi _{h}+\gamma _{s}+\delta _{s}}$ and ${\textstyle h_{2}=\gamma _{h}+\delta _{h}}$.

The basic reproduction number, which measures the secondary cases of COVID-19 caused by an infected individual in a susceptible population during your infectious period, is given by ${\textstyle R_{0}=\rho (FV^{-1})}$, where ${\textstyle \rho }$ denotes the spectral radius. That is

Note that ${\textstyle R_{0}}$ is the sum of three reproductive numbers :

${\textstyle R_{I}}$ is given by the product of: the infection rate of susceptible individuals by interaction with symptomatic individuals ${\textstyle \beta _{s}}$, the proportion of individuals who develop symptoms ${\textstyle (1-\sigma )}$ and the average duration of an individual in the infectious class ${\textstyle {\frac {1}{h_{1}}}}$. ${\textstyle R_{A}}$ is given by the product of: the infection rate of susceptible individuals by interaction with asymptomatic individuals ${\textstyle \beta _{a}}$, the proportion of individuals who do not develop symptoms after the incubation, ${\textstyle \sigma }$, and the average duration of an individual in the asymptomatic class ${\textstyle {\frac {1}{\mu }}.}$ Similarly, ${\textstyle R_{H}}$ is given by the product of: the infection rate of susceptible individuals by interaction with hospitalized individuals ${\textstyle \beta _{h}}$, the proportion of individuals who develop symptoms and require hospitalization, ${\textstyle \phi _{h}(1-\sigma )}$, and the average duration of an individual to travel from the first contact with the virus to leave the class of hospitalized ${\textstyle {\frac {1}{h_{1}h_{2}}}.}$

Theorem 1: If ${\textstyle R_{0}<1}$ the disease-free equilibrium ${\textstyle E_{0}}$ is locally asymptotically stable.

Proof 1: Noting that the two last equations of system (1) are uncoupled to the remaining equations of the system, we can easily obtain, by direct integration, the following analytical results:

In the same way as the total population N is constant, then

${\displaystyle S(t)=N-[E(t)-I(t)-A(t)-H(t)-R(t)-F(t)].}$Thus, the local stability of model (1) can be studied through the remaining coupled system of state variables, namely, the variables ${\textstyle E}$, ${\textstyle I}$, ${\textstyle A}$ and ${\textstyle H}$. The Jacobian matrix of (1) restricted to these variables and evaluated at the equilibrium point ${\textstyle E_{0}=(N,0,0,0,0,0,0)}$ is given by

The characteristic polynomial of ${\textstyle J(E_{0})}$ is

where

Using the Liénard-Chipard test [15], all the roots of ${\textstyle P(\lambda )}$ (the eigenvalues of ${\textstyle J(E_{0})}$) are negative or have negative real part if, and only if, the following conditions are satisfied:

• ${\textstyle a_{i}>0}$, ${\textstyle i=1,2,3,4}$.
• ${\textstyle a_{1}a_{2}>a_{3}}$.

In order to check these conditions of the Liénard-Chipard test, we rewrite the coefficients ${\textstyle a_{1},a_{2},a_{3}}$, and ${\textstyle a_{4}}$ of the characteristic polynomial in terms of the basic reproduction number ${\textstyle R_{0}}$:

In this way, the first Liénard-Chipard's condition is evident. On the other hand,

where

${\displaystyle \Omega _{1}=k\left(h_{1}^{2}+\left(\mu +h_{2}\right)\left(k+\mu +h_{2}\right)+h_{1}\left(k+\mu +2h_{2}\right)\right)}$.

${\textstyle \Omega _{2}=k\phi \beta _{h}+kh_{2}\beta _{s}+{\frac {k^{2}\mu \phi \beta _{h}+k\mu ^{2}\phi \beta _{h}+k\mu \phi h_{2}\beta _{h}}{h_{1}h_{2}}}}$.

${\textstyle {\begin{array}{l}\Omega _{3}=k\phi \left(k+\mu +2h_{2}\right)\beta _{h}\\+h_{1}\left(h_{2}^{3}+(k+\mu )h_{2}\left(\mu +k\sigma \beta _{a}\right)+2h_{2}^{2}\left(\mu +k\sigma \beta _{a}\right)+k\phi \beta _{h}\right).\end{array}}}$

This shows that Liénard-Chipard's second condition is fulfilled. Therefore, the proof is complete.

5 Numerical results (National Healthy Distance Period)

In this section, we show numerical results for the model (1), using Matlab and Mathematica. To solve SIR model we used the algorithms shown in [16]. The simulations were runed with the official data of the cities of México, CDMX and Tabasco, that was reported by the Healthy Secretariat during the National Healthy Distance Period, from March 20th to May 31st, 2020. The Table 1 shows the parameters used in the numerical simulations. Many of them were obtained from official data from the Mexican government. The infection rate of susceptibles by contact with symptomatic individuals, ${\textstyle \beta _{s}}$, was estimated using least squares in the SIR model with reported data.

The estimated basic reproductive number ${\textstyle R_{0}}$ during the National Period of Healthy Distance, for the epidemic in the CDMX using the SIR model is ${\textstyle R_{0}=1.21}$, while for the model (1) with the parameters of the table 1 is ${\textstyle R_{0}=1.79}$. On the other hand, the estimated ${\textstyle R_{0}}$ for the epidemic in Tabasco, that was obtained with the SIR model is ${\textstyle R_{0}=1.51}$ and the estimation using the model (1) is ${\textstyle R_{0}=1.79}$. This means that the model (1) predicts that the epidemic will be stronger with respect to the prediction of the SIR model.

(a) Estimated daily cases for the CDMX, using the SIR model.	(b) Estimated daily cases for CDMX, using the model (1).
Figure 2: Daily infectious cases for CDMX.

Figures 2a and 2b show the estimations for infections per day using the SIR model and the model (1) for the CDMX. It can be seen that both estimations adjust well the known data. However, with the model (1) the maximum daily infection rate is around 1450 infected, which is reached in mid-June, This is consistent with the reported data. This rate is aproximated better using our nodel instead the SIR model.

(a) Estimated daily new cases for Tabasco using the SIR model.	(b) Estimated daily cases for Tabasco, using the model (1).
Figure 3: Daily infectious cases for Tabasco.

Figures 3a and 3b show the estimates of infections per day using the SIR model and model (1) for Tabasco. Like in the CDMX, the model (1) predicts a maximum of daily cases around 380, that is reached by the end of June 2020, when the confinement was implemented.

(a) Estimated cumulative infected cases in CDMX using the SIR model.	(b) Estimated cumulative infected cases in CDMX, using the model (1).
Figure 4: Estimated cumulative infected cases for CDMX.

The figures 4a and 4b show the cumulative case estimates for the CDMX, using the SIR model and the model (1), respectively. We can note that, both models adjust appropriately the official data. However, they differ in the total number of people who were infected during the pandemic. The SIR model predicted around 48,000 total cases and the model (1) predicted 105,000 total cases. On the other hand, the Healthy Secretariat reported a total of 99,564 infections at September 1st, 2020.

(a) Estimation of cumulative infected cases in Tabasco, using the SIR model.	(b) Estimation of cumulative infected cases in Tabasco, using the model (1).
Figure 5: Estimated cumulative infected cases for Tabasco.

Fgures 5a and 5b show the estimated cumulative case for Tabasco using the SIR model and the model (1), respectively. The SIR model estimates around 12,000 total infected cases while the model (1) estimates around 29000 cumulative infected cases. In September 1st, 2020, the Healthy Secretariat reported 28,471 confirmed infected cases.

6 Numerical results (New Normality Data)

On June 1st, the Mexican government activated the so-called new normality strategy in order to open up activities which are considered essential to reactivate the country's economy. This situation caused a minor increase in almost the entire national territory, furthermore, it caused two plateaus to be observed in the epidemiological curve. If we consider the data on the state of Tabasco after confinement until August 21st, 2020, the parameters were adjusted to this normality and we can see that the model predicts the maximum number of daily cases by mid-July as seen in figure 6. The official reports indicated that the maximum number of daily cases in Tabasco ocurrred in July 12th.

Figure 6: Estimation of daily cases using the model (1) for Tabasco considering the new normality.

For the case of Tabasco, in the new normality, the model (1) predicts a maximum of 640 hospitalized individuals at the beginning of August and around 4,100 deaths are estimated at the end of the pandemic. In the same way, the model estimates the end of the pandemic by the end of November or the beginning of December as can be seen in Figures 7 and 8 (and 6).

Furthermore, also, we have obtained satisfactory respective estimations for the CDMX, which are not included here to obtain a document not so extensive.

Figure 7: Estimation of hospitalized cases in Tabasco using the model (1).

Figure 8: Estimated accumulated deaths in Tabasco using the model (1).

7 Sensitivity analysis of the parameters

In this section we describe the sensitivity analysis of the parameters of model (1), this information is useful to determine the robustness of model predictions to parameter values. Furthermore, it is used to know the parameters that have a high impact on the threshold ${\textstyle R_{0}}$ and these must be handled carefully . More specifically, sensitivity indices allows us to measure the relative change in a variable when a parameter changes. In order to do this analysis, we introduce the following definition.

Definition 1: The normalized forward sensitivity index of ${\textstyle R_{0}}$, which is differentiable with respect to a given parameter ${\textstyle \theta }$, is defined by

The values of the sensitivity indices for the parameters values of CDMX and Tabasco, of Table 1, are showed in Table 2. These parameters have been taken from official data. We have seen the mathematical model (1) describes well the real data, giving rise to numerical results showed in the previous section.

Note that the sensitivity index may depend on several parameters of the system, but also can be constant, independent of any parameter. For example, ${\textstyle Y_{\theta }^{R_{0}}=1}$ means that increasing by a given percentage increases always ${\textstyle R_{0}}$ by that same percentage. For the mathematical model (1), the estimation of a sensitive parameter should be carefully done, since a small perturbation in such parameter leads to relevant quantitative changes. On the other hand, the estimation of a parameter with a rather small value for the sensitivity index does not require as much attention to estimate, because a small perturbation in that parameter leads to small changes.

From Table 2, we see that the most sensitive parameters to the basic reproduction number ${\textstyle R_{0}}$ of the COVID-19 model (1) are ${\textstyle \beta _{a}}$ and ${\textstyle \mu }$. However, these parameters were chosen from report of Healthy Secretary. In contrast, the lest sensitive parameter is ${\textstyle \delta _{h}}$.

8 Conclusion

We have introduced a compartmental mathematical model to understand the dynamic of the COVID 19 pandemic. This model considers a constant population splitted into susceptible class ${\textstyle S}$, exposed class ${\textstyle E}$, infectious with symptoms class ${\textstyle I}$, asymptomatic infectious class ${\textstyle I_{a}}$, hospitalized class ${\textstyle H}$, recovery class ${\textstyle R}$, and death class ${\textstyle F}$. The used data correspond to two entities: Mexico city and the State of Tabasco, both considered during the pandemic as “red spotlights”. First of all we showed the existence and stability of the free illnes equilibrium and then we considedered two stages for the numerical study of model: The National Healthy Distance period and the "New normality" period. For comparison reasons, our simulations were done using our model an the classical SIR model. It must be mentioned that for all the numerical simulations, most of the parameters were taken from official government data such as the CONACYT and the updated reference literature linked to the presented models up to 2020.

As a first stage, the study is done during the lockdown (National Healthy Distance) strategy implemented from March 20th to May 31st. For both models, ours and SIR, the obtained numerical estimations fit in an excellent way the reported data, as the accumulated number of infectious disease and the daily infectious cases as well as for the daily deaths (see Figures 4a-5). However, it can be noticed that our model predicts better the infected accumulated cases, the accumulated deaths and the duration of the pandemic (see Figures 5 and 8). With respect to the reproductive number, for the model (1) and using data and parameters for the CDMX we obtain ${\textstyle R_{0}=1.79}$, while for the SIR model the value is ${\textstyle R_{0}=1.21}$. This means that our model is more realistic with respect to pandemic transmision speed, which can be corroborated with the monthly data until May 31, 2020.

In a second stage, simulations were done using the updated data until August 22th 2020, corresponding to the “New Normality” period when the Mexican government open the activities in sectors which they considered important for the economy, but when the peak's pandemic was not reached yet. The obtained numerical results with the model (1) fit in a satisfactory way and the predictions improve with respect to the obtained with the classic SIR Model (see Figure 6-8).

Finally it is important to mention that the inclusion of classes in the population, such as hospitaized and deaths cases, the models provide valuable information to the goverments during the management of pandemics.

BIBLIOGRAPHY

[1] B. Tang and X Wang and Q. Li and NL. Bragazzi and S. Tang and Y. Xiao and J. Wu. (2020) "Estimation of the transmission risk of the 2019-ncov and its implication for public health interventions.", Volume 9(2). Journal of Clinical Medicine 462–532.

[2]. (2020) "Wuhan seafood market pneumonia virus isolate Wuhan-Hu-1". Complete genome 462–532.

[3] L.S. Hung. (2003) "The SARS epidemic in Hong Kong: what lessons have we learned?", Volume 96(8). J R Soc Med 374-378.

[4] K. K. To and O. T. Tsang and K. H Chan and T. C. Wu and J. M. Chan and W. S. Leung and T. S. Chik and C. Y. Choi and D. H. Kandamby. (2020) "Consistent detection of 2019 novel coronavirus in saliva", Volume 71(15). Clin Infect Dis Off Publ Infect Dis, Soc Am 841–843.

[5] E. Zowalaty and E. Mohamed and J. D. Jarhult. (2020) "From SARS to COVID-19: A previously unknown SARS- related coronavirus (SARS-CoV-2) of pandemic potential infecting humans – Call for a One Health approach", Volume 9. One Health.

[6] J. Djordjevic and C. J. Silva and D. F. R. Torres. (2018) "A stochastic SICA epidemic model for HIV transmission", Volume 84. Applied Mathematics Letters 168-175.

[7] A. Rachah and D. F. M Torres. (2016) "Dynamics and optimal control of Ebola transmission", Volume 10. Math Comput Sci 331-342.

[8] F. Brauer and C. Castillo-Chavez and Z. Feng. (2019) "Mathematical models in epidemiology". New York: Springer-Verlag.

[9] F. Gompertz. (1825) "On the nature of the function expressive of the law of h man mortality, and on a new mode of determining the value of life contingencies. In a letter to Francis Baily, Esq. FRS &c", Volume 115. Philosophical transactions of the Royal Society of London 513-583.

[10] P.F. Verhulst. (1838) "Notice sur la loi que la population poursuit dans son accroisse- ment. Corresp", Volume 10. Math Phys 113–121.

[11] O. Torrealba-Rodríguez and R. A. Conde-Gutiérrez and A. L. Hernández-Javier. (2020) "Modeling and prediction of COVID-19 in Mexico applying mathematical and computational models", Volume 138. Chaos, Solitons and Fractals 109946.

[12] F. Ndaïrou and I. Area and J. J. Nieto and D. F. M. Torres "Mathematical modeling of COVID-19 transmission dynamics with a case study of Wuhan". Chaos, Solitons and Fractals 109846.

[13] A. G. M. Neves and G. Guerrero. (2020) "Predicting the evolution of the COVID-19 epidemic with the A-SIR model: Lombardy, Italy and So Paulo state,Brazi", Volume 413. Nonlinear Phenomena Physica D 132693.

[14] P. Van den Driessche and J. Watmough. (2002) "Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission", Volume 180. Math Biosci 29–48.

[15] A. Liénart and H. Chipart. (1914) "Sur le signe de la partie réelle des racines d'une équation algébrique", Volume 10. J Math Pures Appl (6 éme série) 291-346.

[16] Batista, M. (2020) "Estimation of the final size of the COVID-19 epidemic". medRxiv.

[17]. (2020) "COVID-19, México. https://coronavirus.gob.mx/datos/.". Conahcyt. Mexico.

[18] N. M. Ferguson and D. Laydon and G. Nedjati-Gilani and N. Imai and K. Ainslie and M. Baguelin and S. Bhatia and A. Boonyasiri and Z. Cucunuba and G. Cuomo-Dannenburg. (2020) "Impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand", Volume. Imperial College COVID-19 Response Team.