Mathematical Analysis of a Cholera Model with Vaccination

Cui, Jing'an; Wu, Zhanmin; Zhou, Xueyong

doi:https://doi.org/10.1155/2014/324767

Journal of Applied Mathematics

On this page

Abstract Introduction Conclusion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2014 | Article ID 324767 | https://doi.org/10.1155/2014/324767

Mathematical Analysis of a Cholera Model with Vaccination

Jing'an Cui,¹Zhanmin Wu,¹and Xueyong Zhou²

Academic Editor: Tin-Tai Chow

Received14 Mar 2013

Revised18 May 2013

Accepted26 Sept 2013

Published13 Feb 2014

Abstract

We consider a SVR-B cholera model with imperfect vaccination. By analyzing the corresponding characteristic equations, the local stability of a disease-free equilibrium and an endemic equilibrium is established. We calculate the certain threshold known as the control reproduction number . If , we obtain sufficient conditions for the global asymptotic stability of the disease-free equilibrium; the diseases will be eliminated from the community. By comparison of arguments, it is proved that if , the disease persists and the unique endemic equilibrium is globally asymptotically stable, which is obtained by the second compound matrix techniques and autonomous convergence theorems. We perform sensitivity analysis of on the parameters in order to determine their relative importance to disease transmission and show that an imperfect vaccine is always beneficial in reducing disease spread within the community.

1. Introduction

Cholera is an acute intestinal infection caused by the ingestion of food or water contaminated with the bacterium Vibrio cholera. Among the 200 serogroups of Vibrio cholera, it is only Vibrio cholera o1 and o139 that are known to be the cause of the cholera disease [1]. The etiological agent, Vibrio cholera o1 (and more recently vibrio cholera o139), passes through and survives the gastric acid barrier of the stomach and then penetrates the mucus lining that coats the intestinal epithelial [2]. Once they colonise the intestinal gut, then produce enterotoxin (which stimulates water and electrolyte secretion by the endothelial cells of the small intestine) that leads to copious, painless, and watery diarrhoea that can quickly lead to severe dehydration and death if treatment is not promptly given. Vomiting also occurs in most patients. In human volunteer studies, the infection was determined to be 10²–10³ [3]. Cholera can either be transmitted through interaction between humans (i.e., fecal-oral) or through interaction between humans and their environment (i.e., ingestion of contaminated water and food from the environment). To come on urgent, transmission fast, sweeping range widely are the characteristics of cholera which is one of the international quarantine infectious diseases as stipulated by the International Health Regulations (IHR), as well as one of a class of infectious diseases as stipulated by law for infectious diseases prevention and control of China.

Globally, cholera incidence has increased steadily since 2005 with cholera outbreaks affecting several continents (see Figure 1). Cholera continues to pose a serious public health problem among developing world populations which have no access to adequate water and sanitation resources. In 2011, 32% of cases were reported from Africa whereas between 2001 and 2009, 93% to 98% of total cases worldwide were reported from that continent [4]. In 2011, 61.2% of cases were reported from Americas where a large outbreak that started in Haiti at the end of October 2010 also affected the Dominican Republic. The outbreak was still ongoing at the end of 2011 with 522335 cases including 7001 deaths that were reported by 25 December in Haiti [4]. So, cholera remains a global threat and is one of the key indicators of social development. The history and reality have warned: we are facing the serious threat of cholera, the importance for the study of cholera’s pathogenesis, regular transmission and prevention and control strategy have become increasingly prominent, which has also become a major problem that needs to be solved. Up till now, a number of mathematical models have been used to study the transmission dynamics of cholera. Capasso and Serio [5] introduced an incidence rate in the form of (with human-to-human transmission model only) in 1973 [6]. Codeço [7] proposed an incidence form of (with environment-to-human transmission model only) in 2001 which, in the first time, explicitly incorporated the pathogen concentration into cholera modeling. Mukandavire et al. [8] included both transmission pathways in the form of . In 2012, Liao and Wang [9] generalised Codeco’s model [7], incorporating the theory of Volterra-Lyapunov stable matrices into the classic method of Lyapunov functions; they studied the global dynamics of the mathematical model.

Vaccination is a major factor in the resurgence and epidemic outbreaks of some infectious diseases. Since the pioneering work of Edward Jenner on smallpox [10], vaccination has been a commonly used method for diseases control [11–13] and works by reducing the number of susceptible individuals in a population. In modern times, vaccination has a large impact on the incidence and persistence of children infections, such as measles and whooping cough [14]. Hethcote [11] investigated a pertussis infectious model and showed that the vaccination can make the infection undergo fluctuation. Although vaccination offers a very powerful tool for disease control, generally, vaccines are not 100% affective and sometimes they only provide limited immunity due to the natural waning of immunity in the host or antigenic variation in the pathogen [15].

Here, we develop a cholera model with an additional equation for the vaccinated individuals in the population. Since Koch found Vibrio cholera in 1883, the research for cholera vaccine had been going on for over one hundred years. People have developed a variety of vaccines. However, these vaccines were parenteral, which have short effective protection and big side effects. In 1973, the World Health Organization canceled the vaccine inoculation which attracted a major concern to oral vaccines. At present, there are three kinds of oral vaccines (i.e., WC/BS vaccine, WC/rBS vaccine, and CVD₁₀₃-HgR vaccine) have been proved to be safe, effective, and immunogenic, which were approved to apply in some countries [16].

In this paper, according to the natural history of cholera, we improve the model of [9] in the following two aspects. Firstly, if the cholera persists for a long time, it will cause death [17], especially in the area where water and sanitation resources are not adequate [4]; a parameter is added to describe the rate of disease-related death. Secondly, we propose a proportion of the vaccination in susceptible individuals as shown in the following differential equations:

The flow diagram of the model is depicted in Figure 2. Since the first three and last equations in (1) are independent of the variable , it suffices to consider the following reduced model:

Here, , , , and refer to the susceptible individuals, infected individuals, vaccinated individuals, and recovered individuals, respectively. The pathogen population at time is given by . The parameter denotes the natural human birth and death rate, denotes the rate of recovery from the disease, represents the rate of human contribution to the growth of the pathogen, and represents the death rate of the pathogen in the environment. The coefficients and represent the contact rates for the human-environment and human-human interactions, respectively. Constants and adjust the appropriate form of the incidence which determines the rate of new infection. If , the corresponding incidence is reduced to the standard bilinear form based on the mass action law, which is most common in epidemiological models. If , then the corresponding incidence represents a consequence of saturation effects: when the infected number is high, the incidence rate will respond more slowly than linearly to the increase in . Similar meanings stand for . The rate at which the susceptible population is vaccinated is , and the rate at which the vaccine wears off is . All parameters are assumed nonnegative.

The organization of this paper is as follows: the positivity and boundedness of solutions are obtained in Section 2. In Section 3, we obtain the existence of the endemic equilibrium. We get the local and global stability of the disease-free equilibrium in Section 4. In Section 5, we present the persistence of the system. In Section 6, we show the local and global stability of the endemic equilibrium. We analyze the sensitivity of on the parameters, and we present the numerical simulation in Section 7. The paper ends with a conclusion in Section 8.

2. Positivity and Boundedness of Solutions

In the following, we show that the solutions of system (2) are positive with the nonnegative initial conditions.

Theorem 1. The solutions of model (2) are nonnegative for all with the non-negative initial conditions.

Proof. System (2) can be put into the matrix form where and is given by We have Therefore, Due to Lemma 2 in [18], any solution of system (2) is such that for all . This completes the proof of Theorem 1.

Theorem 2. All solutions of model (2) are bounded.

Proof. System (2) is split into two parts, the human population (i.e., , , and ) and pathogen population (i.e., ). It follows from the first three equations of system (2) that then it follows that . From the first equation, we can get Thus , as . It is easy to obtain Thus , as . From the last equation, we can obtain Hence, , when . Therefore, all solutions of model (2) are bounded.

From above discussion, we can see that the feasible region of human population for system (2) is and the feasible region of pathogen population for system (2) is

Define . Let denote the interior of . It is easy to verify that the region is a positively invariant region (i.e., the solutions with initial conditions in remain in ) with respect to system (2). Hence, we will consider the global stability of (2) in region .

3. Equilibria

In this section, we investigate the existence of equilibria of system (2). Solving the right hand side of model system (2) by equating it to zero, we obtain the following biologically relevant equilibria.

It is easy to see that model (2) always has a disease-free equilibrium (the absence of infection, i.e., ) , where and .

In the following, we will discuss the existence and uniqueness of the endemic equilibrium. The components of the endemic equilibrium satisfy

Substituting (13a), (13b), and (13d) into (13c), we obtain a single equation for : After dropping the solution , we obtain where Note that represents a straight line with a negative slope and a vertical intercept . Meanwhile, we have We see that is increasing for , and , where is the control reproduction number of infection. When . Hence, there is one and only one intersection between the curves of and ; that is, there is a unique solution to the equation . Consequently, , , and are uniquely determined by .

Theorem 3. System (2) has a unique endemic equilibrium when and no positive endemic equilibrium when .

4. Stability of Disease-Free Equilibrium

Now, we will discuss the local and global stability of the disease-free equilibrium.

Theorem 4. The disease-free equilibrium is locally asymptotically stable for and unstable for .

Proof. The Jacobian matrix of system (2) at is The characteristic polynomial of the matrix is given by where If , then Further, So , , , , , and (see Appendix A). Thus, using the Routh-Hurwitz criterion, all eigenvalues of have negative real part; is local asymptotically stable for system (2). If , then and we show that has at least one eigenvalue with nonnegative real parts. Consequently, is not asymptotically stable.

Theorem 5. When , the disease-free equilibrium is globally asymptotically stable.

We will prove the global stability of the disease-free equilibrium using Lemma 6.

Lemma 6 (see [19]). If a model system can be written in the form where denotes (its components) the number of uninfected individuals and denotes (its components) the number of infected individuals including latent and so forth, denotes the disease-free equilibrium of the system.
Assume that(H1)for , is globally asymptotically stable;(H2), for , where the Jacobian matrix is an Metzler matrix (the off-diagonal elements of are nonnegative) and is the region where the model makes biological sense. Then the fixed point is a globally asymptotically stable equilibrium of cholera model system (2) provided that .
We begin by showing condition (H1) as For the equilibrium , the system reduces to The characteristic polynomial of the system is given by There are two negative characteristic foots: . Hence, is always globally asymptotically stable.
Next, applying Lemma 6 to the cholera model system (2) gives So is a Metzler matrix. Meanwhile, we find . Hence, the disease-free equilibrium is globally asymptotically stable.

5. Persistence

Persistence is an important property of dynamical systems and of the systems in ecology, epidemics, and so forth, they are modeling. Biologically, persistence means the survival of all populations in future time. Mathematically, persistence of a system means that strictly positive solutions do not have any omega limit points on the boundary of the nonnegative cone [20]. In this section, we will present the persistence of system (2).

For various definitions of persistence [21, 22], we utilize the definitions of persistence developed by Freedman et al. [23]. System (2) can be defined to be uniformly persistent if for some for all initial points in .

A uniform persistence result given in [23] requires the following hypothesis (H) to be satisfied.

We denote that is a closed positively invariant subset of on which a continuous flow is defined and is the maximal invariant set of on . Suppose is a closed invariant set and there exists a cover of , where is a nonempty index set; , , and are pairwise disjoint closed invariant sets. Furthermore, we propose the following hypothesis.

Hypothesis (H):(a)all are isolated invariant sets of the flow ;(b) is acyclic; that is, any finite subset of does not form a cycle [24];(c)any compact subset of contains, at most, finitely many sets of .

Lemma 7 (see [24]). Let be a closed positively invariant subset of on which a continuous flow is defined. Suppose there is a constant such that is point dissipative on and the assumption (H) holds. Then the flow is uniformly persistent if and only if for any , where .

Now, we can obtain the following result.

Theorem 8. System (2) is uniformly persistent in if .

Proof. Suppose . We show that system (2) satisfies all the conditions of Lemma 7. Choose and . The vector field of system (2) is transversal to the boundary of on its faces except the S-axis and V-axis, which are invariant with respect to system (2) and on the S-axis and V-axis the equations for and are and , which implies that and as . Therefore, is the only -limit point on the boundary of . As the maximal invariant set on the boundary of is the singleton and is isolated when , thus the hypothesis (H) holds for system (2). The flow induced by is point dissipative by the positive invariance of . Because of , where is the omega limit set of , when , we have that is contained in the set and for , . Therefore, the uniform persistence of system (2) is equivalent to being unstable, and the theorem is proved.

Remark 9. Theorems 3 and 8 show that is a threshold parameter for the model; that is, when , its epidemiological implication is that the infected fraction of the population vanishes, so the cholera dies out; when , the disease is endemic and the infected fraction remains above a certain positive level for sufficiently large time.

6. Stability of the Endemic Equilibrium

6.1. Local Stability of the Endemic Equilibrium

Now we consider the case with . The stability of the endemic equilibrium is established as follows.

Theorem 10. If , is locally asymptotically stable.

Proof. Let The Jacobian matrix at is The characteristic polynomial of the matrix is given by where Based on (13c) and (13d), we have . It is then easy to observe that Further, So , , , , , and (see Appendix B). Using the well-known Routh-Hurwitz criterion, the proof is thus complete.

6.2. Global Stability of the Endemic Equilibrium

Theorem 11. When , is globally asymptotically stable in if and .

This approach to global dynamics is developed in the papers of Smith [25] and Li and Muldowney [26–28]. Let , where ,,, and represent the right-hand sides of system (2), respectively. Furthermore, let . Then, the Jacobian matrix for system (2) is

The second additive compound matrix (see Appendix C for details) of iswhere

Let Then we havewhere In (41), is the directional derivative of in the direction of the vector field in system (2).

Here, we will use the theorem in [28] to give a sufficient condition on the parameters, which when satisfied, implies that the endemic equilibrium is globally asymptotically stable.

Lemma 12 (see [28]). If is a compact absorbing subset in the interior of and there exists and a Lozinskiǐ measure for all , then every omega point of system (2) in is an equilibrium in .

For , the disease-free equilibrium is repelling towards the interior. In fact, there is a compact absorbing set in int which attracts all orbits that intersect int. This gives the following results.

Corollary 13. If and there exists a Lozinskiǐ measure such that for all , then each orbit of system (2) which intersects limits to the endemic equilibrium.

For a norm on , the Lozinskiǐ measure associated with can be evaluated for an matrix as where is the right-hand derivative [29]. Hence, if we can find a norm on for which the associated Lozinskiǐ measure satisfied for all then the endemic equilibrium is globally asymptotically stable for .

We now define a norm on [30] for which the definition varies from one orthant to another. Let where , with components and

Theorem 14. Assume that Then there exits , such that for all .

Proof. We should show that

Case 1 ( and ). In this case, .

Subcase 1 (). Then and . Taking the right-hand derivative of , we obtain Since , we have Hence,
It is easy to see that (50) also holds for and when , which can be obtained by linearity.

Subcase 2 (). Then and .
Since , taking the right-hand derivative of , we obtain Therefore, It is easy to see that (52) also holds for and when , which can be obtained by linearity.

Case 2 ( and ). In this case, .

Subcase 3 (). Then and . Taking the right-hand derivative of , we obtain Therefore, It is easy to see that (54) also holds for and when , which can be obtained by linearity.

Subcase 4 (). Then and . Taking the right-hand derivative of , since , we obtain Therefore,

It is easy to see that (56) also holds for and when , which can be obtained by linearity.

Case 3 ( and ). In this case, .

Subcase 5 (). Then and . Taking the right-hand derivative of , we obtain Recalling that yields It is easy to see that (58) also holds for and when , which can be obtained by linearity.

Subcase 6 (). Then and . Taking the right-hand derivative of , we obtain Therefore, It is easy to see that (60) also holds for and when , which can be obtained by linearity.

Combing the results of the six cases presented here in (50)–(60), as well as the remaining 10 cases, we obtain the result when . Hence, . Therefore, in int if and . Thus, Theorem 14 is satisfied. From Corollary 13, we can obtain that all solutions intersect the interior of and the endemic equilibrium is globally asymptotically stable, completing the proof of Theorem 11.

7. Sensitivity Analysis of

To facilitate the interpretation of the sensitivity of , we now present some numerical simulations by using the set of parameters values in Table 1.

Now, we regard the vaccinated rate and the waning rate as the control parameters, while the other parameters are fixed. From Figures 3 and 4, the effects of various parameters, that is, and on the control reproduction number have been shown. It is noted that as the parameter increases, decreases and as decreases, decreases. In fact, we can obtain the critical values of and that reduce to 1: In Figure 3, we select , , , , corresponding , ,,, respectively. We can see that when the wanning rate has a greater value, then there is no such that . Similarly, in Figure 4, we select , , , , , corresponding , , , , , respectively. We can see that when is smaller, then there is no such that . Thus, the control reproduction number cannot reduce below unity only by increasing or decreasing . The critical values and are important in regulating the infection magnitude. In order to reduce to 1, a greater vaccinated rate than and a smaller wanning rate than have to be achieved simultaneously. We will deduce below 1 by using both and at the same time, which can control cholera. Otherwise, the cholera persists (see Figure 5).

8. Conclusion

In this paper, we have conducted global stability analysis of a SVIR-B cholera model. Based on the imperfective vaccine, with the environment component incorporated and multiple transmission pathways coupled, the cholera models distinguish themselves from regular SIR epidemiological models. The mathematical results show that the control reproduction number satisfies a threshold property with threshold value 1. When , it has been proved that the disease-free equilibrium is globally asymptotically stable under some sufficient conditions. And, when , the unique endemic equilibrium is globally asymptotically stable. This shows that cholera can be eliminated from the community if the imperfect vaccine brings to a value less than unity.

Now, if we consider that the case there is no vaccination. Model (1) can be written as follows:

According to Theorem 2 in [31], the basic reproduction number of model is . We can express the control reproduction number as . Note that with equality only if . That is, despite being imperfect, the vaccine will always reduce the reproduction number of the disease, and in the absence of the vaccination, the disease transmission will be high. Further simulations in Section 7 also show that the vaccinated rate () and the wanning rate () play equal important roles in reducing the control reproduction number. The disease can be controlled if and only if the reproduction number is reduced to values less than unity if the vaccinated rate () exceed the threshold and the wanning rate () less than the threshold simultaneously. Therefore, vaccination is a good way to control cholera.

However, there are inherent disadvantages towards the vaccination modeling. For cholera with incubation period, it is hard to rapidly identify those with ambiguous symptoms [4]. Moreover, the vaccination does not always work well due to the limitations of medical development level and financial budget (some vaccines are very expensive and some portions of people cannot be covered) [32]. Nonetheless, [33] indicated that during cholera outbreaks (periodic control) vaccination campaigns can be a good strategy to control cholera epidemics. Besides, they pointed out that vaccination and improvement in the sanitation system and food/personal hygiene are the most efficient control strategies to prevent cholera transmission and outbreaks. Hence, incorporating some other control strategies, we may consider the more realistic ordinary differential equation model. The theoretical study of cholera models has been in progress and is an exciting area of future research.

Appendices

A. The Proof of

The proof of = + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + .

From Section 4, we know that . After some algebraic manipulations, we have . Thus, , when .

B. The Proof of

The proof of = + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + .

From Section 6, we know that . After some algebraic manipulations, we have . Thus, , when .

C. The Second Additive Compound Matrix

If is an matrix, its second additive compound is matrix defined as follows. For any integer , let such that . Then the element in the -row and -column of is In the special case when , we have . In general, each entry of is a linear expression of those of . For instance, when , the second additive compound matrix of is when , the second additive compound matrix of is

For detailed discussions of compound matrices and their properties, readers can refer to [34, 35]. A comprehension survey on compound matrices and their relations to differential equations is given in [34].

Conflict of Interests

The authors declare that they have no financial and personal relationships with other people or organizations that can inappropriately influence their work. There is no professional or other personal interest of any nature or kind in any product, service, or company that could be construed as influencing the position presented in or the review of the paper.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (no. 11071011), Funding Project for Academic Human Resources Development in Institutions of Higher Learning under the Jurisdiction of Beijing Municipality (no. PHR201107123), and Basic and Frontier Technology Research Program of Henan Province (nos. 132300410025 and 1323000410364).

References

A. K. T. Kirschner, J. Schlesinger, A. H. Farnleitner et al., “Rapid growth of planktonic Vibrio cholerae non-O1/non-O139 strains in a large alkaline lake in Austria: Dependence on temperature and dissolved organic carbon quality,” Applied and Environmental Microbiology, vol. 74, no. 7, pp. 2004–2015, 2008.
View at: Publisher Site | Google Scholar
J. Reidl and K. E. Klose, “Vibrio cholerae and cholera: out of the water and into the host,” FEMS Microbiology Reviews, vol. 26, no. 2, pp. 125–139, 2002.
View at: Publisher Site | Google Scholar
D. M. Hartley, J. G. Morris Jr., and D. L. Smith, “Hyperinfectivity: a critical element in the ability of V. cholerae to cause epidemics?” PLoS Medicine, vol. 3, no. 1, pp. 63–69, 2006.
View at: Publisher Site | Google Scholar
World Health Organization, 2012, http://www.who.int/gho/epidemic_diseases/cholera/en/index.html.
V. Capasso and G. Serio, “A generalization of the Kermack-McKendrick deterministic epidemic model,” Mathematical Biosciences, vol. 42, no. 1-2, pp. 43–61, 1978.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
V. Capasso and S. L. Paveri-Fontana, “Mathematical model for the 1973 cholera epidemic in the European Mediterranean region,” Revue d'Epidemiologie et de Sante Publique, vol. 27, no. 2, pp. 121–132, 1979.
View at: Google Scholar
C. T. Codeço, “Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir,” BMC Infectious Diseases, vol. 1, article 1, 2001.
View at: Publisher Site | Google Scholar
Z. Mukandavire, S. Liao, J. Wang, H. Gaff, D. L. Smith, and J. G. Morris Jr., “Estimating the reproductive numbers for the 2008-2009 cholera outbreaks in Zimbabwe,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 21, pp. 8767–8772, 2011.
View at: Publisher Site | Google Scholar
S. Liao and J. Wang, “Global stability analysis of epidemiological models based on Volterra-Lyapunov stable matrices,” Chaos, Solitons and Fractals, vol. 45, no. 7, pp. 966–977, 2012.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
F. Fenner, D. A. Henderson, I. Arita, Z. Jezek, and I. D. Ladnyi, “Smallpox and its eradition,” Journal of Community Health, vol. 13, no. 4, pp. 245–592, 1988.
View at: Google Scholar
H. W. Hethcote, “Oscillations in an endemic model for pertussis,” The Canadian Applied Mathematics Quarterly, vol. 6, no. 1, pp. 61–88, 1998.
View at: Google Scholar | Zentralblatt MATH | MathSciNet
H. W. Hethcote, “The mathematics of infectious diseases,” SIAM Review, vol. 42, no. 4, pp. 599–653, 2000.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
X. H. Tian, “Stability analysis of a delayed SIRS epidemic model with vaccination and nonlinear incidence,” International Journal of Biomathematics, vol. 5, no. 6, 18 pages, 2012.
View at: Publisher Site | Google Scholar | MathSciNet
M.-L. Garly and P. Aaby, “The challenge of improving the efficacy of measles vaccine,” Acta Tropica, vol. 85, no. 1, pp. 1–17, 2003.
View at: Publisher Site | Google Scholar
J. Arino, C. C. Mccluskey, and P. V. Van Den Driessche, “Global results for an epidemic model with vaccination that exhibits backward bifurcation,” SIAM Journal on Applied Mathematics, vol. 64, no. 1, pp. 260–276, 2003.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
Y. Kang and H. Z. Zhang, “The study status of oral cholera vaccine,” Chinese Frontier Health Quarantine, vol. 28, no. s1, pp. 91–92, 2005.
View at: Google Scholar
S. Singh, P. Chandra, and J. B. Shukla, “Modeling and analysis of the spread of carrier dependent infectious diseases with environmental effects,” Journal of Biological Systems, vol. 11, no. 3, pp. 325–335, 2003.
View at: Publisher Site | Google Scholar
X. Yang, L. Chen, and J. Chen, “Permanence and positive periodic solution for the single-species nonautonomous delay diffusive models,” Computers and Mathematics with Applications, vol. 32, no. 4, pp. 109–116, 1996.
View at: Publisher Site | Google Scholar | MathSciNet
C. Castillo-Chavez, Z. Feng, and W. Huang, “On the computation of $R_{0}$ and its role on global stability,” in Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An introduction, vol. 65 of The IMA Volumes in Mathematics and Its Applications, pp. 229–250, Springer, New York, NY, USA, 2003.
View at: Google Scholar
D. Mukherjee, “Uniform persistence in a generalized prey-predator system with parasitic infection,” BioSystems, vol. 47, no. 3, pp. 149–155, 1998.
View at: Publisher Site | Google Scholar
H. I. Freedman and P. Waltman, “Persistence in a model of three competitive populations,” Mathematical Biosciences, vol. 73, no. 1, pp. 89–101, 1985.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
V. Hutson and G. T. Vickers, “A criterion for permanent coexistence of species, with an application to a two-prey one-predator system,” Mathematical Biosciences, vol. 63, no. 2, pp. 253–269, 1983.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
H. I. Freedman, S. Ruan, and M. Tang, “Uniform persistence and flows near a closed positively invariant set,” Journal of Dynamics and Differential Equations, vol. 6, no. 4, pp. 583–600, 1994.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
G. Butler and P. Waltman, “Persistence in dynamical systems,” Journal of Differential Equations, vol. 63, no. 2, pp. 255–263, 1986.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
R. A. Smith, “Some applications of Hausdorff dimension inequalities for ordinary differential equations,” Proceedings of the Royal Society of Edinburgh Section A, vol. 104, no. 3-4, pp. 235–259, 1986.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Y. Li and J. S. Muldowney, “A geometric approach to global-stability problems,” SIAM Journal on Mathematical Analysis, vol. 27, no. 4, pp. 1070–1083, 1996.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Y. Li and J. S. Muldowney, “On Bendixson's criterion,” Journal of Differential Equations, vol. 106, no. 1, pp. 27–39, 1993.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Y. Li and J. S. Muldowney, “On R. A. Smith's autonomous convergence theorem,” The Rocky Mountain Journal of Mathematics, vol. 25, no. 1, pp. 365–378, 1995.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
R. H. Martin Jr., “Logarithmic norms and projections applied to linear differential systems,” Journal of Mathematical Analysis and Applications, vol. 45, no. 2, pp. 432–454, 1974.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
A. B. Gumel, C. C. McCluskey, and J. Watmough, “An sveir model for assessing potential impact of an imperfect anti-SARS vaccine,” Mathematical Biosciences and Engineering, vol. 3, no. 3, pp. 485–512, 2006.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
P. Van Den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Chen, “The opinion of WHO about cholera vaccine,” International Journal of Biologicals, vol. 34, no. 4, pp. 214–218, 2011.
View at: Google Scholar
R. P. Sanches, C. P. Ferreira, and R. A. Kraenkel, “The role of immunity and seasonality in cholera epidemics,” Bulletin of Mathematical Biology, vol. 73, no. 12, pp. 2916–2931, 2011.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
J. S. Muldowney, “Compound matrices and ordinary differential equations,” The Rocky Mountain Journal of Mathematics, vol. 20, no. 4, pp. 857–872, 1990.
View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Fiedler, “Additive compound matrices and an inequality for eigenvalues of symmetric stochastic matrices,” Czechoslovak Mathematical Journal, vol. 24, no. 99, pp. 392–402, 1974.
View at: Google Scholar | Zentralblatt MATH | MathSciNet

Copyright

Copyright © 2014 Jing'an Cui 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.

PDF Download Citation

Download other formats

Order printed copies

Views

4461

Downloads

2457

Citations