Stochastic Modelling of the Transmission Dynamics of Measles with Vaccination Control
Kitengeso Raymond E.
Independent Scholar, Dar es Salaam, Tanzania
Email address:
To cite this article:
Kitengeso Raymond E. Stochastic Modelling of the Transmission Dynamics of Measles with Vaccination Control. International Journal of Theoretical and Applied Mathematics. Vol. 2, No. 2, 2016, pp. 6073. doi: 10.11648/j.ijtam.20160202.16
Received: September 6, 2016; Accepted: November 4, 2016; Published:December 5, 2016
Abstract: Measles is still endemic in many parts of the world including developed nations, despite the availability of the infectious disease vaccine since 1963. Elimination of measles requires maintaining the effective reproduction number by achieving and maintaining low levels of susceptibility R_{0} <1. In this project, we concentrate on the stochastic modelling of the transmission dynamics of measles with vaccination control. We have obtained the stochastic differential equations model from the deterministic model. Simulation of the stochastic differential equations model have been performed as well as the deterministic model. The stochastic differential equations model has described the transmission dynamics of measles with more information compared to the deterministic counterpart. Mathematical technique used in the simulation of the stochastic differential equations model is EulerMaruyama numerical scheme and discussions of the model.
Keywords: Measles, Vaccination, Immunity, Stochastic Modeling
1. Introduction
In this chapter we discuss mathematical modelling, types of models and use of models.
Modelling is the process of producing a model; a model is the presentation of the construction and working of some system of interest. It is similar to but simpler than the system it represents. Mathematical modelling is the process of expressing real world phenomena or system using mathematical principles and formula. Mathematical Modelling can classified into deterministic models and stochastic models based on certainty. In this study we will study stochastic model of dynamics of measles with vaccination by simulation of the transformed deterministic model into stochastic differential equations model.
Measles (also called rubeola) is a highly contagious viral infection that can be found around the world through persontoperson transmission mode, with over 90% attack rate among susceptible persons. The measles virus is a paramyxovirus, genus morbillivirus. Even though an effective vaccine is available and widely used, measles continues to occur even in developed countries. Children under five years are most at risk.
The main symptoms of measles are fever, runny nose, cough and a rash all over the body, it also produces characteristicsred rash and can lead to serious and fatal complications including pneumonia, diarrhea and encephalitis. Many infected children subsequently suffer blindness, deafness or impaired vision. Measles confer lifelong immunity from further attacks [15].
Measles is a highly contagious virus that lives in the nose and throat mucus of an infected person. It can spread to others through coughing and sneezing. Also, measles virus can live for up to two hours in an airspace where the infected person coughed or sneezed. If other people breathe the contaminated air or touch the infected surface, then touch their eyes, noses, or mouths, they can become infected. Measles is so contagious that if one person has it, 90% of the people close to that person who are not immune will also become infected. Infected people can spread measles to others from four days before through four days after the rash appears. Measles is a disease of humans; measles virus is not spread by any other animal species.
There is no specific treatment for measles. People with measles need bed rest, fluids, and control of fever. Patients with complications may need treatment specific to their problem.
In the work [17] performed a study on mathematical modeling on the control of measles by vaccination. In their study SEIR model was used to show control of measles by vaccination. Their study recommended introduction of mass vaccination programmer and improvement in early detection of measles cases to the minimize transmission.
In the work [18] performed a mathematical model of measles with vaccination and two phases of infectiousness. They followed the SIR modeling approach hence they partitioned the total population is into Susceptible, Infectious and Recovered compartments. Their study realized that the disease will certainly be eliminated if all susceptible are vaccinated. Achieving a 100% vaccination coverage is impractical but if the goal is set to 100% then we just might hit the ≥94% vaccine coverage which is the herd immunity for measles. Since measles is predominantly found among children aged 5 years and below, they therefore suggested that the measles vaccine should be made compulsory such that no child is allowed to enter school without evidence of at least two dose measles vaccination.
In the work [12] performed a deterministic mathematical model for transmission dynamics of measles it shows that vaccination is an important control strategy for the transmission of this disease.
Despite the availability of the measles vaccine since 1963, the infectious disease is still endemic in many parts of the world including developed nations. The disease has continued causing both economic and health problems to large population worldwide mostly affecting children. Due to these impacts, this study aims to develop a stochastic model for control and elimination of the transmission dynamics of measles. This research in tends to achieve the following specific objectives; to formulate a stochastic model for control and elimination of the transmission dynamics of measles, to compare the stochastic and deterministic mathematical models for measles and to perform simulation and analysis of the stochastic mode. The significance of this research is to show the strength of stochastic methods in the analysis compared to deterministic methods and to emphasize on the importance of vaccination in the control of transmission dynamics of measles.
2. Stochastic Differential Equations
In this chapter we will discuss random variables, stochastic process, Markov process, transition pdf, time homogeneous, discrete random walk, diffusion process, wiener process and Ito formula ([13] and [2]).
Suppose is a collection of random variables defined on a probability space, a stochastic process that is continuous in time the associated pdf is denoted by such that
If we assume is a stochastic process which is continuous in time and in state, then is a Markov Process if for any given sequence of times =
The transition pdf for a continuous time and state Markov process is the density function for a transition from state at time to state at time ,, it is denoted by
The transition pdf is said to be homogeneous or time homogeneous if so the transitions only depend on the length of time between states and the transition pdf is denoted as
If we consider a random walk on the set , let be the probabilities of moving to the right and to the left respectively such that , let be DTMC for this random walk where , it follows that by using Taylor series expansion about the point we obtain:
Letting approach zero, then represents the pdf of the continuous time and state process that satisfy the following partial differential equation which is also known as the forward kolmogorov equation:
where;
When the random walk is unbiased or symmetric, the limiting stochastic process is known as Brownian motion where since so that:
The standard Brownian motion with is also known as Wiener Process.
The assumptions on the limits in the random walk model were necessary to obtain the diffusion equation with drift. These assumptions are very important in obtaining the Kolmogorov Differential equationsm which are related to infinitesimal mean and variance of the process.
Let be a Markov Process with state space having continuous sample paths and transition pdf given by , then is a diffusion process if its pdf satisfy the following three assumptions in terms of the expectation:
where;
drift coefficient and diffusion coefficient
The backward Kolmogorov Differential Equation for a time homogeneous process is
The forward Kolmogorov Differential Equation for a time homogeneous process is
The pdf with is a solution of forward Kolmogorov equation and therefore, we can replace by , it follows from the forward Kolmogorov differential equation
A solution of a stochastic differential equation is a sample path of a diffusion process, if satisfies:
Then, the sample path of the process is a solution of the It SDE in the form
, where is denoted by Wiener Process.
Then the differential equation is equivalent to the It stochastic integral equation
where, the first integral is a Rieman integral and the second integral is an It stochastic integral.
If the drift and diffusion coefficient are and , then we obtain the following diffusion equation
for
which is the solution of the pdf in wiener process
which is the pdf of normal distribution with mean zero and variance t.
Suppose is a solution of the following It SDE:, if is a real valued function defined for and with continuous partial derivatives, then
where;
3. Model Formulation and Analysis
In this chapter we discuss the modified SEIR model with control strategy of vaccination [12] in deterministic and its corresponding SDEs for measles transmission. Mathematical Models of Epidemics can be broadly classified into two main categories: deterministic models and stochastic models. These models often result into nonlinear systems observed through partial noise data. After a deterministic system of ODEs has been formulated for the population dynamics it is possible to transform the system into various stochastic models such as DTMC, CTMC, MCMC, SDEs and others that take into account the uncertainties in the population.
There are many arguments for transforming deterministic to stochastic models. The following table shows the difference between deterministic and stochastic models [13]
DETERMINISTIC  STOCHASTIC  
1  They are simple to analyse but give less information  They are difficult to analyse though they are preferred when the analysis is possible 
2  Depends on the law of large numbers  There are cases where the law of large numbers is not satisfied 
3  Difficult to analyse extinction  Extinction analysis is possible since it occurs when epidemic deviate from mean level 
4  Hard to perform estimation because it is not probabilistic  The uncertainties makes it easier to perform estimations 
5  Preferable for new phenomena  Not preferable for new phenomena 
6  Repeated simulation from identical model result to one realization  Repeated simulation from identical model result to different realizations 
7  Depends on the assumption of mass action to describe spread of the disease  Depends on the natural way to describe spread of the disease 
3.1. Formulation of the Deterministic Mathematical Model
In this section we formulate a deterministic, compartmental mathematical model to describe the transmission dynamics of measles. We assume that the population is homogeneously mixing and reflects increasing dynamics such as birth and immigration, Per Capita birth rate is time constant, Per Capita natural mortality rate is time constant, individual can be infected through direct contact with an infectious individual, on recovery the individual obtains permanent infectionacquired immunity that is an individual cannot be infected again by measles and Individual who has attended first and second dose of vaccine consecutively receive permanent immunity to measles.
The total population is divided into the following epidemiological classes Susceptible, S (Individuals who may get the disease); Exposed or Latent, E (Individuals who are exposed to the disease); Infected, I (Individuals who have the disease and are able to transfer it to others); Recovered, R (Individuals who have permanent infectionacquired immunity and those who received the second dose of vaccine) and Vaccinated, V (Individuals who have received first dose of vaccine). In this study we assume that newborns and immigrants who received first dose of vaccine join the Vaccinated class, V at rate ϕπN and ρΛ respectively. On the other hand newborns and immigrants who have not received first dose of vaccine join the Susceptible class, S at rate (1ϕ)πN and (1ρ)Λ respectively. The Susceptible individuals who received the first dose of vaccine may join the Vaccinated class, V at rate εS. When there is an adequate contact of a Susceptible individual with an Infective individual so that transmission occurs, then the susceptible individual may join the Exposed class, E at the rate λS. After Latent period ends, exposed individuals may progress to the Infectious class, I at rate σE. When the infectious period ends, the individuals may join the recovery class, R at rate ηI otherwise they may die at rate (δ+μ)I due to the disease or naturally. The recovery class, R consists of those with permanent infectionacquired immunity and those who received the second dose of vaccine at rate ωθV. The Vaccinated individuals who did not receive second dose of vaccine may return to Susceptible class, S at rate (1θ)αV as the first dose of vaccine is waning at rate α.
3.1.1. Description of Variables and Parameters
The following tables describe the variables and parameters used in our model:
VARIABLE  DESCRIPTION 
S  The number of Susceptible individuals at time t 
E  The number of Exposed individuals at time t 
I  The number of Infected individuals at time t 
R  The number of Recovered individuals at time t 
V  The number of Vaccinated individuals at time t 
N  The total population at time t 
PARAMETER  DESCRIPTION 
π  Per Capita birth rate 
Λ  Constant Immigration rate 
ϕ  Proportions of newborns who are vaccinated 
1ϕ  Proportion of newborns who are not vaccinated 
ρ  Proportions of immigrants who are vaccinated 
1ρ  Proportion of immigrants who are not vaccinated 
a  Arrival rate 
c  Per Capita contact rate 
δ  Death due to disease 
μ  Per Capita natural mortality rate 
β  Probability of one infected individual to become infectious 
λ  Force Of Infection, λ= 
σ  Progression rate from latent to infectious 
η  Recovery rate of treated infectious individuals 
α  The rate of waning of first dose of vaccine 
ω  The rate of receiving second dose of vaccine 
ε  Proportion of individuals who received a first dose of vaccine 
θ  Proportion of individuals who received a second dose of vaccine at rate α 
1θ  Proportion of individuals who are not vaccinated for a second time and return to Susceptible class 
3.1.2. Compartmental Diagram
The description of measles dynamics can be summarized by compartmental diagram below:
3.1.3. Ordinary Differential Equations
From the above explanation and compartmental diagram Figure 1, the transition between compartments can now be expressed by the following differential equations:
(1)
(2)
(3)
(4)
(5)
The above system of equations (1)(5) can be simplified into:
(6)
(7)
(8)
(9)
(10)
where is the force of infection given by
(11)
The total population size is:
Where by adding the system of equations (6)(10) we get:
(12)
3.1.4. The Basic Reproduction Number, R0
Epidemiologists have always been interested in finding the basic reproduction number of an emerging disease because this threshold parameter can tell whether a disease will die out or persist in a population. Denoted by, this parameter is arguably the most important quantity in infectious disease epidemiology. It is defined as the average number of new cases (infections) produced by a single infective when introduced into a susceptible population. It is one of the first quantities estimated for emerging infectious diseases in outbreak situations (Diekmann et al, 2009). It is a key epidemiological quantity, because it determines the size and duration of epidemics and is an important factor in determining targets for vaccination coverage (Grais et al, 2006). The basic reproduction number is sought after principally because:
If < 1, then throughout the infectious period, each infective will produce less than one new infective on the average. This in turn implies that the disease will die out as the DFE is stable.
If > 1, then throughout the infectious period, each infective will produce more than one new infective on the average. This in turn implies that the disease will persist as the DFE is unstable. In other words, there will be an outbreak.
If can be determined, then the transmission parameters which will force to be less than or greater than 1 can easily be identified and control measures effectively designed.
Next, we shall find the Basic Reproduction Number of the system using the next generation method [6] and we obtain the basic reproduction number denoted by equation (13) below:
(13)
3.2. Formulation of Stochastic Mathematical Model
In this section we will transform the deterministic model to stochastic model. The corresponding Stochastic Differential Equations model for control of the transmission dynamics of measles with vaccination will be formulated by construction of equivalent Stochastic Differential Equations Models [13].
Let denote the number of individuals for respectively
Thus we can have
(14)
(15)
(16)
(17)
(18)
Define , to transform the above deterministic model we need to compute the expectations and . The table below contains the possible changes of the process together with their associated transition probabilities.
POSSIBLE CHANGES  PROBABILITIES 


























From the table above we formulate three systems of SDEs that have the same Forward Kolmogorov Differential equation and sample paths:
First SDE model:
is a vector of five independent wiener processes
The drift part and the diffusion part are computed by the following neglecting higher orders of (Yuan and Allen, 2011)
If we define
and
where; and
Thus computing from the table we obtain:
This is the same as the RHS of the deterministic model.
The resulting SDE is a multidimensional and multiplicative, hence it is hard to simulate. Also, the diffusion coefficient of the formed SDE is a square root of a matrix, that is
Second SDE model:
is a vector of thirteen (13) independent wiener processes
Since matrix is hermitian (symmetric) and positive definite, cholesky decomposition can be applied obtain matrix such that; =V
The diffusion matrix G is a 5×13 matrix of the form:
Third SDE model:
is a vector of five independent wiener processes
Column reduction can be performed on the matrix ∆X to obtain matrix H such that;
The diffusion matrix is a 5×5 matrix of the form:
4. Simulation and Discussion
In this chapter we simulate the deterministic and stochastic models formulated.
For simplicity this study will perform simulation the Third SDE model at thus we have the following system of nonlinear equations:
is a vector of five independent wiener processes
PARAMETER  DESCRIPTION  VALUE/RANGE  SOURCE 
π  Per Capita birth rate  0.02755 per year  [9] 
ϕ  Proportions of newborns who are vaccinated  0.5 also varies with scenario (0.0 – 1.0)  [10] 
ρ  Proportions of immigrants who are vaccinated  0.7 also varies with scenario (0.0 – 1.0)  [10] 
a  Arrival rate  0.02755 per year  [9] 
c  Per Capita contact rate  0.09091 per year  [9] 
δ  Death due to disease  0.125 per year  [4] 
μ  Per Capita natural mortality rate  0.00875 per year  [9] 
β  Probability of one infected individual to become infectious  varies with scenario (0.08 – 0.7)  [5] 
λ  Force Of Infection, λ=  0.096 per year  [10] 
σ  Progression rate from latent to infectious  0.125 per year  [9] 
η  Recovery rate of treated infectious individuals  0.14286 per year varies with scenario (0.0 – 1.0)  [9] 
α  The rate of waning of first dose of vaccine  0.167 per year  [10] 
ω  The rate of receiving second dose of vaccine  0.8 per year  [10] 
ε  Proportion of individuals who received a first dose of vaccine  0.7 per year varies with scenario (0.0 – 1.0)  [10] 
θ  Proportion of individuals who received a second dose of vaccine  0.5 varies with scenario (0.0 – 1.0)  [10] 
From Figure 2 we observe that individuals who have received first and second dose of vaccine and those who have been treated which comprise the recovered population in the sense that they can’t be infected is increasing hence the infectious disease is controlled in the population.
The Figure 3 shows that Susceptible population is decreasing while the exposed and infected ones are increasing from the beginning till a certain time, t where it starts decreasing to be zero. This means that vaccination is controlling the infectious disease and eradicated to die out from the population hence the recovered population is increasing.
Figure 4 shows that the stochastic model describe the dynamics of measles transmission more effectively.
Figure 5 shows that the Susceptible population varies as time goes on but eventually decreases hence the disease is controlled in the population.
Figure 6 shows that the exposed population is decreasing since vaccination is controlling the spread of the disease.
Figure 7 shows that the infectious population increasing and decreasing in the population till a certain time where it starts decreasing to be zero because of vaccination control strategy.
Figure 8 shows the variation of the recovered population which shows the vaccination strategy is controlling the disease, since as time goes on the recovered population is increasing.
Figure 9 shows that the variation of the vaccinated population
5. Conclusion, Recommendation and Future Work
In this chapter we provide conclusion and recommendation of the project.
5.1. Conclusion
The study has shown the effectiveness of stochastic analysis in studying the dynamics of measles compared to deterministic analysis.
5.2. Recommendation
This study recommend the use of stochastic analysis in studying dynamics of infectious diseases.
5.3. Future Work
Based on the model of this study, it is proposed that future work should consider the following:
i. Application of other stochastic approaches such as Monte Carlo Markov Chains
ii. Use of more accurate numerical schemes like Rungekutta and others
List of Abbreviations and Symbols
AIDS  Acquire Immune Deficiency Syndrome 
CFRs  Case Fatality Ratios 
CTMC  Continuous Time Markov Chain 
DTMC  Discrete Time Markov Chain 
HIV  Human Immunodeficiency Virus 
MCC  Measles Control Campaign 
MCMC  Monte Carlo Markov Chain 
MCV1  First Dose of Measles Containing Vaccine 
MCV2  Second Dose of Measles Containing Vaccine 
MMR  Measles MumpsRubella 
ODEs  Ordinary Differential Equations 
Probability Density Function  
SEIR  Susceptible Latent Period Infectious Recovery 
SIA  Supplemental Immunization Activity 
SDEs  Stochastic Differential Equations 
SIR  Susceptible Infectious Recovery 
U.S  United States 
UNICEF  United Nation Children’s Fund 
WHO  World Health Organization 
References