Global Dynamics of Avian Influenza Epidemic Models with Psychological Effect

Liu, Sanhong; Pang, Liuyong; Ruan, Shigui; Zhang, Xinan

doi:https://doi.org/10.1155/2015/913726

Computational and Mathematical Methods in Medicine

On this page

Abstract Introduction Discussion Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2015 | Article ID 913726 | https://doi.org/10.1155/2015/913726

Global Dynamics of Avian Influenza Epidemic Models with Psychological Effect

Sanhong Liu,¹Liuyong Pang,^1,2Shigui Ruan,³and Xinan Zhang¹

Academic Editor: Chung-Min Liao

Received17 Dec 2014

Revised09 Feb 2015

Accepted12 Feb 2015

Published12 Mar 2015

Abstract

Cross-sectional surveys conducted in Thailand and China after the outbreaks of the avian influenza A H5N1 and H7N9 viruses show a high degree of awareness of human avian influenza in both urban and rural populations, a higher level of proper hygienic practice among urban residents, and in particular a dramatically reduced number of visits to live markets in urban population after the influenza A H7N9 outbreak in China in 2013. In this paper, taking into account the psychological effect toward avian influenza in the human population, a bird-to-human transmission model in which the avian population exhibits saturation effect is constructed. The dynamical behavior of the model is studied by using the basic reproduction number. The results demonstrate that the saturation effect within avian population and the psychological effect in human population cannot change the stability of equilibria but can affect the number of infected humans if the disease is prevalent. Numerical simulations are given to support the theoretical results and sensitivity analyses of the basic reproduction number in terms of model parameters that are performed to seek for effective control measures for avian influenza.

1. Introduction

There are three types of influenza viruses: A, B, and C. Influenza A viruses are divided into subtypes based on two proteins on the surface of the virus: the hemagglutinin (H) and the neuraminidase (N) (CDC [1]). Avian influenza viruses with all 16 haemagglutinin H1–H16 and all 9 neuraminidase N1–N9 influenza A subtypes in the majority of possible combinations have been isolated from avian species (Alexander [2]). These viruses occur naturally among wild aquatic birds worldwide and can infect domestic poultry and other bird and animal species. Avian flu viruses do not normally infect humans. However, sporadic human infections with avian flu viruses have occurred (CDC [3]).

Avian influenza A viruses are classified into two categories: low pathogenic avian influenza A (LPAI) and highly pathogenic avian influenza A (HPAI). Three prominent subtypes of avian influenza A viruses known to infect both birds and people are influenza A H9, influenza A H5, and influenza A H7. All H9 viruses identified worldwide in wild birds and poultry are LPAI viruses. Rare, sporadic H9N2 virus infections of humans have been reported to cause generally mild upper respiratory tract illness. Most H5 and H7 viruses identified worldwide in wild birds and poultry are LPAI viruses. HPAI H5N1 virus infections in humans have been reported from 15 countries, often resulting in severe pneumonia with approximately 60 percent of mortality worldwide. H7 virus infection in humans is uncommon, but it has been documented in persons who have direct contact with infected birds, especially during outbreaks of H7 virus among poultry. In humans, LPAI (H7N2, H7N3, and H7N7) virus infections have caused mild to moderate illness; HPAI (H7N3, H7N7) virus infections have caused mild to severe and fatal illness. In March 2013, the first known human cases of infection with avian influenza H7N9 viruses were reported. They were associated with severe respiratory illness and death (CDC [4]).

In order to understand infectious diseases, mathematical models have been expensively used to analyze the epidemiological characteristics and then to provide useful control and prevention measures (Anderson and May [5], Keeling and Rohani [6]). In 2007, Iwami et al. [7] proposed ordinary differential equation (ODE) models to characterize the dynamical behavior of avian influenza between human and avian populations. Since then, various models have been developed to investigate different aspects of avian influenza transmitted by the H5N1 virus (see Lucchetti et al. [8], Iwami et al. [9], Jung et al. [10], Iwami et al. [11], Gumel [12], Agusto [13], Ma and Wang [14], Bourouiba et al. [15], Gourley et al. [16], Tuncer and Martcheva [17], etc.). Most recently, mathematical models have also been proposed to study the transmission of the avian influenza A H7N9 virus; see Zhang et al. [18], Xiao et al. [19], and Hsieh et al. [20].

Human behavior and social response play a very important role in the transmission of infectious diseases (Bauch and Galvani [21], Ferguson [22], and Funk et al. [23]). Understanding the impacts of human behavior and social response on the spread of infectious diseases are the key to improve control efforts (Funk et al. [24]). Toward the avian influenza, the general public reacts behaviourally, psychologically, and socially. After the outbreaks of the avian influenza A H5N1 virus, to determine the knowledge, attitudes, and practices relating to avian influenza in the general populations, cross-sectional surveys conducted in Thailand (Olsen et al. [25]) and China (Xiang et al. [26]) show a high degree of awareness of human avian influenza in both urban and rural populations and a higher level of proper hygienic practice among urban residents. To investigate human exposure to live-poultry and changes in risk perception and behavior after the March–May 2013 influenza A H7N9 outbreak in China, Wang et al. [27] surveyed 2,504 urban residents in 5 cities and 1,227 rural residents in 4 provinces and found that most (77%) urban respondents reported that they visited live markets less often after influenza A H7N9 outbreak. These certainly helped in controlling the further spread of avian influenza and and provide scientific support to assist the governments in developing strategies and health-education campaigns to prevent avian influenza infection among the general population.

Though behavioral and social responses to the spread of infectious diseases have been reported frequently, there is very little systematic study on their effect on the spread of infectious diseases (Funk et al. [24], Liu et al. [28]). Incorporating the dynamics of human behavior and social response in infectious diseases faces various challenges (Funk et al. [29]), for example, how to incorporate human behavior and social response in models of infectious disease dynamics and how to parameterize and measure the relevant human behavior and social response in such models.

It is well known that the incidence rate is an important factor in the transmission of infectious disease. In almost all recent models related to avian influenza, the incidence rate between susceptible avians and infective avians (and between susceptible humans and infective avians) takes the mass-action form with bilinear interactions, which is increasing and unbounded. As the surveys by Olsen et al. [25] and Xiang et al. [26] indicate, when more and more infected cases are reported via media, the general population would reduce their chances to visit the markets, so the incidence rate will in fact decrease. The nonmonotone incidence function (increasing at first, then decreasing when the number of infectives reaches a critical value, and bounding) proposed by Liu et al. [30] is suitable to model such a phenomenon; see also Ruan and Wang [31], Tang et al. [32], and Xiao and Ruan [33]. Similarly, the poultry farmers would take extra cautions and protection measures when the number of infective birds becomes larger and larger so that the incidence rate from infective birds to susceptible birds may saturate. This can be described by a saturated incidence function (increasing and bounding) used in Capasso and Serio [34].

In this paper, we construct an avian-human epidemic model with a nonlinear incidence rate for the spread of avian influenza virus from infective birds to susceptible birds describing the saturation effect within avian population and a nonmonotone incidence rate for the spread of avian influenza virus from infective birds to susceptible humans characterizing the psychological effect within humans. The paper is organized as follows. We describe the epidemic model and analyze the existence of the equilibria in Section 2. In Section 3, we study the global property of avian-only subsystem. The complete system is analyzed in Section 4. In Section 5, we present some numerical simulations to illustrate the theoretical results. A brief discussion is given in the last section.

2. The Avian-Human Influenza Epidemic Model and Equilibria

In our paper, we always assume that the avian influenza virus does not spread from person to person and mutate, and the domestic birds are the main infection source. The avian population is classified into two subclasses: susceptible and infective, denoted by and , respectively, and the human population is classified into three subclasses: susceptible, infective, and recovered, denoted by , , and , respectively. In order to construct the corresponding model, we make the following assumptions.

(1) All new recruitments and newborns of the avian population and the human population are susceptible; the rate is denoted by and , respectively.

(2) The avian influenza virus is contagious from an infective bird to a susceptible bird and from an infected bird to a susceptible human.

(3) When the number of the infective avians becomes larger, the incidence rate between infective avians and susceptible avians tends to a saturation level due to the protection measures taken by the poultry farmers or the crowding of the infected birds, which can be expressed as follows (Capasso and Serio [34]): where is the transmission coefficient so that measures the infection force of the infective birds and is the parameter which measures the inhibitory effort so that describes the saturation due to the protection measures of the poultry farmers or the crowding of infected birds when the number of infective birds increases. Similarly, as the number of the infective human individuals increases, the susceptible human population may tend to reduce the number of contacts with infective avian population per unit time due to the psychological effect, so we use a nonmonotone incidence function (Liu et al. [30], Xiao and Ruan [33]) to describe the transmission of the virus from infective birds to susceptible humans; that is, where measures the infection force of the disease and describes the psychological effect from the behavioral change of the susceptible humans when the number of infective individuals becomes larger.

(4) An infected avian remains in the state of disease and cannot recover, but an infected human may recover and a recovered human has permanent immunity.

Based on the above assumptions, we have the following SI-SIR avian influenza model:where () is the natural death rate of the avian population (the human population) respectively; () is the disease-related death rate of the infected avian (the infected human), respectively; is the natural recovery rate of the infective human. We assume that , are nonnegative constants and all the other parameters are positive.

Notice that the model proposed in Iwami et al. [7] is a special case of model (3) with being and if we ignore the virus mutation and spreading from human to human.

In the following, we will analyze the global dynamics of the above system. We first analyze the existence of the equilibria. It is clear that system (3) has a unique solution satisfying the initial conditions in which is a positively invariant set. The equilibrium of system (3) must satisfy the following equations:We consider two cases.

Case 1. If , then , , , and , which shows that system (3) always has a unique disease-free equilibrium given by with and .
We adopt the notation and method of van den Driessche and Watmough [35] to calculate the basic reproduction number. According to [35], system (3) can be written as whereThe infective compartments are and . Then,Hence, following [35], we obtain the basic reproduction number

Case 2. If , then from the second equation of (4). Substituting for in the first equation of system (4), we obtain that when . Hence, , , and satisfy the following equations:From the equations in (10), we have , , and satisfies the following equation:Obviously, , as , and which indicates that is strictly monotonically increasing. Hence, equation has only one positive root .

In summary, we have the following results.

Lemma 1. (i) If , system (3) has only a disease-free equilibrium given by . (ii) If , system (3) has two equilibria: a disease-free equilibrium and an endemic equilibrium , where , , , , , , and satisfies (11).

3. Analysis of the Avian-Only Submodel

Note that the first two equations are independent in system (3), so we study the dynamical behavior of the avian-only system firstly, which takes the following form:Based on Lemma 1, system (12) always has a disease-free equilibrium and a unique endemic equilibrium if , where the expressions of , , , and have been given in Section 2. Note that is also the basic reproduction of system (12). We will discuss the properties of system (12) in the positively invariant set .

Lemma 2. (i) The disease-free equilibrium is locally asymptotically stable if . (ii) The endemic equilibrium is locally asymptotically stable if .

Proof. (i) The characteristic equation of the Jacobian matrix around the equilibrium is Clearly, the eigenvalues of the above characteristic equation are , . Hence, the equilibrium is locally asymptotically stable if but unstable if .
(ii) If , system (12) has a unique positive equilibrium . Similarly, we can obtain the characteristic equation of the Jacobian matrix around the equilibrium as follows:whereSince , if , then , , so the two eigenvalues of characteristic equation have negative real parts. Hence, the positive equilibrium is locally asymptotically stable.

Lemma 3. (i) The disease-free equilibrium is globally asymptotically stable if . (ii) The endemic equilibrium is globally asymptotically stable if .

Proof. (i) If , we can choose a Liapunov function Then, Since , according to LaSalle’s invariance principle [36], is globally asymptotically stable.
(ii) In order to prove the global stability of the endemic equilibrium in , we use Bendixson-Dulac criteria [37, 38]. For the sake of convenience, let In order to use the Bendixson-Dulac criteria, we choose a positive smooth function . Then, , , and are continuously differentiable functions in the region , and It is obvious that the sign of the above expression is negative, which shows the nonexistence of any closed orbit. Therefore, is globally asymptotically stable if .

4. Analysis of the Avian-Human Influenza Epidemic Model

Since the first four equations of system (3) are independent of the variable , we only need to analyze the dynamical behavior of the following equivalent system:We will discuss the dynamical behavior of system (20) in the positive invariant set .

4.1. Local Stability

According to Lemma 1, system (20) always has a disease-free equilibrium given by ; if , system (20) also has a unique endemic equilibrium , where , , , , , , and have been given in Section 2. And is also the basic reproduction of system (20).

Similarly, we have the following properties.

Lemma 4. (i) The disease-free equilibrium is locally asymptotically stable for positive trajectories if but unstable if . (ii) The endemic equilibrium is locally asymptotically stable for positive trajectories if .

Proof. The Jacobian matrix of system (20) around an arbitrary equilibrium is
(i) If , the eigenvalues are , , , and .
If , all the eigenvalues are negative; hence, the equilibrium is locally asymptotically stable. But if , due to the eigenvalue , then the equilibrium is unstable.
(ii) If , system (20) exists a unique endemic equilibrium . The characteristic equation of the Jacobian matrix at the endemic equilibrium iswhere and satisfy (15) andAccording to Lemma 2, the eigenvalues of equation have negative real parts, so the stability of equilibrium is decided by the eigenvalues of equation .
Since , when , then , . Hence, both eigenvalues of equation have negative real parts. Therefore, the endemic equilibrium is locally asymptotically stable.

4.2. Global Stability

Theorem 5. (i) If , then the disease-free equilibrium of system (20) is globally asymptotically stable. (ii) If , the endemic equilibrium of system (20) is globally asymptotically stable.

Proof. (i) According to Lemma 3, the disease-free equilibrium of the avian-only system (12) is globally asymptotically stable if . To prove the global stability of , we only need to consider system (20) with the avian components already at the disease-free steady state, given by Clearly, we can obtain that , if . Hence, the disease-free equilibrium is globally asymptotically stable if .
(ii) Similarly, by Lemma 3, the endemic equilibrium of subsystem (12) is globally asymptotically stable in the if . To prove the global stability of the equilibrium in the region , we only need to consider system (20) with the avian components already at the endemic steady state, given byWe can easily deduce that system (25) has a unique positive equilibrium which is locally asymptotically stable. To prove the global stability of the equilibrium in the region , we adopt the Bendixson-Dulac criteria [37, 38].
For convenience, we note that In order to use the Bendixson-Dulac criteria, we choose a positive smooth function and then , , and are continuously differentiable functions in the region and According to the Bendixson-Dulac Theorem [37, 38], there is no closed orbit in the region . Therefore, is globally asymptotically stable. Hence, the endemic equilibrium is globally asymptotically stable if .

Corollary 6. (i) The disease-free equilibrium of system (3) is globally asymptotically stable for positive trajectories if but unstable if . (ii) The unique endemic equilibrium of system (3) is globally asymptotically stable if .

Remark 7. The basic reproduction number is a threshold value which determines whether the avian influenza disappears or not and parameter , the transmission rate from infective birds to susceptible birds, is a key parameter which influences the basic reproduction number .

Remark 8. Although our results are similar to those in Iwami et al. [7] when the virus mutation and spreading from human to human are not considered, parameters and influence the number of infected humans dramatically. We will show this point in the next section via numerical simulations.

5. Numerical Simulations and Sensitivity Analysis

In this section, for the sake of convenience and simplicity, we always keep some parameters fixed as follows: , , , , , , and . Generally speaking, avian influenza mainly outbreaks in a specific location. We estimate that the number of susceptible avian population is between 100000 and 1000000 and the number of infective avian population is between 0 and 100 every day in the region. Then, we choose initial values as , , , , , , , , .

5.1. Numerical Simulations

It should be noted that the transmission rate from infective birds to susceptible birds is a key parameter which influences the size of the basic reproduction number , the parameter describes the saturation effect within the avian population, and represents the psychological effect of the general public toward the infective individuals within the human population. Moreover, is also an important parameter which influences the number of infective human individuals. In this subsection, we investigate the influence of parameters , , , and on the number of infected humans by performing some numerical simulations.

Firstly, we study the influence of parameters , on the number of infective individuals. When parameters , , and are fixed as the above, the threshold value , so that . If , the disease will die out, and the solution is asymptotically stable and converges to the disease-free state value. If , the disease is prevalent and the solution is asymptotically stable and converges to the endemic state value (see Figure 1(a)). When we fix , for instance, , then . We can see that the solution is approaching the endemic state value and increases with the increasing of (see Figure 1(b)).

(a)

(b)

Secondly, we investigate the influence of parameters and on the number of infected humans. If parameter , then , the disease is endemic, and the solution is asymptotically stable and converges to the endemic state value.

For comparison, when parameters , are chosen as , and , , respectively, we observe that subtle changes on parameters decrease the number of infected humans dramatically (see Figure 2). When parameter is fixed and the other parameter takes different values, we can see that the number of infected humans decreases obviously when parameter increases (see Figure 3(a)). When parameter is fixed and parameter takes different values, we have a similar result (see Figure 3(b)).

(a)

(b)

(a)

(b)

If parameter , then , the disease disappears, and the solution is asymptotically stable and converges to the disease-free state value.

Similarly, when parameter is fixed, we can observe that the number of infected humans decreases with the increase of parameter (see Figure 4(a)). When we fix parameter and let parameter change, we obtain a similar result (see Figure 4(b)). When both parameters and change simultaneously, we can see that the number of infected humans decreases obviously with the increase of parameters and (see Figure 5).

(a)

(b)

5.2. Sensitivity Analysis

In order to find better control strategies for avian influenza infections, we firstly investigate which parameters can reduce the basic reproduction number . From Figure 6, we can observe that decreases if or decreases or if parameters or increase. If the contact rate is sufficiently small, then avian influenza could be eliminated even if or . But it is difficult to control . is hard to control too because the disease-related death rate is a relative fixed constant. Even if is not small, could be less than 1 as long as is large enough and is small enough. It is feasible if we shorten the lifetime of domestic birds or reduce all new recruitments and newborns of domestic birds. The optimal control strategy will be a combination of shortening lifetime of domestic birds and reducing all new recruitments and newborns of domestic birds.

(a)

(b)

(c)

Secondly, we investigate the influence of parameters , , and on the endemic state value of the infective human individuals . We can observe that decreases with the increase of parameters or and will tend to 0 when parameters or increase to infinity from Figures 7(a) and 7(b). We can also observe that increases with the increase of parameter from Figure 7(c). It suggests that we can increase parameters and or reduce parameter to control the disease to a lower level. Similarly, it is difficult to control parameter . However, we can control parameters and . For example, we can enhance the intensity of media coverage to enhance the psychological effect of the human population; we can reduce through susceptible human individuals by avoiding as far as possible the contact with the infective avian population. From Figure 7(c), we can also see that if , avian influenza will be eliminated which implies that closing the retail live-poultry markets is an effective control measure.

(a)

(b)

(c)

6. Discussion

In this paper, to study the transmission dynamics of avian influenza from birds to humans, we constructed a nonlinear ordinary differential equation model considering the saturation effect within the avian population and the psychological effect of the general public toward the outbreaks of avian influenza. We obtained a threshold value for the prevalence of avian influenza and discussed the local and global asymptotical stability of each equilibrium of the nonlinear system. Our results indicate that the asymptotic dynamics of the model are completely determined by the threshold value : the disease-free equilibrium exists and is globally asymptotically stable if ; the disease-free equilibrium becomes unstable and the endemic equilibrium exists and is globally asymptotically stable if . In other words, the avian influenza disappears if but is prevalent and becomes endemic if . Our numerical simulations (see Figures 1–5) support the theoretical results.

Our results demonstrate that transmission dynamics of the avian influenza are greatly determined by the infection force of the disease in birds. It should be noted that the saturation effect within the avian population , the psychological effect within the human population , and the transmission coefficient from infective birds to susceptible humans do not affect the stability of the equilibria and thus the transmission and outbreaks of the disease, since infected humans do not spread the virus any further. However, the theoretical analyses and numerical simulations indicate that increasing parameters and and the reducing parameter could reduce the number of the infected humans and may help to control the disease.

Numerical simulation suggests that the optimal control strategy should be a combination of shortening the lifetime of domestic birds, reducing all new recruitments and newborns of domestic birds, and reducing contacts between susceptible birds and infective birds. If there is an outbreak of avian influenza, the effect control measures should be a combination of enhancing the intensity of media coverage and avoiding contact with the infective avian population. Moreover, closing the retail live-poultry markets is the fastest control measure in controlling the disease.

The roles of wild birds and domestic birds in the transmission of the H5N1 avian influenza are different and mathematical models have been proposed to include both types of birds (Bourouiba et al. [15], Gourley et al. [16], Lucchetti et al. [8], and Tuncer and Martcheva [17]). It will be very interesting to include both wild birds and domestic birds in modeling the bird-to-human transmission of the H7N9 avian influenza. We leave this for future consideration.

Conflict of Interests

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

Acknowledgments

The authors are very grateful to the reviewers for their constructive suggestions and comments which helped us to improve the paper significantly. This work was partially supported by the National Natural Science Foundation of China (nos. 11371161 and 11228104) and the National Science Foundation (DMS-1412454).

References

Centers for Disease Control and Prevention (CDC), “Types of influenza virus,” http://www.cdc.gov/flu/about/viruses/types.htm.
View at: Google Scholar
D. J. Alexander, “An overview of the epidemiology of avian influenza,” Vaccine, vol. 25, no. 30, pp. 5637–5644, 2007.
View at: Publisher Site | Google Scholar
Centers for Disease Control and Prevention (CDC), “Information on Avian Influenza,” http://www.cdc.gov/flu/avianflu/index.htm.
View at: Google Scholar
Centers for Disease Control and Prevention (CDC), “Influenza Type A Viruses,” http://www.cdc.gov/flu/avianflu/influenza-a-virus-subtypes.htm.
View at: Google Scholar
R. M. Anderson and R. M. May, Infectious Diseases of Humans: Dynamics and Control, Oxford University Press, Oxford, UK, 1991.
M. J. Keeling and P. Rohani, Modeling Infectious Diseases in Humans and Animals, Princeton University Press, Princeton, NJ, USA, 2008.
View at: MathSciNet
S. Iwami, Y. Takeuchi, and X. Liu, “Avian-human influenza epidemic model,” Mathematical Biosciences, vol. 207, no. 1, pp. 1–25, 2007.
View at: Publisher Site | Google Scholar | MathSciNet
J. Lucchetti, M. Roy, and M. Martcheva, “An avian influenza model and its fit to human avian influenza cases,” in Advances in Disease Epidemiology, J. M. Tchuenche and Z. Mukandavire, Eds., pp. 1–30, Nova Science Publishers, New York, NY, USA, 2009.
View at: Google Scholar
S. Iwami, Y. Takeuchi, and X. Liu, “Avian flu pandemic: can we prevent it?” Journal of Theoretical Biology, vol. 257, no. 1, pp. 181–190, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
E. Jung, S. Iwami, Y. Takeuchi, and T.-C. Jo, “Optimal control strategy for prevention of avian influenza pandemic,” Journal of Theoretical Biology, vol. 260, no. 2, pp. 220–229, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
S. Iwami, Y. Takeuchi, X. Liu, and S. Nakaoka, “A geographical spread of vaccine-resistance in avian influenza epidemics,” Journal of Theoretical Biology, vol. 259, no. 2, pp. 219–228, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
A. B. Gumel, “Global dynamics of a two-strain avian influenza model,” International Journal of Computer Mathematics, vol. 86, no. 1, pp. 85–108, 2009.
View at: Publisher Site | Google Scholar | MathSciNet
F. B. Agusto, “Optimal isolation control strategies and cost-effectiveness analysis of a two-strain avian influenza model,” BioSystems, vol. 113, no. 3, pp. 155–164, 2013.
View at: Publisher Site | Google Scholar
X. Ma and W. Wang, “A discrete model of avian influenza with seasonal reproduction and transmission,” Journal of Biological Dynamics, vol. 4, no. 3, pp. 296–314, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
L. Bourouiba, S. A. Gourley, R. Liu, and J. Wu, “The interaction of migratory birds and domestic poultry and its role in sustaining avian influenza,” SIAM Journal on Applied Mathematics, vol. 71, no. 2, pp. 487–516, 2011.
View at: Publisher Site | Google Scholar | MathSciNet
S. A. Gourley, R. Liu, and J. Wu, “Spatiotemporal distributions of migratory birds: patchy models with delay,” SIAM Journal on Applied Dynamical Systems, vol. 9, no. 2, pp. 589–610, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
N. Tuncer and M. Martcheva, “Modeling seasonality in avian influenza H5N1,” Journal of Biological Systems, vol. 21, no. 4, Article ID 1340004, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
J. Zhang, Z. Jin, G.-Q. Sun, X.-D. Sun, Y.-M. Wang, and B. Huang, “Determination of original infection source of H7N9 avian influenza by dynamical model,” Scientific Reports, vol. 4, article 4846, 2014.
View at: Publisher Site | Google Scholar
Y. Xiao, X. Sun, S. Tang, and J. Wu, “Transmission potential of the novel avian influenza A(H7N9) infection in mainland China,” Journal of Theoretical Biology, vol. 352, pp. 1–5, 2014.
View at: Publisher Site | Google Scholar | MathSciNet
Y.-H. Hsieh, J. Wu, J. Fang, Y. Yang, J. Lou, and G.-Q. Sun, “Quantification of bird-to-bird and bird-to-human infections during 2013 novel H7N9 avian influenza outbreak in China,” PLoS ONE, vol. 9, no. 12, Article ID e111834, 2014.
View at: Publisher Site | Google Scholar
C. T. Bauch and A. P. Galvani, “Social factors in epidemiology,” Science, vol. 342, no. 6154, pp. 47–49, 2013.
View at: Publisher Site | Google Scholar
N. Ferguson, “Capturing human behaviour,” Nature, vol. 446, article 733, 2007.
View at: Publisher Site | Google Scholar
S. Funk, E. Gilad, C. Watkins, and V. A. A. Jansen, “The spread of awareness and its impact on epidemic outbreaks,” Proceedings of the National Academy of Sciences of the United States of America, vol. 106, no. 16, pp. 6872–6877, 2009.
View at: Publisher Site | Google Scholar
S. Funk, M. Salathé, and V. A. A. Jansen, “Modelling the influence of human behaviour on the spread of infectious diseases: a review,” Journal of the Royal Society Interface, vol. 7, no. 50, pp. 1247–1256, 2010.
View at: Publisher Site | Google Scholar
S. J. Olsen, Y. Laosiritaworn, S. Pattanasin, P. Prapasiri, and S. F. Dowell, “Poultry-handling practices during avian influenza outbreak, Thailand,” Emerging Infectious Diseases, vol. 11, no. 10, pp. 1601–1603, 2005.
View at: Publisher Site | Google Scholar
N. Xiang, Y. Shi, J. Wu et al., “Knowledge, attitudes and practices (KAP) relating to avian influenza in urban and rural areas of China,” BMC Infectious Diseases, vol. 10, article 34, 2010.
View at: Publisher Site | Google Scholar
L. Wang, B. J. Cowling1, and P. Wu, “Human exposure to live poultry and psychological and behavioral responses to influenza A(H7N9), China,” Emerging Infectious Diseases journal, vol. 20, pp. 1296–1305, 2014.
View at: Google Scholar
R. Liu, J. Wu, and H. Zhu, “Media/psychological impact on multiple outbreaks of emerging infectious diseases,” Computational and Mathematical Methods in Medicine, vol. 8, no. 3, pp. 153–164, 2007.
View at: Publisher Site | Google Scholar | MathSciNet
S. Funk, S. Bansal, C. T. Bauch et al., “Nine challenges in incorporating the dynamics of behaviour in infectious diseases models,” Epidemics, 2014.
View at: Publisher Site | Google Scholar
W. M. Liu, S. A. Levin, and Y. Iwasa, “Influence of nonlinear incidence rates upon the behavior of SIRS epidemiological models,” Journal of Mathematical Biology, vol. 23, no. 2, pp. 187–204, 1986.
View at: Publisher Site | Google Scholar | MathSciNet
S. Ruan and W. Wang, “Dynamical behavior of an epidemic model with a nonlinear incidence rate,” Journal of Differential Equations, vol. 188, no. 1, pp. 135–163, 2003.
View at: Publisher Site | Google Scholar | MathSciNet
Y. Tang, D. Huang, S. Ruan, and W. Zhang, “Coexistence of limit cycles and homoclinic loops in a SIRS model with a nonlinear incidence rate,” SIAM Journal on Applied Mathematics, vol. 69, no. 2, pp. 621–639, 2008.
View at: Publisher Site | Google Scholar | MathSciNet
D. Xiao and S. Ruan, “Global analysis of an epidemic model with nonmonotone incidence rate,” Mathematical Biosciences, vol. 208, no. 2, pp. 419–429, 2007.
View at: Publisher Site | Google Scholar | MathSciNet
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 | 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 | MathSciNet
J. K. Hale, Ordinary Differential Equations, Wiley-Interscience, New York, NY, USA, 1969.
View at: MathSciNet
J. Zhang and B. Feng, The Geometric Theory and Bifurcation Problems of Ordinary Differential Equations, Peking University Press, Beijing, China, 1997, (Chinese).
Z. Zhang, T. Ding, W. Huang, and Z. Dong, Qualitative Theory of Differential Equations, Science Press, Beijing, China, 1985, (Chinese), English Edition, Translations of Mathematical Monographs vol. 101, American Mathematical Society, Providence, RI, USA, 1992.

Copyright

Copyright © 2015 Sanhong Liu 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

2772

Downloads

2065

Citations