Next Article in Journal
A Game—Theoretic Model for a Stochastic Linear Quadratic Tracking Problem
Previous Article in Journal
EEG-Based Person Identification and Authentication Using Deep Convolutional Neural Network
Previous Article in Special Issue
A New Emergency-Risk-Evaluation Approach under Spherical Fuzzy-Information Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The SQEIRP Mathematical Model for the COVID-19 Epidemic in Thailand

by
Sowwanee Jitsinchayakul
,
Usa Wannasingha Humphries
* and
Amir Khan
Department of Mathematics, Faculty of Science, King Mongkut’s University of Technology Thonburi, 126, Pracha-Uthit Road, Bang Mod, Thung Khru, Bangkok 10140, Thailand
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(1), 75; https://doi.org/10.3390/axioms12010075
Submission received: 9 December 2022 / Revised: 6 January 2023 / Accepted: 6 January 2023 / Published: 11 January 2023
(This article belongs to the Special Issue Mathematics of the COVID-19)

Abstract

:
The spread of COVID-19 started in late December 2019 and is still ongoing. Many countries around the world have faced an outbreak of COVID-19, including Thailand, which must keep an eye on the spread and find a way to deal with this extreme outbreak. Of course, we are unable to determine the number of people who will contract this disease in the future. Therefore, if there is a tool that helps to predict the outbreak and the number of people infected, it will be able to find preventive measures in time. This paper aims to develop a mathematical model suitable for the lifestyle of the Thai population facing the COVID-19 situation. It has been established that after close contact with an infected person, a group of individuals will be quarantined and non-quarantined. If they contract COVID-19, they will enter the incubation period of the infection. The incubation period is divided into the quarantine class and the exposed class. Afterwards, both classes will move to the hospitalized infected class and the infected class, wherein the infected class is able to spread the disease to the surrounding environment. This study describes both classes in the SQEIRP model based on the population segmentation that was previously discussed. After that, the positive and bounded solutions of the model are examined, and we consider the equilibrium point, as well as the global stability of the disease-free point according to the Castillo-Chavez method. The SQEIRP model is then numerically analyzed using MATLAB software version R2022a. The cumulative percentage of hospitalized and non-hospitalized infections after 7 days after the commencement of the infection was determined to be 11 and 34 percent of the entire population, respectively. The Next-Generation Matrix approach was used to calculate the Basic Reproduction Numbers ( R 0 ). The SQEIRP model’s R 0 was 3.78, indicating that one infected individual can result in approximately three additional infections. The results of this SQEIRP model provide a preliminary guide to identifying trends in population dynamics in each class.

1. Introduction

Humans are social animals. The existence of human society may lead to the emergence of contagions [1]. Some contagions can cause epidemics that spread rapidly among many people in the community at the same time, such as Spanish Flu, SARS, H7N9, Ebola, Zika, MERS, and others [2].
In December 2019, the coronavirus disease 2019 (COVID-19) was discovered. COVID-19 is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), previously provisionally known as the 2019 novel coronavirus (2019-nCoV) [3,4,5,6]. Subsequently, this virus spread rapidly from human to human and spread to many countries around the world, causing the World Health Organization (WHO) to declare a global health emergency in late January 2020 [7,8,9] due to the rapid increase in the number of infected people. Consequently, the number of medical personnel, as well as the amount of medication and medical equipment, were insufficient to cope with the increasing number of infected people [10,11]. Many countries are interested in the spread of this virus and are urgently trying to find ways to support the growing number of patients to prevent increased mortality. According to the WHO report on 8 December 2022, the total number of cases worldwide has reached 642,924,560 and there have been over 6,625,029 deaths. In Thailand, there were 4,711,528 cases in total, with 33,285 deaths [12]. It is anticipated that the number of people infected and deaths will undoubtedly increase in the future if we are unable to control this virus.
Transmission of SARS-CoV-2 occurs when infected people cough and sneeze. After that, contamination of the virus with these aerosols can cause infection in the air or environment around the infected individual [13,14,15]. For this reason, if an infected person has contact with or lives with others, the infection will easily spread to others. Thus, there have been many SARS-CoV-2 infections. The diagnosis of COVID-19 requires either laboratory testing or the use of specialized detection equipment [6,16,17,18]. Typical symptoms of people who contract COVID-19 include fever, coughing, sneezing, shortness of breath, sore throat, loss of sense of smell or taste, as well as complications such as pneumonia. Moreover, the virus can affect the respiratory system of the infected person as well [19]. The incubation period after this virus enters a patient’s body is approximately 2 to 14 days after infection, but there is an incubation period of 5 to 6 days in most cases [9,20]. As I mentioned earlier, the incubation period of COVID-19 and specialized diagnostic tests is quite long. Therefore, it is difficult to detect the infected person in a timely manner.
The spread of the virus not only affects the medical field. Other areas, such as transportation, industries, public health, consumption, agriculture, pandemic vulnerability index (PVI), and other fields are also affected [21]. Many countries around the world have been affected in this way. The outbreak illustrates the urgent need for public health response and preparedness to help mitigate the impact of COVID-19. In each country, it is imperative to assess the situation and find appropriate ways to control the outbreak in a timely manner. However, it is difficult to predict future situations in order to mitigate their potential impact. As such, if there is a tool that helps to assess the situation in advance, it will be possible to find ways to prevent the outbreak in a timely manner.
For this reason, one of the interesting tools for assessing the epidemic situation is a mathematical model. Mathematical models are applications of mathematics to describe natural problems such as biology, chemistry, environmental engineering, environmental ecology, physics, social sciences, statistics, wildlife management statistics, etc. In this epidemic situation, mathematical models are tools for predicting future epidemics. The result will guide the development of strategies for the prevention and elimination of epidemics [22,23]. For instance, what is the trend of increasing the number of infected people, so that each country finds a suitable way to prevent the number of infected people, leading to a decrease in the number of deaths? Therefore, the model is valuable for assessing the epidemic situation, providing guidelines for finding appropriate measures to control the epidemic. Mathematical models are abstract structures that use mathematical concepts to describe and develop behavioral systems [6,22]. The mathematical tools convert real situations into problems that mimic real situations, which solve problems mathematically and then interpret them to predict events in real situations [24,25,26].
Therefore, this paper aims to develop a mathematical model for the COVID-19 epidemic in Thailand. This model is based on the SEIR epidemic model, together with the quarantined population and the pathogens that are transmitted from the infected population to the environment. The model is illustrated in the Materials and Methods section.

2. Materials and Methods

The aim of this paper is to develop mathematical models for the COVID-19 epidemic in Thailand. The methodology is as follows. First, we study the epidemic of COVID-19 and the epidemic model. Next, we develop a mathematical model for the COVID-19 epidemic. After that, we perform a stability analysis to verify the stability of the model and find the basic reproduction number ( R 0 ) to determine the degree of infection. Finally, we consider finding solutions from the model to assess the COVID-19 outbreak.

2.1. Model Development for COVID-19

Several models have been used to study the COVID-19 epidemic such as SIR, SIRS, SEIR, and others. This research studies the spread of COVID-19 by dividing the human population into six main groups: susceptible class (S), quarantine class (Q), exposed class (E), hospitalized infected class ( I h ), infected class (I), and recovered class (R), in conjunction with pathogens (P) in the environment spread by an infected class as shown in Figure 1. The model name is SQEIRP, and the model was presented in Equation (1) and Equations (3)–(8). The idea for this model comes from Samuel Mwalili and Ebrahem A. Algehyne, who studied the spread of COVID-19 using the SEIRP model in 2020 and SQIR model in 2021, respectively.
In the first section, consider the susceptible class as representing a healthy or disease-free population. The susceptible class is transferred to the incubation period after exposure to an infected population or to pathogens transmitted from an infected population to the environment. Another inevitable factor is that population dynamics in the susceptible class depend on the natural birth and death of the population. The dynamic of the susceptible class is shown in Equation (1).
d S d t = A + m R σ α S P 1 + ω 1 P σ β S ( I + I h ) 1 + ω 2 ( I + I h ) ( 1 σ ) α S P 1 + ω 1 P ( 1 σ ) β S ( I + I h ) 1 + ω 2 ( I + I h ) d S
d S d t = A + m R α S P 1 + ω 1 P β S ( I + I h ) 1 + ω 2 ( I + I h ) d S
Second, consider the incubation period population. The population during the incubation period is divided into two classes: quarantine and exposed classes. Both classes represent the populations in quarantine and not in quarantine during the incubation period, respectively. After that, the quarantine and exposed classes are transformed into the hospitalized infected and infected classes, respectively, but if symptoms are severe, some exposure classes may require hospitalization. In addition, consider that the population of quarantine and exposed classes has decreased due to natural mortality. The dynamic of the quarantine and exposed class are shown in Equations (3) and (4), respectively.
d Q d t = σ α S P 1 + ω 1 P + σ β S ( I + I h ) 1 + ω 2 ( I + I h ) γ 1 Q d Q
d E d t = ( 1 σ ) α S P 1 + ω 1 P + ( 1 σ ) β S ( I + I h ) 1 + ω 2 ( I + I h ) θ γ 2 E ( 1 θ ) γ 2 E d E
Third, consider the dynamics of the infected population, comprising the hospitalized infected and the infected classes. The increase in both classes as a result of the incubation period converts to the infected class, and then they are converted to the recovered class. Population dynamics in the hospitalized infected and infected classes are also associated with natural mortality and mortality caused by COVID-19. In addition, the infected class can also spread the pathogens to the environment. The dynamic of the hospitalized infected and the infected classes are shown in Equations (5) and (6), respectively.
d I h d t = γ 1 Q + ( 1 θ ) γ 2 E η 1 I h ( d + δ ) I h
d I d t = θ γ 2 E η 2 I ( d + δ ) I
Next, the part of the recovered class depends on the recovery of the population in the infected class. In this population, if there is a loss of immunity, it becomes the susceptible class. After recovering from COVID-19, permanent immunity is not possible. The dynamic of the infected class is shown in Equation (7).
d R d t = η 1 I h + η 2 I m R d R
Finally, the part of the pathogens class determines the increase in pathogens in the environment caused by coughing or sneezing of the infected class. The pathogens also have a natural death rate. The dynamic of the pathogens class is shown in Equation (8).
d P d t = μ I d p P
Here, this research considers the total human population as N = S + Q + E + I h + I + R and the number of pathogens present in the environment as P. The initial conditions are as follows: S ( 0 ) > 0 , Q ( 0 ) > 0 , E ( 0 ) > 0 , I h ( 0 ) > 0 , I ( 0 ) > 0 , R ( 0 ) > 0 and P ( 0 ) > 0 . An overview of the population for the spread of COVID-19 is shown in Figure 2.
The parameters and parameter descriptions of the SQEIRP model are shown in Table 1.

2.2. The Preconditions and Boundedness of the SQEIRP Model Solutions

The consideration of positive solutions and the bounded solutions of SQEIRP are considered in Theorems 1 and 2, respectively.
Theorem 1. 
The solutions of the mathematical model SQEIRP are positive for all t > 0 .
Proof of Theorem 1. 
Let S ( 0 ) = S 0 > 0 , Q ( 0 ) = Q 0 > 0 , E ( 0 ) = E 0 > 0 , I h ( 0 ) = I h 0 > 0 , I ( 0 ) = I 0 > 0 , R ( 0 ) = R 0 > 0 and P ( 0 ) = P 0 > 0 . Considering Equation (2), we get
d S d t α S P 1 + ω 1 P β S ( I + I h ) 1 + ω 2 ( I + I h ) d S ,
d S d t α P 1 + ω 1 P β ( I + I h ) 1 + ω 2 ( I + I h ) d S ,
1 S d S α P 1 + ω 1 P β ( I + I h ) 1 + ω 2 ( I + I h ) d d t ,
1 S d S α P 1 + ω 1 P β ( I + I h ) 1 + ω 2 ( I + I h ) d t d d t ,
ln | S | α P 1 + ω 1 P + β ( I + I h ) 1 + ω 2 ( I + I h ) t d t + C ,
S t e α P 1 + ω 1 P + β ( I + I h ) 1 + ω 2 ( I + I h ) t d t + C ,
S ( t ) S 0 e α P 1 + ω 1 P + β ( I + I h ) 1 + ω 2 ( I + I h ) t d t > 0 .
As for Q ( t ) , E ( t ) , I h ( t ) , I ( t ) , R ( t ) and P ( t ) , considering Equations (3)–(8), respectively, results in the following
Q ( t ) Q 0 e ( γ 1 + d ) t > 0 ,
E ( t ) E 0 e ( γ 2 + d ) t > 0 ,
I h ( t ) I h 0 e ( η 1 + d + δ ) t > 0 ,
I ( t ) I 0 e ( η 2 + d + δ ) t > 0 ,
R ( t ) R 0 e ( m + d ) t > 0 ,
and
P ( t ) P 0 e d p t > 0 .
Then, S ( t ) , Q ( t ) , E ( t ) , I h ( t ) , I ( t ) , R ( t ) and P ( t ) of the SQEIRP mathematical model are non-negative for all t > 0 . □
Theorem 2. 
The solution for the SQEIRP model is within the boundary
Ω = S , Q , E , I h , I , R , P R 7 : 0 N ( t ) A d , 0 P ( t ) μ A d d p .
Proof of Theorem 2 
For the human population, let N ( t ) = S ( t ) + Q ( t ) + E ( t ) + I h ( t ) + I ( t ) + R ( t ) and N ( 0 ) = N 0 .>
Consider
d N d t = d S d t + d Q d t + d E d t + d I h d t + d I d t + d R d t ,
d N d t = A d S d Q d E ( d + δ ) I h ( d + δ ) I d R ,
d N d t = A d ( S + Q + E + I h + I + R ) δ ( I h + I ) ,
d N d t = A d N δ ( I h + I ) ; N = S + Q + E + I h + I + R ,
d N d t A d N .
Therefore,
N ( t ) A d ( 1 e d t ) + N 0 e d t ,
and get,
lim sup t N ( t ) lim sup t A d ( 1 e d t ) + N 0 e d t = A d .
For the pathogen population P(t), let P ( 0 ) = P 0 .
Consider
d P d t = μ I d p P ,
d P d t μ A d d p P .
Therefore,
P ( t ) μ A d p d ( 1 e d p t ) + P 0 e d p t ,
and get,
lim sup t P ( t ) lim sup t μ A d p d ( 1 e d p t ) + P 0 e d p t = μ A d p d .

3. Stability Analysis of Development for COVID-19

This research presents the SQEIRP model for the COVID-19 epidemic, as shown in Equation (1) and Equation (3) through (8). This research examines the behavior of this model using stability analysis. First, we find the equilibrium point of the system of equations of the developed COVID-19 model. Then, the stability by the roots of the characteristic equation is determined by considering | J λ I | = 0.

3.1. Equilibrium Point of the SQEIRP Model

The equilibrium point of the SQEIRP model for the COVID-19 epidemic can be found by assigning the derivative of the model equal to 0.

3.1.1. Disease-Free Equilibrium Point of the SQEIRP Model

The disease-free point of the SQEIRP model is
X 0 = S * , Q * , E * , I h * , I * , R * , P * = A d , 0 , 0 , 0 , 0 , 0 , 0 .

3.1.2. Epidemic Equilibrium Point of the SQEIRP Model

The epidemic equilibrium point of the SQEIRP model is
X * = S * , Q * , E * , I h * , I * , R * , P * ,
where
S * = A d + m d R * γ 1 + d d Q * γ 2 + d d E * , Q * = η 1 + d + δ γ 1 I h * 1 θ γ 2 γ 1 E * , E * = η 2 + d + δ θ γ 2 I * ,
I h * = A d N * δ I * , I * = A 2 + A 2 2 4 A 1 A 3 2 A 1 , R * = η 1 I h * + η 2 I * m + d , P * = μ I * d p ,
and
A 1 = 1 σ α μ + y 1 ω 1 μ y 2 + y 3 y 4 ω 1 μ ,
A 2 = 1 σ α μ + y 1 ω 1 μ A + y 2 A d N * δ y 1 d p y 2 + y 3 y 4 d p ,
A 3 = y 1 d p A + y 2 A d N * δ ,
with
y 1 = 1 σ β A d N * δ + ω 2 A d N * , y 2 = m η 1 m + d γ 1 + d η 1 + d + δ γ 1 ,
y 3 = m η 2 m + d + η 2 + d + δ θ γ 2 1 θ γ 2 γ 1 + d γ 1 γ 2 + d , y 4 = d γ 2 + d η 2 + d + δ θ γ 2 .

3.2. Basic Reproduction Numbers ( R 0 ) of the SQEIRP Model

The value of R 0 is calculated from the next-generation matrix (G) method, as follows. First, consider
d X d t = F i X V i X ,
where X = Q , E , I h , I , P t , F i X is the matrix of the rate of appearance of new infections in compartment i and V i X is the matrix of rate of other transitions between compartment i and other infected compartments. Thus, we get
F i X = σ α S P 1 + ω 1 P + σ β S ( I + I h ) 1 + ω 2 ( I + I h ) ( 1 σ ) α S P 1 + ω 1 P + ( 1 σ ) β S ( I + I h ) 1 + ω 2 ( I + I h ) 0 0 0 , V i X = γ 1 Q + d Q θ γ 2 E + ( 1 θ ) γ 2 E + d E η 1 I h + ( d + δ ) I h γ 1 Q ( 1 θ ) γ 2 E η 2 I + ( d + δ ) I θ γ 2 E μ I + d p P .
Next, the Jacobian of F i X at the disease-free equilibrium is
F = 0 0 σ β A d σ β A d σ α A d 0 0 ( 1 σ ) β A d ( 1 σ ) β A d ( 1 σ ) α A d 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
and the Jacobian of V i X at the disease-free equilibrium is
V = γ 1 + d 0 0 0 0 0 γ 2 + d 0 0 0 γ 1 ( 1 θ ) γ 2 η 1 + d + δ 0 0 0 θ γ 2 0 η 2 + d + δ 0 0 0 0 μ d p .
Then, the next-generation matrix is
F V 1 = σ β A γ 1 d k 1 k 3 σ β A 1 θ γ 2 d k 2 k 3 + σ β A θ γ 2 d k 2 k 4 + σ α A θ γ 2 μ d k 2 k 4 d p σ β A d k 3 σ β A d k 4 + σ α A μ d k 4 d p σ α A d d p 1 σ β A γ 1 d k 1 k 3 1 σ β A 1 θ γ 2 d k 2 k 3 + 1 σ β A θ γ 2 d k 2 k 4 + 1 σ α A θ γ 2 μ d k 2 k 4 d p 1 σ β A d k 3 1 σ β A d k 4 + 1 σ α A μ d k 4 d p 1 σ α A d d p 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ,
with k 1 = γ 1 + d , k 2 = γ 2 + d , k 3 = η 1 + d + δ and k 4 = η 2 + d + δ .
Therefore, the basic reproductive number given by the spectral radius of the matrix is
R 0 = ρ F V 1 = 1 σ β A 1 θ γ 2 d k 2 k 3 + 1 σ β A θ γ 2 d k 2 k 4 + 1 σ α A θ γ 2 μ d k 2 k 4 d p + σ β A γ 1 d k 1 k 3 .

Global Stability Analysis of Disease-free Equilibrium Point of the SQEIRP Model When R 0 < 1

Next, we consider the global stability of the disease-free equilibrium point of the SQEIRP model with the method of Castillo-Chavez et al. [29]. First, the model is divided into two parts, as shown in Equations (17) and (18).
d Y 1 d t = G Y 1 , Y 2
d Y 2 d t = H Y 1 , Y 2
where Y 1 represents the uninfected populations; Y 1 = S , R R 2 ,and Y 2 represent the infected populations; Y 2 = Q , E , I h , I , P R 5 .
Next, the global stability of the disease-free equilibrium was determined according to the following conditions, with the disease-free equilibrium X 0 = Y 1 0 , 0 .
(i) For d Y 1 d t = G Y 1 , 0 , Y 1 0 is globally asymptotically stable.
(ii) H Y 1 , Y 2 = B Y 2 H ^ Y 1 , Y 2 where h a t H Y 1 , Y 2 0 for all Y 1 , Y 2 Ω .
Theorem 3. 
The disease-free equilibrium point X 0 of the SQEIRP model is globally asymptotically stable.
Proof of Theorem 3 
Let Y 1 = S , R R 2 , Y 2 = Q , E , I h , I , P R 5 and X 0 = Y 1 0 , 0 ; Y 1 0 = A d , 0 , and get
d Y 1 d t = G Y 1 , Y 2 = A + m R α S P 1 + ω 1 P β S ( I + I h ) 1 + ω 2 ( I + I h ) d S η 1 I h + η 2 I m R d R ,
d Y 2 d t = H Y 1 , Y 2 = σ α S P 1 + ω 1 P + σ β S ( I + I h ) 1 + ω 2 ( I + I h ) γ 1 Q d Q ( 1 σ ) α S P 1 + ω 1 P + ( 1 σ ) β S ( I + I h ) 1 + ω 2 ( I + I h ) θ γ 2 E d E γ 1 Q + ( 1 θ ) γ 2 E η 1 I h ( d + δ ) I h θ γ 2 E η 2 I ( d + δ ) I μ I d p P .
First, consider d Y 1 d t = G Y 1 , 0 by defining S = S 0 , R = R 0 and G Y 1 , 0 = 0 , and get
G Y 1 , 0 = A + m R d S m R d R = 0 .
According to Equation (21), then Y 1 Y 1 0 as t . So Y 1 = Y 1 0 is globally asymptotically stable.
Next, consider H Y 1 , Y 2 = B Y 2 H ^ Y 1 , Y 2 by defining
B = γ 1 d 0 σ β A d σ β A d σ α A d 0 γ 2 d ( 1 σ ) β A d ( 1 σ ) β A d ( 1 σ ) α A d γ 1 1 θ γ 2 η 1 d δ 0 0 0 θ γ 2 0 η 2 d δ 0 0 0 0 μ d p ,
where B is an M-matrix (the off-diagonal elements are non-negative) and
H ^ Y 1 , Y 2 = σ β A I h + I d + σ α A P d σ α S P 1 + ω 1 P σ β S ( I + I h ) 1 + ω 2 ( I + I h ) 1 σ β A I h + I d + 1 σ α A P d ( 1 σ ) α S P 1 + ω 1 P ( 1 σ ) β S ( I + I h ) 1 + ω 2 ( I + I h ) 0 0 0 ,
and get
σ β A I h + I d + σ α A P d > σ α S P 1 + ω 1 P + σ β S ( I + I h ) 1 + ω 2 ( I + I h ) ,
and
1 σ β A I h + I d + 1 σ α A P d > ( 1 σ ) α S P 1 + ω 1 P + ( 1 σ ) β S ( I + I h ) 1 + ω 2 ( I + I h ) .
This implies that H ^ Y 1 , Y 2 0 for all Y 1 , Y 2 Ω .
Therefore, by proving the above two conditions, the disease-free equilibrium point X 0 of the SQEIRP model is globally asymptotically stable.

4. Results

The development of a mathematical model for the outbreak of COVID-19 in Thailand was achieved by dividing the population into two groups, such as the human population (N) and the pathogen population (P) transmitted from infected people to the environment. This article considers the model of the COVID-19 epidemic in Thailand, dividing the incubation period and infected population into two lines. The first line represents the quarantined population. This population has contracted the virus. After division, the population in the incubation period is quarantined in hospitals or government-provided residences and eventually infected. The second line represents the mild or asymptomatic population. This population is incubated but not quarantined. Some infected people are not admitted to hospitals or government-provided residences, so these people can spread the disease into the environment. The non-quarantined population will be hospitalized after they present severe symptoms.
The parameter values and the source parameters of this model are shown in Table 1. The natural birth and death rates of the population in Thailand in 2021 are 10.25 per 1000 people and 7.66 per 1000 people, respectively, with A = 10.25 356 × 1000 days 1 and d = 7.66 356 × 1000 days 1 [27]. In addition, those returning to a susceptible class can be re-infected, considering the loss of immunity after recovery of approximately 3 months or 90 days [30]. Therefore, the rate of the recovered population converting to the susceptible class is equal to 1 / 90 days 1 . The model’s hospital infection estimate used the parameters shown in Table 1 compared with actual data in Thailand in January 2022, where the COVID-19 variant omicron outbreak was detected [31]. Figure 3 shows the model solution based on actual data from 1 January 2022 to 17 January 2022. Based on this time period, the model’s solution compared to the actual data shows that the trend of the graph is in the same direction.

4.1. The Numerical Analysis of the SQEIRP Model

The numerical analysis of the SQEIRP model using MATLAB software version R2022a with the initial values of S ( 0 ) , Q ( 0 ) , E ( 0 ) , I h ( 0 ) , I ( 0 ) , R ( 0 ) , and P ( 0 ) are equal to 66,170,000; 12; 6,331; 0; 0; 0; and 500, respectively [23,27,32]. Therefore, the numerical analysis of the SQEIRP model is shown as follows. Figure 4 shows an overview of the dynamics of the human population and the pathogen population from 0–90 days. The dynamics of the human population infected with the virus but during the incubation period and the dynamics of the infected population are shown in Figure 5 and Figure 6, respectively.
The dynamic of the susceptible population declines rapidly over a period of 0–3 days, with the initial population set at 66,170,000. This decline in population shifts to populations that are in the incubation period after exposure to the virus. The incubation period population dynamics are divided into two types: quarantine class, and exposed class, as shown in Figure 5a,b respectively. The dynamics of the quarantine and exposed class increase rapidly over the period of 0–1 day. The maximum population of the quarantine class is approximately 2,200,000, which is less than the maximum population of the exposed class of approximately 53,000,000. After day 1, the dynamics of both populations decline rapidly.
The dynamics of the hospitalized infected and infected population are related to the dynamics of the quarantine and exposed population, respectively. As a result of the incubation period population changes to the infected population, the dynamics of hospitalized infected population and infected population increase rapidly over 0–7 days, with the peak population estimated at 7,500,000 and 22,700,000, respectively. After that, it declines rapidly. The dynamics of the hospitalized infected and infected population are shown in Figure 6a,b respectively.
The final part of the human population dynamics is expressed in the recovered population dynamics, increasing gradually over 0–30 days, and after that, it begins to stabilize. Part of the pathogens population in the environment is caused by the infectious population. The initial population of pathogens was set at 500. At 0–12 days, the number of pathogens increases sharply and then decreases. The dynamics of the recovered population and the pathogen population are shown in Figure 4.

4.2. The Value of Basic Reproduction Numbers ( R 0 )

The dynamic simulation of each population is shown in Figure 4. This paper considers the basic reproduction numbers ( R 0 ) to determine how many infections will originate from one infected person. Where R 0 is greater than 1, it means that one infected person can increase the number of infected people to more than one, which could cause an outbreak. Similarly, a value of R 0 less than 1 would not result in an outbreak. In this work, the R 0 values of the SQEIRP model were considered, together with the parameters assigned to be as relevant as possible for the outbreak in Thailand. The R 0 value of the SQEIRP model is 3.78, meaning that one patient can lead to approximately three additional infections.

5. Discussions

The concept of pathogens being transmitted to the environment around infected people is supported by Samuel Mwalili’s research. The concept of quarantined and non-quarantined population segments is supported by Ebrahem A. Algehyne’s research, as well as the COVID-19 situation report provided by the Ministry of Health in Thailand. The parameters used in the model to estimate the epidemic situation in Thailand are presented in Table 1, which refers to the research of Jankhonkhan et al., who studied the spread of COVID-19 in Thailand from 1 November 2020 to 19 August 2021. The parameters also related to the pathogen research by Samuel Mwalili, as well as the report on the COVID-19 epidemic in Thailand. However, this article considers adjusting parameters that correspond to an increase in the number of infected people based on the outbreak of the COVID-19 variant omicron that began to spread in Thailand in early January 2022. Here, we consider the tendency to change the solution of patients admitted to the hospitalization compared to the actual data from 1 January 2022, to 17 January 2022.
The solution of this model only applies to the beginning of the outbreak of the COVID-19 variant omicron, during which time people were very careful. Thus, there was strict prevention of the disease. As a result, the number of infected people was low. Afterward, people reduced epidemic prevention, and new COVID-19 variants occurred; as a result, the number of infected people increased. Therefore, the parameter representing the number of infected people had to be adjusted. The result of this model assigns all populations in the country as a susceptible class. Then, the number of infected people is quite high. Therefore, it is possible to analyze the susceptible class more accurately.

6. Conclusions

This article aimed to develop a mathematical model suitable for Thai people’s lifestyle during the COVID-19 epidemic in Thailand. The SEIR model was developed alongside the SQEIRP model by adding parts of the quarantine population and the pathogens that spread in the environment. Then, we found parameters suitable for the epidemic situation. The results of the SQEIRP model showed that approximately 7 days after the onset of the infection, the number of hospitalized infected and infected cases was as high as 7,500,000 and 22,700,000 or 11 and 34 percent of the total population, respectively. Moreover, the maximum number of infected people provides a guideline for the government and people to prevent the spread of COVID-19 such that the number of patients does not increase to exceed the estimated number of the model. Therefore, we suggest adjusting parameters according to infected tendencies, including estimating the number of susceptible classes accurately. Thus, the model results will be able to predict the epidemic situation better.
Finally, the researcher expects the SQEIRP model will adapt to the Thai people’s lifestyle and the transmission characteristics of COVID-19. This will lead to results used to predict the spread of COVID-19 in Thailand in order to provide guidelines for assessing the transmission situation epidemic and finding ways to prevent the epidemic in a timely manner next time.

Author Contributions

Writing—original draft preparation, S.J.; writing—review and editing, U.W.H.; writing—review and editing, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Petchara Pra Jom Klao Doctoral Scholarship, King Mongkut’s University of Technology Thonburi (KMUTT), Thailand grant number 68/2563.

Data Availability Statement

Not applicable.

Acknowledgments

We gratefully acknowledge the support from the Petchra Pra Jom Klao Doctoral Scholarship, King Mongkut’s University of Technology Thonburi (KMUTT), Thailand.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Saladino, V.; Algeri, D.; Auriemma, V. The Psychological and Social Impact of Covid-19: New Perspectives of Well-Being. Front. Psychol. 2020, 11, 1–6. [Google Scholar] [CrossRef] [PubMed]
  2. Reperant, L.A.; Osterhaus, A.D.M.E. AIDS, Avian flu, SARS, MERS, Ebola, Zika what next? Vaccine 2017, 35, 4470–4474. [Google Scholar] [CrossRef] [PubMed]
  3. Li, H.; Liu, S.M.; Yu, X.H.; Tang, S.L.; Tang, C.K. Coronavirus disease 2019 (COVID-19): Current status and future perspectives. Int. J. Antimicrob. Agents 2020, 55, 1–8. [Google Scholar] [CrossRef] [PubMed]
  4. Jha, N.K.; Jeyaraman, M.; Rachamalla, M.; Ojha, S.; Dua, K.; Chellappan, D.K.; Muthu, S.; Sharma, A.; Jha, S.K.; Jain, R.; et al. Current Understanding of Novel Coronavirus: Molecular Pathogenesis, Diagnosis, and Treatment Approaches. Immuno 2021, 1, 30–66. [Google Scholar] [CrossRef]
  5. Thompson, H.A.; Mousa, A.; Dighe, A.; Fu, H.; Arnedo-Pena, A.; Barrett, P.; Bellido-Blasco, J.; Bi, Q.; Caputi, A.; Chaw, L.; et al. Severe Acute Respiratory Syndrome Coronavirus 2 (SARSCoV-2) Setting-specific Transmission Rates: A Systematic Review and Meta-analysis. Clin. Infected Dis. 2021, 73, 754–764. [Google Scholar] [CrossRef]
  6. Arenas, A.; Cota, W.; Gómez-Gardeñes, J.; Gómez, S.; Granell, C; Matamalas, J. T.; Soriano-Paños, D.; Steinegger, B. Modeling the Spatiotemporal Epidemic Spreading of COVID-19 and the Impact of Mobility and Social Distancing Interventions. Phys. Rev. X 2020, 10, 1–21. [Google Scholar]
  7. Algehyne, E.A.; Ud Din, R. On global dynamics of COVID-19 by using SQIR type model under non-linear saturated incidence rate. Alex. Eng. J. 2021, 60, 393–399. [Google Scholar] [CrossRef]
  8. Yuen, K.S.; Ye, Z.W.; Fung, S.Y.; Chan, C.P.; Dong, Y.J. SARS-CoV-2 and COVID-19: The most important research questions. Cell Biosci. 2020, 10, 1–5. [Google Scholar] [CrossRef] [Green Version]
  9. Hafeez, A.; Ahmad, S.; Siddqui, S.A.; Ahmad, M.; Mishra, S. A Review of COVID-19 (Coronavirus Disease-2019) Diagnosis, Treatments and Prevention. Eurasian J. Med. Oncol. 2020, 4, 116–125. [Google Scholar]
  10. Aljadeed, R.; AlRuthia, Y.; Balkhi, B.; Sales, I.; Alwhaibi, M.; Almohammed, O.; Alotaibi, A.J.; Alrumaih, A.M.; Asiri, Y. The Impact of COVID-19 on Essential Medicines and Personal Protective Equipment Availability and Prices in Saudi Arabia. Healthcare 2021, 9, 1–14. [Google Scholar] [CrossRef]
  11. Ayati, N.; Saiyarsarai, P.; Nikfar, S. Short and long term impacts of COVID-19 on the pharmaceutical sector. DARU J. Pharm. Sci. 2020, 28, 799–805. [Google Scholar] [CrossRef]
  12. World Health Organization. WHO Coronavirus (COVID-19) Dashboard. Available online: https://covid19.who.int/ (accessed on 8 December 2022).
  13. World Health Organization. Coronavirus Disease (COVID-19): How Is It Transmitted? Available online: https://www.who.int/news-room/questions-and-answers/item/coronavirus-disease-covid-19-how-is-it-transmitted (accessed on 14 August 2021).
  14. Cheng, V.C.C.; Wong, S.K.; Chen, J.H.K.; Yip, C.C.Y.; Chuang, V.W.M.; Tsang, O.T.Y.; Sridhar, S.; Chan, J.F.W.; Ho, P.L.; Yuen, K.Y. Escalating infection control response to the rapidly evolving epidemiology of the coronavirus disease 2019 (COVID-19) due to SARS-CoV-2 in Hong Kong. Soc. Healthc. Epidemiol. Am. 2020, 41, 493–498. [Google Scholar] [CrossRef]
  15. Ong, S.W.X.; Tan, Y.K.; Chia, P.Y.; Lee, T.H.; Ng, O.T.; Wong, M.S.Y.; Marimuthu, K. Air, surface environmental, and personal protective equipment contamination by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) from a symptomatic patient. JAMA 2020, 323, 1610–1612. [Google Scholar] [CrossRef] [Green Version]
  16. Biancolella, M.; Colona, V.L.; Mehrian-Shai, R.; Watt, J.L.; Luzzatto, L.; Novelli, G.; Reichardt, J.K.V. COVID-19 2022 update: Transition of the pandemic to the endemic phase. Hum. Genom. 2022, 16, 1–12. [Google Scholar] [CrossRef]
  17. World Health Organization. COVID-19—Tests. Available online: https://www.who.int/westernpacific/emergencies/covid-19/information/covid-19-testing (accessed on 14 August 2021).
  18. De Walque, D.; Friedman, J.; Gatti, R.; Aaditya, M. How Two Tests Can Help Contain COVID-19 and Revive the Economy; World Bank: Washington, DC, USA, 2020; pp. 1–5. [Google Scholar]
  19. Elliott, J.; Whitaker, M.; Bodinier, B.; Eales, O.; Riley, S.; Ward, H.; Cooke, G.; Darzi, A.; Chadeau-Hyam, M.; Elliott, P. Predictive symptoms for COVID-19 in the community: REACT-1 study of over 1 million people. PLoS Med. 2021, 18, 1–14. [Google Scholar] [CrossRef]
  20. Özdemir, Ö. Coronavirus Disease 2019 (COVID-19): Diagnosis and Management. Erciyes Med. J. 2020, 42, 242–247. [Google Scholar]
  21. Shrestha, N.; Muhammad, Y.S.; Ulvi, O.; Khan, M.H.; Karamehic-Muratovic, A.; Nguyen, U.S.D.T.; Baghbanzadeh, M.; Wardrup, R.; Aghamohammadi, N.; Cervantes, D.; et al. The impact of COVID-19 on globalization. One Health 2020, 11, 100180. [Google Scholar] [CrossRef]
  22. Banerjee, S. Mathematical Modeling; Chapman and Hall Book: London, UK; New York, NY, USA, 2014; p. 1. [Google Scholar]
  23. Mwalili, S.; Kimathi, M.; Ojiambo, V.; Gathungu, D.; Mbogo, R. SEIR model for COVID-19 dynamics incorporating the environment and social distancing. BMC Res. Notes 2020, 13, 1–5. [Google Scholar] [CrossRef]
  24. Kerr, C.C.; Stuart, R.M.; Mistry, D.; Abeysuriya, R.G.; Rosenfeld, K.; Hart, G.R.; Núñez, R.C.; Cohen, J.A.; Selvaraj, P.; Hagedorn, B.; et al. Covasim: An agent-based model of COVID-19 dynamics and interventions. PLoS Comput. Biol. 2021, 17, e1009149. [Google Scholar] [CrossRef]
  25. Lozano, M.A.; Garibo i Orts, O.; Piñol, E.; Rebollo, M.; Polotskaya, K.; Garcia-March, M.A.; Conejero, J.A.; Escolano, F.; Oliver, N. Machine Learning and Knowledge Discovery in Databases, Applied Data Science TrackOpen Data Science to Fight COVID-19: Winning the 500k XPRIZE Pandemic Response Challenge; Springer: Berlin/Heidelberg, Germany, 2021; Volume 12978, pp. 384–399. [Google Scholar]
  26. LANL COVID-19 Cases and Deaths Forecasts, Los Alamos. Available online: https://covid-19.bsvgateway.org/ (accessed on 30 December 2022).
  27. BOI: The Board of Investment of Thailand, Demographic. Available online: https://www.boi.go.th/ (accessed on 14 October 2022).
  28. Jankhonkhan, J.; Sawangtong, S. Model Predictive Control of COVID-19 Pandemic with Social Isolation and Vaccination Policies in Thailand. Axioms 2021, 10, 1–17. [Google Scholar] [CrossRef]
  29. Brauer, F.; Castillo-Chavez, C. A Global Asymptotic Stability Result. In Mathematical Models in Population Biology and Epidemiology; Springer: London, UK; New York, NY, USA, 2012; p. 403. [Google Scholar]
  30. Buskermolen, M.; Te Paske, K.; van Beek, J.; Kortbeek, T.; Götz, H.; Fanoy, E.; Feenstra, S.; Richardus, J.H.; Vollaard, A. Relapse in the first 8 weeks after onset of COVID-19 disease in outpatients: Viral reactivation or inflammatory rebound? J. Infect. 2021, 83, e6–e8. [Google Scholar] [CrossRef] [PubMed]
  31. Department Of Disease Control. COVID-19 Infected Situation Reports. Available online: https://ddc.moph.go.th/covid19-daily-dashboard/ (accessed on 28 October 2022).
  32. Department Of Disease Control. COVID-19 Infected Situation Reports. Available online: https://ddc.moph.go.th/viralpneumonia/eng/index.php (accessed on 28 October 2022).
Figure 1. The population class of SQEIRP model used to determine the COVID-19 epidemic.
Figure 1. The population class of SQEIRP model used to determine the COVID-19 epidemic.
Axioms 12 00075 g001
Figure 2. The SQEIRP mathematical model for the COVID-19 epidemic.
Figure 2. The SQEIRP mathematical model for the COVID-19 epidemic.
Axioms 12 00075 g002
Figure 3. The simulation of the hospitalized infected population compared to the actual data from 1 January 2022 to 17 January 2022.
Figure 3. The simulation of the hospitalized infected population compared to the actual data from 1 January 2022 to 17 January 2022.
Axioms 12 00075 g003
Figure 4. The simulation of human and pathogen population from 0 to 90 days.
Figure 4. The simulation of human and pathogen population from 0 to 90 days.
Axioms 12 00075 g004
Figure 5. The simulated of quarantine (a) and exposed population (b) from 0 to 90 days.
Figure 5. The simulated of quarantine (a) and exposed population (b) from 0 to 90 days.
Axioms 12 00075 g005
Figure 6. The simulated of hospitalized infected (a) and infected population (b) from 0 to 90 days.
Figure 6. The simulated of hospitalized infected (a) and infected population (b) from 0 to 90 days.
Axioms 12 00075 g006
Table 1. The value of the parameter of the SQEIRP model.
Table 1. The value of the parameter of the SQEIRP model.
Parameter DescriptionSymbolValueSource
Birth rate of the human populationA0.000028 days 1 [27]
Natural human death rated0.000021 days 1 [27]
Probability of the population of σ 0.0397[28]
the susceptible class that has been quarantined
Probability of the exposed population θ 0.8[28]
changing to the infected population
Rate of transmission from susceptible class α 0.00414 days 1 [23]
to exposed class due to the pathogens
Rate of transmission from susceptible class β 0.333333 days 1 [Assume]
to exposed class due to the infected
Proportion of interaction with an infected environment ω 1 0.1[23]
Proportion of interaction with an infected individual ω 2 0.1[23]
Rate of transmission from quarantine class γ 1 1/5.2 days 1 [28]
to hospitalized infected class
Rate of transmission from exposed class to infected class γ 2 1/5.2 days 1 [28]
Rate of recovery of the hospitalized infected population η 1 1/10 days 1 [28]
Rate of recovery of the infected population η 2 1/8 days 1 [28]
Death rate due to the coronavirus δ 0.00011 days 1 [28]
Rate of virus spread to environment by infected class μ 0.1 days 1 [23]
Natural death rate of pathogens in the environment d p 0.1724 days 1 [23]
Rate of recovered population transferring into a susceptible classm1/90 days 1 [Assume]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jitsinchayakul, S.; Humphries, U.W.; Khan, A. The SQEIRP Mathematical Model for the COVID-19 Epidemic in Thailand. Axioms 2023, 12, 75. https://doi.org/10.3390/axioms12010075

AMA Style

Jitsinchayakul S, Humphries UW, Khan A. The SQEIRP Mathematical Model for the COVID-19 Epidemic in Thailand. Axioms. 2023; 12(1):75. https://doi.org/10.3390/axioms12010075

Chicago/Turabian Style

Jitsinchayakul, Sowwanee, Usa Wannasingha Humphries, and Amir Khan. 2023. "The SQEIRP Mathematical Model for the COVID-19 Epidemic in Thailand" Axioms 12, no. 1: 75. https://doi.org/10.3390/axioms12010075

APA Style

Jitsinchayakul, S., Humphries, U. W., & Khan, A. (2023). The SQEIRP Mathematical Model for the COVID-19 Epidemic in Thailand. Axioms, 12(1), 75. https://doi.org/10.3390/axioms12010075

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop