International Journal of Stochastic Analysis

Volume 2017 (2017), Article ID 6313620, 7 pages

https://doi.org/10.1155/2017/6313620

## Global Stability of Nonlinear Stochastic SEI Epidemic Model with Fluctuations in Transmission Rate of Disease

Department of Mathematics, Marshall University, One John Marshall Drive, Huntington, WV, USA

Correspondence should be addressed to Olusegun Michael Otunuga

Received 23 October 2016; Accepted 4 January 2017; Published 23 January 2017

Academic Editor: Manuel Osvaldo Cáceres

Copyright © 2017 Olusegun Michael Otunuga. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We derive and analyze the dynamic of a stochastic SEI epidemic model for disease spread. Fluctuations in the transmission rate of the disease bring about stochasticity in model. We discuss the asymptotic stability of the infection-free equilibrium by first deriving the closed form deterministic () and stochastic () basic reproductive number. Contrary to some author’s remark that different diffusion rates have no effect on the stability of the disease-free equilibrium, we showed that even if no epidemic invasion occurs with respect to the deterministic version of the SEI model (i.e., ), epidemic can still grow initially (if ) because of the presence of noise in the stochastic version of the model. That is, diffusion rates can have effect on the stability by causing a transient epidemic advance. A threshold criterion for epidemic invasion was derived in the presence of external noise.

#### 1. Introduction

Many mathematical models have been developed in order to understand disease transmissions and behavior of epidemics. Among these model is the SEI susceptible-exposed-infectious model. This model is used by some author in studying disease transmission of the Severe Acute Respiratory Syndrome (SARS) disease. Several authors [1–3] have studied other models like SEIR and SEIRS to analyze the spread of the disease. Guihua and Zhen [3] studied the deterministic SEI model by providing conditions for the global asymptotic stability of the infection-free and epidemic equilibrium using the higher dimensional Poincare-Bendixson property.

In this paper, we are interested in studying the effect of stochastic fluctuations in the disease transmission rates in the susceptible-exposed-infected epidemic model. We assume that a susceptible individual first goes through latent period after infection before becoming infectious. We consider a case where the disease is infectious in the latent period and the infected period. We study the role of external noise in the transmission rate. We assume the external noise is a Gaussian white noise. According to Méndez et al. [4], Langevin equations that describe system with real noise should be interpreted as a Stratonovich equation, instead of Ito equation. Due to this reason, we develop a Stratonovich stochastic dynamic SEI model by introducing noise in the transmission rates.

The paper is organized as follows.

In Section 2, we present a Stratonovich stochastic SEI model by allowing the transmission rates to fluctuate around a mean value. The Stratonovich model is now converted into its Ito version. In Section 3, we show that the solution of the stochastic SEI model discussed in Section 2 exists and it is positive. By linearizing the Ito version of the stochastic SEI model around the infection-free equilibrium, we give a closed form expectation of the susceptible, exposed and infected. In Section 4, the closed form value for the stochastic reproductive number is given. This is used to discuss and analyze the stability of the infection-free equilibrium. In Section 5, a numerical simulation is presented to verify our claim. The conclusion of the work is given in Section 6.

#### 2. Stochastic SEI Model

We consider the SEI model for description of the population dynamics for SARS and other similar diseases. The host population is partitioned into three compartments: the susceptible, exposed (latent), and infectious, with sizes denoted by , , and , respectively. The total population . The SEI model is described by the following system of differential equation:where is the recruitment constant, and are rates of efficient contact in the latent period and infected period, respectively, is the natural death rate, is the transfer rate from the exposed to the infectious compartment, , are rates of disease-caused death, and is the rate coefficient of segregation after disease. From (1) and the fact that the total population , we have satisfying the equationIt follows from (2) that the population size may vary with time and . Hence, we consider model (1) in the feasible region:Here, denotes nonnegative real number. It can be shown that is positively invariant with respect to (1).

The system has two equilibriums: the infection-free equilibrium and the endemic equilibrium . The infection-free equilibrium exists on the boundary, , of while the endemic equilibrium exists in the interior of with

By setting , we make the sizes , , and into percentages. This reduces the feasible region to

Let If , is the only equilibrium in . If , the unique endemic equilibrium exists in . Note that model (1) is similar to the model considered by Guihua and Zhen in [3]. They showed using LaSalle’s invariance principle that the disease-free equilibrium, , is globally asymptotically stable in if and unstable if . We define in (6) as the deterministic basic reproductive number.

By allowing the transmission rates and to fluctuate around a mean value, we introduce external fluctuations in the model as follows: where , , are independent noise term with zero mean, and , are noise intensity, a measure of the amplitude of fluctuation with respect to , , respectively. Substituting (7) into (1), the dynamic model (1) reduces to the Stratonovich stochastic model:where , , are independent standard Wiener processes defined on stochastic basis ; denotes the Stratonovich integral [4, 5].

We convert (8) to its Ito’s equivalent using the Stratonovich-Ito conversion (Bernardi et al. [6]) given below.

Theorem 1. *The Ito Stochastic differential equationhaving the same solution as the -dimensional Stratonovich SDE with an -dimensional Wiener process has drift coefficient that is defined in terms of , componentwise, by *

Using Theorem 1, the Ito equivalent of (8) is given by

#### 3. Existence and Uniqueness of Positive Solution of (12)

In this section, we show that not only does the stochastic model (12) have a unique global solution but also the solution will remain within whenever it starts from there.

Following Theorem of Khasminskii [7], we use Theorem 4 below to show the existence and uniqueness of positive solution of (12).

*Definition 2. *Let be the family of nonnegative functions defined on such that they are continuously differentiable with respect to and twice continuously differentiable with respect to .

*Definition 3. *Define the domain by . One defines the differential operator for a function corresponding to a stochastic differential equation with drift and diffusion coefficients and , respectively, bywhere and . Define to be set of positive real numbers. Using the substitution , one reduces (12) to model governing only and and applies Theorem 4 to show existence and uniqueness of positive solution of the reduced model.

Theorem 4. *Suppose that (12) satisfies the classical existence and uniqueness theorem in every cylinder and, moreover, that there exists a nonnegative function such that for some constant Then, for every random variable independent of the processes , there exists a solution of the reduced stochastic differential equation (12) (with ) which is an almost surely continuous stochastic process and is unique up to equivalence.*

*Proof. *It is easy to show that (12) satisfies the classical existence and uniqueness theorem in every cylinder . Define byNote that, for any positive real number , we have and . It then follows that . Also, using the fact that , , we havewhere .

It can be easily shown that as . The result follows. The existence of follows immediately.

The fact that the solution of (12) remains in follows from Corollary of Khasminskii [7].

##### 3.1. Closed Form Expectation of Susceptible, Exposed, and Infected Population Near

We study the condition under which system (12) evolves into an endemic state by analyzing the endemic behavior of the linearized version of (12) around the infection-free equilibrium .

Using the transformation we rewrite (12) to get the nonlinear version: where

The linearization of (12) around the infection-free equilibrium is equivalent to the linearization of (18) around its trivial solution (), given bywhere

It follows that the expected value, , of the solution of (20) is given bywhere , , , , , , , .

#### 4. Stability Analysis of Infection-Free Equilibrium

Using Theorem of Tornatore et al. [2], we show the global stability of the nonlinear stochastic system (12). Notice that, in order to avoid epidemic invasion, we must have , . This is equivalent to Also (23) implies that .

Define

*Remark 5. *It follows that (23) is equivalent to . We can rewrite in terms of by combining (6) and (24) as follows: We call the constant defined in (25) the stochastic basic reproductive number.

Using the following theorem, we get conditions for stochastic asymptotic stability of the infection-free equilibrium.

*Definition 6. *The solution of system (20) is said to be(1)-stable () for , if(2)asymptotically -stable, if it is -stable and moreover as ;(3)exponentially -stable, if, for some positive constants and , (4)when , one says stability in the mean and for , one says stability in mean square.

Theorem 7. *Ifholds, then the solution of (20) is globally asymptotically stable. Furthermore, if , the solution is unstable.*

*Proof. *If (28) holds, it follows from (25) that and (22) implies . If , then and as .

*Remark 8. *Contrary to Remark of Kim and Lin [1], if then epidemic can grow initially, leading to transient epidemic advance. Note that condition (29) implies and . Hence, , , and . The transient epidemic advance is caused by the noise intensity, , in the rate of efficient contact in the latent period, .

*Remark 9. *Note that the global stability of the trivial solution of (20) is equivalent to the global stability of the infection-free equilibrium .

The following theorem shows that the expected value in (22) does not always converge to the trivial solution of (20) if . In fact, if , then as if .

Theorem 10. *If , then in (22) converge to the trivial solution of (20) if .*

*Proof. *From (23) and (24), is equivalent to . If , then , , , , and as , where , . If , then as .

We show the global stability of the infection-free equilibrium and trivial solution of the nonlinear stochastic systems (12) and (18), respectively.

Theorem 11. *If the trivial solution for linear system (20) with drift and diffusion coefficients and , respectively, is asymptotically stable and the drift and diffusion coefficients and , respectively, of nonlinear system (18) satisfy the inequalityin a sufficiently small neighbourhood of , with a sufficiently small constant , then the trivial solution of system (18) is globally asymptotically stable.*

*Proof. *In a sufficiently small neighbourhood of , choose sufficiently small so that . We have where .

#### 5. Numerical Simulations

Figure 1 verifies the global stability criteria discussed in Theorems 7, 10, and 11 and Remark 8.