Abstract and Applied Analysis

Volume 2014 (2014), Article ID 131257, 7 pages

http://dx.doi.org/10.1155/2014/131257

Research Article

An SIR Epidemic Model with Time Delay and General Nonlinear Incidence Rate

Mingming Li and Xianning Liu

Key Laboratory of Eco-Environments in Three Gorges Reservoir Region (Ministry of Education), School of Mathematics and Statistics, Southwest University, Chongqing 400715, China

Received 26 December 2013; Accepted 6 January 2014; Published 20 February 2014

Academic Editor: Weiming Wang

Copyright © 2014 Mingming Li and Xianning Liu. 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.

An SIR epidemic model with nonlinear incidence rate and time delay is investigated. The disease transmission function and the rate that infected individuals recovered from the infected compartment are assumed to be governed by general functions and , respectively. By constructing Lyapunov functionals and using the Lyapunov-LaSalle invariance principle, the global asymptotic stability of the disease-free equilibrium and the endemic equilibrium is obtained. It is shown that the global properties of the system depend on both the properties of these general functions and the basic reproductive number .

1. Introduction

The mechanism of transmission is usually qualitatively known for most diseases from epidemiological point of view. For modeling the spread process of infectious diseases mathematically and quantitatively, many classical epidemic models have been proposed and studied, such as SIR, SIS, SEIR, and SIRS models. Recently, considerable attention have been paid to study the dynamics of epidemic models with epidemiological meaningful time delays.

The fundamental assumption in epidemic models is that the population can be divided into distinct groups. The most common groups are the susceptible which contains individuals that may be infected by the disease; the infected which contains individuals that are already infected and can spread the disease to susceptible individuals; the removed which contains individuals that have the immunity and cannot be infected. Therefore such models are referred to SIR models. The simplest forms of these models are ordinary differential equations (ODEs).

It is well known that the disease transmission progress plays an important role in the epidemic dynamics; that is, applying different incidence rates can potentially change the behaviors of the system. In many epidemic models, following incidence functions with delay or without delay are widely used in different epidemiological backgrounds.

(1)The bilinear incidence rate (e.g., [1–8]), where is the average number of contacts per infected individual per day.(2)The standard incidence rate ([9–12]), where .(3)The Holling type incidence rate of the form ([13–15]), where is a positive constant.(4)The saturated incidence rate of the form ([16–21]), where is a positive constant.(5)The saturated incidence rate of the form ([22–25]), where , is a positive constant.

The bilinear incidence rate in (1) is based on the law of mass action, which is more appropriate for communicable diseases, such as influenza, but not suitable for sexually transmitted diseases. It has been pointed out that standard incidence rate in (2) may be a good approximation when the number of available partners is large enough and everybody could not make more contacts than that is practically feasible. In fact, the infection probability per contact is likely influenced by the number of infective individuals, because more infective individuals can increase the infection risk.

In the incidence rates in (3) and (4), measures the infection force of the disease and , measure the inhibition effect from the behavioral changes of the susceptible individuals when their number increases or from the crowding effect of the infective individuals. In these incidence rates, the number of effective contacts between infectived and susceptible individuals may saturate at high infective levels. These incidence rates seem more reasonable than the bilinear incidence rate , because they include the behavioral changes of susceptible individuals and crowding effect of the infective individuals which prevent the unboundedness of the contact rate by choosing suitable parameters.

Obviously, the incidence rate in (5) includes the former three incidence rates: the bilinear incidence rate (when , ), the Holling type incidence rate (when ), and the saturated incidence rate (when ).

The incidence rate can also be modeled by many other kinds of more general functions. It is interesting that whether the functional form of the incidence rate can change the epidemic dynamics or not. Korobeinikov studied the global properties for epidemiological models with various nonlinear incidence rates, such as in [26], in [27–29]. By constructing Lyapunov functions, [27, 28] established the global stability for ordinary differential equations models of epidemiological dynamics with nonlinear incidence rate under some conditions.

These models have not included time delays, which are usually used to model the fact that an individual may not be infectious until some time after becoming infected. In the context of epidemiology, delays can be caused by a variety of factors. The most common reasons for a delay are (i) the latency of the infection in a vector and (ii) the latency of the infection in an infected host [30]. In these cases, some time should elapse before the level of infection in the infected host or the vector reaching a sufficiently high level to transmit the infection further.

Motivated by all the above, we present a model described by delay differential equations (DDEs) with two general nonlinear terms as follows:

where , , and , as mentioned above, represent the population of the susceptible, the infected, and the removed at time , respectively. The parameters in the equations are explained as below. The positive is the recruitment rate of the population, is the natural death rate of the population, is the death rate due to disease, all is the latent period. The term represents the survival rate of population and the time they take to become infectious is . We assume that the force of infection at any time is given by the general function , and the recovered infected individuals at any time is given by the function .

Since does not appear in equations for and , it is sufficient to analyze the behaviors of solutions of (1) by the following system of DDEs:

The initial conditions for system (2) take the form

where . Here denotes the Banach space of continuous functions mapping the interval into , equipped with the supremum norm. The nonnegative cone of is defined as .

The organization of this paper is as follows. In Section 2, we study the existence of a positive equilibrium. In Section 3, we show that the global asymptotic stability of the disease-free equilibrium and the endemic equilibrium of model (2) depend only on the basic reproductive number under some hypotheses. A brief discussion is given in the last section to conclude this paper.

2. The Existence of Positive Equilibrium

In this section, we prove the existence of a positive equilibrium. We assume that and are always positive, continuously differentiable, and monotonically increasing for all and . That is, they satisfy the following conditions:

(H1), , for and .(H2), , for and .(H3). for .

Global behaviors of system (2) may depend on the basic reproduction number , which is the average number of secondary cases produced by a single infective individual introduced into an entirely susceptible population. The basic reproductive number for system (2) can be computed as

where . Usually, implies that the number of infected individuals will tend to zero and implies that the number will increase.

The epidemiologically natural condition ensures that system (2) always has a disease-free equilibrium . And it may also admit an endemic equilibrium which depends on . And , satisfy the following equations:

We will show that under certain epidemiologically reasonable conditions, the existence of the positive equilibrium state is ensured. We have the following theorem.

Theorem 1. Assume that satisfies (H1) and (H2), and satisfies (H3). Then system (2) has a positive equilibrium state if .

Proof. Let the right-hand sides of the three equations in system (2) equal zero; we have that

Substituting the expression of by , we obtain the following equation for :

It is obvious that , and we can compute that there exists a positive such that . Hence

And when , since is continuously differentiable, we have

Thus, ensures that . And is continuous on ; then there exists some such that . Since is strictly monotonically increasing, we have . Therefore , and we have proved the existence of the endemic equilibrium for system (2) under condition . This completes the proof.

3. Global Dynamics of the Model

In this section, we will analyze the global dynamics of system (2) and show the global asymptotic stability of the disease-free equilibrium and the endemic equilibrium.

3.1. Stability of the Disease-Free Equilibrium

In this subsection, we will study the global stability of the disease-free equilibrium of system (2). We propose the following conditions:

(H4) is increasing with respect to .(H5) with respect to .(H6) with respect to .

By (H4), the following inequalities hold true:

Under these conditions, we have the following theorems.

Theorem 2. Suppose that conditions (H1)–(H3) are satisfied. Then the disease-free equilibrium of system (2) is locally asymptotically stable for any if ; is unstable if .

Proof. The characteristic equation of system (2) at is

It has a negative real root . Moreover, it has a root of

In (12), if , becomes ; one can see that . Hence the is locally asymptotically stable. In (12), if ,  ,  ,   is a root of (12), then we have

Further we have

which is a contradiction. Hence the is locally asymptotically stable.

If , let ; then we have

And when , since is continuously, then has at least one positive root. Hence is unstable. This completes the proof.

Theorem 3. Suppose that conditions (H1)–(H6) are satisfied. Then the disease-free equilibrium of system (2) is globally asymptotically stable for any if .

Proof. Define a Lyapunov functional

where

By (H1)–(H6), it is obvious that is defined and continuously differentiable for all , , and at . The system (2) at has . The time derivative of along the solutions of system (2) is given by

Further, we have

Thus

By (10), we have

Note that , for , , and by and (H4), the equality in (21) holds if and only if . Furthermore, (H5) and (H6) imply that

Therefore, ensures that for all , , where holds only for . It is easy to verify that the disease-free equilibrium is the only fixed point of the systems on the plane and hence it is easy to show that the largest invariant set in is the singleton . By the Lyapunov-LaSalle asymptotic stability theorem in [31], is globally asymptotically stable for any . This completes the proof.

3.2. Global Stability of the Endemic Equilibrium

In this subsection, we will study the global stability of the endemic equilibrium of system (2) by the Lyapunov direct method. We propose the following hypotheses:

(H7) for ,   for .(H8) for , for .

Based on these, we have the following theorem.

Theorem 4. Suppose that conditions (H1)–(H3) and (H7)-(H8) are satisfied. Then the endemic equilibrium of system (2) is globally asymptotically stable for any if .

Proof. Define a Lyapunov functional

where

By (H1)–(H6), is defined and continuously differentiable for all , . And at . At , system (2) has

The time derivative of along the solutions of system (2) is given by

Further, we have

Then we have

The function is monotonically increasing for any ; hence the following inequality holds:

And by the properties of the function , (), we note that has its global maximum . Hence when and the following inequalities hold true:

Furthermore, by (H7) the following inequality holds:

And by (H8) we have the following inequality:

By (29)–(32), we see that for all , . It is easy to verify that the largest invariant set in is the singleton . By the Lyapunov-LaSalle asymptotic stability theorem in [31], is globally asymptotically stable for any . This completes the proof.

4. Discussion and Conclusion

In this paper, we formulated an SIR epidemic model with delay and two general functions, one is which represents the incidence rate, and the other is which represents the recovered infected individuals from the infected compartment. We studied the global asymptotic stability of disease-free equilibrium and endemic equilibrium of system (2), respectively. We showed that in Theorem 2 the disease-free equilibrium is locally asymptotically for any if the basic reproduction number and are unstable if ; in Theorem 3 the disease-free equilibrium is globally asymptotically for any if , while in Theorem 4, the endemic equilibrium is globally asymptotically for any if .

In order to obtain the global properties of the system (2), we proposed assumptions (H1)–(H8) for functions of and . Conditions (H1)–(H3) are some basic assumptions; for example, (H1) implies that the function is a nonnegative differentiable function on nonnegative quadrant and is positive if and only if both arguments are positive. We used (H4)–(H6), (H7)-(H8) to establish the global asymptotic stability of disease-free equilibrium and endemic equilibrium of system (2), respectively. These hypotheses seem to be mathematical techniques; however, they may be obviously true for many concrete forms of the functions of and in previous studies.

A special case of system (2) is that when , is the recovery rate of the infective individuals. System (2) becomes the following DDEs:

The basic reproductive number for system (33) is presented as

Using Theorems 3 and 4, we can easily obtain the global asymptotic stability of the disease-free equilibrium and the endemic equilibrium of system (33). Regarding to system (33), we now give examples of incidence function that satisfies the required hypotheses obviously.

Example 1. Without delay: let , and . Then system (33) becomes to the SIR model studied in [27].

Example 2. Holling type II incidence rate: let for some constant . Then the hypotheses on are satisfied and the global properties are determined by the basic reproductive number. This model was introduced by [32, 33] for considering delays in the standard bacterial growth model in a chemostat. Its global dynamics were first proved by [32] by the fluctuation lemma and a different proof was given in [34, Theorem 5.16] by comparison.

Example 3. Saturate incidence rate: let for some constant , . Then hypotheses about are also satisfied and the global properties are determined by the basic reproductive number. The behaviors of this model were previously studied in [23, 24]. In [23], the local stability of disease-free equilibrium and endemic equilibrium was obtained. And in [24], the global stability of disease-free equilibrium and endemic equilibrium was studied.

From these Examples 1–3, we can see that system (2) is reasonably established and it can contain many classical epidemic models and imply their global dynamics as special cases.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Nature Science Foundation of China (11271303 and 10971168) and the Fundamental Research Funds for the Center Universities (XDJK2009B012).

References

M. Gabriela, M. Gomes, L. J. White, and G. F. Medley, “The reinfection threshold,” Journal of Theoretical Biology, vol. 236, no. 1, pp. 111–113, 2005.

Y. Zhou and H. Liu, “Stability of periodic solutions for an SIS model with pulse vaccination,” Mathematical and Computer Modelling, vol. 38, no. 3-4, pp. 299–308, 2003.

A. Gray, D. Greenhalgh, L. Hu, X. Mao, and J. Pan, “A stochastic differential equation SIS epidemic model,” SIAM Journal on Applied Mathematics, vol. 71, no. 3, pp. 876–902, 2011.

M. Song, W. Ma, and Y. Takeuchi, “Permanence of a delayed SIR epidemic model with density dependent birth rate,” Journal of Computational and Applied Mathematics, vol. 201, no. 2, pp. 389–394, 2007.

W. Ma, M. Song, and Y. Takeuchi, “Global stability of an SIR epidemic model with time delay,” Applied Mathematics Letters, vol. 17, no. 10, pp. 1141–1145, 2004.

F. Zhang, Z. Li, and F. Zhang, “Global stability of an SIR epidemic model with constant infectious period,” Applied Mathematics and Computation, vol. 199, no. 1, pp. 285–291, 2008.

X. Liu, Y. Takeuchi, and S. Iwami, “SVIR epidemic models with vaccination strategies,” Journal of Theoretical Biology, vol. 253, no. 1, pp. 1–11, 2008.

J. Wang, J. Zhang, and Z. Jin, “Analysis of an SIR model with bilinear incidence rate,” Nonlinear Analysis, vol. 11, no. 4, pp. 2390–2402, 2010.

A. Korobeinikov and G. C. Wake, “Lyapunov functions and global stability for SIR, SIRS, and SIS epidemiological models,” Applied Mathematics Letters, vol. 15, no. 8, pp. 955–960, 2002.

S. M. O'Regan, T. C. Kelly, A. Korobeinikov, M. J. A. O'Callaghan, and A. V. Pokrovskii, “Lyapunov functions for SIR and SIRS epidemic models,” Applied Mathematics Letters, vol. 23, no. 4, pp. 446–448, 2010.

L. Liu, X. Q. Zhao, and Y. Zhou, “A tuberculosis model with seasonality,” Bulletin of Mathematical Biology, vol. 72, no. 4, pp. 931–952, 2010.

N. Yoshida and T. Hara, “Global stability of a delayed SIR epidemic model with density dependent birth and death rates,” Journal of Computational and Applied Mathematics, vol. 201, no. 2, pp. 339–347, 2007.

R. M. Anderson and R. M. May, “Regulation and stability of host-parasite population interactions. I. Regulatory processes,” Journal of Animal Ecology, vol. 47, no. 1, pp. 219–267, 1978.

C. Wei and L. Chen, “A delayed epidemic model with pulse vaccination,” Discrete Dynamics in Nature and Society, vol. 2008, Article ID 746951, 12 pages, 2008.

J. Zhang, Z. Jin, Q. Liu, and Z. Zhang, “Analysis of a delayed SIR model with nonlinear incidence rate,” Discrete Dynamics in Nature and Society, vol. 2008, Article ID 636153, 16 pages, 2008.

Z. Jiang and J. Wei, “Stability and bifurcation analysis in a delayed SIR model,” Chaos, Solitons and Fractals, vol. 35, no. 3, pp. 609–619, 2008.

R. Xu and Z. Ma, “Stability of a delayed SIRS epidemic model with a nonlinear incidence rate,” Chaos, Solitons and Fractals, vol. 41, no. 5, pp. 2319–2325, 2009.

R. Xu, Z. Ma, and Z. Wang, “Global stability of a delayed SIRS epidemic model with saturation incidence and temporary immunity,” Computers & Mathematics with Applications, vol. 59, no. 9, pp. 3211–3221, 2010.

R. Xu and Z. Ma, “Global stability of a SIR epidemic model with nonlinear incidence rate and time delay,” Nonlinear Analysis, vol. 10, no. 5, pp. 3175–3189, 2009.

X. Zhang and X. Liu, “Backward bifurcation of an epidemic model with saturated treatment function,” Journal of Mathematical Analysis and Applications, vol. 348, no. 1, pp. 433–443, 2008.

C. C. McCluskey, “Global stability for an SIR epidemic model with delay and nonlinear incidence,” Nonlinear Analysis, vol. 11, no. 4, pp. 3106–3109, 2010.

A. Kaddar, “On the dynamics of a delayed SIR epidemic model with a modified saturated incidence rate,” Electronic Journal of Differential Equations, vol. 2009, no. 133, pp. 1–7, 2009.

A. Kaddar, A. Abta, and H. T. Alaoui, “A comparison of delayed SIR and SEIR epidemic models,” Nonlinear Analysis, vol. 16, no. 2, pp. 181–190, 2011.

A. Abta, A. Kaddar, and H. T. Alaoui, “Global stability for delay SIR and SEIR epidemic models with saturated incidence rates,” Electronic Journal of Differential Equations, vol. 2012, no. 23, pp. 1–13, 2012.

Z. Liu, “Dynamics of positive solutions to SIR and SEIR epidemic models with saturated incidence rates,” Nonlinear Analysis, vol. 14, no. 3, pp. 1286–1299, 2013.

A. Korobeinikov and K. Philip Maini, “Nonliear incidence and stability of infectious diseasemodels,” Mathematical Medicine and Biology, vol. 22, pp. 113–128, 2005.

A. Korobeinikov, “Lyapunov functions and global stability for SIR and SIRS epidemiological models with non-linear transmission,” Bulletin of Mathematical Biology, vol. 68, no. 3, pp. 615–626, 2006.

A. Korobeinikov, “Global properties of infectious disease models with nonlinear incidence,” Bulletin of Mathematical Biology, vol. 69, no. 6, pp. 1871–1886, 2007.

A. Korobeinikov, “Stability of ecosystem: global properties of a general predator-prey model,” Mathematical Medicine and Biology, vol. 26, no. 4, pp. 309–321, 2009.

G. Huang, Y. Takeuchi, W. Ma, and D. Wei, “Global stability for delay SIR and SEIR epidemic models with nonlinear incidence rate,” Bulletin of Mathematical Biology, vol. 72, no. 5, pp. 1192–1207, 2010.

J. K. Hale and S. M. Verduyn Lunel, Introduction to Functional Differential Equations, Springer, New York, NY, USA, 1993.

S. F. Ellermeyer, “Competition in the chemostat: global asymptotic behavior of a model with delayed response in growth,” SIAM Journal on Applied Mathematics, vol. 54, no. 2, pp. 456–465, 1994.

S. F. Ellermeyer, J. Hendrix, and N. Ghoochan, “A theoretical and empirical investigation of delayed growth response in the continuous culture of bacteria,” Journal of Theoretical Biology, vol. 222, no. 4, pp. 485–494, 2003.

H. Smith, An Introduction to Delay Differential Equations with Applications to the Life Sciences, Springer, New York, NY, USA, 2011.