Stochastic differential equation model of Covid-19: Case study of Pakistan.

In recent years, the world lived a horrible nightmare named the Covid-19 pandemic. It's changing lifestyle caused the closure of critical establishments like schools, sports halls, companies, etc., as well as causing damage and menace for humanity. Mathematical models are helpful tools for modeling and analysing the infectious disease transmission This article presents a stochastic model of the Covid-19 epidemic for a population with five compartments. We give a numerical analysis of the proposed stochastic model. Also, we compare it with results of the corresponding deterministic model.

Introduction

Recently, the coronavirus pandemic has changed the healthy balance of humanity; the world has known a rapid increase in the number of fatalities. Moreover, the novel coronavirus is a respiratory disease that generally transmits through direct contact with an infected person. However, it also can spread by indirect contact with a virus- infected environment. The hardness degree caused by the infection coronavirus differs from one to another according to immune power. Generally, the syndromes of this disease are chills, weakness, stiffness and loss of smell or taste. According to the report done on 12 October 2021 by World Health Organization, the rate of infection by coronavirus is around 237 655 302, and the number of deaths is 4 846 981 worldwide. The first time when announced the outbreak of a new disease named coronavirus was in December 2019. After a new breed of coronavirus SARS-CoV-2 was discovered in Janvier 2020, the town of Wuhan. It’s spread fastly in all countries and provides a great menace to humanity, caused bringing public life to a halt. Because of Covid 19, universities, sports parks, and high-end establishments shuttered, affecting the countries economies and causing economic harm. Furthermore, increasing multimedia awareness programs about Covid-19’s various can be reduced and controlled. Also, the awareness campaigns play a significant role in urging individuals to take the required measures against the disease, as well as help to refute the fake news about Covid-19 and stressing preventative measures like hand washing and mask use.

In literature, many scholars have adopted mathematic approaches to model the transmission of Covid-19 in different countries. For example, in [1], the authors used the classical deterministic SIR epidemic model to formulate a compatible model for the transmission of Covid-19. In addition, they modeled the process of the propagation of the disease by the convex incidence rate. Khan et al. [2] have presented a fractional-order Covid-19 model by the fractional-order differential equation type. Moreover, they gave a detailed mathematical analysis of the dynamics of the Covid-19 with a case study of Saudi Arabia. A new approach based on piecewise differential and integral operators developed by Atangana and Seda in [3] and prolonger in [4] on the epidemiological fields precisely to modeled the third waves of Covid-19 in Turkey, Spain, and Czechia. In [5], the authors study the dynamics of Covid-19 disease with the fuzzy Caputo approach. In the above-cited works, the authors used deterministic ways to formulate mathematical models of the transmission of Covid-19. Whereas these models have some limitations to describe the reality compared with the stochastic models [6] because the parameters of the system are not constant, but it fluctuates around an average value. In addition, the transmission of epidemics is compatible with randomness because in real-life environmental fluctuations influence the propagation of epidemics. Many works on the application of stochastic calculus in epidemiology have been developed. For example, El Koufi et al. in [7] proposed a stochastic epidemic model with vaccination and nonlinear incidence rate and investigated the asymptotic propriety of the system. In [8], the authors have used a system with Markovian switching to study the change of regimes of a human population divided into three classes Susceptible, Infected, and Recovered. In addition, many works integrated the Lévy processes in deterministic systems to give some degree of realism to this model (see Refs. [6], [9], [10], [11]) and to express the effect of certain phenomena such as earthquakes, volcanoes, tsunami, and pandemic on the dynamics of the system. In [12], the authors have proposed an influential work on the Lotka–Volterra population-driven with Lévy noise. Gao and Wang [13] have studied a stochastic mutualism model integrating the Lévy and Markov processes.

The purpose of this paper is to investigate the effect of white noise on the dynamics of the Covid-19 epidemic. Then, we introduce a stochastic model of Covid-19 that includes the randomness effect, and we discuss numerical conditions for the extinction of the Covid-19 epidemic. The rest of this paper is organized as follows. In Section “Deterministic model”, we presented the deterministic system and its analytical results. The limitations of the deterministic model and the presentation of the stochastic differential equation model existing in Section “Stochastic model”. In Section “Numerical investigation”, numerical simulations of our stochastic model with a brief discussion of the obtained results. Finally, we give a conclusion of the paper and futur investigation in Section “Conclusion”.

Deterministic model

In [14], Li et al. investigated an epidemic model of Covid-19 virus transmission and have divided the population into five groups of individuals, namely, Susceptible , Exposed , Infected , Asymptomatic infected , Recovered . They presented their model by the following ordinary differential equations: where is the recruitment rate from the population, represents the disease transmission rate, denotes the transmissibility multiple of asymptomatic infection. The parameter is the natural mortality rate. However, The terms represent the infection rate of the exposed individuals before a successful incubation period. During the test of individuals, the asymptomatic infections observed for the individuals and its recovery is shown by the rate , the recovery rate from the symptomatic infection is . The parameters , and represent the disease death rates in the exposed, infected, and asymptomatic infected groups, respectively. The theoretical results are shown in [14] by the following:

System (1) has a basic reproduction number defined by where , , .

The system (1) has a disease-free equilibrium defined by If , there exists a unique endemic equilibrium expressed as In addition, they have shown the following analytic results

If , then the disease-free equilibrium of Covid-19 model (1) is locally asymptotically stable.

If , the disease-free equilibrium of Covid-19 model (2) is globally asymptotically stable.

Stochastic model

In the system (1), the parameters are all deterministic and do not consider the effects of environmental fluctuation and randomness. Whereas, in a biological environment, it is more suitable and reasonable to consider how environmental noise affects the spread of the Covid-19 epidemic. Numerous investigations have shown that environmental fluctuations affect the propagation of an epidemic on the population (see, [7], [8], [11], [12], [13], [15], [16]). Therefore, there are many approaches to consider the random fluctuations in the deterministic systems one of them suppose that the epidemic is subject to some small and standard random fluctuations that can be expressed by white noise. In this paper, in order to examine the effect of environmental noise on the system (1), we suppose that the Covid-19 model (1) is subject to environmental fluctuations and the stochastic perturbations are of white noise type and proportional to the variable. Then, we introduce a stochastic Covid-19 model defined by the following system of stochastic differential equation where , are a standard Brownian motion (see Definition 1) defined on a complete probability space be a complete probability space with a filtration satisfying the usual conditions (i.e. it is increasing and right continuous while contains all -null sets), , represent the intensity of noise.

[17]

The positive real-time numbers process is called the Wiener process if

, .

, are independent of .

, .

The n-dimensional stochastic differential equation (briefly, SDE) is defined by: where is a n-dimensional standard Wiener process, is the drift coefficient, and is the diffusion coefficient.

For any initial value , there is a unique solution to Eq. (3) on and the solution will remain in with probability one, namely for all almost surely (briefly a.s.).

By using the same steps of proof of Theorem 2.1. in [11] we can prove the above theorem. In the next Section, we will show numerical simulations for the model (3) to show the effect of white noise on the dynamics of the Covid-19 epidemic.

Numerical investigation

To study the dynamics of the model (3), we present numerical simulations by using the Milstein numerical approximation method for the SDE presented in [18]. Then, we can write the system (3) by the following discretized form: where are sequences of random numbers uniformly distributed in and is a time step. The initial conditions are given as follows: , , , . In this paper, we choose the parameters of the third wave in Pakistan between March 06, 2021, till April 30, 2021 [14], [19]. According to [19], the total population in Pakistan at is , and the recruitment number in model (3) is per day. The natural mortality rate is  [19]. In addition, we regrouped the rest of the parameters in Table 1. Then, by simple computation, we have , then there exists an endemic equilibrium to the system (1) which is globally asymptotically stable.

Table 1

The values for parameters given in model (3).

Parameter	Value
Λ	14
β	0.5
η	0.05
μ	0.1
δ	0.02
ϖ	0.1
μ_0	0.8
μ_1	0.01
μ_2	0.02
γ_1	0.15
γ_2	0.15

Choosing the value of transmission rate of Covid-19 equal to 0.792, we obtain that , then the disease-free equilibrium of system (1) is globally asymptotically stable Fig. 3 confirms this remark. In addition, all values less than 0.792 lead to obtaining a value less than one of the basic reproduction numbers of the system (1). Thus, the coefficient can likewise decrease if we reduce the contact between the infected and susceptible groups. Therefore, we keep the same parameters values used with , and to study the effect of the environmental fluctuation on the dynamics of Covid-19, choosing the different white noise values as: , , , . The integration of white noise in the Covid-19 system leads to the extinction of the epidemic. Then, increasing the white noise value leads to the disappearance of disease (see, Fig. 4, Fig. 5). In addition, the , and in system (3) going extinct rapidly comparing with , and in a deterministic system (1). On the other hand, when (see, Fig. 1, Fig. 2). We compare in this situation the dynamics on the stochastic system (3) with the corresponding deterministic system (1) for this taking the noise valued on the system (3) as follows: , , , , . So, the result of the computer running presented in Fig. 7, Fig. 8, Fig. 9 shows that the integration of the random noise in the system leads to the extinction and the reduction of the Infected, Exposed, and Asymptomatic infected individuals. Moreover, by taking the following different noise values in the system (3): , , , (i 1,2,3,4,5) we obtain the graph shown in Fig. 10, Fig. 11, Fig. 12. that show the impact of noise intensity in the spread of the Covid-19 epidemic. Then, a large value of noise led to control the propagation of the epidemic (see, Fig. 10, Fig. 11, Fig. 12).

Fig. 3

Paths of the deterministic systems (1).

Fig. 4

Comparison between the paths of stochastic (3) and deterministic (1) systems with different noise values.

Fig. 5

Comparison between the paths of stochastic (3) and deterministic (1) systems with different noise values.

Fig. 1

Paths of the deterministic systems (1).

Fig. 2

Paths of the deterministic systems (1).

Fig. 7

Comparison between the paths of stochastic (3) and deterministic (1) systems.

Fig. 8

Comparison between the paths of stochastic (3) and deterministic (1) systems.

Fig. 9

Comparison between the paths of stochastic (3) and deterministic (1) systems.

Fig. 10

Comparison between the paths of stochastic (3) and deterministic (1) systems with different noise values.

Fig. 11

Comparison between the paths of stochastic (3) and deterministic (1) systems with different noise values.

Fig. 12

Comparison between the paths of stochastic (3) and deterministic (1) systems with different noise values.

Conclusion

In this investigation, we presented a stochastic differential equation mathematical model to study the dynamic of Covid-19 in Pakistan. The proposed model considers five classes; Susceptible, Exposed, Infected, and Asymptomatic infected, Recovered. We have affected an environmental perturbation of white noises type in the model (1). We employed the Milstein numerical approximation method to give a comparison between the model (1), (3). Then, with the help of MATLAB and the data presented in [14], we presented graphs that confirm the power of randomness to reduce the number of infected people of Covid-19. We know that the parameters of the epidemic model can depend on the nature of each environmental regime, so to make the system (3) more realistic, we will integrate the telegraph noise into the (3) system in our future work.

CRediT authorship contribution statement

All the authors have equal contribution in this work.

Uncited references

Fig. 6

Fig. 6

Comparison between the paths of stochastic (3) and deterministic (1) systems with different noise values.

Declaration of Competing Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.