Measles is a highly contagious viral infection that manifests as a rash associated with signs of respiratory infections. It is caused by a virus of the paramyxovirus family whose reservoir is exclusively human  and is transmitted by direct contact with secretions from the nose, throat and through the air  . The virus primarily infects the respiratory tract. Upon infection, the patient passes through a latent period of 6 to 9 days, followed by 6 to 7 day infective period . The infection results in either death or full recovery of the host. In the last case, the host develops lifelong immunity. However, immunity can also be acquired by vaccination before infection, hence its essential role in any measles control initiative.
The children under 5 years remain the most affected. 90% who die have less than 5 years. In developing countries, like Niger where children under one year old represent 4.32%, those under 5 years old 19.73% and those under 15 years 51.18%. The measles remains one of the main causes of infant mortality    .
Our stochastic model is a stochastic SVIR (Susceptible, Vaccinated, Infected, Recovered) model for the measles , where the process is a continuous-time Markov chain resulting from a set of transient states which evolves until it escapes to a set of absorbing states corresponding to disease-free equilibrium. denote respectively the number of susceptible and infected. When the process reaches the set of absorbing states, it remains there permanently. However, before the instant of absorption (which is relatively long), the process passes through a quasi-stationary state.
To understand this phenomenon, we study the long time behavior of the process conditioned on non extinction, which leads us to consider the quasi-stationary distribution introduced by Danoch and Seneta in biology. It allows to describe the steady state behaviour of population models which exhibit discernible stationarity before to become extinct . The term quasi-stationarity refers to the distribution of the Markov chain by conditioning on the event that absorption has not occurred yet . It gives a good measure of the behavior before absorption when the absorption time is very long. But this measure has a number of flaws. Indeed if the set of transient states is finite and irreducible, it is well known that the quasi-stationary distribution exists . But if this set is infinite, the existence of a quasi-stationary distribution is not guaranteed, and even if it does exist, it is practically impossible to determine it explicitly.
For the continuous time SVIR model under some conditions on the effective reproduction number , the quasi-stationary distribution of the number of infected exists and can be closely approximated by geometric distribution.
The main results are stated in theorem 4 and theorem 5. Precisely, let be the effective reproduction number. In theorem 4, we prove that, if , the quasi-stationary distribution of the number of infected can be closely approximated by geometric distribution with parameter . However in theorem 5, we note that if this latter distribution is approximatively geometric with parameter , where , and stands respectively, for the disease transmission coefficient, the natural death rate and the endemic equilibrium point, for the number of infected.
The rest of the paper is organized as follows: The Section 2 describes the deterministic SVIR model by a system of differential equations. The equilibrium points of the system of differential equations are also given. In Section 3, we use the continuous time Markov chains model to form our stochastic SVIR model . Section 4 is devoted to the study of the quasi-stationary distribution of the stochastic SVIR model, followed by numerical simulations in the fifth section. Finally, in the last section, we discuss our stochastic approach and scientific conclusions.
2. The Deterministic SVIR Model
In what follows, denote respectively the number of susceptible, infected and immunized (susceptible vaccinated and recovered patients) at time t.
In this model, the new susceptible (newborns) are introduced at a constant rate n. A fraction, pn, of newborns has acquired immunity by vaccination. The other fraction remains susceptible. p is the probability that a newborn will acquire immunity after being vaccinated. In addition, we assume that:
● the natural death rate is for each compartment.
● infectious patients recover at the rate of .
● infectious patients have an additional death rate from measles.
● we consider the standard incidence , is the disease transmission coefficient. is the average probability of an adequate contact (contact sufficient for transmission) between an infected and a susceptible per unit of time.
The compartment diagram of the transitions in the SVIR model is in Figure 1.
The dynamics of a well-mixed population can be described by the following system of differential equations:
Figure 1. Compartment diagram of model SVIR.
1) In the case of equilibrium without disease, the system (1) admits an equilibrium point ( ) with
is the basic reproduction number  and the effective
reproduction number is defined by . Recall that is defined as average number of secondary infections produced when one infected individual is introduced into a host population where everyone is susceptible  . If this equilibrium point is asymptotically stable . In addition, we have and if and only if . We say that is the critical vaccination coverage of newborns.
2) If , an endemic equilibrium point appears ( ) asymptotically stable , where
3. The Continuous Stochastic SVIR Model
Let be a continuous-time homogeneous Markov chain on the denumerable state space . First, assume that can be chosen sufficiently small such that at most one change in state occurs during the time interval . In particular, there can be either a new infection, a birth, a death, or a recovery. From of state , only the following states are accessible:
corresponding to the possible transitions starting from the state (see Figure 2). has an absorbing set corresponding to disease-free equilibrium states .
Let be the set of neighbors of state :
Setting , the transition rates are defined by:
The transition probabilities of are defined by
Figure 2. States transition.
We have ,
The distribution of is if or and if . Therefore, the marginal distributions are given by:
From the Equations (5), we obtain the Kolmogorov Forward equations, for all and
Hence the system of differential equations verified by the mathematical expectations:
In the previous article, see , we established the following result:
Theorem 1. Let with . Then, for all , and .
Theorem 2. Let with and ( , , ).
If , then (1) and (2) .
Theorem 3. Let , and ( , , )
If , then (1) and (2)
Unlike the deterministic approach, we note that the epidemic is extinguished independently of the threshold with a probability equal to 1. More precisely, if extinction occurs in a time of finite mean, and if the disease eventually disappears in a time of infinite mean. However, before the instant of absorption (which is relatively long) the process passes through a quasi-stationary state. To understand this phenomenon, we study the long time behavior of the process conditioned on non extinction, which leads us to consider the quasi-stationary distribution introduced by Danoch and Seneta in biology.
4. Quasi-Stationary Distribution
The term quasi-stationarity refers to the distribution of the Markov chain by conditioning on the event that absorption has not occurred yet . It gives a good measure of the behavior before absorption when the absorption time is very long. In all the following refer to the probability measure conditional on and for any probability measure . and are the corresponding expectations
Definition 4.1. A probability distribution on the set of transient states is called a quasi-stationary distribution for the process if for all and any measurable set we have .
1) Equivalently, is the unique limiting conditional probability distribution such that , independently of initial distribution .
2) If the set of transient states is finite and irreducible, it is well known that the quasi-stationary distribution exist. But if this set is infinite the existence of quasi-stationary distribution is not guaranteed, furthermore even if it exist, it is typically impossible to evaluate it explicitly. One is therefore lead to consider iterative methods  or asymptotic solutions by diffusion processes   for the quasi-stationary distribution.
Theorem 4. Let be the quasi-stationary distribution of the process and the marginal distribution of the number of infected in a quasi-stationary regime. If , for all ,
Proof. For all , setting
and according to the process definition ,
we deduce that
As in the case of disease-free equilibrium, , thanks to the Equation (11), we have
thus asymptotically the process is a linear birth-death process with infinitesimal generator:
In this case, under the condition , it is well known  there is a unique quasi-stationary distribution for the process which follows the geometric law with parameter . Hence if for all , we obtain:
The proof is completed for the theorem 4. □
Remark. Under the condition , the irreducible Markov chain is positive recurrent. Then an unique invariant probability measure exists and
Thus, the theorem 4 simply states that for all , if .
Theorem 5. Let be the quasi-stationary distribution of the process and the marginal distribution of the number of infected in a quasi-stationary regime. If , for all ,
Proof. If , we have .
The approximation of the process by a birth-death process does not lead to a satisfactory result. In fact we will use the recursive method of Nåssell  and estimate the characteristics of quasi-stationary distribution of the process by those of a diffusion in stationary regime.
Let the conditional distribution of given that the epidemic has not extinguished:
denotes the marginal distribution of , conditional on non-extinction. From Kolmogorov forward equations (Equation (6)) we obtain the following system:
The quasi-stationary distribution is the stationary probability distribution satisfying:
the recursive method of Nåsell  gives the following recursion relationship for :
where is the conditional expectation of given that . Let denote the expectation of , by summing the Equation (18) over we obtain:
It follows that the marginal distribution satisfies the recursion relationship
Let denote the endemic equilibrium point in Equation (3) and . We approximate the process by a diffusion process with drift matrix and covariance matrix
its stationary distribution is approximately bivariate normal with mean 0 and covariance matrix :
from Equation (18) we approximate by and by , thus we obtain .
In other hand since the stochastic means are close to equilibrium points in stationary regime we have set where c is a constant, and from the Equation (18) necessarily . It follows that the marginal distribution satisfies the recursion relationship
The only solution of which is , pour tout . The condition
of the theorem ensures that so the probability distribution is well defined. The proof of theorem 4 is complete. □
Two sample paths of in Figure 3 are simulated over the interval ,
Figure 3. Two sample paths of in quasi-stationary regime for the parameter values ; ; ; ; ; ; ; . . All simulations started from the quasi-stationary distribution (Equation (23)).
Figure 4. approximations and (dashed curve), for the parameter values ; ; ; ; ; . .
with initial distributions the quasi-stationary distribution estimated at Equation (23). In the case , the approximation of the quasi-stationary distribution of I by a diffusion process gives a normal distribution whit mean and variance  .
since , the approximation of this distribution is
where and denote the normal density function and the normal cumulative distribution function respectively. Figure 4 gives and , for the values of .
6. Conclusion and Discussion
The asymptotic analysis of our model identifies three regions in parameter space with qualitatively different behaviors of the quasi-stationary distribution. is significantly greater or less than the deterministic threshold value 1 in the first region or the third region respectively and that is in a second region can be the transition region close to the deterministic threshold value 1 namely
or . This result is analogous to that of Nåsell  for models SIS and SIR with finite population size N where the transition region is determined by a parameter . The fact that a given population belongs to one of these three regions is governed by the size N of the population and the number . In our study, this is a function of . For example, for parameters and the transition region corresponds to , and if in addition we have . As Nåsell  noted, the approximation of the quasi-stationary distribution in the transition region remains a relatively complicated problem.
Note that is average life length in the target population. For measles the target group remains the children of 0 to 5 years. The mortality rate of children before five years of age in Niger is 280% . This means that one in three children does not reach the age of five years and therefore we have years. The parameter in our model hardly exceeds 1. We see by this estimate that the transition region for is in the interval regardless of the size of the population.
Another approach for understanding the dynamics of the system before absorption is the ratio of means approach. Given the initial state, the ratio of expectations distribution (RE) is defined as a ration between the time that the process spends at each transient state and the expected time to absorption (provided that the expected time to absorption is finite). Precisely let be the time that the process , starting from state j, spending in state i before absorption, the
RE is defined by : , where T denotes the absorption time. If is a birth and death process with state space and infinitesimal generator
we have :
as and is a linear birth-death process with infinitesimal generator given in Equation (13), It is easy to see that
Indeed, the RE is a another natural measure of the behavior of absorbing Markov chains before absorption, but the approximation is good only if the convergence to quasi-stationarity is relatively fast . This measure is useful in studying the steady state of outbreaks measles epidemic.
 Edward, S., Kitengeso, R.E., Kiria, G.T., Felician, N., Mwema, G.G., et al. (2015) A Mathematical Model for Control and Elimination of the Transmission Dynamics of Measles. Applied and Computational Mathematics, 4, 396-408.
 Alkassoum, S., et al. (2016) Surveillance épidémiologique de la rougeole au Niger: Analyse de la base de données des maladies à déclaration obligatoire(mdo) de 2003 à 2015. International Journal of Innovation and Scientific Research, 17, 264-274.
 Mitku, S.N. and Koya, P.R. (2017) Mathematical Modeling and Simulation Study for the Control and Transmission Dynamics of Measles. American Journal of Applied Mathematics, 5, 99-107.
 Darroch, J.N. and Seneta, E. (1967) On Quasi-Stationary Distributions in Absorbing Continuous-Time Finite Markov Chains. Journal of Applied Probability, 4, 192-196.
 Artalejo, J.R. and Lopez-Henero, M.J. (2010) Quasi-Stationary and Ratio of Expectations Distributions: A Comparative Study. Journal of Theoretical Biology, 266, 264-274.
 Jacquez, J.A. and O’Neill, P. (1991) Reproduction Numbers and Thresholds in Stochastic Epidemic Models I. Homogeneous Populations. Mathematical Biosciences, 107, 161-186.
 Lahrouz, A., Omari, L. and Kiouach, D. (2011) Global Analysis of a Deterministic and Stochastic Nonlinear SIRS Epidemic Model. Nonlinear Analysis: Modelling and Control, 16, 59-76.
 van der Driessche, P. and Watmough, J. (2002) Reproduction Numbers and Sub-Threshold Endemic Equilibria for Compartmental Models of Disease Transmission. Mathematical Biosciences, 180, 29-48.
 Nasell, I. (1999) On the Time to Extinction in Recurrent Epidemics. Journal of the Royal Statistical Society: Series B (Statistical Methodology), 61, 309-330.