key: cord-1046908-dd78rhf9
authors: Hussain, Sultan; Zeb, Anwar; Rasheed, Akhter; Saeed, Tareq
title: Stochastic mathematical model for the spread and control of Corona virus
date: 2020-10-14
journal: Adv Differ Equ
DOI: 10.1186/s13662-020-03029-6
sha: a85f339910e078de18e61089a2765a180c53acb9
doc_id: 1046908
cord_uid: dd78rhf9

This work is devoted to a stochastic model on the spread and control of corona virus (COVID-19), in which the total population of a corona infected area is divided into susceptible, infected, and recovered classes. In reality, the number of individuals who get disease, the number of deaths due to corona virus, and the number of recovered are stochastic, because nobody can tell the exact value of these numbers in the future. The models containing these terms must be stochastic. Such numbers are estimated and counted by a random process called a Poisson process (or birth process). We construct an SIR-type model in which the above numbers are stochastic and counted by a Poisson process. To understand the spread and control of corona virus in a better way, we first study the stability of the corresponding deterministic model, investigate the unique nonnegative strong solution and an inequality managing of which leads to control of the virus. After this, we pass to the stochastic model and show the existence of a unique strong solution. Next, we use the supermartingale approach to investigate a bound managing of which also leads to decrease of the number of infected individuals. Finally, we use the data of the COVOD-19 in USA to calculate the intensity of Poisson processes and verify our results.

control techniques of such diseases. It can also be used to predict the expected number of patients in the future under any controlling strategy and to set their goals. Fundamental work is done by researchers to model viral diseases and used by policy makers to control them. For example, in 2009, Pang et al. [1] considered dynamical behavior of a Hepatitis B virus transmission model with vaccination, and Zou et al. [2] discussed modeling the transmission dynamics and control of Hepatitis B virus in China. These models were implemented in China to study the number of HBV patients and to control HBV virus. To understand the spread and control of coronavirus, Chen et al. [3] proposed a mathematical model to understand the transmissibility of coronavirus. Zhou et al. [4] studied pneumonia outbreak associated with coronavirus, whereas Li et al. [5] discussed the early transmission dynamics in Wuhan, China. Huang et al. [6] provided clinical features of the patients infected with coronavirus, whereas Chan et al. [7] discussed familywise transmission of the novel coronavirus. Wu et al. [8] provided newscasting and forecasting of the COVID-19 outbreak in Wuhan, China, through modeling. Zhao et al. [9] have used the model approach to estimate the unreported number of COVID-19 in China in the first half of January 2020. Chen et al. [3] proposed a mathematical model for simulating the transmission of this virus in Wuhan. Ivorra et al. [10] developed and implemented a mathematical model on the spread of COVID-19 using data from China. Sameni [11] studied the epidemic patterns of coronavirus through a mathematical model. Kochańczyk et al. [12] studied constant and time-dependent contact rates of COVID-19 pandemic. For more detail on mathematical modeling of such viral disease, we refer the readers to [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] . Besides the integer-order models, fractional calculus and stochastic differential equations play an important role in the epidemic models; see [23] [24] [25] [26] . A Poisson process is a random process that counts the number of occurrences of certain events that happen at certain rate called the intensity of the Poisson process. For more detail on this process, we refer the readers to [27] .

In this work, we study a stochastic model on the spread and control of coronavirus in which the total population of an infected area is divided into susceptible, infected, and recovered classes. Generally, the exact number of individuals who get the disease, the number of deaths due to coronavirus, and the number of recovered are random and unknown, and thus these numbers are stochastic. Any model containing these terms must be stochastic. These numbers are estimated and counted by a Poisson process. We construct an SIR type model in which the above numbers are stochastic and counted by a Poisson process. To understand the flow and prevent of coronavirus in a better way, we first study the stability of the corresponding deterministic model, calculate the unique nonnegative strong solution, and investigate an inequality managing of which leads to control of the virus. After this, we pass to the stochastic model and show the existence of a unique strong solution. Next, we use the supermartingale approach to investigate a bound managing of which also leads to decrease of the number of infected individuals. Finally, we use COVID-19 data from USA (April 1 to July 19, 2020) to calculate the intensity of Poisson processes and verify our results.

The rest of the paper is organized as follows. In Sect. 2, we give a deterministic model, discus the stability, calculate solution of the model, and investigate an inequality managing of which leads to control the virus. In Sect. 3, we pass to a stochastic model and show the existence and uniqueness of a nonnegative strong solution. Next, we use COVID-19 data from USA, calculate the parameters of the Poison process, verify our results, and 

In this section, we consider an SIR deterministic model on the spread and control of coronavirus, discuss the local and global stability of the model and the reproductive number, calculate the unique strong solution, which can be used to calculate the future number of infected individuals, and investigate an inequality depending on some parameters of the model and the number of susceptible individuals, which would be managed to control the virus and to set the future goals.

Let us divide the total population in three categories, susceptible, infected, and recovered individuals, and denote by S = S(t), t ≥ 0, the number of susceptible individuals, I = I(t) the number of infected individuals, and by R = R(t) the number of recovered individuals. The flow and control are shown in the Fig. 1 . The mathematical model on these classes is expressed by the following system of differential equations:

where a is the rate of new born and migrated individuals, b is the transmission rate from susceptible to infected, c is the contact rate of susceptible with infected, μ is the natural death rate, k is the recovery rate, and λ is the death rate of infected class due to virus.

In the following result, we discuss the local and global stability of the model and calculate the reproductive number R 0 .

Proof The free virus equilibrium point is P 1 ( a μ , 0, 0), whereas the positive virus point is

Putting the point P 1 in the Jacobean matrix J gives the eigenvalues λ 1 = λ 2 = -μ and λ 3 = abc μ -(μ + k + λ). Putting λ 3 < 0, we get the reproductive number R 0 = abc μ(μ+k+λ) . Thus all the eigenvalues are negative if and only if R 0 < 1, and hence the model is locally asymptotically stable.

Similarly, at point P 2 , all the eigenvalues are negative if R 0 > 1. Hence the proposed model is globally stable if R 0 > 1.

Next, we come to the nonnegative unique strong solution of the model, formula for calculating the future number of infected individuals, and an inequality that would be managed to control the virus.

Proof As coefficients in the model are locally Lipschitz, the solution of the model is strong and unique (see Arnold [28] ). Using integrating factors, we can express the solution of the proposed model as

From solution (2) we first observe that I(t) is nonnegative, and using this, we find that S(t) and R(t), are also nonnegative for all t.

The future value of the infected number can be calculated from the expression of I(t).

For the second part, we put the exponent bcSμkλ into I(t) less than zero and get the answer.

Note that the product bc can be easily calculated by using past information in the expression of I(t).

In this section, we pass from deterministic to stochastic model as the number of infected population, deaths due to virus, and recovered individuals are generally not deterministic. These numbers are random and can be counted by a Poisson process. To do this, let us consider a sample space, F a sigma algebra, F t a filtration of F, and P a probability measure. On a filtered probability space ( , F, F t , P) t≥0 , we consider three classes, the susceptible class S = (S(t)) t≥0 , the infected individuals I = (I(t)) t≥0 , the recovered class R = (R(t)) t≥0 , and three independent Poisson processes N t , M t , and K t with intensities λ, γ , and η, respectively. Let M t be the number of individuals who got disease in the time interval (0, t] due to I, N t be the number of deaths, and K t be the number of recovered in this interval. We replace the transmission term The flow in the model is shown in Fig. 2 . 

Lapeyre [27] ), and therefore the increment in a Poisson random process can be given as

0 if there is no jump, and the derivative of the Poisson process coincides with the increment (see Lamberton and Lapeyre [27] ). Thus the Poisson random processes M t , N t , and K t must satisfy

Since μ is the natural death rate, using the later expressions, the differential form of system (3) is (see Shreve [30] in case of assets modeling)

Since most of the data is available daywise, so if t is time in days, then if X i , Y i , and Z i , i = 1, 2, . . . , denote the numbers of infected individuals, deaths, and recovered individuals, respectively, at the ith day since the starting day t 0 = 0, then the Poisson processes M t , N t , and K t satisfy

where the random variables X i , Y i , and Z i are independent and identically distributed with uniform distribution.

If t is time in days, then M t , N t , and K t satisfy

where

Similarly, 

The expected values of the above expressions can be expressed (see Lamberton and Lapeyre [27] and Karatzas and Shreve [31] ) as

with the following quadratic variations (see Karatzas and Shreve [31] , p. 31)

Moreover, the random processes M tγ t, N tλt, and K tηt are martingales (see Shreve [30] ).

Here, we show the stability of the stochastic model and investigate an inequality managing of which leads to decrease of the number of COVID-19 patients. (3), there exists a unique nonnegative solution (S(t), I(t), R(t)), t > 0, with probability one.

Proof As all the coefficients in stochastic system (4) are locally Lipschitz continuous (see Mao [32] ), there exists a unique maximum local solution (S(t),

where T e is the explosion time.

To show that the solution is global, we will show that the explosion time T e = ∞ a.s. Otherwise, we suppose that there is some bounded time such that the solution (S(t), I(t), R(t)) cannot explode to infinity. Let N be a large positive number such that the initial values S(0), I(0), and R(0) belong to the interval [ 1 N , N] . For each integer n ≥ N , let us define the sequence of stopping times

and let inf(∅) = ∞. Note that the above sequence increases and denote T ∞ = lim n→∞ T n ; then T ∞ ≤ T e a.s. If T ∞ = ∞ a.s., then also T e = ∞ a.s. Thus all the values (S(t), I(t), and R(t)) are positive a.s. for all t. Otherwise, suppose there exists a pair (t, ε), t > 0 and 0 < ε < 1, such that the probability P{T ∞ ≤ t} > ε. Thus there exists an integer n 1 ≥ N such that the probability P{T ∞ ≤ t} ≥ ε for all integers n 1 ≥ N .

Next, let us define the smooth function W :

Using Itô's formula (see Karatzas and Shreve [31] ), we calculate

where W · and W ·· are the first-and second-order partial derivatives of W , with respect to the space variable respectively, and S , I , and R are the quadratic variations. Using (3), (4), and (9), we calculate

Using the martingale property of M tγ t, N tλt,and K tηt, this stochastic equation becomes

where

with q = a + 2μ + λ. Further, as the positive space variables S, I, and R represent the numbers of individuals, we have I S ≤ 1, 1 I ≤ 1, and 1 R ≤ 1. Using this, we write LW (S, I, R) ≤ a + 2μ + 3γ + 2λ + 4η.

Thus by Theorem 2.2 in [33] we deduce that T ∞ = ∞, which completes the proof.

Now we come to the inequality among some parameters of the stochastic model, managing of which leads to decrease of the number I(t) of the infected COVID-19 patients.

Applying mathematical expectation to both sides and simplifying, we obtain

and differentiating with respect to t, we find

Using inequality (15) , we get the result.

Note that in this work the parameter γ is calculated from M t by using equations (5) and (8) , which is replaced by transmission and contact rate parameters. These parameters decrease by applying lock down, using a sanitizer, washing hands, and keeping social distance and awareness campaign. Moreover, η is calculated from K t replaced by the recovery term, and we know that recovery terms increase through strong immunity system of the body and good treatment. Decrease of γ and increase of η strick inequality (15) . The strictness of (15) leads to control of disease COVID-19.

A controlling strategy is optimal if it decreases the value of I(t) and strick inequality (15) . Finally, we use relations (5) in I(t) and graph I(t) using the COVID-19 data of USA from first April to 19th July, 2020 taken from [34] . We also calculate the values of γ , λ, and η to verify Theorem 4.

The following table shows the values of γ , λ, and η for April, May, June, and July, 2020, in case of COVID-19 data of USA [34]. In Table 1 the parameter γ is calculated from

Similarly, the parameters λ and η are calculated from 

In this work, we analyze a stochastic model on the spread and control of coronavirus. First,

we considered an SIR model to understand the spread and control of current infectious disease, then investigated the stability of the proposed model. Further, we investigated an inequality managing of which leads to decrease of the number of infected individuals. After this, we presented the nonnegative unique strong solution of the stochastic model and showed an inequality among some parameters of the model through martingale theory, managing of which leads to minimizing the number of COVID-19 patience. We found that decreasing the number of susceptible individuals, transmission, and contact parameters and increasing the recover parameter strick the inequalities. We found that applying lock down, use of sanitizer, washing hands, keeping social distance and awareness campaign decrease the value of susceptible individuals, transmission, and contact parameters, whereas a strong immunity system of the body and good treatment increase the value of recover parameter. We concluded that managing these terms leads to dropping the number of COVID-19 patients. The results are verified through the data of USA and shown in Fig. 3 , and the managing parameters are calculated in the Table 1 . Analysis of the results leads to understanding the flow and control of coronavirus in a clear way.

Dynamical behavior of a Hepatitis B virus transmission model with vaccination

Modeling the transmission dynamics and control of Hepatitis B virus in China

A mathematical model for simulating the phase-based transmissibility of a novel coronavirus

A pneumonia outbreak associated with a new coronavirus of probable bat origin

Early transmission dynamics in Wuhan, China, of novel coronavirus-infected pneumonia

Clinical features of patients infected with 2019 novel coronavirus in Wuhan China

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person to person transmission: a study of a family cluster

Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study

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

Mathematical modeling of the spread of the coronavirus disease 2019 (COVID-19) considering its particular characteristics. The case of China

Mathematical modeling of epidemic diseases; A case study of the COVID-19 coronavirus

Dynamics of COVID-19 pandemic at constant and time-dependent contact rates

COVID-19 control in China during mass population movements at new year

Mathematical formulation and validation of the Be-FAST model for classical swine fever virus spread between and within farms

Early dynamics of transmission and control of COVID-19: a mathematical modelling study

Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study

The reproductive number of COVID-19 is higher compared to SARS coronavirus

The outbreak of Coronavirus Disease 2019 (COVID-19)-an emerging global health threat

A novel covid-19 mathematical model with fractional derivatives: singular and nonsingular kernels

Mathematical model for Coronavirus Disease 2019 (COVID-19) containing isolation class

Statistical analysis of forecasting COVID-19 for upcoming month in Pakistan

Qualitative analysis of a mathematical model in the time of COVID-19

Blind in a commutative world: simple illustrations with functions and chaotic attractors

Fractional discretization: The African's tortoise walk

Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system

Some new edge detecting techniques based on fractional derivatives with non-local and non-singular kernals

Stochastic Calculus Applied to Finance

Stochastic Differential Equations: Theory and Applications

Ruin Probabilities

Stochastic Calculus for Finance

Brownian Motion and Stochastic Calculus

Stochastic Differential Equations and Applications

Environmental Brownian noise suppresses explosions in population dynamics

The authors thank to the editor and anonymous referees for useful comments, which led to a great improvement of the paper and also thank to the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, Saudi Arabia for funded this project, under grant no. (FP-73-42).

The authors confirm that the data supporting the findings of this study are available within the paper cited therein.

The authors declare that there is no conflict of interest regarding the publication of this paper.

Authors have equally contributed in preparing this manuscript. All authors read and approved the final manuscript. 

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.Received: 2 August 2020 Accepted: 4 October 2020