ABSTRACT
In general, media coverage would not be implemented unless the number of infected cases reaches some critical number. To reflect this feature, we incorporate the media effect and a critical number of infected cases into the disease transmission rate and consider an susceptible-infected-susceptible epidemic model with logistic growth. Our model analysis shows that early media alert and strong media effects are preferable to decrease the numbers of infected cases at endemic equilibria. Furthermore, we noticed that the model may have up to three endemic equilibria and bi-stability can occur in a threshold interval for the critical number. Note that the interval depends on parameters for the focal disease and the media effect. It is possible to roughly estimate the interval for re-emerging diseases in a given region. Therefore, the result could be useful to health policymakers. Global stability is also obtained when the model admits a unique endemic equilibrium.
1. Introduction
Media has been utilized as a disease control measure, especially for epidemics associated with emerging and re-emerging infectious diseases [Citation19] such as HIV/AIDS, SARS, H1N1, Ebola virus disease (EVD), Middle East Respiratory Syndrome. During the outbreak of the influenza A (H1N1) in 2009, mass media was extensively used by the Centers for Disease Control and Prevention of United States and WHO to keep the public aware of information related to the pandemic [Citation6]. It is believed that media use contributed to the control of the pandemic. WHO also indicated that media played an important role in controlling the spread of H7N9 in China in 2013 [Citation31]. Media does not only alert the general public on the hazard from the infectious diseases but also informs the public of the requisite preventive measures like wearing protective masks [Citation25], vaccination, voluntary quarantine, avoidance of congregated places, etc. Therefore, the extensive use of media may bring in changes in public behaviour and reduce the frequency and probability of contacts with infected individuals so that the severity of a disease outbreak would be diminished [Citation4, Citation9, Citation10, Citation13, Citation14, Citation21, Citation24].
In order to study the impact of media-like control measures on disease transmission dynamics, several types of media function forms have been proposed to describe reduced disease transmission rates due to media use and compartmental models with these rates have been analysed (e.g. [Citation9, Citation10, Citation13–17, Citation23, Citation24]). The deduction in the transmission rate was described by the form of with the parameter m>0 reflecting how strongly media coverage can affect contact infection [Citation9]. With the rate, the analysis of a susceptible–exposed–infected model (SEI) shows that the model may exhibit periodic oscillations for weak media effects while it may have three endemic equilibria for strong media effects [Citation9]. The form of was also used as the transmission rate with the deduction due to media use [Citation10, Citation13, Citation24]. A threshold dynamics was obtained for an SIS epidemic model. It is also shown that media coverage can lower infection and delay the arrival of the infection peak [Citation10, Citation13]. However, the study of a susceptible–vaccinated–infected–recovered epidemic model with a vaccinated class indicates that media effects could be complicated and simplified understandings may even make the disease worse due to possible public panic [Citation24]. The third function type with psychological/media effects is of the form , identified by Collinson and Heffernan (see [Citation8, Citation33] and the references therein). Using a simple susceptible-exposed-infected-recovered model, Collinson and Heffernan [Citation8] found that important measurements of an epidemic outbreak (such as peak magnitude of infection, peak time of infection peak and end of the outbreak) depend on the chosen media function. Their sensitivity analysis also showed such dependence for the sensitivities of model parameters. This makes it difficult to identify effective disease control strategy and calls for more study on the effects of mass media on disease transmission dynamics.
Theoretical studies usually assume that media coverage affects disease transmission during the whole time period of the disease spread (e.g. [Citation9, Citation10, Citation14, Citation24]). In the reality, media coverage generally does not occur in the beginning stage of a disease spread. For instance, the early suspected cases of EVD died in December 2013 while the first notification by WHO on the outbreak [Citation32] was not issued until 21 March 2014. In general, media delivers alerts and timely reports infected cases only when certain number of infected cases is reached [Citation28, Citation34, Citation38]. To include such feature of media/psychological effect, Xiao et al. [Citation34] introduced a critical number of infected cases and proposed a piecewise and discontinuous control function [Citation7, Citation26] for disease control strategy (sliding mode control). Wang and Xiao [Citation28] further constructed a Filippov SIR epidemic model to describe media effects using the following transmission rate (1) (1) Their analysis shows that the model system can stabilize at either one of the equilibria for the resulted subsystems or the new endemic state induced by the on-off media effect, depending on the critical level . They also demonstrated that proper combinations of critical levels and control intensities can lead to the desired case number. The transmission rate (Equation1(1) (1) ) was then generalized in a time-dependent way to study an influenza outbreak in Shannxi, China [Citation35].
It was demonstrated that compartmental models can exhibit distinct dynamics, depending on the chosen incidence rate (e.g. [Citation20, Citation33]). To explore the media effect on disease transmission dynamics we here propose a new transmission rate. Following the idea of the critical number of infected cases , we consider the following non-smooth but continuous transmission rate: (2) (2) where represents the intensity of the media effect on contact infection. If p=0, is equal to the background transmission rate β, implying that media coverage does not occur. With this rate, we shall consider an SIS endemic model.
Classical compartmental models with media effects assume either a constant size of the total population or constant recruitment rate for the susceptible class. The assumption of varying total populations may be more reasonable for a relatively long-lasting disease or for a disease with high mortality rates. In fact, varying total populations were discussed before (e.g. [Citation1, Citation3, Citation5, Citation9, Citation11, Citation22, Citation29, Citation36]). Here, we assume that the population of a community follows the logistic growth. For the sake of mathematical simplicity, we assume that newborns directly enter into the susceptible class and infected persons do not contribute to births and deaths in the susceptible class. Following works in [Citation2, Citation9, Citation28], our SIS model reads (3) (3) where r is the intrinsic growth rate of the susceptible population, a denotes the carrying capacity of the community in the absence of infection, d is natural death rate, γ represents the recovered rate, and ε is the disease-induced death rate. The analysis of model (Equation3(3) (3) ) with the transmission rate (Equation2(2) (2) ) shows that there exists a threshold interval Γ for the critical number in which the model may be stabilized at one of two stable equilibria with different levels of infected cases. This implies that the policymaker may have to choose the critical number according to the focal disease in order to minimize infected cases and also avoid unnecessary public panic. The global stability was also obtained for not in the threshold interval. In the following, the analysis of existence of equilibria is presented in Section 2 while the local and global stabilities are given in Sections 3 and 4, respectively. A discussion section comes to the end of the work.
2. Existence of equilibria
For our convenience, denote . Model (Equation3(3) (3) ) with the transmission rate (Equation2(2) (2) ) can be decomposed into two sub-systems (4) (4) and (5) (5) The origin and the disease-free equilibrium always exist. The basic reproductive number can be easily calculated, given by from which one can see that media coverage does not change the basic reproduction number (e.g. [Citation9, Citation10, Citation13, Citation15–17, Citation23, Citation24]). To obtain the existence of endemic equilibria, denote and consider two cases: and , separately.
In the case of , the sub-system (Equation4(4) (4) ) has a unique positive equilibrium if and only if , that is, , where .
In the case of , the I component of a positive equilibrium for the sub-system (Equation5(5) (5) ) satisfies the following equation (6) (6) From , positive solutions to Equation (Equation6(6) (6) ) satisfy (7) (7) That is, the positive solutions of Equation (Equation6(6) (6) ) must be in the interval . Also, . Therefore, we must have due to . Next, we discuss the existence of positive solutions to Equation (Equation6(6) (6) ) in in two cases, and p>1.
Case I. . In this case, is strictly increasing. In fact, the derivative of is given by (8) (8) and we can show that for . If , one clearly derives . If , we can calculate the zero point of as (9) (9) which, obviously, is the unique maximum point of . Since (10) (10) we have for , and hence on . To sum up, we always have if . Meanwhile, note that (11) (11) may be positive, negative or zero. If , which is equivalent to , then Equation (Equation6(6) (6) ) has no positive real root in . If , which is equivalent to , then Equation (Equation6(6) (6) ) admits a unique positive root, denoted by , satisfying . That is, in the case of , the sub-system (Equation5(5) (5) ) has a unique positive equilibrium, denoted by , if and only if , where and .
The following lemma is needed to discuss the case of p>1.
Lemma 2.1:
Let then always holds, where
Proof:
Consider the concave function . Assign Note that , and is a strictly concave function. This implies that , and we have The monotonicity of leads to . The proof is completed.
Case II. p>1. By Lemma 2.1, holds. It follow from the formula (Equation10(10) (10) ) that . Note that the derivative given by Equation (Equation8(8) (8) ) is positive. is the unique minimum point of . Clearly, . Next, we want to compare the sizes of and and determine the signs of and . The following six cases are considered. Let us denote
Case 1. and . It follows from that Since , we have , and . Furthermore, from Equation (Equation9(9) (9) ) one can get . Thus . Direct calculation shows that holds. If , we have . Thus, Equation (Equation6(6) (6) ) has a unique positive root , satisfying . If , the conclusion is still true. Accordingly, the sub-system (Equation5(5) (5) ) has a unique positive equilibrium .
Case 2. and . It follows from that that is, . Thus, we have from Equation (Equation9(9) (9) ) and because of . Moreover, one can show that . Therefore, when , Equation (Equation6(6) (6) ) has two different positive roots in , denoted by , , where , . That is, the sub-system (Equation5(5) (5) ) admits two different positive equilibria , where , i=2,3.
Case 3. and . Similar to the argument in Case 2, it follows that , and hence we have and . Note that the equivalent relationship . We must have . This suggests that Equation (Equation6(6) (6) ) only has one positive solution . Accordingly, for the sub-system (Equation5(5) (5) ), there exists exactly one positive equilibrium , where .
Case 4. and . Recall that and is the unique minimum point of in . One immediately deduces that Equation (Equation6(6) (6) ) has no positive root in no matter or . Hence, the sub-system (Equation5(5) (5) ) has no positive equilibrium.
Case 5. and . From , one deduces that and . Hence, Equation (Equation6(6) (6) ) only has one positive solution no matter or , where . That is, the sub-system (Equation5(5) (5) ) has only one positive equilibrium .
Case 6. and . Clearly, one can have from and . Hence, . Notice that is an increasing function on the interval , and . Equation (Equation6(6) (6) ) has no positive root in the interval . Namely, the sub-system (Equation5(5) (5) ) has no positive equilibrium.
Now, one summarizes the existence of the equilibria of model (Equation3(3) (3) ) as follows:
Theorem 2.2:
Model (Equation3(3) (3) ) always admits an equilibrium and a disease-free equilibrium . If i.e. then model (Equation3(3) (3) ) has no endemic equilibrium. If i.e. we have the following conclusions.
Assume that .
If then is a unique endemic equilibrium;
If then is a unique endemic equilibrium.
Assume that p>1 and .
If is a unique endemic equilibrium;
If there are two endemic equilibrium, and
If there exist three endemic equilibria, and
If both and exist;
If is a unique endemic equilibrium.
Assume that p>1 and .
If is a unique endemic equilibrium;
If is a unique endemic equilibrium.
Remark 1:
From (ii) of (1), (v) of (2) and (iii) of (3) in Theorem 2.2 one can deduce that is a unique endemic equilibrium if .
By Theorem 2.2, if , model (3) only has one endemic equilibrium, either (in the case of ), or (in the case of ). The existence of positive equilibria for the model in the case of p>1 is illustrated in Figure . Here, we define the following different curves and regions for parameters a and as follows: , , , , , (where ), and .
Note that the condition with in Theorem 2.2 is equivalent to , where In the case of and p>1, multiple endemic equilibria exist if the critical number is in the threshold interval (12) (12) The existence region in Figure is given by the union of , and .
3. Local stability of equilibria
This section focuses on the local stability of equilibria. Corresponding to the equilibria , and , the Jacobian matrix of the sub-system (Equation4(4) (4) ) reads (13) (13) At the determinant . Thus is a saddle point. At it can be shown that (14) (14) If , then , which means that is a saddle point. If , then and . Hence is locally asymptotically stable. Since (15) (15) is a stable node or critical node or degenerate node. If , it follows from , that is a saddle-node (see Theorem 7.1 in [Citation37] or Theorem 2.11.1 in [Citation18]).
The Jacobian matrix at can be written as (16) (16) And hence, (17) (17) Recall that holds when exists. We have , . Thus, is asymptotically stable.
In the following, we discuss the stability of the equilibria , . The Jacobian matrix of the sub-system (Equation5(5) (5) ) at is given by (18) (18) By Equation (Equation6(6) (6) ), we have (19) (19) Hence, (20) (20) If , then always holds. If , from , it follows that (21) (21) Therefore, we always have .
From formula (Equation19(19) (19) ), one can have (22) (22) Let us first determine the sign of at the equilibrium . If , both terms in the bracket of formula (Equation22(22) (22) ) are positive and hence . If , from Case I in Section 2 we know . Therefore, (23) (23)
If p>1, from the discussion of Cases 1 and 5 in Section 2 it follows that , implying that inequality (Equation23(23) (23) ) still holds. In summary, we have and . Therefore, is asymptotically stable. The stability of the equilibria , and can be determined similarly. At , one can have , implying that is a degenerate node. At , it follows from that . Thus, is a saddle point. At , we have (since ). Hence, is stable.
To sum up, one obtains the following result.
Theorem 3.1:
For model (Equation3(3) (3) ), the local stability of the equilibria is stated as follows:
The origin is a saddle point;
If i.e. then the disease-free equilibrium is locally stable. If then is a saddle-node. If i.e. then is a saddle point;
are locally asymptotically stable whenever they exist. is a saddle point, and is a degenerate node.
To illustrate the existence and stabilities of multiple equilibria for model (Equation3(3) (3) ), we utilize some parameter values estimated from the influenza A (H1N1). Set the recovered rate year−1 [Citation12], the disease-induced death rate year−1 [Citation25], and fix the natural death rate year−1, and . The basic reproductive number of influenza A was estimated as 1.5−3.1 [Citation30]. For ease of demonstration, we naively set p=3, and a=600, and hence the basic reproductive number is equal to 3.0. We shall consider two examples. Both examples show the occurrence of bi-stability, in which solutions may converge to one of two stable equilibria, depending on initial conditions.
Example 3.1:
Set . We can calculate and . Therefore, , and . In this case (see Theorem 2.2(2)), model (Equation3(3) (3) ) admits two stable endemic equilibria and , that is, bi-stability occurs (see the line in Figure ). Some solutions to the model are illustrated in Figure . Note that is equivalent to . The endemic equilibrium lies on the horizontal line .
Example 3.2:
Set . In this case, , and . The conditions of (iii) of (2) in Theorem 2.2 are satisfied. Hence, model (Equation3(3) (3) ) has three positive equilibria , and . Figure illustrates some numerical solutions to the model. The stable manifolds of the saddle split the phase plane into two regions. In the lower region, solutions approach to while in the upper region solutions approach to (see Figure ).
4. Global stability analysis
In this section, we study the global stability of model (Equation3(3) (3) ). It can be shown that the state variables of model (Equation3(3) (3) ) remain non-negative for non-negative initial conditions. Consider the biologically feasible region Choosing the straight line S+I−Ab=0, similar to the proof of Corollary in [Citation36], we have the following two results.
Lemma 4.1:
The closed set Π is a positively invariant set for model (Equation3(3) (3) ).
Theorem 4.2:
is globally asymptotically stable if and unstable if .
To the best of our knowledge, by precluding the existence of a limit cycle we are able to prove global stability of the unique endemic equilibrium of model (Equation3(3) (3) ).
Theorem 4.3:
There exist no limit cycles for model (Equation3(3) (3) ).
Proof:
The following two steps are considered to achieve our conclusion.
Step 1. We shall prove that there are no limit cycles in the region below the line in the feasible region Π and in the region above the line . Denote these two regions by and .
Take the Dulac function . In the case of , by the transformation (24) (24) one can transfer sub-system (Equation4(4) (4) ) into (25) (25) and . Let , be the right-hand side functions of (Equation25(25) (25) ). Then (26) (26) Therefore, there are no limit cycles in the region below the line .
In the case of , we can extend subsystem (Equation5(5) (5) ) to the case of because of the continuity of the transmission function . That is, subsystem (Equation5(5) (5) ) can be rewritten into (27) (27) Consider two cases p=1 and . If p=1, one can set x=S and y=I. With and being the right sides of system (Equation27(27) (27) ), direct calculations show that (28) (28) If , using the transformation (29) (29) one obtains (30) (30) Therefore, (31) (31) Consequently, for all p>0, one can have (32) (32) Hence, there is no limit cycle in , the region above the line . We should point out that inequalities (Equation26(26) (26) ) and (Equation32(32) (32) ) hold for .
Step 2. We are now ready to show that model (Equation3(3) (3) ) has no limit cycle crossing the line . The idea is similar to that in [Citation27, Citation28, Citation34].
Assume that Γ is a limit cycle across the line . Let be the part of the cycle below of the cycle Γ, and be the part above , with the direction designated in Figure . Let both and include two intersection points , of Γ with the line . The region enclosed by and the segment is denoted by , and the region enclosed by and the segment is denoted by .
Let us choose two the directed-paths and , as shown in Figure . It can be seen that (33) (33) Meanwhile, it is easy to obtain that (34) (34) Hence, from Green's Theorem it follows that (35) (35) Obviously, one can have (36) (36) From Step 1, however, we know that there is no limit cycle in the regions or , and (37) (37) A contradiction to Equation (Equation35(35) (35) ). Therefore, there are no limit cycles crossing the line .
To sum up, model (Equation3(3) (3) ) has no limit cycles.
From Theorems 2.2 and 4.3, we immediately have
Corollary 4.4:
If model (Equation3(3) (3) ) admits a unique endemic equilibrium either of or then it is globally asymptotically stable.
We can illustrate the global stability by numerical simulations. Set and the remaining parameters take the same values as in Figure . Then . From Remark 1 and Corollary 4.4, the unique endemic equilibrium (see the region in Figure ) of model (Equation3(3) (3) ) is globally asymptotically stable. Figure shows that the number I of infections is stabilized at the level . In the figure, it is also shown that the number I of infected cases is stabilized at decreased levels as either the intensity p of the media effect increases or the critical number decreases. That is, stronger media effects and/or lower critical numbers lead to a decreased number of infections at the endemic equilibrium . In fact, such decreasing effects are also true for all endemic equilibria except .
5. Discussion
During a disease spread, only when the number of infected cases and/or the severity of infection are high enough to draw the attention of media and public health organizations, alerts are issued and all related information is brought to the public through media coverage. The information gradually changes public behaviour, which reduces the chance of potential contact infection and eventually helps to curb the disease spread. Fast and dramatic changes in public behaviour can also occur subject to some intensive control measures such as closing school and distancing certain groups of persons. Previous studies have introduced a critical number to be the level for media and health organizations to take action. Functions with jump-discontinuity at were used to describe the changing transmission rate due to the media effect [Citation28, Citation34, Citation35]. With the rate functions, media effects on disease outbreaks were studied through compartmental epidemic models. It turns out that proper use of media coverage can curb disease outbreaks [Citation28, Citation34]. In this paper, following the idea of the critical number , we proposed the non-smooth function (2) to describe the transmission rate which is continuous at . With this description, change in public behaviour is continuous. Using a susceptible–infected–susceptible model with a logistic growth in the susceptible class, we studied the media effect on the transmission dynamics of an infectious disease in a given region.
Our model analysis shows that without the media effect or with relatively weak effect (), the endemic equilibrium or (depending on the chosen ) is globally asymptotically stable. With relatively strong media effect p>1, the model may have up to three endemic equilibria for the chosen critical number in the threshold interval . Otherwise, the model admits a unique endemic equilibrium (which is globally asymptotically stable, see Figure ). In the former case, solutions to the model can converge to either one of two stable endemic equilibria (see Figures and ), depending on initial conditions. That is, bi-stability can occur. To avoid such uncertainty in practice, it is necessary to choose the critical number below the value . It is worthwhile to point out that the critical values , depend only on the basic reproduction number and the intensity p of the media effect and hence they are prescribed by the focal disease and the population in the given region. Therefore, it is possible for policymaker to roughly estimate the threshold interval for re-emerged diseases and then choose a reasonable critical number to initiate media coverage.
From the point of view of disease control, early media alerts (i.e. setting a small critical number for ) and strong media effects (i.e. p>1) are definitely preferable. Our analysis shows that if , equivalently , the endemic equilibrium is globally asymptotically stable. However, if , model solutions approach the unique equilibrium , causing that media coverage loses its impact on disease transmission. This should be avoided by health policymakers.
Early media alerts and strong effect can decrease the numbers of infected cases at endemic equilibria , . For example, as the critical number decreases, or the intensity p of the media effect increases, the number of infected cases can be stabilized at with decreased numbers of cases (see Figure ). Therefore, properly choosing the critical number and strengthening the media effect can reduce disease prevalence. The analysis of an SEI model also implies that media coverage can reduce the number of infected cases at endemic equilibria [Citation9].
The existence of multiple endemic equilibria was obtained in some previous studies (see e.g. [Citation9, Citation20]). Using the transmission rate with media parameter m, Cui et al. [Citation9] found that the model may exhibit periodic oscillations for sufficiently small media effects (small m) and may have multiple endemic equilibria for strong media effects (large m). Unfortunately, the stabilities of the equilibria are not available and hence the solution behaviour is unknown for strong media effects. Media/psychological effects may be also included in the incidence rate in a nonlinear and saturation way such as and simple compartmental models can exhibit complex dynamics. For instance, saddle-node bifurcation, Hopf bifurcation and homocyclic bifurcation can occur in a simple SIR epidemic model with the rate [Citation20]. Therefore, the impact of media coverage on a disease transmission dynamics could be complicated and simplified understandings may even make the disease worse (e.g. [Citation24]). It needs to be further studied from distinct aspects [Citation8].
In order to avoid the complexity of mathematical analysis we assumed that the logistic growth in the susceptible class depends only on the number of the susceptible, instead of the total population size. This is a significant simplification. As an end of the paper, we would like to point out that this simplification is unlikely to change our major qualitative results about the media effect, such as the existence of the threshold interval Γ and bi-stability.
Acknowledgments
The authors thank the editor and the anonymous reviewers for their constructive comments that help to improve the early version of this paper.
Disclosure statement
No potential conflict of interest was reported by the authors.
Additional information
Funding
References
- R.M. Anderson, H.C. Jackson, R.M. May, and A.D.M. Smith, Population dynamics of fox rabies in Europe, Nature 289 (1981), pp. 765–771. doi: 10.1038/289765a0
- F. Brauer, Models for the spread of universally fatal diseases, J. Math. Biol. 28 (1990), pp. 451–462. doi: 10.1007/BF00178328
- F. Brauer and C. Castillo-Chavez, Mathematical Models in Population Biology and Epidemiology, Springer, New York, 2001.
- S. Broder and R.C. Gallo, A pathogenic retrovirus (HTLV-III) linked to AIDS, New Engl. J. Med. 311 (1984), pp. 1292–1297. doi: 10.1056/NEJM198411153112006
- S. Busenberg and P. van den Driessche, Analysis of a disease transmission model in a population with varying size, J. Math. Biol. 28 (1990), pp. 257–270. doi: 10.1007/BF00178776
- CDC, H1N1 Flu, Center for Disease Control and Prevention Website. Available at http://www.cdc.gov/h1n1flu/.
- N.S. Chong and R.J. Smith, Modeling avian influenza using Filippov systems to determine culling of infected birds and quarantine, Nonlinear Anal. Real World Appl. 24 (2015), pp. 196–218. doi: 10.1016/j.nonrwa.2015.02.007
- S. Collinson and J.M. Heffernan, Modelling the effects of media during an influenza epidemic, BMC Public Health 14 (2014), pp. 1–10. doi: 10.1186/1471-2458-14-376
- J. Cui, Y. Sun, and H. Zhu, The impact of media on the control of infectious diseases, J. Dynam. Differential Equations 20 (2008), pp. 31–53. doi: 10.1007/s10884-007-9075-0
- J-A. Cui, X. Tao, and H. Zhu, An SIS infection model incorporating media coverage, Rocky Mountain J. Math. 38 (2008), pp. 1323–1334. doi: 10.1216/RMJ-2008-38-5-1323
- H.W. Hethcote, A thousand and one epidemic models, in Frontiers in Mathematical Biology, S.A. Levin, ed., Lecture Notes in Biomathematics, Vol. 100, Springer, Berlin, 1994, pp. 504–515.
- S. Leekha, N.L. Zitterkopf, M.J. Espy, T.F. Smith, R.L. Thompson, and P. Sampathkumar, Duration of influenza: A virus shedding in hospitalized patients and implications for infection control, Infect. Control Hosp. Epidemiol. 28 (2007), pp. 1071–1076. doi: 10.1086/520101
- Y. Liu and J-A. Cui, The impact of media coverage on the dynamics of infectious disease, Int. J. Biomath. 1 (2008), pp. 65–74. doi: 10.1142/S1793524508000023
- R. Liu, J. Wu, and H. Zhu, Media/psychological impact on multiple outbreaks of emerging infectious diseases, Comput. Math. Methods Med. 8 (2007), pp. 153–164. doi: 10.1080/17486700701425870
- Y. Li and J. Cui, The effect of constant and pulse vaccination on SIS epidemic models incorporating media coverage, Commun. Nonlinear Sci. Numer. Simul. 14 (2009), pp. 2353–2365. doi: 10.1016/j.cnsns.2008.06.024
- Y. Li, C. Ma, and J. Cui, The effect of constant and mixed impulsive vaccination on SIS epidemic models incorporating media coverage, Rocky Mountain J. Math. 38 (2008), pp. 1437–1455. doi: 10.1216/RMJ-2008-38-5-1437
- J. Pang and J. Cui, An SIRS epidemiological model with nonlinear incidence rate incorporating media coverage, Second International Conference on Information and Computing Science, USA, IEEE, 2009, pp. 116–119.
- L. Perko, Differential Equations and Dynamical Systems, Springer-Verlag, New York, 1996.
- M.S. Rahman and M.L. Rahman, Media and education play a tremendous role in mounting AIDS awareness among married couples in Bangladesh, AIDS Res. Therapy 4 (2007), pp. 1–7. doi: 10.1016/0005-7967(66)90037-4
- S. Ruan and W. Wang, Dynamical behavior of an epidemic model with a nonlinear incidence rate, J. Differential Equations 188 (2003), pp. 135–163. doi: 10.1016/S0022-0396(02)00089-X
- F. Schweitzer and R. Mach, The epidemics of donations: Logistic growth and power-laws, PLoS One 3 (2008), p. e1458. doi: 10.1371/journal.pone.0001458
- A.M. Spagnuolo, M. Shillor, L. Kingsland, A. Thatcher, M. Toeniskoetter, and B. Wood, A logistic delay differential equation model for Chagas disease with interrupted spraying schedules, J. Biol. Dynam. 6 (2012), pp. 377–394. doi: 10.1080/17513758.2011.587896
- C. Sun, W. Yang, J. Arino, and K. Khan, Effect of media-induced social distancing on disease transmission in a two patch setting, Math. Biosci. 230 (2011), pp. 87–95. doi: 10.1016/j.mbs.2011.01.005
- J.M. Tchuenche, N. Dube, C.P. Bhunu, R.J. Smith, and C.T. Bauch, The impact of media coverage on the transmission dynamics of human influenza, BMC Public Health 11 (2011), pp. 1–16. doi: 10.1186/1471-2458-11-1
- S.M. Tracht, S.Y. Del Valle, J.M. Hyman, and D.A. Carter, Mathematical modeling of the effectiveness of facemasks in reducing the spread of novel influenza A (H1N1), PLoS One 5 (2010), p. e9018. doi: 10.1371/journal.pone.0009018
- V.I. Utkin, Sliding Modes in Control and Optimization, Springer-Verlag, Berlin, 1992.
- W. Wang, Backward bifurcation of an epidemic model with treatment, Math. Biosci. 201 (2006), pp. 58–71. doi: 10.1016/j.mbs.2005.12.022
- A. Wang and Y. Xiao, A Filippov system describing media effects on the spread of infectious diseases, Nonlinear Anal. Hybr. Syst. 11 (2014), pp. 84–97. doi: 10.1016/j.nahs.2013.06.005
- J-J. Wang, J-Z. Zhang, and Z. Jin, Analysis of an SIR model with bilinear incidence rate, Nonlinear Anal. Real World Appl. 11 (2010), pp. 2390–2402. doi: 10.1016/j.nonrwa.2009.07.012
- L.F. White, J. Wallinga, L. Finelli, C. Reed, S. Riley, M. Lipsitch, and M. Pagano, Estimation of the reproductive number and the serial interval in early phase of the 2009 influenza A/H1N1 pandemic in the USA, Influenza Other Respir. Viruses 3 (2009), pp. 267–276. doi: 10.1111/j.1750-2659.2009.00106.x
- WHO, Avian influenza A(H7N9) overview, 2013. Available at http://www.who.int/influenza/human_animal_interface/influenza_h7n9/WHA_H7N9_update_KeijiFukuda_21May13.pdf?ua=1.
- WHO, Situation Report 1 Ebola virus disease, Guinea, 2014. Available at http://www.afro.who.int/en/clusters-a-programmes/dpc/epidemic-a-pandemic-alert-and-response/sitreps/4070-sitrep-1-ebola-guinea-28-march-2014.html.
- D. Xiao and S. Ruan, Global analysis of an epidemic model with nonmonotone incidence rate, Math. Biosci. 208 (2007), pp. 419–429. doi: 10.1016/j.mbs.2006.09.025
- Y. Xiao, X. Xu, and S. Tang, Sliding mode control of outbreaks of emerging infectious diseases, Bull. Math. Biol. 74 (2012), pp. 2403–2422. doi: 10.1007/s11538-012-9758-5
- Y. Xiao, S. Tang, and J. Wu, Media impact switching surface during an infectious disease outbreak, Scientific Reports, 5 (2015), article 7838.
- X. Zhang and L. Chen, The periodic solution of a class of epidemic models, Comput. Math. Appl. 38 (1999), pp. 61–71. doi: 10.1016/S0898-1221(99)00206-0
- Z. Zhang, T. Ding, W. Huang, and Z. Dong, Qualitative Theory of Differential Equations, Vol. 101, Amer. Math. Soc., Providence, 1992.
- J. Zu and L. Wang, Periodic solutions for a seasonally forced SIR model with impact of media coverage, Adv. Difference Equ. 2015 (2015), pp. 1–10. doi:10.1186/s13662-015-0477-8.