International Journal of Differential Equations
Volume 2019 (2019), Article ID 9275051, 9 pages
https://doi.org/10.1155/2019/9275051

Research Article

Analysis of a Stochastic SIR Model with Vaccination and Nonlinear Incidence Rate

Amine El Koufi, Jihad Adnani, Abdelkrim Bennar, and Noura Yousfi

Laboratory of Analysis Modeling and Simulation (LAMS), Faculty of Sciences Ben M’sik, Hassan II University, P.O. Box 7955, Sidi Othman, Casablanca, Morocco

Correspondence should be addressed to Amine El Koufi; elkoufiamine1@gmail.com

Received 1 March 2019; Accepted 1 August 2019; Published 21 August 2019

Academic Editor: Elena Kaikina

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.

Abstract.

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.

1. Introduction

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 [10]. 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 [11]. Meng and Chen [12] 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 [13]. 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 [17](3)If , becomes a Beddington–DeAngelis functional response [18](4)If , becomes the Crowley–Martin functional response presented in [19]

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.

Denote by the family of all nonnegative functions defined on    such that they are continuously twice differentiable in x and once in t. The differential operator [16] associated with (4) is defined by

If the differential operator acts on a function , then

where , , 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 τ by

with 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 expression

Using Itô’s formula, we get

Since , and , , we get,

where . Then, we have

Note that some components of equal 0; thus, .
Letting in (12), we obtain

which contradicts our assumption. Then, a.s. This completes the proof of theorem.

3. Extinction

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 to

Integrating both sides of (16) from 0 to , we get

where , which is a local continuous martingale, and . Moreover, its quadratic variation is

By the large number theorem for martingales [21], we obtain

If condition (i) is satisfied, dividing by t and taking the limit superior of both sides of (17), we get

If the condition (ii) is satisfied, note that

It follows that

Taking the superior limit of both sides of (22), we obtain

which 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.

4. Persistence

Here, we investigate the condition for the persistence of the disease. The basic reproduction number (see [20]), 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 if

We define

Theorem 3. If then the solution of system (3) with initial initial value is persistent in mean. Moreover, we have

where

Proof 3. We have

Integrating from 0 to t and dividing by the third equation of system (3), we get

From (28), one can get

where From Itô’s formula and (3), we obtain

Integrating from 0 to t and dividing by on both sides of yields

We can rewrite the inequality (32) as

where We can see that and Thus, one has and Taking the inferior limit of both sides of (33) yields

This 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.

(a)

(b)

Figure 1: Simulations of the path S(t), I(t), and R(t) for the corresponding deterministic system (1).

Example 2. We choose the parameters in system (3) as follows:

A simple computation shows that

It follows that the condition of Theorem 2 is satisfied. We conclude that the disease dies out; Figure 2 illustrates this result.

(a)

(b)

(c)

Figure 2: Simulations of the path (a) S(t), (b) I(t), and (c) R(t) for the deterministic system (1) and the corresponding stochastic system (3).

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).

(a)

(b)

(c)

Figure 3: Simulations of the path (a) S(t), (b) I(t), and (c) R(t) for the deterministic system (1) and the corresponding stochastic system (3).

6. Conclusion

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.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

References

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.