Research Article | Open Access
Analysis of a Stochastic SIR Model with Vaccination and Nonlinear Incidence Rate
We expand an SIR epidemic model with vertical and nonlinear incidence rates from a deterministic frame to a stochastic one. The existence of a positive global analytical solution of the proposed stochastic model is shown, and conditions for the extinction and persistence of the disease are established. The presented results are demonstrated by numerical simulations.
Mathematical models play an important role in the analysis and control of infectious diseases, and thus effective measures can be taken to reduce its transmission as much as possible. The study of mathematical models in epidemiology has received much attention from many scientists, and some novel results are obtained [1–4]. Recently, stochastic analysis has widely been applied in mathematical modeling in biology [5–9].
For compartmental mathematical models, the total population is divided into three classes, namely, susceptible population S(t), infected population I(t), and recover population R(t). For more details, see . For some diseases, such as AIDS, rubella, varicella, hepatitis B, hepatitis, syphilis, and mumps, it is one of the main transmission modes that infected mothers infect their unborn or newborn offsprings, called vertical transmission . Meng and Chen  proposed an SIR epidemic model with vaccination and vertical transmission mode as follows:where b is the mortality rate in the susceptible and the recovered individuals and d is the mortality rate in the infective individuals. The constants p and q () are vertical transmission rates, namely, p and q are, respectively, the proportion of the offspring of infective parents that are susceptible individuals and the rest are born infected. The arrival of newborns constitutes a recruitment rate of into the susceptible individuals and into the infectious individuals. is the successful vaccination proportion to the newborn from S and R, r is the recover rate in the infective individuals into recovered individuals, and β is the contact rate. System (1) has a basic reproduction number defined by
On the contrary, environmental fluctuations have great influence on all aspects of real life. The aim of this work is to study the effect of these environmental fluctuations on the transmission rate β. We assume that the stochastic perturbations are of white noise type, that is, , where is a Brownian motion and σ is the intensity. Then, the stochastic version corresponding to the deterministic model (1) with general incidence rate is as follows:where is the independent standard Brownian motions defined on a complete probability space with a filtration satisfying the usual conditions (i.e., it is increasing and right continuous while contains all null sets).
is the incidence rate, where . In addition, the incidence rate is the number of new infected situations by population in a determined time period. To model the disease transmission process, several authors employ the following bilinear incidence rate , where β is a positive constant . Yet, there exist many forms of nonlinear incidence rate and every form presents some advantages [14–16]. It is very important to note that is a general form which represents mutual interference between S and I. in particular cases:(1)If , becomes a bilinear incidence rate(2)If or , becomes a saturated incidence rate (3)If , becomes a Beddington–DeAngelis functional response (4)If , becomes the Crowley–Martin functional response presented in 
Next, we consider the d-dimensional stochastic system:where is a function in defined in and is a matrix, and f and are locally Lipschitz functions in x. is a d-dimensional standard Wiener process defined on the above probability space.
If the differential operator acts on a function , thenwhere , , and .
The organization of this paper is as follows. In Section 2, we show the existence of unique positive global solution to the given SDE system. Extinction and persistence in mean results are explored in Section 3 and Section 4, respectively. In Section 5, the analytical results are illustrated with the support of numerical examples. Finally, we close the paper with conclusion and future directions.
2. Existence and Uniqueness of the Nonnegative Solution
As we are dealing with the population model, the positive solution of the model is of our interest. The coefficients of (3) are locally Lipschitz continuous and do not satisfy the linear growth condition, so the solution of (3) may explode at a finite time. The following theorem shows that the solution is positive and will not explode at a finite time.
We define a subset of as follows:
Theorem 1. For any given initial value , there is a unique positive solution of (3) on , and the solution will remain in with probability 1, namely, for all almost surely (briefly a.s.).
Proof 1. Since the coefficients of system (3) are locally Lipschitz continuous, for any initial value , there is a unique local solution on , where is the explosion time. To show this solution is global, we need to show that For this, we define the stopping time τ bywith the traditional setting (as usual denotes the empty set). We have . If we can show that a.s., then a.s. and for all . In other words, to complete the proof, all we only need to show is that . Assume that this statement is false, then there exists a such that . Define a function U, by the expressionUsing Itô’s formula, we getSince , and , , we get,where . Then, we haveNote that some components of equal 0; thus, .
Letting in (12), we obtainwhich contradicts our assumption. Then, a.s. This completes the proof of theorem.
In the following, we give a condition for the extinction of the disease. Let
Theorem 2. Let be the solution of system (3) with initial value . Assume that(i) or(ii)Then,
Namely, tends to zero exponentially a.s., i.e., the disease dies out with probability 1.
Proof 2. Applying Itô’s formula to system (3) leads toIntegrating both sides of (16) from 0 to , we getwhere , which is a local continuous martingale, and . Moreover, its quadratic variation isBy the large number theorem for martingales , we obtainIf condition (i) is satisfied, dividing by t and taking the limit superior of both sides of (17), we getIf the condition (ii) is satisfied, note thatIt follows thatTaking the superior limit of both sides of (22), we obtainwhich implies that .
Remark 1. From Theorem 1, we can get that the disease will die out if , and the white noise is not large such that , while if the white noise is large enough such that condition (i) is satisfied, then the infectious disease of system (3) goes to extinction almost surely.
Here, we investigate the condition for the persistence of the disease. The basic reproduction number (see ), is the threshold between disease extinction and persistence, with extinction for and persistence for in the deterministic model. In the stochastic model, we define the threshold of persistence for disease as
Definition 1. System (3) is said to be persistent in the mean ifWe define
Theorem 3. If then the solution of system (3) with initial initial value is persistent in mean. Moreover, we havewhere
Proof 3. We haveIntegrating from 0 to t and dividing by the third equation of system (3), we getFrom (28), one can getwhere From Itô’s formula and (3), we obtainIntegrating from 0 to t and dividing by on both sides of yieldsWe can rewrite the inequality (32) aswhere We can see that and Thus, one has and Taking the inferior limit of both sides of (33) yieldsThis completes the proof of theorem.
5. Numerical Simulations
In order to illustrate our theoretical results, we give some numerical simulations. The values of m, , and β will be varied over the different examples.
Example 1. We choose the parameters in system (3) as follows:By calculation, we have ; in this case, the disease dies out as shown in Figure 1(a). By choosing , we obtain , and we deduce that the disease persists in the population. Figure 1(b) illustrates this result.
Example 2. We choose the parameters in system (3) as follows:A simple computation shows thatIt follows that the condition of Theorem 2 is satisfied. We conclude that the disease dies out; Figure 2 illustrates this result.
Example 3. On the contrary, we choose and σ = 0.1 by simple calculation, and it can be found that which implies that the disease persists (Figure 3).
This article discusses a stochastic SIR epidemic model with vertical transmission and vaccination and nonlinear incidence rate. We have shown that when the noise is so small such that , the extinction of the disease can be determined by the value of , i.e., if , the disease dies out. Moreover, the disease dies out when the white noise is large enough such that .
The persistence of the disease is determined by , i.e., if , the disease persists. We presented some numerical simulations to illustrate the obtained analytical results. To go further in this study, we can give a new dimension to the stochastic SIR epidemic model (3) by introducing a different type of noise which is the random telegraph noise (RTN), modeled by the Markov chain. We will investigate this case in our future works.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
- J. Zhang and Z. Ma, “Global dynamics of an SEIR epidemic model with saturating contact rate,” Mathematical Biosciences, vol. 185, no. 1, pp. 15–32, 2003.
- Z. Feng, W. Huang, and C. Castillo-Chavez, “Global behavior of a multi-group SIS epidemic model with age structure,” Journal of Differential Equations, vol. 218, no. 2, pp. 292–324, 2005.
- J. M. Tchuenche, A. Nwagwo, and R. Levins, “Global behaviour of an SIR epidemic model with time delay,” Mathematical Methods in the Applied Sciences, vol. 30, no. 6, pp. 733–749, 2007.
- K. Hattaf, N. Yousfi, and A. Tridane, “Mathematical analysis of a virus dynamics model with general incidence rate and cure rate,” Nonlinear Analysis: Real World Applications, vol. 13, no. 4, pp. 1866–1872, 2012.
- Y. Cai, Y. Kang, and W. Wang, “A stochastic SIRS epidemic model with nonlinear incidence rate,” Applied Mathematics and Computation, vol. 305, pp. 221–240, 2017.
- D. Jiang, J. Yu, C. Ji, and N. Shi, “Asymptotic behavior of global positive solution to a stochastic SIR model,” Mathematical and Computer Modelling, vol. 54, no. 1-2, pp. 221–232, 2011.
- J. Adnani, K. Hattaf, and N. Yousfi, “Stability analysis of a stochastic SIR epidemic model with specific nonlinear incidence rate,” International Journal of Stochastic Analysis, vol. 2013, Article ID 431257, 4 pages, 2013.
- A. El Koufi, J. Adnani, A. Bennar, and N. Yousfi, “Dynamical behaviors of a stochastic SIRS epidemic model,” Journal of Mathematical and Computational Science, vol. 8, no. 3, pp. 421–436, 2018.
- A. El Koufi, M. Naim, A. Bennar, and N. Yousfi, “Stability analysis of a stochastic SIS model with double epidemic hypothesis and specific nonlinear incidence rate,” Communications in Mathematical Biology and Neuroscience, vol. 2018, 2018.
- M. Zhien, Y. Zhou, and J. Wu, Modeling and Dynamics of Infectious Diseases, vol. 11, World Scientific, Singapore, 2009.
- S. Busenberg and K. Cooke, Vertically Transmitted Diseases: Models and Dynamics, Springer Science & Business Media, Berlin, Germany, 2012.
- X. Meng and L. Chen, “The dynamics of a new SIR epidemic model concerning pulse vaccination strategy,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 582–597, 2008.
- M. G. M. Gomes, L. J. White, and G. F. Medley, “The reinfection threshold,” Journal of Theoretical Biology, vol. 236, no. 1, pp. 111–113, 2005.
- A. Kaddar, “Stability analysis in a delayed SIR epidemic model with a saturated incidence rate,” Nonlinear Analysis: Modelling and Control, vol. 15, no. 3, pp. 299–306, 2010.
- J. Cui, Y. Sun, and H. Zhu, “The impact of media on the control of infectious diseases,” Journal of Dynamics and Differential Equations, vol. 20, no. 1, pp. 31–53, 2008.
- P. Das, D. Mukherjee, and A. K. Sarkar, “Study of an SI epidemic model with nonlinear incidence rate: discrete and stochastic version,” Applied Mathematics and Computation, vol. 218, no. 6, pp. 2509–2515, 2011.
- X. Liu and L. Yang, “Stability analysis of an SEIQV epidemic model with saturated incidence rate,” Nonlinear Analysis: Real World Applications, vol. 13, no. 6, pp. 2671–2679, 2012.
- J. R. Beddington, “Mutual interference between parasites or predators and its effect on searching efficiency,” The Journal of Animal Ecology, vol. 44, no. 1, pp. 331–340, 1975.
- P. H. Crowley and E. K. Martin, “Functional responses and interference within and between year classes of a dragonfly population,” Journal of the North American Benthological Society, vol. 8, no. 3, p. 211, 1989.
- P. V. D. Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, no. 1-2, pp. 29–48, 2002.
- X. Mao, Stochastic Differential Equations and Applications, Elsevier, Amsterdam, Netherlands, 2007.
Copyright © 2019 Amine El Koufi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.