Stability Analysis and Stochastic SI Modelling of Endemic Diseases

In this paper, we study a stochastic epidemic model in Meta-population setting. The stochastic model is obtained from the deterministic model by set up random perturbations about the endemic equilibrium state. The outcome of random perturbations on the stability actions of endemic equilibrium is discussed. Stability of the two equilibriums is studied using the Lyapunov function.


Introduction
Modeling of infectious diseases with stochastic differential equation (SDE) has increased foundation lately due to its extensive variety of applications and its aptitude to reflect actuality in epidemiology [1].The diseases outbreaks in a population of susceptibles rationally go behind stochastic processes [2].
Stochastic process occurs naturally in lots of physical applications where randomness is to be incorporated in the mathematical model [3] [4].In recent years, main studies on stochastic model that have been published by researchers have recognized the growing significance of study the stability of stochastic positive equilibrium, as well as the global stability of the endemic equilibrium [5] [6] [7].In this paper we approached by using deterministic and stochastic model.
Briefly deterministic models are model processes which are often described by differential equations, with a unique input leading to a unique output for well-defined linear models and with multiple outputs possible for non-linear models.Throughout this paper, let ( ) , , Ω   be a complete probability space ( ) = for all 0 t ≥ , so (1.0) has the solution ( ) 0 x t = .This solution is called the trivial solution.
Definition 1.1.The trivial solution ( ) 0 x t = of (1.0) is said to be as follows: 1) stable in probability if for all 0 >  , ( ) 2) asymptotically stable if it is stable in probability and, moreover, From the proposed schematics of the compartment model shown in [8], we will extract a metapopulation model for HIV dynamics among the youth coupled with awareness/education i.e., we extended the single patch disease model to include multiple patches (see Figure 1).ϕ , rate at which children grow to become uneducated male youths, 4 ϕ , rate at which children grow to become educated male youth, ρ , awareness/education rate, α , rate at which infected of the two equilibria for system (1.0).We extend our stability analysis to the stochastic system (5.0), which is obtained by random perturbation of the deterministic system (1.0) and find the stability of its positive equilibrium.
Finally, numerical examples which shows the dynamics of systems (1.0) and (5.0) are given, which gives the explicit difference in the dynamics of the models.

Basic Properties of the Model
In this section, the basic properties of model system (1.1) which are useful in the proofs of stability are studied.These are the invariant region and positivity of solutions.The former describes the region in which the solutions of system (1.1) makes biological sense while the latter describes non-negativity of solutions of system (1.1).The model under consideration monitors a human population and as such, we need to have that all the parameters and the variables of the model are positive for all 0 t ≥ .

Positivity of Solutions
The theory of ordinary differential equations requires that, for every set of initial conditions ( )

Invariant Region
Note that  ( ) where A is a constant.Now, applying the initial condition 3), we get Making N the subject in (2.5) we have, As t → ∞ in (2.6) above, the population size N, approaches Therefore, the feasible solutions set of system (2.7) enters the region In this case, whenever

Basic Reproduction Number
The basic reproduction ratio ( 0 R ) is defined as an infections originating from an infected individual that invades a population originally of susceptible individuals.
The above system can be represented in matrix form as .

I fI vI = +
where f is the matrix of the infection rates and v is the matrix of the transition rates.
The spectral radius of the Metzler Matrix, ( ) , is defined as the largest eigenvalue of the Metzler Matrix.Thus: , , Sub-Population j infects on average more than one other person; the infection could therefore establish itself in the population and become endemic.An SIR epidemic model, where the presence or absence of an epidemic wave is characterized by the value of 0 j R .

The Global Stability of the Endemic Equilibrium
In this part, we analyse the global stability of the endemic equilibrium point * E by construction a appropriate Lyapunov function.For simplicity, we consider the reduced model system (6) to prove for global stability.We use the come up to of [8] as it is used for several complicated epidemiological models.We consider the Lyapunov function of the form ( ) ( ) is a properly chosen positive constant in the given region i P . is a population of compartment i and * i P is the equilibrium level.So we define the Lyapunov function as  The time derivative of L is ( ) 1 , F is non-positive by following the approach of [9].Thus, < and is zero when < is the singleton * E which is our endemic equilibrium point.By LaSalles invariant principle [9] we conclude that * E is globally asymptotically stable (g.a.s).Thus, we establish the following theory.

The Stochastic Model
Stochastic perturbations were bring in some of the major parameters involved in the model equations.
Here, we bring in stochastic perturbations in the major parameters of the deterministic model (1.1).Thus we permit stochastic perturbations of the variable

S S S S S I I I I I T
around their values at positive equilibrium * E .Hence, we assume that the white noise of the stochastic perturbations of the variable around values of are proportional to the distances of

Stochastic Stability of the Positive Equilibrium
It can be shown clearly that, the deterministic model (1.1) has one disease-free equilibrium , ,0,0,0,0,0 which is globally asymptotically stable when 0 1 R ≤ .However, when 0 1 R > , the disease-free equilibrium 0 E is unstable.Obviously, there is also a unique positive endemic equilibrium , , This equilibrium is globally asymptotically stable.The stochastic system (5.0) has the similar equilibria as the deterministic system (1.1).Assuming that   ( ) Clearly, the endemic equilibrium * E corresponds to the trivial solution ( ) 0 x t = in (5.2).We denote L to be the differential operator associated with (5.2), defined for the family of nonnegative functions ( ) ( ) such that it is continuously differentiable with respect to t and twice with respect to x.
Theorem 5.1: Suppose a function ( , , and 0 p > .Then the trivial solution of (5.2) is p th moment exponentially stable.Again, given that p = 2 the trivial solution is supposed to be exponentially stable in mean square and the equilibrium x = 0 is globally asymptotically stable.From theorem 5.1, the conditions for stochastic asymptotic stability of trivial solution of (5.0) are given theorem 5.2.Theorem 5.2: Suppose and hold, then the zero solution of (5.0) is asymptotically mean square stable.
Proof: We consider the Lyapunov ( ) ( ) i w i = non-negative constants that will be chosen in the course of the proof.It can be easily ascertained that inequality (5.4) hold true when 2 p = Applying the operator L on ( ) , , , , , According to Theorem 5.1, we conclude that the trivial solution of stochastic system is globally asymptotically stable.
Hence, according to theorem 5.1, the proof is completed.

Conclusion
In this article, the dynamics of deterministic epidemic model and its stochastic variant are presented.The stability analyses of the deterministic model were investigated.Suitable Lyapunov functions were constructed for the global stability of the two equilibria.Mathematical analysis was done and it was established that in the absence of the disease a disease free equilibrium will always exist if 0 1 j R ≤ for j = 1, 2, 3, 4,.We also established that the endemic equilibrium exists in the presence of the disease that is when 0 1 j R > for j = 1, 2, 3, 4, with the infectious population greater than zero.Reducing the infection in the vector population reduces 0 j R for j = 1, 2, 3, 4, greatly.Thus the best methods of controlling HIV transmission is to target the Infected uneducated female youth, Infected educated female youth, Infected uneducated male youth, Infected educated male youth.0 j R is a threshold that completely determines the global dynamics of disease transmission.Our major purpose of the study was to examine the asymptotic stability behavior of the endemic equilibrium of the stochastic version of the deterministic epidemic model in Metapopulation setting.
stable in the large if it is stable in probability and, moreover, In this paper, we consider the epidemic model in Meta-population setting.

Figure 1 .
Figure 1.A schematic of the metapopulation model.A schematic of the metapopulation model for HIV transmission in the youths coupled with awareness/education in each patch i, i = 1, •••, n.
youth take the ART.The solid lines represent movement between classes and the dashed lines represent rate at which a susceptible individual moves in to infected class.The model answers one important underlying research subjects; the determination of the existence of the threshold parameter which hints on the spreading or dying out of an invading epidemic into a population of susceptible.In this research article, we first study the positivity and boundedness of the D. P. Kelkile DOI: 10.4236/apm.2018.85030519 Advances in Pure Mathematics system (1.0).The basic reproduction ratio is determined.Applying the hypothetical theorem of the Lyapunov functional, we find out the global stability

2 , 3 , 4 , 1 jR
then each infectious individual in Sub-Population j infects on average less than one other person and the disease is likely to die out Otherwise, If 0 > for j = 1, 2, 3, 4, then each infectious individual in ,3,4,5,6,7,8,9,10,11 i B = are independent wiener processes.We examine the asymptotic stability behavior of the equilibrium * E of the stochastic equation (5.0) and contrast results with the deterministic model (1.1).

0 1 R
≤ , we examine the stability of the endemic equilibrium * E of (5.0).The stochastic differential equation (5.0) can be centered at its positive equilibrium * E by the change of variables d , death rate of the youth and children, δ , disease induced rate in the youth and children before ART, 1 δ , disease induced deaths in the youth and children after ART, λ , birth rate, θ , rate of vertical transmission, 1 YUF I , infected uneducated female youth, YEF I , infected educated female youth, YUM I , infected uneducated male youth, YEM I , infected educated male youth, and Y T , antiretroviral therapy treatment among the youths.Parameters: µ From the assumptions of the theorem, we deduce that 0, 1,2,3,4,5,6,7,8,9,10,11 D > .Hence, D is a symmetric positive definite matrix.Let m λ denote the minimum of its eleven positive eigenvalues