Next Article in Journal
A New Load Frequency Control Technique for Hybrid Maritime Microgrids: Sophisticated Structure of Fractional-Order PIDA Controller
Previous Article in Journal
Fractional-Order Model-Free Predictive Control for Voltage Source Inverters
Previous Article in Special Issue
First-Passage Times and Optimal Control of Integrated Jump-Diffusion Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Lévy Noise with Infinite Activity on the Dynamics of Measles Epidemics

1
School of Science, Hunan University of Technology, Zhuzhou 412007, China
2
School of Statistics and Mathematics, Guangdong University of Finance and Economics, Guangzhou 510320, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(6), 434; https://doi.org/10.3390/fractalfract7060434
Submission received: 17 April 2023 / Revised: 21 May 2023 / Accepted: 26 May 2023 / Published: 27 May 2023
(This article belongs to the Special Issue Stochastic Modeling in Biological System)

Abstract

:
This research article investigates the application of Lévy noise to understand the dynamic aspects of measles epidemic modeling and seeks to explain the impact of vaccines on the spread of the disease. After model formulation, the study utilises uniqueness and existence techniques to derive a positive solution to the underlying stochastic model. The Lyapunov function is used to investigate the stability results associated with the proposed stochastic model. The model’s dynamic characteristics are analyzed in the vicinity of the infection-free and endemic states of the associated ODEs model. The stochastic threshold R s that ensures disease’s extinction whenever R s < 1 is calculated. We utilized data from Pakistan in 2019 to estimate the parameters of the model and conducted simulations to forecast the future behavior of the disease. The results were compared to actual data using standard curve fitting tools.

1. Introduction

The measles infection, also called rubella, is extremely contagious and spreads throughout the globe. It is caused by the virus Morbilli, related to the community of Paramyxoviridae [1,2]. Vaccination is a powerful tool in reducing the impact of the disease over time. However, the force of mortality continues to affect young children under the age of five [3]. Measles infects ten million children each year, and millions of them die due to factors such as the unavailability of a balanced diet, weak digestion, and pneumonia [4]. Disease spread is caused by personal contact, sneezing, coughing, and contact with airborne and nasal droplets. The virus stays active for approximately two hours, and remains highly contagious during this period. The initial symptoms of the disease include nasal discharge, sore throat, and the presence of small white spots on the tongue. In later stages, coughing may occur. The typical incubation period for this infection is about four days before the appearance of the rash, and lasts for nearly five days after the rash appears. The mean incubation time is about fourteen or fifteen days, changing with external factors such as the environment and weather of the surroundings [5]. In real practice, vaccinated people may suffer side-effects, or it may be the case that vaccines are not available. On the other hand, vaccination acts as a preventive measure against measles, reducing the incidence rate by up to 73% over the past eighteen years. According to the World Health Organisation, Measles is present in many developed countries around the world, particularly in Asia and Africa. More than one hundred and fifty thousand people died due to measles in 2018. In the meantime, the death to infection ratio has decreased by about 85 percent [6,7]. According to reports from the WHO, nearly 110,000 people died in 2017 due to measles, with a majority of them being children under six years old. This highlights the need for effective and safe vaccination [8]. The discovery of vaccines for infectious diseases has lowered the ratio of death and infection to a large extent. This process saves about three hundred thousand individuals from death each year. Vaccines produce antibodies that act against reinfection and strengthen the immune system [9]. Specifically, measles can be controlled by the MMR vaccines. These vaccines are considered very safe for both children and the elderly, and can reduce the severity of measles infection to a great extent. According to studies, one dose of the MMR vaccine can prevent disease in around 91% of cases, while two doses can provide approximately 95% protection. Other diseases such as the mumps can be controlled by this vaccine as well [10]. Nigeria experiences periodic outbreaks of measles, as the disease is endemic in the country. The disease spreads in Nigeria during all seasons of the year, and is highly contiguous during dry weather.
In addition, Pakistan is among the countries in the WHO’s Eastern Mediterranean Region most highly burdened by measles [11]. Over the past 9 to 12 years, there has been a sharp increase in measles outbreaks throughout Pakistan; some 2845 cases of measles were reported in Pakistan in 2016. This number increased to 6791 in 2017 and 33,007 in 2018, which represents approximately 43%, 21%, and 50% of all cases reported in the Eastern Mediterranean region comprising 23 countries [11]. Nearly 129 children died from measles in 2017, and in 2018 the number increased to about 300 [12].
Mathematical modeling is regarded as a significant and powerful approach for examining and forecasting the dynamic patterns of epidemics [7,12,13,14]. Most of the models used to date employ the integer order derivative, and their analysis is related to classical theory. The foundation of a mathematical modeling relies heavily on the availability of biological information and data regarding the epidemic under consideration. Mathematical models for the measles have wide application for understanding its spreading dynamics and controlling it. White noise plays a vital role in explaining the dynamics of physical and biological problems. The impact of external environmental factors such as white noise on the dynamics of measles epidemics is significant [15]. Due to the non-homogeneity and nonlinearity of population interactions and other complexities, epidemic prediction cannot be depicted using traditional modeling approaches. Moreover, the dynamics and control of various epidemics may be affected by global environmental factors.
The human population is subject to many complex and random variations in the real world. Therefore, stochastic models are a more appropriate technique for modelling epidemics. It has been demonstrated that stochastic models are more realistic compared to deterministic models. Several researchers have recently focused on perturbations that can capture the true dynamics of the epidemic based on stochastic modelling [13,15,16]. The Lévy white noise is important at various velocities for the spacing of threshold parameters. Adding Lévy noise terms can lead to more realistic and accurate results in stochastically analyzed models for various infectious diseases. The stochastic version of the underlying deterministic mode can be obtained by including this noise in the aforementioned systems. Typically, two types of noise are utilized in such models, namely, Lévy noise and Gaussian noise. The Lévy noise is more suitable than the Gaussian noise because the resulting system is able to model systems with a higher degree of complexity [16,17,18,19]. Due to fluctuations in diffusion problems, disturbances cannot be described by the continuous stochastic model; thus, it is important to model these phenomena using jump processes. For this reason, we consider the problem using Lévy noise.
The remainder of this manuscript is structured as follows. In Section 2, we propose a model based on random processes for the transmission dynamics of measles. In Section 3, we present the dynamic features of the globalized positive model’s solution. In Section 4 and Section 5, we obtain necessary conditions related to disease elimination and persistence of the proposed stochastic model. We optimize the proposed problem using the cases of measles infection in the population of Pakistan for January to October 2019 in Section 6. In Section 7, we verify the theory behind the obtained results qualitatively and quantitatively, and provide numerical simulations. Finally, we conclude our analysis in Section 8 with remarks and recommendations for future research work.

2. Model Formulation

Olumuyiwa et al. [20] recently developed a mathematical model for describing the infection dynamics of the epidemic measles by using a deterministic approach. The authors employed the epidemiological concept of population and divided the entire human population into six classes, including susceptible, vaccinated, exposed, infectious, hospitalised, and recovered individuals, respectively denoted by S ( t ) , V ( t ) , E ( t ) , I ( t ) , H ( t ) , and R ( t ) . The addition to the susceptible population per day is provided by the rate ϕ . The vulnerable class is subject to a vaccination rate τ and loses immunity against the disease at a rate ω , which is commonly known as the “waning rate of vaccination”. The rate of infection for the susceptible class is denoted by α . Therefore, the term α S I represents the total infection rate per time unit. Additionally, the transition from the exposed class to the infected class is represented by β . Individuals in the infected population are hospitalized at a rate ρ and then recover from measles at a rate γ . The model assumes a constant natural death rate, denoted by μ , for all population classes. In addition, the infection-related death rate is provided by the parameter δ . Here, we do not assume any natural recovery from measles infection. The model chart is illustrated in Figure 1 [21]. The mathematical formulation of the above discussion can be converted into a model of ordinary differential equations (ODEs), represented in symbolic form as follows:
d S ( t ) t = α I ( t ) S ( t ) + ϕ ( μ + τ ) S ( t ) + ω V ( t ) , d V ( t ) t = ( ω + μ ) V ( t ) + τ S ( t ) , d E ( t ) t = ( β + μ ) E ( t ) + α I ( t ) S ( t ) , d I ( t ) t = ( μ + ρ + δ ) I ( t ) + β E ( t ) , d H ( t ) t = ( μ + γ + δ ) H ( t ) + ρ I ( t ) , d R ( t ) t = μ R ( t ) + γ H ( t ) .
The threshold parameter for the model (1) is obtained using standard techniques, and has the following form:
R 0 D = ( μ + ω ) ϕ β α ( μ + β ) ( μ + δ + ρ ) ( μ + ω + τ ) μ .
The main theme of the present manuscript is to modify model (1) by including white and Lévy noise as well as the incidence rate of nonlinear shapes; the white noise is used for the continuity part and the Lévy noise for the jumping part. The extended form of the deterministic system can be written in stochastic form as follows:
d S = ϕ α S ( t ) I ( t ) N ( t ) + ω V ( t ) ( τ + μ ) S ( t ) d t + ξ 1 S ( t ) d W 1 ( t ) + Y X 1 ( y ) S t N ˜ ( d t , d y ) , d V = τ S ( t ) ( μ + ω ) V ( t ) d t + ξ 2 V d W 2 ( t ) + Y X 2 ( y ) V t N ˜ ( d t , d y ) , d E = α S ( t ) I ( t ) N ( t ) ( μ + β ) E ( t ) d t + ξ 3 E ( t ) d W 3 ( t ) + Y X 3 ( y ) E t N ˜ ( d t , d y ) , d I = β E ( t ) ( ρ + δ + μ ) I ( t ) d t + ξ 4 I ( t ) d W 4 ( t ) + Y X 4 ( y ) I t N ˜ ( d t , d y ) , d H = ρ I ( t ) ( γ + δ + μ ) H ( t ) d t + ξ 5 I ( t ) d W 5 ( t ) + Y X 5 ( y ) H t N ˜ ( d t , d y ) , d R = γ H ( t ) μ R ( t ) d t + ξ 6 R ( t ) d W 6 ( t ) + Y X 6 ( y ) R t N ˜ ( d t , d y ) ,
where W i ( t ) for i = 1 , , 6 is standard Brownian motion defined in a complete probability space ( Ω , F , P ) with filtration F t t 0 , satisfying the usual condition ξ 1 , ξ 2 , ξ 3 , ξ 4 , ξ 5 , with ξ 6 representing the intensity of noise, S t , V t , E t , I t , H t and R t being the left limit of S , V , E , I , H and R , respectively, N ˜ = N ( d t , d y ) v ( d y ) d t and N ( d y , d t ) a Poisson counting measure with characteristic measure ν on the measurable subset Y of [ 0 , ) , and with ν ( Y ) < and X i : Z × Ω R + , ( i = 1 , 2 , 3 , 4 , 5 , 6 ) representing the effect of random jumps, which are assumed to be bounded and continuous with respect to ν and to be B ( Y ) × F t -measurable.
With regard to the model (3), in the present work we are particularly interested in answering the following questions:
Q1:
Does the Lévy noise influence the dynamic properties of measles outbreaks?
Q2:
How can contaminated vaccinations contribute to the spread of measles, and what measures are in place to prevent such incidents?
Q3:
What criterion is used to determine the extinction of a disease?
Q4:
What are the criteria that indicate the persistence of the system?

3. The Existence of a Positive Solution and Its Uniqueness

In this section, we intend to apply the techniques from [16] to establish a proof of the existence of a global and positive solution to the proposed stochastic model. To prove the existence and uniqueness of a solution of model (3), it is important to consider the following two conditions:
(C1).
For every Q > 0 L Q > 0 ;
Y A i x 1 , y A i x 2 , y 2 v ( d y ) L Q x 1 x 2 2 , i = 1 , , 6 ,
for y 1 y 2 M ; here,
A 1 ( x , y ) = X 1 ( y ) x at x = S t , A 2 ( x , y ) = X 2 ( y ) x at x = V t , A 3 ( x , y ) = X 3 ( y ) x at x = E t , A 4 ( x , y ) = X 4 ( y ) x at x = I t , A 5 ( x , y ) = X 5 ( y ) x at x = H t , A 6 ( x , y ) = X 6 ( y ) x at x = R t .
(C2).
log X i ( x ) C for X i ( x ) > 1 , i = 1 , , 6 , C is a positive constant.
Theorem 1.
For a given initial value ( S , V , E , I , H , R ) ( 0 ) R + 6 , system (3) has one global root ( S , V , E , I , H , R ) ( t ) R + 6 for all t 0 a.s.
Proof. 
Using ( C 1 ) , the results due to drifted and the diffused terms are local Lipschitzian, for every initial conditions ( S , V , E , I , H , R ) ( 0 ) R + 6 , ∃ a unique localized root ( S ( t ) , V ( t ) , E ( t ) , I ( t ) , H ( t ) , R ( t ) ) on t [ 0 , τ e ) , τ e is the explosion time. To show that the solution is global, we have to prove that τ e = a.s. In the first attempt, we prove the condition that ( S ( t ) , V ( t ) , E ( t ) , I ( t ) , H ( t ) , R ( t ) ) are not approaching infinity in a finite duration of time. Let k 0 > 0 is too much large such that ( S ( 0 ) , V ( 0 ) , E ( 0 ) , I ( 0 ) , H ( 0 ) , R ( 0 ) ) lies in [ 1 k 0 , k 0 ] . For every integer value k k 0 , consider the stopping time in the form of
τ k = inf t 0 , τ e / ( S ( t ) , V ( t ) , E ( t ) , I ( t ) , H ( t ) , R ( t ) ) 1 k , k .
Let inf = ; from this, we can see that τ + τ e , which implies that τ + = a.s showing τ e = + a.s. Considering that τ + is less than , there must exist a number T > 0 such that 0 < P τ + < T .
Next, let us consider the operator F : R + 6 R + from the C 2 set in the form of
F ( S , V , E , I , H , R ) = H + V + S + I + E + R 6 ( log S + log V + log E + log I + log H + log R ) .
Applying I t o ^ formula to F for all t [ 0 , τ + ] , we have
d F = L F ( S , V , E , I , H , R ) d t + ξ 1 ( S 1 ) d W 1 ( t ) + x i 2 ( V 1 ) d W 2 ( t ) + ξ 3 ( E 1 ) d W 3 ( t ) + ξ 4 ( I 1 ) d W 4 ( t ) + ξ 5 ( H 1 ) d W 5 ( t ) + ξ 6 ( R 1 ) d W 6 ( t ) + Υ X 1 ( y ) S l o g X 1 ( y ) + 1 N ˜ ( d χ ) + Υ X 2 ( y ) V l o g 1 + X 2 ( y ) N ˜ ( d χ ) + Υ X 3 ( y ) E log 1 + X 3 ( y ) N ˜ ( d χ ) + Υ X 4 ( y ) I log 1 + X 4 ( y ) N ˜ ( d χ ) , + Υ X 5 ( y ) H log 1 + X 5 ( y ) N ˜ ( d χ ) + Υ X 6 ( y ) R log 1 + X 6 ( y ) N ˜ ( d χ ) .
In Equation (7), L F : R + 6 R + is given by using the assumption C 2 ; thus,
L F ϕ + 6 μ + τ + ω + β + ρ + γ + 2 δ + ξ 1 2 + ξ 2 2 + ξ 3 2 + ξ 4 2 + ξ 5 2 + ξ 6 2 2 + Y X 1 ( y ) log 1 + X 1 ( y ) ν ( d y ) + Y X 2 ( y ) log 1 + X 2 ( y ) ν ( d y ) + Y X 3 ( y ) log 1 + X 3 ( y ) ν ( d y ) + Y X 4 ( y ) log 1 + X 4 ( y ) ν ( d y ) + Y X 5 ( y ) log 1 + X 5 ( y ) ν ( d y ) + Y X 6 ( y ) log 1 + X 6 ( y ) ν ( d y ) .
We can refer to Theorem 2.1 of Fatin et al. [16] for the remaining proof, and therefore we skip it here. □

4. Extinction for System (3)

Minimizing the infection effect on any community depends mostly on time and certain useful conditions taken from analysis of the dynamics of the disease in question. In this section, the investigation of the main condition for the vanishing of the infection is discussed through stochastic modeling. To proceed further, we first define the threshold parameter for the deterministic model (3) as follows:
R 0 D = ( μ + ω ) ϕ β α ( μ + β ) ( μ + δ + ρ ) ( μ + ω + τ ) μ .
Based on the standard techniques for stochastic systems, we are now ready to calculate the threshold number for the stochastic model, which is provided by
R s = α ( μ + β ) ( μ + δ + ρ ) + ξ 3 2 2 + ξ 4 2 2 + y ( X 3 ( y ) + X 4 ( y ) ) ln 1 + ( X 3 ( t ) + X 4 ( y ) ) ν ( d y ) .
Subsequently, we introduce the following concept, which is beneficial for the forthcoming discourse.
X ( t ) = 1 t 0 t X ( w ) d w .
Lemma 1
((Strong Law) [21,22]). If an operator fulfills the condition of continuity of process Z = { Z } 0 t within a local Martingale as t 0 , it vanishes; then,
lim t Z , Z t = , a . s . , lim t Z t Z , Z t = 0 , a . s . lim t sup Z , Z t t < 0 , a . s . , lim t Z t t = 0 , a . s .
Theorem 2.
Consider a solution ( S , V , E , I , H , R ) ( t ) of model (3) with initial condition ( S , V , E , I , H , R ) ( 0 ) R 6 . Next, for q > ( ξ 1 2 ξ 2 2 ξ 3 2 ξ 4 2 ξ 5 2 ξ 6 2 ) 2 and R s < 1 , we have
lim t log E ( t ) t < 0 , a n d lim t log I ( t ) t < 0 , a . s .
The above inequality implies that the compartments E ( t ) and I ( t ) approaching 0 a.s., showing the elimination of the disease with unit probability.
Furthermore,
lim t S ( t ) = ( μ + ω ) ϕ ( μ + ω + τ ) μ , lim t V ( t ) = τ ϕ μ ( τ + ω + μ ) , lim t E ( t ) = 0 , lim t I ( t ) = 0 , lim t H ( t ) = 0 , lim t R ( t ) = 0 .
Proof. 
Let us consider a solution ( S , V , E , I , H , R ) ( t ) of the proposed model (3) associated with the initial condition ( S , V , E , I , H , R ) ( 0 ) in the positive cone of R + 6 . Further, let us define
G 1 ( t ) = ( β + μ ) I ( t ) + β E ( t ) .
If we differentiate relation (13) and then follow the Ito’ formula, we have
d ( l n G 1 ( t ) ) = 1 G 1 α β S I N ( β + μ ) ( ρ + δ + μ ) I β 2 E 2 ξ 3 2 + ( μ + β ) 2 ξ 4 2 I 2 2 ( G 1 ) 2 + β ξ 3 [ β E ( t ) + ( μ + β ) I ] E d W 3 ( t ) + ( μ + β ) ξ 4 [ E + ( μ + β ) I ] I d W 4 ( t ) + y ln 1 + β X 3 ( t ) E + ( μ + β ) X 4 ( y ) I G 1 β X 3 ( y ) E + ( μ + β ) X 4 ( y ) I G 1 ν ( d y ) + y ln 1 + β X 3 ( y ) E ( t ) + ( μ + β ) X 4 I ( t ) G 1 ( t ) N ˜ ( d χ ) 1 G 1 α β I ( β + μ ) ( ρ + δ + μ ) I β 2 E 2 ξ 3 2 2 ( G 1 ) 2 ( μ + β ) 2 ξ 4 2 I 2 2 ( G 1 ) 2 y β X 3 ( y ) E + ( μ + β ) X 4 ( y ) I G 1 ln 1 + β X 3 ( t ) E + ( μ + β ) X 4 ( y ) I G 1 ν ( d y ) + β ξ 3 [ β E ( t ) + ( μ + β ) I ] E d W 3 ( t ) + ( μ + β ) ξ 4 [ E + ( μ + β ) I ] I d W 4 ( t ) + y ln 1 + β X 3 ( y ) E ( t ) + ( μ + β ) X 4 I ( t ) G 1 ( t ) N ˜ ( d χ ) , [ S N ] 1 ( β + μ ) ( μ + ρ + δ ) ( β + μ ) + α ξ 3 2 2 ξ 4 2 2 y ( X 3 ( y ) + X 4 ( y ) ) ln 1 + ( X 3 ( t ) + X 4 ( y ) ) ν ( d y ) + β ξ 3 [ β E ( t ) + ( β + μ ) I ] E d W 3 ( t ) + ( μ + β ) ξ 4 [ E + ( μ + β ) I ] I d W 4 ( t ) + y ln 1 + β X 3 ( y ) E ( t ) + ( μ + β ) X 4 I ( t ) G 1 ( t ) N ˜ ( d χ ) . [ I I + β E ( μ + β ) ]
We obtain the following result if we integrate both sides of the previous inequality over the interval [ 0 , t ] :
l n G 1 ( t ) 1 ( β + μ ) { α [ ( μ + β ) ( ρ + δ + μ ) + ξ 3 2 2 + ξ 4 2 2 + y ( X 3 ( y ) + X 4 ( y ) ) ln 1 + ( X 3 ( t ) + X 4 ( y ) ) ν ( d y ) ] } + 0 t E ξ 3 β d W 3 ( t ) [ ( β + μ + ) I + β E ( t ) ] + 0 t ( β + μ ) I ξ 4 d W 4 ( t ) [ ( β + μ ) I + E ] + 0 t y ln 1 + β X 3 ( y ) E ( t ) + ( μ + β ) X 4 I ( t ) G 1 ( t ) N ˜ ( d χ ) , 1 ( β + μ ) { α [ ( μ + ρ + δ ) ( + β + μ ) + ξ 3 2 2 + ξ 4 2 2 + y ( X 3 ( y ) + X 4 ( y ) ) ln 1 + ( X 3 ( t ) + X 4 ( y ) ) ν ( d y ) ] } + 0 t E ξ 3 β d W 3 ( t ) [ ( β + μ + ) I + β E ( t ) ] + 0 t ( β + μ ) I ξ 4 d W 4 ( t ) [ ( β + μ ) I + E ] + 0 t y ln 1 + β X 3 ( y ) E ( t ) + ( μ + β ) X 4 I ( t ) G 1 ( t ) N ˜ ( d χ ) , ( μ + β ) ( μ + ρ + δ ) + ξ 3 2 2 + ξ 4 2 2 + y ( X 3 ( y ) + X 4 ( y ) ) ln 1 + ( X 3 ( t ) + X 4 ( y ) ) ν ( d y ) ( μ + β ) R s 1 + 0 t E ξ 3 β d W 3 ( t ) [ ( β + μ ) I + β E ( t ) ] + 0 t ( β + μ ) I ξ 4 d W 4 ( t ) [ ( β + μ ) I + E ] + 0 t y ln 1 + β X 3 ( y ) E ( t ) + ( μ + β ) X 4 I ( t ) G 1 ( t ) N ˜ ( d χ ) .
By taking the superior limit as t after dividing Equation (15) by t and using Lemma 1, we obtain
lim sup t ( l n G 1 ( t ) ) ( μ + ρ + δ ) ( β + μ ) + ξ 3 2 2 + ξ 4 2 2 + y ( X 3 ( y ) + X 4 ( y ) ) ln 1 + ( X 3 ( t ) + X 4 ( y ) ) ν ( d y ) ( μ + β ) R s 1 .
If 1 > R s , then lim t G 1 = 0 , a . s whenever R s < 1 . As μ + β and β are both positive, per relation (13) we have lim t [ ( β + μ ) I ( t ) + β E ( t ) ] = 0 lim t E = lim t I = 0 ; thus, we reach the conclusion. □

5. Persistence in Mean

The purpose of this section is to perform an analysis of the persistence of the disease and examine the long-term behaviour of the infection. First, we present the mean persistence, as can be seen in [16].
Definition 1
([19]). Under the following assumption, model (3) shows the persistence of the infection
lim inf t 1 t 0 t F ( r ) d r > 0 a . s .
For more details on disease persistence, interested readers are referred to the results provided in [16,17].
Lemma 2.
Let g C ( R 0 × Ω , R > 0 ) and G C ( R 0 × Ω , R ) such that lim t G ( t ) t = 0 a.s. If we assume the following relation to be true for all positive values of t
log g ( t ) λ 0 t + G ( t ) λ 0 t g ( s ) d s , a . s .
then
lim inf t g ( t ) λ 0 λ a . s . ,
where λ 0 is non-negative and λ is a positive real number.
Next, we provide several mathematical assumptions for the mean persistency of system (3), where the conclusion of the said part is provided by the following theorem.
Theorem 3.
If R 0 s > 1 , then for initial approximations ( S , V , E , I , H , R ) ( 0 ) R + 6 the disease class I ( t ) has the following property:
lim inf t I ( t ) 3 ϕ R 0 s 1 C 1 α , a . s . ,
where C 1 = ϕ ( τ + μ + ξ 1 2 2 + Y X 1 ( y ) + log 1 + X 1 ( y ) ν ( d y ) ) implies that infection is present in the community.
Let us now reproduce the threshold for the stochastic system as
R 0 s = α β a b c ,
where
a = τ + μ + ξ 1 2 2 + Y X 1 ( y ) + log 1 + X 1 ( y ) ν ( d y ) ) , b = β + μ + ξ 3 2 2 + Y X 3 ( y ) + log 1 + X 3 ( y ) ν ( d y ) , c = μ + δ + ρ + ξ 4 2 2 + Y X 4 ( y ) + log 1 + X 4 ( y ) ν ( d y ) .
Proof. 
Let
G 1 = C 1 l n S C 2 l n E C 3 l n I ,
here, C 1 , C 1 and C 3 are constants, and will be calculated later. Using the I t o ^ formula with Equation (21), we obtain
d G 1 = L G 1 C 1 ξ 1 d W 1 ( t ) C 2 ξ 3 d W 3 ( t ) C 3 ξ 4 d W 4 ( t ) Y log 1 + X 1 ( y ) N ˜ ( d χ ) Y log 1 + X 3 ( y ) N ˜ ( d χ ) Y log 1 + X 4 ( y ) N ˜ ( d χ ) ,
where
L G 1 = C 1 ϕ S + C 1 α I N C 1 ω V S + C 1 ( τ + μ ) C 2 α S I E + C 2 ( μ + β ) C 3 β E I + C 3 ( μ + δ + ρ ) + C 1 ξ 1 2 2 + C 2 ξ 3 2 2 + C 3 ξ 4 2 2 + Y C 1 X 1 ( y ) + C 1 log 1 + X 1 ( y ) ν ( d y ) + Y C 2 X 3 ( y ) + C 2 log 1 + X 3 ( y ) ν ( d y ) + Y C 3 X 4 ( y ) + C 2 log 1 + X 4 ( y ) ν ( d y ) , C 1 ϕ S C 2 α S I E C 3 β E I + C 1 ( τ + μ + ξ 1 2 2 + Y X 1 ( y ) + log 1 + X 1 ( y ) ν ( d y ) ) + C 2 ( μ + β + ξ 3 2 2 + Y X 3 ( y ) + log 1 + X 3 ( y ) ν ( d y ) ) + C 3 ( μ + δ + ρ + ξ 4 2 2 + Y X 4 ( y ) + log 1 + X 4 ( y ) ν ( d y ) ) + C 1 α I N .
Let
C 1 ( τ + μ + ξ 1 2 2 + Y X 1 ( y ) + log 1 + X 1 ( y ) ν ( d y ) ) = ϕ , C 2 ( μ + β + ξ 3 2 2 + Y X 3 ( y ) + log 1 + X 3 ( y ) ν ( d y ) ) = ϕ , C 2 ( μ + δ + ρ + ξ 4 2 2 + Y X 4 ( y ) + log 1 + X 4 ( y ) ν ( d y ) ) = ϕ ,
namely,
a = τ + μ + ξ 1 2 2 + Y X 1 ( y ) + log 1 + X 1 ( y ) ν ( d y ) ) , b = μ + β + ξ 3 2 2 + Y X 3 ( y ) + log 1 + X 3 ( y ) ν ( d y ) , c = μ + δ + ρ + ξ 4 2 2 + Y X 4 ( y ) + log 1 + X 4 ( y ) ν ( d y ) .
Then, we can write inequality (23) in the form of
L G 1 3 C 1 ϕ S × C 3 α S I E × C 3 β E I + C 1 ( τ + μ + ξ 1 2 2 + Y X 1 ( y ) + log 1 + X 1 ( y ) ν ( d y ) ) + C 2 ( μ + β + ξ 3 2 2 + Y X 3 ( y ) + log 1 + X 3 ( y ) ν ( d y ) ) + C 3 ( μ + δ + ρ + ξ 4 2 2 + Y X 4 ( y ) + log 1 + X 4 ( y ) ν ( d y ) ) + C 1 α I , = 3 ϕ 3 α β a b c + 3 ϕ + C 1 α I , = 3 ϕ R 0 s 1 + C 1 α I .
By putting Equation (25) into Equation (21) and taking the integral of both sides of the stochastic model (3), we have
G 1 ( S , E , I ) ( t ) G 1 ( S , E , I ) ( 0 ) t 3 ϕ R 0 s 1 + C 1 α I C 1 ξ 1 d W 1 ( t ) C 2 ξ 3 d W 3 C 3 ξ 4 d W 4 Y log 1 + X 1 ( y ) N ˜ ( d χ ) Y log 1 + X 3 ( y ) N ˜ ( d χ ) Y log 1 + X 4 ( y ) N ˜ ( d χ ) .
If we define the notion Ψ ( t ) in the form of
Ψ ( t ) = C 1 ξ 1 d W 1 ( t ) C 2 ξ 3 d W 3 C 3 ξ 4 d W 4 Y log 1 + X 1 ( y ) N ˜ ( d χ ) Y log 1 + X 3 ( y ) N ˜ ( d χ ) Y log 1 + X 4 ( y ) N ˜ ( d χ ) .
Then, from Strong’s law as provided in Lemma 1, we have
lim t Ψ ( t ) = 0 .
Further, from Equation (27), we have
C 1 α I ( t ) 3 ϕ R 0 s 1 Ψ ( t ) + G 1 ( S ( t ) , E ( t ) , I ( t ) ) G 1 ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) t , I ( t ) 3 ϕ R 0 s 1 C 1 α Ψ ( t ) C 1 α 1 C 1 α G 1 ( S ( t ) , E ( 0 ) , I ( t ) ) G 1 ( S ( 0 ) , E ( 0 ) , I ( 0 ) ) t .
From Lemma 1 and Equation (28), the superior limit of Equation (5) takes the form
lim inf t I ( t ) 3 ϕ R 0 s 1 C 1 α 0 , a . s . ,
showing that lim inf t I ( t ) 0 . □
Thus, the proof of Theorem 3 is concluded.

6. Estimation

Utilizing practical observations to obtain insights into certain missing epidemiological factors is a widely used technique in biological systems analysis. The verification of analytical results pertaining to the measles model (1) and determination of the parameters were performed by considering the measles data presented in Table 1; accordingly, the model was fitted against the data. From the WHO reports for 2018, the continuous rate of fatality μ for a Pakistani individual is 66.5 years ( 1.253 × 10 4 per month), and with a population of 207,862,518 , the inflow rate is estimated to be ϕ 25,983 individuals per month. In [12], it is reported that the Measles vaccine has an efficacy rate of approximately 97 % , indicating that the vaccination outcome, denoted by τ , is nearly 0.97 . The effects of the other parameters, such as the interactive rate β , the retrieval rate δ , the rate of clinically tested symptoms α , and the rate of vaccine coverage ω , are presented in Table 2 in relation to the calculated numbers. These parameters were estimated, and the modeled predictions using real data are plotted in Figure 2 using the MATLAB software considering the data for the first ten months of 2019. It can be seen from Figure 2 that model (1) provides a good fit, and the actual measles data are almost covered by the curve predicted by model (1). Using the result 1 10 k = 1 10 ϰ k real ϰ k approximate ϰ k real 1.4685 e 01 , we can find the mean relative error of the fitting procedure.

7. Numerical Results and Discussion

In this section, we present a graphical representation of the model dynamics using the available numerical method(s) and values of the parameters discussed in the previous section. For illustrative purposes, we provide the findings of our scheme and compute the approximate paths for the stochastic model (3) and its deterministic counterpart. The desired time interval is [ 0 , 100 ] and the step size is Δ = 0.3 , whereas the initial size of the population is provided by ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) = ( 0.5 , 0.4 , 0.3 , 0.5 , 0.2 , 0.1 ) .
To simulate the model, we utilized two different types of hypothetical data. Initially, we simulated the model by considering the intensity of the noise and the parameter values from Example 1. The figure presented in Figure 3 demonstrates that as R 0 D < 1 , the solution of system (3) with L’evy jump approaches the disease-free equilibrium point of the corresponding deterministic system, indicating that the disease is in the process of disappearing. Theorem 2 provides the necessary conditions for the extinction of system (3). The graphical representation verifies that the result of Theorem 2 is valid only if R s < 1 .
Example 1.
In the present scenario, we have considered the parameter values as ϕ = 1.0 , τ = 0.3 , α = 0.3 , μ = 0.03 , β = 0.005 , δ = 0.4 , γ = 0.02 , ρ = 0.01 , ω = 0.02 , and the intensities of the white noise as ξ 1 = 0.55 , ξ 2 = 0.40 , ξ 3 = 0.30 , ξ 4 = 0.20 , ξ 5 = 0.10 , ξ 6 = 0.20 , and X i ( y ) = k i y 2 1 + y 2 , with y = 0.5 and k i equal to 0.50 , 0.50 , 0.40 , 0.20 , 0.20 , and 0.30 , respectively, for i = 1 , , 6 . The simulations of the model indicate that the epidemic will become extinct, as predicted by Theorem 2. The disappearance of the infection based on the predictions of the stochastic system are depicted in Figure 3a–f.
The findings of Theorem 3 suggest that the disease will continue to exist in the community with the dynamics predicted by model (3) if the value of R 0 D is greater than 1. The statement R 0 D > 1 indicates that R 0 s will be greater than 1 even with low intensities of the noise. This conclusion is supported by the numerical simulations presented in Figure 3c,d and Figure 4c,d, which validate the expected behavior of the epidemic model (3) as per Theorem 3.
Example 2.
In this example, the numerical values of the parameters are ϕ = 0.9 , τ = 0.75 , α = 0.5 , μ = 0.06 , β = 0.05 , δ = 0.4 , γ = 0.2 , ρ = 0.1 , and ω = 0.2 , and the intensity of the white noise is ξ 1 = 0.45 , ξ 2 = 0.30 , ξ 3 = 0.10 , ξ 4 = 0.20 , ξ 5 = 0.15 , ξ 6 = 0.20 , and X i ( y ) = k i y 2 1 + y 2 , with y = 0.5 and where k i equal to 0.50 , 0.35 , 0.30 , 0.35 , 0.10 , and 0.20 , respectively, for i = 1 , , 6 .
For this set of parameters, it is easy to verify that both the stochastic threshold R 0 s and the basic reproduction number R 0 D are greater than one. The simulation results of the mean shown in Figure 4 confirm the persistence of the disease, which supports the conclusion of Theorem 3. Additionally, Figure 4a–f illustrates that the epidemic will continue to persist in the population in the meantime.

The Impact of Lévy Noise on the I Class

The impact of the intensity of the white noise on class I corresponding to system (3) is shown in Figure 5a–c. These figures suggest that increasing values of ξ i for i = 1 , , 6 lead towards the extinction of the disease. This means that the size of the infected class approaches zero as the intensity value of the noise increases. Further, this indicates that for small values of the noise intensity the infected class oscillates around the endemic steady state I m which confirms the result of Theorem 3. However, if the intensity of the white noise term is sufficiently high, the solution I may not exhibit oscillations near the endemic equilibrium point. This demonstrates that continuous efforts to increase stochastic disruptions through mass recovery of susceptible individuals and the effective treatment and care of the infected persons can significantly lower the spread of the measles virus in the population.

8. Concluding Remarks and Future Directions

In this work, we have considered a stochastic S V E I H R epidemic model with Lévy noise while considering the nonlinear incidence rate. Our study aimed to investigate the impact of environmental noise on human society, and the results can shed light on the crucial role of noise in disease persistence and extinction. The proposed model assumed six different classes: a healthy class ( S ) , vaccinated class ( V ) , infected class ( I ) , exposed class ( E ) , hospitalized class ( H ) , and class of recovered individuals ( R ) . After model formulation, the study utilized the uniqueness and existence techniques to derive a positive solution for the underlying stochastic model. The stability results of the model were investigated using the Lyapunov function. The model’s dynamic characteristics were analyzed in the context of the infection-free and endemic states of the associated ODE model. The stochastic threshold R s was calculated to ensure disease extinction whenever R s < 1 . Using data from Pakistan in 2019, we estimated the model parameters and conducted simulations to forecast future disease behavior. The results were compared to actual data using standard curve fitting tools. In our future research, we plan to investigate the impact of regime switching and temporary immunity on system (3).

Author Contributions

Methodology, Y.S.; Software, Y.S.; Formal analysis, P.L.; Investigation, P.L.; Writing—original draft, Y.S.; Writing—review and editing, P.L.; Project administration, Y.S.; Funding acquisition, P.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the project of “Research on Key parameters of structure design of petal vortex shaft torsional section” (No. JG523001), Development of Optimization Algorithm for 3D Data of Light Field (No. 3630019001), and the National Natural Science Foundation of China (No. 11901114), and Guangzhou Science and Technology Innovation general project (No. 201904010010), young innovative talents project of the Guangdong Provincial Department of Education (No. 2018KQNCX087).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Mose, O.F.; Sigey, J.K.; Okello, J.A.; Okwoyo, J.M.; Giterere, J. Mathematical Modeling on the Control of Measles by Vaccination: Case Study of Kisii County, Kenya; Kenya Education Research Database: Kenya, 2014; Available online: https://kerd.ku.ac.ke/ (accessed on 16 April 2023).
  2. Garba, S.M.; Safi, M.A.; Usaini, S. Mathematical model for assessing the impact of vaccination and treatment on measles transmission dynamics. Math. Methods Appl. Sci. 2017, 40, 6371–6388. [Google Scholar] [CrossRef]
  3. Roberts, M.G.; Tobias, M.I. Predicting and preventing measles epidemic in New Zealand: Application of mathematical model. Epidem. Infect. 2000, 124, 279–287. [Google Scholar] [CrossRef] [PubMed]
  4. Qureshi, S.; Jan, R. Modeling of measles epidemic with optimized fractional order under Caputo differential operator. Chaos Solitons Fractals 2021, 145, 110766. [Google Scholar] [CrossRef]
  5. Perry, R.T.; Halsey, N.A. The clinical significance of measles: A review. J. Infect. Dis. 2004, 189 (Suppl. S1), S4–S16. [Google Scholar] [PubMed]
  6. Guo, B.; Khan, A.; Din, A. Numerical Simulation of Nonlinear Stochastic Analysis for Measles Transmission: A Case Study of a Measles Epidemic in Pakistan. Fractal Fract. 2023, 7, 130. [Google Scholar] [CrossRef]
  7. Mossong, J.; Muller, C.P. Modelling measles re-emergence as a result of waning of immunity in vaccinated populations. Vaccine 2003, 21, 4597–4603. [Google Scholar] [CrossRef] [PubMed]
  8. Bolarin, G. On the dynamical analysis of a new model for measles infection. Int. J. Math. Trends Technol. 2014, 2, 144–155. [Google Scholar] [CrossRef]
  9. Taiwo, L.; Idris, S.; Abubakar, A.; Nguku, P.; Nsubuga, P.; Gidado, S. Factors affecting access to information on routine immunization among mothers of under 5 children in Kaduna state Nigeria, 2015. Pan Afr. Med. J. 2017, 27, 1–8. [Google Scholar] [CrossRef] [PubMed]
  10. Mere, M.O.; Goodson, J.L.; Chandio, A.K.; Rana, M.S.; Hasan, Q.; Teleb, N.; Alexander, J.P., Jr. Progress toward measles elimination—Pakistan, 2000–2018. Morb. Mortal. Wkly. Rep. 2019, 68, 505. [Google Scholar] [CrossRef] [PubMed]
  11. World Health Organization. Eastern Mediterranean Vaccine Action Plan 2016–2020: A Framework for Implementation of the Global Vaccine Action Plan (No. WHO-EM/EPI/353/E); World Health Organization, Regional Office for the Eastern Mediterranean: Geneva, Switzerland, 2019.
  12. Memon, Z.; Qureshi, S.; Memon, B.R. Mathematical analysis for a new nonlinear measles epidemiological system using real incidence data from Pakistan. Eur. Phys. J. Plus 2020, 135, 378. [Google Scholar] [CrossRef] [PubMed]
  13. Din, A.; Ain, Q.T. Stochastic optimal control analysis of a mathematical model: Theory and application to non-singular kernels. Fractal Fract. 2022, 6, 279. [Google Scholar] [CrossRef]
  14. Khan, F.M.; Khan, Z.U.; Lv, Y.P.; Yusuf, A.; Din, A. Investigating of fractional order dengue epidemic model with ABC operator. Results Phys. 2021, 24, 104075. [Google Scholar] [CrossRef]
  15. Liu, P.; Ikram, R.; Khan, A.; Din, A. The measles epidemic model assessment under real statistics: An application of stochastic optimal control theory. Comput. Methods Biomech. Biomed. Eng. 2022, 26, 138–159. [Google Scholar] [CrossRef] [PubMed]
  16. El Fatini, M.; Sekkak, I. Lévy noise impact on a stochastic delayed epidemic model with Crowly–Martin incidence and crowding effect. Phys. Stat. Mech. Its Appl. 2020, 541, 123315. [Google Scholar] [CrossRef]
  17. Sabbar, Y.; Khan, A.; Din, A.; Tilioua, M. New Method to Investigate the Impact of Independent Quadratic Stable Poisson Jumps on the Dynamics of a Disease under Vaccination Strategy. Fractal Fract. 2023, 7, 226. [Google Scholar] [CrossRef]
  18. Berrhazi, B.E.; El Fatini, M.; Caraballo Garrido, T.; Pettersson, R. A stochastic SIRI epidemic model with Lévy noise. Discret. Contin. Dyn.-Syst.-Ser. B 2018, 23, 3645–3661. [Google Scholar] [CrossRef]
  19. Li, X.P.; Din, A.; Zeb, A.; Kumar, S.; Saeed, T. The impact of Lévy noise on a stochastic and fractal-fractional Atangana–Baleanu order hepatitis B model under real statistical data. Chaos Solitons Fractals 2022, 154, 111623. [Google Scholar] [CrossRef]
  20. James Peter, O.; Ojo, M.M.; Viriyapong, R.; Abiodun Oguntolu, F. Mathematical model of measles transmission dynamics using real data from Nigeria. J. Differ. Equ. Appl. 2022, 28, 753–770. [Google Scholar] [CrossRef]
  21. Zhao, Y.; Jiang, D. The threshold of a stochastic SIRS epidemic model with saturated incidence. Appl. Math. Lett. 2014, 34, 90–93. [Google Scholar] [CrossRef]
  22. Zhao, Y.; Jiang, D. The threshold of a stochastic SIS epidemic model with vaccination. Appl. Math. Comput. 2014, 243, 718–727. [Google Scholar] [CrossRef]
Figure 1. Moments of individuals among the classes in the measles model (1) [21].
Figure 1. Moments of individuals among the classes in the measles model (1) [21].
Fractalfract 07 00434 g001
Figure 2. Fitting of the model (1) against reported measles cases from Table 1.
Figure 2. Fitting of the model (1) against reported measles cases from Table 1.
Fractalfract 07 00434 g002
Figure 3. The plot shows the dynamics of the stochastic system (3) and the associated deterministic system (1) subject to jumps and without jumps. (a) The plot shows the S ( t ) —curve with and without jumps. (b) The plot shows the dynamics of V ( t ) with and without jumps. (c) Three different scenarios for class I ( t ) . (d) The dynamics of the exposed class with and without jumps. (e) Time evolution of the number of hospitalized people with and without jumps. (f) Three different scenarios for class R ( t ) .
Figure 3. The plot shows the dynamics of the stochastic system (3) and the associated deterministic system (1) subject to jumps and without jumps. (a) The plot shows the S ( t ) —curve with and without jumps. (b) The plot shows the dynamics of V ( t ) with and without jumps. (c) Three different scenarios for class I ( t ) . (d) The dynamics of the exposed class with and without jumps. (e) Time evolution of the number of hospitalized people with and without jumps. (f) Three different scenarios for class R ( t ) .
Fractalfract 07 00434 g003aFractalfract 07 00434 g003b
Figure 4. Simulation of the stochastic system (3) and deterministic system (1) with and without jumps for the case when the respective thresholds are greater than one. (a) The dynamics of the susceptible class with and without jumps. (b) Time evolution of the number of vaccinated people both with and without jumps. (c) The plot shows the E ( t ) —curve with and without jumps. (d) The dynamics of the infected class with and without jumps. (e) The plot shows the H ( t ) —curve with and without jumps. (f) Time evolution of the number of recovered people when the threshold exceeds one.
Figure 4. Simulation of the stochastic system (3) and deterministic system (1) with and without jumps for the case when the respective thresholds are greater than one. (a) The dynamics of the susceptible class with and without jumps. (b) Time evolution of the number of vaccinated people both with and without jumps. (c) The plot shows the E ( t ) —curve with and without jumps. (d) The dynamics of the infected class with and without jumps. (e) The plot shows the H ( t ) —curve with and without jumps. (f) Time evolution of the number of recovered people when the threshold exceeds one.
Fractalfract 07 00434 g004aFractalfract 07 00434 g004b
Figure 5. Simulations of I ( t ) based on the stochastic and deterministic systems, showing the effect of the intensities on class I when ϕ = 0.5 , τ = 0.90 , α = 0.5 , μ = 0.06 , β = 0.05 , δ = 0.4 , γ = 0.2 , ρ = 0.2 , ω = 0.3 , and the noise intensity is assumed as ξ 1 = 0.50 , ξ 2 = 0.40 , ξ 3 = 0.10 , ξ 4 = 0.20 , ξ 5 = 0.15 , ξ 6 = 0.20 , and X i ( y ) = k i y 2 1 + y 2 , with y = 0.5 and where k i equals 0.60 , 0.34 , 0.30 , 0.35 , 0.10 , 0.20 , respectively, for i = 1 , , 6 and ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) = ( 0.5 , 0.4 , 0.5 , 0.3 , 0.2 , 0.1 ) . (a) The dynamics of I ( t ) subject to sufficiently large values of noise intensity. (b) The plot shows the behavior of the infected class subject to moderate values of noise intensity. (c) The dynamic behavior of the infected class subject to low values of noise intensity.
Figure 5. Simulations of I ( t ) based on the stochastic and deterministic systems, showing the effect of the intensities on class I when ϕ = 0.5 , τ = 0.90 , α = 0.5 , μ = 0.06 , β = 0.05 , δ = 0.4 , γ = 0.2 , ρ = 0.2 , ω = 0.3 , and the noise intensity is assumed as ξ 1 = 0.50 , ξ 2 = 0.40 , ξ 3 = 0.10 , ξ 4 = 0.20 , ξ 5 = 0.15 , ξ 6 = 0.20 , and X i ( y ) = k i y 2 1 + y 2 , with y = 0.5 and where k i equals 0.60 , 0.34 , 0.30 , 0.35 , 0.10 , 0.20 , respectively, for i = 1 , , 6 and ( S 0 , V 0 , E 0 , I 0 , H 0 , R 0 ) = ( 0.5 , 0.4 , 0.5 , 0.3 , 0.2 , 0.1 ) . (a) The dynamics of I ( t ) subject to sufficiently large values of noise intensity. (b) The plot shows the behavior of the infected class subject to moderate values of noise intensity. (c) The dynamic behavior of the infected class subject to low values of noise intensity.
Fractalfract 07 00434 g005
Table 1. Reported measles cases in Pakistan during the first ten months of 2019 [12].
Table 1. Reported measles cases in Pakistan during the first ten months of 2019 [12].
JanFebMarAprMayJuneJulyAugSepOct
23825339839827716971292418
Table 2. Justification and values of the parameters used for simulating the model (1).
Table 2. Justification and values of the parameters used for simulating the model (1).
ParameterDescriptionSource
ϕ 260,479Estimated
α 1.253133 × 10 3 Estimated
ω 0.97Estimated
τ 1.60056 × 10 7 Fitted
μ 9.3408Fitted
β 9.2373 × 10 1 Fitted
δ 5.8306 × 10 1 Fitted
ρ 0Estimated
γ 5.8306 × 10 1 Fitted
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Song, Y.; Liu, P. Impact of Lévy Noise with Infinite Activity on the Dynamics of Measles Epidemics. Fractal Fract. 2023, 7, 434. https://doi.org/10.3390/fractalfract7060434

AMA Style

Song Y, Liu P. Impact of Lévy Noise with Infinite Activity on the Dynamics of Measles Epidemics. Fractal and Fractional. 2023; 7(6):434. https://doi.org/10.3390/fractalfract7060434

Chicago/Turabian Style

Song, Yuqin, and Peijiang Liu. 2023. "Impact of Lévy Noise with Infinite Activity on the Dynamics of Measles Epidemics" Fractal and Fractional 7, no. 6: 434. https://doi.org/10.3390/fractalfract7060434

APA Style

Song, Y., & Liu, P. (2023). Impact of Lévy Noise with Infinite Activity on the Dynamics of Measles Epidemics. Fractal and Fractional, 7(6), 434. https://doi.org/10.3390/fractalfract7060434

Article Metrics

Back to TopTop