1. Introduction
Measles is still a major worldwide health issue, particularly in the less developed countries. Measles (also known as Rubella or morbilli) is an extremely contagious illness caused mostly by the Morbillivirus genus in the Paramyxovirus [
1,
2]. Although efficient vaccines against this severe illness are commonly accessible, still measles is a leading cause of death among children below five of years of age [
3]. The disease infecting hundreds of millions of children each year and resulting in a high mortality and morbidity in the child population, owing primarily to complicated conditions that exist side by side with the disease like malnutrition, diarrhea, and pneumonia [
4]. Sneezing and coughing, touching the nasal or aerosol fluids, or close physical touch with an infected person are all ways to spread measles. It can stay extremely contagious for a maximum of two hours in the atmosphere and on the surfaces. Clinical manifestations include soar throat, cough, nasal congestion, blurred vision, and small white patches in the mouth; these signs and symptoms often develop within 10–12 days post-infection. Subsequently, a rash appears, extending downward out of the nose. The period of peak infectivity (virus shed) starts four days before that and four days just after commencement of the rashes. The usual incubation time is 14 days; however, it can range from 7–18 days [
5]. In reality, even vaccinated people may still be susceptible if the immunization fails or existing immunity from the vaccine wears off. Despite the fact that vaccination has cut worldwide measles fatalities by 73% between 2000–2018, measles continues to remain a widespread in many underdeveloped nations, particularly in Asia and Africa. Around 140,000 individuals died from the measles in 2018. During 2000–2018, worldwide measles immunization results in an 85 percent drop in measles transmission [
6,
7]. Despite the abundance of a safe and effective vaccination in 2017, around 110,000 deaths occurred from measles, primarily children under the age of six, as per the report of the World Health Organization (WHO) [
8]. Vaccination is amongst the most successful health promotion strategy for reducing death rates and the spreading of epidemics; it has been demonstrated that vaccination saves over 3 million people only in Nigeria every year. The vaccination will work with the immune system of the body to establish the body’s natural defenses, reducing the likelihood of relapse [
9]. The MMR vaccination can protect against measles, which is a vaccine-preventable illness. The MMR vaccines are highly effective at protecting both adults and children from the epidemic measles. Only one dose of the MMR vaccines is roughly 92% successful in suppressing the measles, whereas two doses are around 95% effective. The MMR vaccine also protects against rubella and mumps [
10]. This disease is an endemic one in Nigeria, with epidemics occurring at regular intervals. Measles is present across Nigeria at all seasons; however, it is more widespread in the summer months.
Pakistan is one of the most measles-affected nations in Asia [
11]. Every 8–10 years, the nation has a recurring measles epidemic. In fact, 2845 identified measles infections were reported in Pakistan during the year 2016. This figure increased to 6791 in 2017 and, in the year 2018, 33,007 cases were reported. These results represent about 44, 20, and 51 percent of the total number of cases recorded in the East Mediterranean, which includes 22 nations. In 2017, over 130 children lost their lives due to this infection, with the figure rising to nearly 300 in 2018 [
12].
It is strongly advised to employ the methods of mathematical modeling to examine the transmission process and prevention of an epidemic disease [
7,
13,
14,
15], modeling with fractional differential equations also have several applications in all fields of science [
16,
17,
18]. While depicting the natural history of an infectious disease, the tools of modeling can create a balance among the data and its real biology. Thus far, models responsible for describing the dynamics of measles both from population to outbreak levels have demonstrated a wide variety of disease patterns. External noise is usually the main significant feature of physical processes and bio-systems. It has been discovered that environmental variables have a significant impact on the dynamics of measles disease spread [
19]. Because of the uncertainty of person-to-person interactions or inherent characteristics of the population, outbreak onset and propagation are fundamentally unpredictable. As a result, the condition of the disease is influenced by the environment’s heterogeneity and uncertainty.
Changes in the environment likewise have a significant impact on the parasites’ persistence and distribution. Because the stochasticity of parameters and states depicts the exact dynamic behavior of an infectious disease, it is regarded as an essential part in epidemiological studies. Even though the perturbations are varying, these should be autocorrelated in a positive way. Furthermore, the perturbations may be estimated theoretically using the linked problem’s probability density function [
20,
21,
22,
23]. There are two main techniques to epidemic modeling: the deterministic modeling and stochastic modeling. Stochastic differential equation (SDE) models are recommended over deterministic models for mathematical modeling of biological functions because they may provide a higher level of reality than their deterministic equivalents [
23,
24,
25,
26]. We may choose to use SDEs to generate a distribution of the predicted output(s); for example, the number of infected individuals over time
t. Moreover, when tested numerous times, a stochastic model produces more useful outputs than a deterministic one. A deterministic model, on the other hand, will produce only one outcome irrespective of the number of experiments. To explain the viral evolution of COVID-19, numerous deterministic infectious disease models have been suggested; for example, see [
27,
28].
The rest of the manuscript is organized in the following manner:
Section 2 deals with the formulation of the proposed stochastic model for the spreading of the epidemic measles. The uniqueness and existence problem for obtaining a positive global solution is presented in
Section 3. In
Section 4 and
Section 5, we characterized sufficient criteria for the stationary distribution and extinction of the disease. We optimize the proposed model using data from Pakistan compiled in the first ten months of 2019 in
Section 6. We quantitatively compared the obtained analytical results, and graphical illustrations were presented in
Section 7. We concluded the work in
Section 8 by presenting the future directions and a comprehensive summary.
2. Model Formulation
Olumuyiwa et al. [
29] have recently developed a model of the transmission of rubella disease by using the differential equations. Keeping in view the work of Olumuyiwa et al., here we intend to extend the model to a stochastic model. Furthermore, by considering different stages during the measles epidemic, we stratified the total population into six disjoint classes, namely: susceptible, vaccinated, exposed, infectious, hospitalized and recovered individuals whose sizes in mathematical terms are, respectively,
, and
.
The entrance of new persons through this population is captured by the rate
and will be kept in the susceptible compartment. People in the vulnerable group start receiving a vaccination at a rate
and setback immune function at a vaccine wane rate
. The contact rate of susceptible persons is
, and thus the term force of infection becomes
, with the transition again from exposures to infection stages occurring at a rate of
. People who are infected with the measles seek medical attention at a rate of
and recover from the infection after the successful treatment supplemented at a rate of
. We consider two types of death rates: the natural
(that is constant for all classes) and the disease-induced mortality rate
. This study assumes the recovery from measles that is possible due to the treatment only, that is, the study considering no natural recovery. The movements between the compartments are depicted via a flowchart in
Figure 1. The above assumptions will lead to the following deterministic system:
The threshold parameter has the following expression for model (
1) as
In order to consider the noises due the environment in the study (i.e., model (
1)), we shall take into account the standard Brownian motions
for
with
. Furthermore, system (
1) is modified by considering the incidence rate of bi-linear form
. For the noise intensities, we have taken
as the relative weights. By considering these stochastic terms, the deterministic model (
1) takes the form
Keeping in view system (
3), the authors have a keen interest to find the possible answers to the following questions:
- Q1:
What role do the random noises play in the transmission measles?
- Q2:
What role contaminated vaccination in the spreading of measles disease?
- Q3:
Under what condition(s) will the disease tend to go extinct?
- Q4:
Under what condition(s) will the epidemic measles persist in the population?
4. Extinction
While modeling the dynamical aspects of any epidemic diseases, it is important to investigate the situations under which the disease will become exterminated or tend to become extinct in the community. Within this section, we demonstrate that, when the white noises are large enough, the solution of the associated stochastic model (
3) surely vanishes. Let us define
Lemma 1. (Strong Law)[
32,
33]
Let be continuous and real valued along with a local martingale, which vanishes as , then Lemma 2. Assume that corresponds to initial data in the space and is a solution of model (3). Then, Furthermore, if , and , then Then, the solution of (3) To prove Lemma 2, we follow almost the same techniques as performed in proving Lemmas 2.1 and 2.2 carried out in the work of Zhao [
32], and thus the readers can simply prove the results.
Next, to develop the extinction theory related to model (
3), we first define the threshold quantity
for the proposed stochastic model which can be written in the form of
Theorem 2. The exposed population and infected population of system ((3) tends to zero exponentially almost surely if , where is given by Equation (12). Proof. Let
correspond to initial data
in the space
, being a solution of model (
3).
Differentiating Equation (
13) following Ito’s formula, one can obtain
By taking the integral of either sides of the above inequality over the interval
, we have
By assuming the
as
and multiplying both sides of relation (
15) by
while considering Lemma 2, we obtain
If , then if . Again, and we assert that —hence the result. □
5. The Stationary Distribution of the Disease
We understand that there are no endemic states in stochastic models. As a result, the stability of the system cannot be used to predict how long an illness would endure. As a result, one must concentrate on the existence/uniqueness assumption of the stationary distribution. In certain aspects, this assists the illness with survival. For this purpose, we employ a well-known method due to Hasminskii [
34].
Assume a homogeneously time-dependent Markov process
in the space
that satisfies the below stochastic model
The diffusion matrix can be demonstrated as:
Lemma 3. [
34]
. The process has a one and only one stationary distribution whenever there is a bounded domain having a regular boundary in such a way that closure and having the below characteristicsThe smallest eigenvalue of the matrix is bounded below from the origin in the open domain U and in its neighborhood.
For the average time period τ (at which a path starts from ϰ reaching the set U) is bounded, and for every compact subset , we have . When is an integrable function with the measure , then we have
Let us define the following parameter for our future use:
Theorem 3. For , then a solution of system (3) is ergodic and has one and only one stationary distribution . Proof. For verifying condition (2) of Lemma (3), we must develop a positive
function
To do so, we first formulate
where
and
are all real and positive constants, and must be calculated later on. By assuming the formula due to Itô’s and keeping in view system (
3), we obtain
□
The preceding calculation indicates that
Furthermore, we obtain
where
is an unknown real number and must be determined later. It is very helpful to state
where
.
In the below steps, we show that actually the function
has the minimum value,
. The partial derivatives of the function
with respect to the state variables are given by:
One can verify very easily that the function
has only one stagnation point.
Furthermore, at
for
, the Hessian matrix is as follows:
The Hessian matrix is obviously positive definite. As a result, the minimum value of
is
. We may assert that
has one and only one minimum value
inside
based on Equation (
20) and the continuity of
. Then, as follows, we delineate a non-negative
function
:
Using
formula and the suggested system, we arrive at
Based on the above result, we have the following assertion:
where
Overall, for the solution to the model, we have the following set
where
are extremely small positive real values to be calculated later on. For the sake of simplicity, the whole set was partitioned into the following subsets:
Then, we show that on , which is the same is displaying the solution on the eight sub-regions.
Case 1. If
, then, by Equation (
24), we obtain
For every , picking returns .
Just as in the proof above, we conclude that for , .
As a result, we arrived to the point that there must be positive constant
in such a way
Consider
, and
is the amount of time it takes for a path to start from
to achieve set
D,
The next relation could be obtained if one integrates the above inequality from 0 to
, considering the expectation and using Dynkin’s formula:
Hence,
is non-negative; then,
We have
as a result of the proof of Theorem (3). This also shows that model (
3) is regular and, consequently, by letting
, almost surely we have
.
Moreover, by utilizing the Fatou’s lemma, we have
where K is the compact set, i.e., a subset of . This proves the condition (2) of Lemma (3) in an alternative approach.
Additionally, the diffusion matrix for the system (
3) is
Taking
we obtain
where
Similarly, this proves condition (1) of Lemma (3). Keeping in view the previous statements, the ergodicity of model (
3) is insured by Lemma (3) and hence proving that the model has stationary distribution.
6. Parameter Estimation
Among the most important processes to carry through out the testing process is the utilization of real data (if available) to acquire findings for certain unidentified biological factors used in the epidemiology system under study. Real-world measles cases, as shown in
Table 1, are used to test the proposed rubella model and to choose the best fitted parameters for numerous unknown biological parameters that emerge in the system. Considering the 2018 statistics of WHO, the natural death rate of a Pakistani individual is
since the life expectancy of a Pakistani is
. In addition, the entire size of the country is
, whereas the recruitment rate is calculated to be
. Additionally, Memon et al. [
13] indicates that the rate of measles vaccination is generally 97 percent effective, implying that the vaccines’ effectiveness
equals
. Some values of the parameters are estimated and others were fitted against the data, and these were presented in
Table 2. In
Figure 2, the results via simulations for measles resurgence cases obtained by adapting the proposed model (
3) with real data from the first ten months of the year 2019 are shown. As shown in
Figure 2, it gives a rather good match, lending veracity to the predictions generated from the proposed measles model (
3). We employed the relation
to measure the associated relative error for fitting the model against the data.
7. Numerical Simulations and Discussion
It is of great concern to simulate each mathematical model against the real data and to verify theoretical results, and this step is very important when dealing with modeling biological problems. The researchers seek to simulate problem (
3) employing classic numerical algorithms that rapidly converged. Almost all of the theoretical conclusions are quantitatively validated by using the well-known RK-4 approach.
where
stands for the standard Gaussian stochastic variables, having distribution
and
is the constant time-step. The terms
reflect the intensities of the white noises.
To quantitatively validate the analytical conclusions, we require the parameters’ value being used in model (
3). In Example 1 and 2, we established two sets of parameters’ values for this reason, and the population levels of each compartment at
were displayed. For every scenario, we simulated the model over the interval
.
We formulated Theorem 2 based on the stochastic theory of stability, which indicates that, under the condition of
, the infection would continue to perish from the community regardless of the levels of the variables at
. Furthermore, the theory demonstrates that the infection will be eradicated from the community with unit probability.
Figure 3 shows that the random curves reach the infection-free state after a limited period of time, confirming the analytical results.
Example 1. The values of the parameter are assumed as: and , where the initial values of the state vector are . Similarly, the intensities of the white noises are: . Using those same model parameters, we estimated , which was found to be lower than 1. As a result, the assumption of Theorem 2 is met, and hence the components of the solution of the model adhere to the following relations: Essentially, these relations explain the elimination of the virus from the community, as seen by Figure 3. As a consequence, the derived research findings on elimination are valid and may be relied on. Likewise, Theorem 3 guarantees the disease’s prevalence in the population given permissible limits. By considering data from Example 2, we calculated the value of
, and it was found that it is greater than unity.
Figure 4 depicts the numerical results based on the theorem’s assumptions. The figure implies that the infection persists inside the population, and that there is persistence of the solution of the proposed model (
3). This verifies the conclusion of Theorem 3 in the case of deterministic model (
1). These results further elaborate that, when the related threshold parameter of the perturbed system exceeds unity, the solution of the model (
3) lies within the neighborhood of endemic equilibrium. Thus, policymakers must provide strong preventative measures against various variations in order to restrict the spreading of multiple strains throughout the community. Moreover, Theorem 3 ensures the existence of an ergodic stationary distribution for model (
3), and it is confirmed by
Figure 5.
Example 2. In this case, the chosen parameter values are of the form: and . Similarly, the initial population sizes of the state variables are , whereas the intensities are given by . Considering the measles data and hence the estimated and fitted parameters, we find that exceeds unity. It is also explored that the model parameters taken in this example satisfy the premise of Theorem 3. We tested the model using this input, and the outcomes are displayed visually in Figure 4. The figure shows that the disease tends to stay inside the community, and the model exhibits a homogeneous stationary distribution in this situation. From Theorem 3, there is a single stationary distribution of system (
3). To numerically illustrate this statistical property, we present in
Figure 4 and
Figure 5, the trajectories and the associated joint density function for each class of the population. For a good visibility, we offer the
and upper view of the aforementioned joint densities in
Figure 6,
Figure 7,
Figure 8,
Figure 9 and
Figure 10. This indicates that the infection is still present in the population over time. We talk here about persistence in the mean of the epidemic. In
Figure 4, we illustrate the continuation of all groups of the studied population. We remark that the stochastic trajectories fluctuate around the deterministic solution with reasonable distances according to magnitude of the noises.
8. Conclusions and Future Research
In this study, we presented a detail analysis of a stochastic model that describes the spreading behavior of the measles virus while accounting for the white noises and the influence of immunizations. It is assumed that the randomness being used in the model is nonlinear, and that a person may lose the resistance after vaccination, implying that vaccination might create temporary protection against the disease. First of all, we formulated a deterministic model and then it was expanded to a stochastic model. It is shown that the stochastic model is both theoretically and practically viable by demonstrating that the model has a global solution which is positive and stochastically bounded. Next, we developed sufficient criteria for the disease’s elimination and permanence. Furthermore, the presence of a stationary distribution is examined by constructing a suitable Lyapunov function, wherein we noticed that the disease will persist for and that the illness will vanish from the community for . We simulated the model against the available data of measles in Pakistan during the first ten months of 2019, by using the conventional curve fitting methods, and the values of the parameters were calculated therein. These values of the parameters were used in simulating the model, and the theoretical findings of the research were evaluated and conclusions were made. Simulations of the study suggest that, in order to fully understand the dynamic behavior of the measles epidemic, time-delay must be included in such analyses, and that advancements in every vaccine campaign are unavoidable in order to stop or minimize the spreading of the disease.
We fit both stochastic and determinism models to known data on rubella in Pakistan during the first ten months of 2019, and the values of parameters were obtained using lsqcurvefit methods. These model parameters are used, and the threshold number, which is around 13, is determined. This demonstrates that measles is extremely harmful and might have a negative impact on this community if adequate control tactics are not implemented in time. It also suggests that, if no appropriate measures are implemented, the cases of the measles may increase in the coming years. To further reduce the transmission rate, health authorities and lawmakers must launch awareness campaigns, mass immunizations, particularly among youngsters, and, most crucially, care and treatment for people having chronic conditions. This was also discovered that serves as the most critical indicator to the threshold parameter; thus, lowering the disease spreading co-efficient (by lowering the household and sexual interactions of chronic patients with vulnerable) for acutely infected individuals who become chronic is also an effective control to reduce the spread of measles.
In the next research studies, the authors intend to use fractional calculus and modify this and the related model while utilizing different definitions of fractional derivatives such as Riesz, Caputo, Atangana–Baleanu, Caputo–Fabrizio, and many others for capturing the real dynamics of such and related diseases.