BMC Public Health. 2011; 11(Suppl 1): S5. Published online Feb 25, 2011. PMCID: PMC3317585
The impact of media coverage on the transmission dynamics of human influenza
Jean M Tchuenche,1 Nothabo Dube,2 Claver P Bhunu,2 Robert J Smith?,3 and Chris T Bauch1
1Department of Mathematics and Statistics, University of Guelph, Guelph, Ontario N1G 2W1, Canada
2Department of Applied Mathematics, National University of Science and Technology, Box AC 939 Ascot, Bulawayo, Zimbabwe
3Department of Mathematics and Faculty of Medicine, The University of Ottawa, 585 King Edward Ave, Ottawa ON K1N 6N5, Canada
Corresponding author.
Jean M Tchuenche: jmt_biomaths/at/yahoo.co.uk; Nothabo Dube: mamthe75/at/gmail.com; Claver P Bhunu: cpbhunu/at/gmail.com; Robert J Smith?: rsmith43/at/uottawa.ca; Chris T Bauch: cbauch/at/uoguelph.ca
##### Supplement
Mathematical Modelling of Influenza
Jean M Tchuenche and Robert Smith?
1471-2458-11-S1-info.pdf
Background
There is an urgent need to understand how the provision of information influences individual risk perception and how this in turn shapes the evolution of epidemics. Individuals are influenced by information in complex and unpredictable ways. Emerging infectious diseases, such as the recent swine flu epidemic, may be particular hotspots for a media-fueled rush to vaccination; conversely, seasonal diseases may receive little media attention, despite their high mortality rate, due to their perceived lack of newness.
Methods
We formulate a deterministic transmission and vaccination model to investigate the effects of media coverage on the transmission dynamics of influenza. The population is subdivided into different classes according to their disease status. The compartmental model includes the effect of media coverage on reporting the number of infections as well as the number of individuals successfully vaccinated.
Results
A threshold parameter (the basic reproductive ratio) is analytically derived and used to discuss the local stability of the disease-free steady state. The impact of costs that can be incurred, which include vaccination, education, implementation and campaigns on media coverage, are also investigated using optimal control theory. A simplified version of the model with pulse vaccination shows that the media can trigger a vaccinating panic if the vaccine is imperfect and simplified messages result in the vaccinated mixing with the infectives without regard to disease risk.
Conclusions
The effects of media on an outbreak are complex. Simplified understandings of disease epidemiology, propogated through media soundbites, may make the disease significantly worse.
Infectious diseases are responsible for a quarter of all deaths in the world annually, the vast majority occurring in low- and middle-income countries [1]. There are diseases such as SARS and flu that exhibit some distinct features such as rapid spatial spread and visible symptoms [2]. These features, associated with the increasing trend of globalization and the development of information technology, are expected to be shared by other emerging/re-emerging infectious diseases. It is therefore important to refine classical mathematical models to reflect these features by adding the dimensions of massive news coverage that have great influence not only on the individual behaviours but also on the formation and implementation of public intervention and control policies [2].
People’s response to the threat of disease is dependent on their perception of risk, which is influenced by public and private information disseminated widely by the media. While government agencies for disease control and prevention may attempt to contain the disease [3], the general information disseminated to the public is often restricted to simply reporting the number of infections and deaths. Mass media are widely acknowleged as key tools in risk communication [4,5], but have been criticised for making risk a spectacle to capitalise on audience anxiety [6,7].
The original interpretation of media effects in communication theory was a “hypodermic needle” or “magic bullet” theory of the mass media. Early communication theorists [8,9] imagined that a particular media message would be directly injected into the minds of media spectators. This theory of media effects, in which the mass media has a direct and rapid influence on everyday understanding, has been substantially revised. Contemporary media studies analyses how media consumers might only partially accept a particular media message [10], how the media is shaped by dominant cultural norms [11,12] and how media consumers resist dominant media messages [13,14]. It follows that media effects may sway people into panic (eg swine flu), especially with a disease where scientific evidence is thin or nonexistent. Conversely, media may have little effect on seasonal diseases (eg regular influenza).
Media reporting plays a key role in the perception, management and even creation of crisis [6]. Since media reports are retrievable and because the messages are widely distributed, they gain authority as an intersubjective anchorage for personal recollection [4]. At times of crisis, non-state-controlled media thrive, while state-controlled media are usually rewarded for creating an illusion of normalcy [6]. Media exposure and attention partially mediate the effects of variables such as demographics and personal experience on risk judgments [5]. The role of media coverage on disease outbreaks is thus crucial and should be given prominence in the study of disease dynamics.
Klein et al., [15] noted that much more research is needed to understand how provision of information influences individual risk perception and how it shapes the evolution of epidemics; for example, individuals may overprotect, which can have additional consequences for the spread of disease. An example of such complex dynamics is the 1994 outbreak of plague in a state in India: after the announcement of the disease, many people fled the state of Surat in an effort to escape the disease, thus carrying the disease to other parts of the country [16]. Even though information on the number of cases and deaths can have an adverse effect, the number of those vaccinated has not been given prominence.
A handful of mathematical models have described the impact of media coverage on the transmission dynamics of infectious diseases. Liu et al.[2] examined the potential for multiple outbreaks and sustained oscillations of emerging infectious diseases due to the psychological impact from reported numbers of infectious and hospitalized individuals. Liu & Cui [3] analysed a compartment model that described the spread and control of an infectious disease under the influence of media coverage. Li & Cui [17] incorporated constant and pulse vaccination in SIS epidemic models with media coverage. Cui et al., [18] showed that when the media impact is sufficiently strong, their model – with incidence rate being of the exponential form capturing the alertness to the disease of each susceptible individual in the population – exhibits multiple positive equilibria (also see [2]) which poses a challenge to the prediction and control of the outbreaks of infectious diseases.
The aim of this study is to investigate the impact of media coverage on the spread and control of an influenza strain when a vaccine is available, and where the media reporting of both disease dynamics and vaccination is high. Vaccination is one of the most effective tools for reducing the burden of infectious diseases [19]. However, despite their public-health benefit, vaccination programs face obstacles. Individuals often refuse or avoid vaccinations they perceive to be risky. Recently, rumours that the polio vaccine could cause sterility and spread HIV have hampered polio eradication in Nigeria [20], while misplaced fears of autism in the developed world have stoked vaccination fears [21]. Reporting the number of individuals who vaccinate may have a positive effect on the disease transmission by increasing the vaccination rate.
Conversely, behavioural interventions can also have an enormous effect on the course of a disease [22,23] Our model considers the same contact rate after a media alert, as proposed by Liu & Cui [3], but there are fundamental differences in both models. They consider the classical SIR type model, while vaccination is included in ours to reflect transmission dynamics of human influenza.
We divide the population (N) into four sub-populations, according to their disease status: susceptible (S), vaccinated (V ), infected (I), and recovered (R). Our model monitors the dynamics of influenza based on a single strain without effective cross-immunity against the strain. The susceptible population is increased by recruitment of individuals (either by birth or immigration), and by the loss of immunity, acquired through previous vaccination or natural infection. This population is reduced through vaccination (moving to class V ), infection (moving to class I) and by natural death or emigration. The population of vaccinated individuals is increased by vaccination of susceptible individuals. Since the vaccine does not confer immunity to all vaccine recipients, vaccinated individuals may become infected, but at a lower rate than unvaccinated. The vaccinated class is thus diminished by this infection (moving to class I) by waning of vaccine-based immunity (moving to class S) and by natural death. The population of infected individuals is increased by infection of susceptibles, including those who remain susceptible despite being vaccinated. It is diminished by natural death, death due to disease and by recovery from the disease (moving to class R). The recovered class is increased by individuals recovering from their infection and is decreased as individuals succumb to natural death. Media coverage is introduced into the model via a saturated incidence function.
A schematic model flowchart is depicted in Figure Figure11.
 Figure 1The model Schematic model flow diagram
Model equations
The transmission model with media coverage is given by the following deterministic system of nonlinear ordinary differential equations:
 (1)
 (2)
 (3)
 (4)
where Λ is the rate at which individuals are recruited into the population (recruitment of infectives is ignored for now); θ is the rate at which susceptible individuals receive the vaccine; µ is the the rate at which people leave the population, through natural death or emigration. We assume this rate to be the same for all sub-populations. β1 is the rate at which susceptibles get infected; ω is the rate at which vaccine-based immunity wanes; γ is the vaccine efficacy; α is the death rate due to the infection; and λ is the recovery rate from infection. The terms and measure the effect of reduction of the contact rate when infectious and vaccinated individuals are reported in the media [2,3,18]. The half-saturation constant mI > 0 reflects the impact of media coverage on the contact transmission. The function is a continuous bounded function which takes into account disease saturation or psychological effects [24]. We note that recovered individuals cannot be vaccinated. Also, a vaccinated individual who gets infected and then recovers will return to the susceptible class with no vaccine protection. This is true even if ω is quite small but σ and λ are large. For example, if vaccination lasts three years, but recovery and loss of immunity takes 6 months, then we are assuming this person is subsequently unvaccinated.
In the Michaelis-Menten functional response, the rate at which information is spread by the media rises as infectives increase, but eventually levels off at a plateau (or asymptote) at which the information (rate) remains constant (i.e. it has reached a maximum number of individuals due to information saturation) regardless of the increase in infections. Such dynamics can easily be observed in the spread of rumours, gossip and jokes (also known as randomized broadcast) [25,26]. This constant coverage is extended by examining more complex effects which involve more than just reducing contacts down the line. The news in particular is extremely fickle so that what is news one day may be forgotten about next week; including the media effects in some more sophisticated way such as by an impulsive pulsing is also investigated. The limited power of the infection due to contact is accounted for by the saturation incidence. The first available information is the reported number of infected individuals when the disease is emerging. We assume that media coverage can slow but not prevent disease spread, so β1β2 and β1β3.
The above model is closely related to those in [27,28] to analyze the transmission dynamics of human influenza, but there are some differences. In [27], the authors consider the inflow of infective immigrants, while in [28] the model includes treatment. Neither of these are considered here. Our model is clearly a crude reflection of the complicated nonlinear phenomena of the transmission dynamics, and it does not incorporate the self-control property due to the change of avoidance patterns of individuals at different stages of the infectious process [2]. News coverage may have a significant impact on avoidance behaviours at both individual and society levels, which may reduce the effective contact between susceptible and infectious individuals; we include this via a saturation incidence functional response.
Since the model monitors human populations, all the variables and parameters of the model are nonnegative. Based on biological considerations the system of equations (1)-(4) will be studied in the following region,
which is positively invariant and attracting (thus, the model is mathematically and epidemiologically well-posed); it is therefore sufficient to consider solutions in Ω. Existence, uniqueness and continuation results for model system (1)-(4) hold in this region and all solutions of this system starting in Ω remain in Ω for all t ≥ 0.
Stability of the equilibrium states
The disease-free equilibrium of the system is given by
The endemic equilibrium of the system is given by
It satisfies and
where h1(I) = m1 + I, h2(I) = (θ + µ + β1I)h1(I) – β2I2, h3(I) = (θ + ω + (1 – γ)β1I)h1(I) – β3(1 – γ)I2. Substituting the above into the second equation at equilibrium will yield the expression for Î after some rearrangement. For illustration, suppose θ = β2 = β3 = 0. Then the endemic equilibrium satisfies
where Î is the positive solution to the quadratic
The basic reproductive ratio, Rv, is defined as the expected number of secondary infections caused by an infective individual upon entering a totally susceptible population [29-31]. This quantity is not only important in describing the infectious power of the disease, but can also can supply information for controlling the spread of the disease [32]. The linear stability of Ev0 is governed by the basic reproductive ratio Rv. Using the next-generation method [31], we have
and
The basic reproductive ratio is the spectral radius ρ(FV–1) which is
 (5)
Local stability of the disease-free equilibrium
Lemma 1The disease-free equilibrium Ev0 is locally asymptotically stable if Rv < 1, and unstable if Rv > 1.
Proof. The Jacobian of the system evaluated at Ev0 is given by
The eigenvalues of JEV0 are
For local stability of the disease-free equilibrium, we require that all the eigenvalues be negative. Three of the eigenvalues satisfy this condition while ς2 < 0 implies that Rv < 1 and, consequently, all the eigenvalues of the Jacobian matrix above have negative real part. Thus, the disease-free equilibrium is locally asymptotically stable.
Global stability of the disease-free equilibrium
We adopt the method of Castillo-Chavez et al,[33] and we rewrite the set of model equations in the form
with G(XG,0) = 0. XG 3 denotes the number of uninfected classes and ZG denotes the number of infected classes. denotes the disease-free equilibrium of the system where For the set of equations in (1)-(4), we set XG = (S, V, R) and ZG = (I). The conditions (H1) and (H2) below must be met for global stability.
(H1) For is globally asymptotically stable.
(H2) G(XG, ZG) = AGZGĜ(XG, ZG), Ĝ(XG, ZG) ≥ 0 for (XG, ZG) Ω where is an M-matrix (the off-diagonal elements of A are nonnegative) and Ω is the region where the model makes biological sense.
If the above two conditions are satisfied, then the following theorem holds.
Theorem 2 (Castillo-Chavez et al,[33]): The fixed pointis a globally stable equilibrium of (2.28) provided that Rv < 1 and that assumptions (H1) and (H2) are satisfied.
Therefore, Ev0 is globally asymptotically stable (GAS) since Ĝ(XG,ZG) > 0. The GAS of Ev0 excludes any possibility of the phenomenon of backward bifurcation. We note that the GAS of the DFE Ev0 when σ = 0 is straightforward.
Our objective in this section is to extend the initial model to include two intervention methods, called controls, represented as functions of time and assigned reasonable upper and lower bounds, each representing a possible method of influenza intervention. Using optimal control theory and numerical simulations, we determine the benefit of vaccination and media coverage when the latter has positive or negative effect on the former.
We will integrate the essential components into one SIVR-type model to accommodate the dynamics of an influenza outbreak determined by population-specific parameters such as the effect of contact reduction when infectious and vaccinated individuals are reported in the media.
Let uv and um be the control variables for vaccination and media coverage respectively. Thus, model (1)-(4) now reads
 (6)
 (7)
 (8)
 (9)
A balance of multiple intervention methods can differ between populations. A successful mitigation scheme is one which reduces influenza-related deaths with minimal cost. A control scheme is assumed to be optimal if it maximizes the objective functional
 (10)
The first two terms represent the benefit of the susceptible and vaccinated populations. The parameters B1 and B2 represent the weight constraints for the infected population and the control, respectively. They can also represent balancing coefficients transforming the integral into dollars expended over a finite time period of T days [34]. The goal is to maximize the populations of susceptible and vaccinated individuals, minimize the population of infectives, maximize the benefits of media coverage and vaccination, while minimizing the systemic costs of both media coverage and vaccination. The value uv(t) = um(t) = 1 represents the maximal control due to vaccination and media coverage, respectively. The terms and represent the maximal cost of education, implementation and campaigns on both vaccination and media coverage. S(t) and V(t) account for the fitness of the susceptible and the vaccinated groups as a result of a reduction in the rate at which the vaccine wanes, and vaccination and treatment efforts are implemented [35]. We thus seek optimal controls and such that
where U = {(uv, um)|uv, um measurable, 0 ≤ a11uvb11 ≤ 1, 0 ≤ a22umb22 ≤ 1, t [0, tf]} is the control set, with t [t0, tf]. The basic framework of this problem is to characterize the optimal control.
Existence of an optimal control
The existence of an optimal control can be obtained by using a result by Joshi [36] and Fister et al.[37].
Theorem 3Consider the control problem with the system of Equations (4.1)-(4.4). There exists an optimal controlsuch that max
Proof. To prove this theorem on the existence of an optimal control, we use a result from Fleming and Rishel [38] (Theorem 4.1 pp. 68-69), where the following properties must be satisfied.
1. The set of controls and corresponding state variables is nonempty.
2. The control set U is closed and convex.
3. The right-hand side of the state system is bounded above by a linear function in the state and control.
4. The integrand of the functional is concave on U and is bounded above by c2c1(|uv|k + |um|k), where c1, c2 > 0 and k > 1.
An existence result in Lukes [39] (Theorem 9.2.1) for the system of equations (6)-(9) for bounded coefficients is used to give the first condition. The control set is closed and convex by definition. The right-hand side of the state system (Equations (4.1)-(4.4)) satisfies Condition 3 since the state solutions are a priori bounded. The integrand in the objective functional, , is concave on U. Furthermore, c1, c2 > 0 and k > 1, so
 (11)
Therefore, the optimal control exists, since the left-hand side of (11) is bounded; consequently, the states are bounded.
Since there exists an optimal control for maximizing the functional (10) subject to equations (6)-(9), we use Pontryagin’s Maximum Principle to derive the necessary conditions for this optimal control. Pontryagin’s Maximum Principle introduces adjoint functions that allow us to attach our state system (of differential equations), to our objective functional. After first showing existence of optimal controls, this principle can be used to obtain the differential equations for the adjoint variables, corresponding boundary conditions and the characterization of an optimal control and . This characterization gives a representation of an optimal control in terms of the state and adjoint functions. Also, this principle converts the problem of minimizing the objective functional subject to the state system into minimizing either the Lagrangian or the Hamiltonian with respect to the controls (bounded measurable functions) at each time t[40]. The Lagrangian is defined as
where w11(t) ≥ 0, w12(t) ≥ 0 are penalty multipliers satisfying w11(t)(a11uv(t)) + w12(t)(uv(t) – b11) at the optimal , and w21(t) ≥ 0, w22(t) ≥ 0 are penalty multipliers satisfying w21(t)(a22um(t)) + w22(t)(um(t) – b22) at the optimal .
Given optimal controls and , and solutions of the corresponding state system (6)-(9), there exist adjoint variables λi, for i = 1, 2, 3, 4 satisfying the following equations
with transversality conditions λi[tf] = 0, for i = 1, 2, 3, 4. To determine the interior maximum of our Lagrangian, we take the partial derivatives of L with respect to uv and um, respectively, and set it to zero. Thus,
To determine an explicit expression for our controls , (without w11, w12, w21, w22), a standard optimality technique is utilized. The following cases are considered to determine a specific characterization of the optimal control.
Case 1: Optimality of
1. On the set . Hence, the optimal control is
2. On the set .We have
or
since w12 ≥ 0.
3. On the set . Hence
or
Combining all the three sub-cases in a compact form gives
 (12)
Case 2: Optimality of
1. On the set . We have
2. On the set . We have
or
since w22 ≥ 0.
3. On the set . Hence
or
Combining all the three sub-cases in a compact form gives
 (13)
The optimal system
The optimality system consists of the state system coupled with the adjoint system, with the initial conditions, the transversality conditions and the characterization of the optimal control:
where and are given by expressions (12) and (13), respectively, with S(0) = S0, I(0) = I0, V(0) = V0, R(0) = R0 and λi[tf] = 0 for i = 1,··· ,4. Due to the a priori boundedness of the state and adjoint functions and the resulting Lipschitz structure of the ODEs, we obtain the uniqueness of the optimal control for small [tf] [36]. The uniqueness of the optimal control follows from the uniqueness of the optimality system.
The state system of differential equations and the adjoint system of differential equations together with the control characterization above form the optimality system solved numerically and depicted in Figures Figures2,2, ,3,3, ,4,4, ,55.
 Figure 2Optimality effect when the weight constraint for the infected population varies and media has a beneficial effect on the vaccine. Graphs of the optimality system when media coverage has a beneficial effect on the vaccination rate and when the weight constraint (more ...)
 Figure 3Optimality effect when the weight constraint for the control varies and media has a beneficial effect on the vaccine. Graphs of the optimality system when media coverage has a beneficial effect on the vaccination rate and when the weight constraint for (more ...)
 Figure 4Optimality effect when the weight constraint for the infected population varies and media has an adverse effect on the vaccine. Graphical representation of the optimality system when media coverage has an adverse effect on the vaccination rate and when (more ...)
 Figure 5Optimality effect when the weight constraint for the control varies and media has an adverse effect on the vaccine. Graphs of the optimality system when media coverage has an adverse effect on the vaccination rate and when the weight constraint for the (more ...)
The model with pulse vaccination
The general model with pulse vaccination is given as
for ttk, where tk is the time of the kth vaccination. We may have tk+1tk either constant or not, as we choose. The impulsive effect is given by
ΔS = –θS
ΔV = θS
when t = tk. Here, is the change in state at the impulse time.
In this model, vaccination occurs at fixed times, not continuously. This is closer to reality, since vaccination centres are only open at certain times, when people may get vaccinated in waves. Similarly, media stories tend to clump together, so that a big news story occurs on one day, which may trigger a short period of intense vaccination. We shall use a simplified version of this model to illustrate the possibility that media may have an adverse effect.
Consider the following scenario. At the onset of the outbreak, the media - and hence the general population - is unaware of the disease and thus nobody gets the vaccine, allowing the disease to spread in its initial stages. At some point, there is a critical number of infected individuals, whereupon people are sufficiently aware of the infection to change their behaviour. We suppose that, initially, new infected people arrive at fixed times.
We further assume that vaccinated people mix more than susceptibles. In this case, people who are vaccinated feel confident enough to mix with the infected, even though they may still have the possibility to contract the virus. This might be the case for health-care workers, for instance, who get vaccinated and then have to tend to the sick.
Mathematically, we have a threshold for the critical number of infectives, Icrit.
For I <Icrit, this model would look like
For I >Icrit, the model becomes
with β4β6 ≥ 0.
However, to illustrate the adverse affect, we shall simplify the model even further. For a short timescale, we can assume recovery is permanent, so σ = 0. Thus, we can ignore the R equation.
For I <Icrit, we assume that there is no mixing, but rather that new infectives arrive impulsively into the system at fixed times tk and in numbers Ii, where Ii Icrit. (If the new infectives arrive at irregular times, then the broad results will be unchanged.)
For I >Icrit, fear of the disease keeps susceptibles from mixing with the infected, but the vaccinated will.
Thus β4 = β6 = 0. Since Ii Icrit, we can assume that, for I >Icrit, the effects of new infectives are negligible.
The model then becomes
 (14)
 (15)
 (16)
 (17)
for I <Icrit and
 (18)
 (19)
 (20)
for I >Icrit.
Thus, the effects of the media are to trigger a vaccinating panic whenever the number of infectives is large enough. We kept the model with impulse vaccination as simple as possible since even this simplified version shows that media reports could have an adverse effect.
Suppose new infectives appear regularly, so that tk+1tk = τ. (If not, the analysis generalizes quite easily.) For tk <t <tk+1, we have
where is the value immediately after the kth impulse. Then, since the period is constant, we have
This is a recursion relation with solution
Consequently,
Thus, if m+ >Icrit, then eventually the system will switch from model (14)-(17) to model (18)-(20). The endemic equilibrium in model (18)-(20) satisfies
The Jacobian is
At the endemic equilibrium, . Thus, we have
The characteristic equation is
It follows that the endemic equilibrium is stable if Î >Icrit. Thus, even in an extremely simplified version of the model, the media may make things significantly worse than if no media effect were included. We kept this model deliberately simple, partly for mathematical tractability and partly to show that the media effects apply even in this idealised scenario.
Note that, in reality, the fluctuations would apply in the upper region as well, making the actual value even
larger. In the lower region, we ignored interaction between susceptibles and infectives (ie we assume β4 = β6 = 0). The effect of including these terms would be to slow the exponential decay between impulses (or possibly cause it to increase). This would only increase the effect seen here.
In summary, a small series of outbreaks that would equilibrate at some maximum level m+ >Icrit will, as a result of the media, instead equilibrate at a much larger value I >m+ >Icrit. The driving factor here is if an imperfect vaccine causes overconfidence, so that people who have been vaccinated mix significantly more with infectives than susceptibles do. If this happens (as would be quite likely; most people who have been vaccinated feel invulnerable, even if the vaccine is not perfect, largely thanks to media oversimplifications), then the media effect is likely to be adverse. A simplified version of the model with pulse vaccination shows that the media can make things worse, if the vaccine is imperfect because the vaccinated mix over-confidently with the infectives.
We now return to model (6)-(9) and illustrate some of the properties discussed in the previous sections. The parameter values that we use for numerical simulations are in Table Table1.1. Initial conditions: S(0) = 200.0, I(0) = 1.0, V(0) = 10.0, R(0) = 0.0. The parameter θ varies between 0.3 and 0.7 with an average of 0.5 [41]. We consider an imperfect flu vaccine for which the waning rate is about 0.15. The relationship between β2 and β3 is not very obvious; consequently, we can either assume equality or that the former is slightly greater than the latter. Transmission dynamics of infectious diseases with and without media coverage have already been carried out in previous studies, but these models do not account for the vaccination coverage. Therefore, we illustrate some numerical results for the model with optimal control when media coverage has (i) a beneficial effect (see Figures Figures22 and and3)3) (ii) and adverse effect on the vaccination rate (see Figures Figures44 and and55).
 Table 1Parameter values
The optimality system is solved using an iterative method with a fourth order Runge-Kutta scheme. Starting with a guess for the adjoint variables, the state equations are solved forward in time. Then the state values obtained are used to solve the adjoint equations backward in time; the iterations continue until convergence. Simulations are carried out to determine how maximizing media coverage enhances vaccination. The effects of costs that can be incurred, which include education, implementation and campaigns on media coverage, are also studied to evaluate how these costs can affect the transmission of human influenza. We increase the value of B2 (the cost weight) in Figure Figure22 to assess how the populations of susceptibles, infectives, vaccinated and recovered individuals are altered. In Figure Figure3,3, we investigate how increasing minimization of infectives through increasing the weight B1 affects the control of human influenza transmission. We do the same in Figures Figures44 and and5,5, respectively, to see how, if media coverage has an adverse effect, the various populations behave. In Figure Figure4,4, we vary the cost weight, while in Figure Figure5,5, we vary the weight of minimizing infectives.
We note from Figure Figure22 that, during the initial days, there is a very sharp drop in the population of infected individuals, while other populations show increases. Increasing minimization of infectives, while keeping costs low, can lead to the disease being controllable. The slight rise and fall, after the initial 20 days, in the population of infectives could be attributed to complacency on the part of some individuals (or may be due to oscillations in the system independent of external factors). We find that people tend to relax after the initial shock of the disease threat. However, we note that this is not for long, and this could be attributed to the fact that vaccination levels continue to rise, so as people continue to receive vaccination, infection is controlled. Thus, if costs are kept minimal, and more people are able to access media and vaccination, then infection can be controlled. Both vaccination and media coverage continue at optimum levels as a result of the low costs and minimization of infectives.
From Figure Figure3,3, as costs are increased, few people have access to media and vaccination; as a result, low numbers get vaccinated against the disease. In the long run, the infection levels rise. The degree of media coverage and vaccination also decrease as a result of the exorbitant costs. With the little available media coverage and the few vaccinated individuals, we find that, due to information filtration, there is a jump in the vaccination levels, though these only last briefly; as the degree of media coverage and vaccination decrease, so do the vaccination levels.
From Figure Figure4,4, even though costs are kept at minimal levels, the negative reports concerning vaccination result in a drastic reduction in the vaccination levels. After some time, we note a slight increase in the vaccination levels; however, these numbers remain very low. This could be due to the fact that, as infection rises, a few will risk getting vaccinated in the hope of being cured. Thus, media coverage can have adverse effects if people’s perception towards the vaccine is negatively influenced by the media.
In Figure Figure5,5, both media coverage and vaccination are eventually withdrawn. Very low numbers get vaccinated. It is only when infection escalates that vaccination levels also increase as some might find it better to try to prevent the infection, despite the negativity towards vaccination in the media. Figure Figure66 illustrates other potential adverse effects that media may have, if the effect is to trigger a vaccinating panic where vaccinated individuals are not fully protected and mix with infected individuals but susceptible individuals do not. In this case, the number of infected individuals may increase sharply as a result of the media. Figure Figure77 illustrates the long-term results of a media-induced vaccination panic. Without media effects, the result is a low-level infection. When the media triggers a vaccinating panic, there is a large outbreak, followed by an endemic level of infected individuals significantly higher than the level of infected individuals without the media effects. Note that these examples assume no post-vaccination mixing of susceptible and infected individuals.
 Figure 6The vaccination panic threshold The effect of the vaccination panic threshold using the simplified model (14)-(20). Without media triggering a vaccinating panic, the number of infected individuals remains low (solid purple curve). However, if the media (more ...)
 Figure 7A media-induced epidemic Long-term dynamics for the simplified model (15)-(20). Without media effects, the maximum number of infected individuals remains low (solid purple curve). When media effects trigger a vaccinating panic where partially protected (more ...)
Figure Figure88 illustrates the cases when post-vaccination mixing of susceptible and infected individuals is maximal (β4 = β5, β6 = 0), 50% (β4 = β5, β6 = β4/2) or zero (β4 = β6 = 0). Thus, if susceptible and infected individuals mix after a vaccinating panic has occurred, the effect is an earlier outbreak and a larger number of infected individuals.
 Figure 8The effect of different post-panic mixing rates Comparison of mixing rates between susceptibles and infected individuals after a vaccinating panic. The more susceptibles mix with infecteds, the earlier the outbreak occurs and the larger the number of (more ...)
Media simplifications can lead to overconfidence in the idea of a vaccine as a cure-all. The result is not just a vaccinating panic and a blow-out epidemic, but a net increase in the endemic equilibrium. Thus, media coverage of an emerging epidemic can fan the flames of fear and also implicitly reinforce an imperfect solution as the only answer.
We have formulated and investigated a simple deterministic vaccination model describing the effects of media coverage on the transmission dynamics of influenza. The media effect due to reporting the number of infections as well as the number of individuals successfully vaccinated is introduced into the compartmental model via a saturated incidence-type function. The impact of costs that can be incurred, which include vaccination, education, implementation and campaigns on media coverage, are also investigated using optimal control theory applied via the Pontryagin’s maximum principle. A simplified version of the model with pulse vaccination shows that the media can have an adverse effect if the vaccine is imperfect and the vaccinated mix over-confidently with the infectives. Numerical simulations are carried out to support the analytical results.
We note, however, that our caricature model is not complete; a more comprehensive study will require interdisciplinary research across traditional boundaries of social, natural, medical sciences and mathematics [2]. Nevertheless, our work provides some insights into the effects of media reporting on the transmission dynamics of infectious diseases for which a vaccine exists. The present study is in no way exhaustive and can be extended in various ways: for example, to investigate the case in which there is media coverage but people ignore it (in which case the vaccination rate is unchanged despite the control). Thus, the effects of media on an outbreak of influenza with a partially effective vaccine may be complicated. While the media may encourage more people to get vaccinated, they may also trigger a vaccinating panic or promote overconfidence in the ability of a vaccine to fully protect against the disease. This may have potentially disastrous consequences in the face of a new pandemic.
Competing interests
The authors declare that they have no competing interests.
Authors contributions
JMT, RJS and CTB developed the model. JMT and CTB designed and formulated the study framework and analyzed the model. CPB and ND carried out the optimal control analysis and numerical simulations for Figures Figures2,2, ,3,3, ,4,4, ,5.5. RJS wrote the section on adverse effects, some of the introduction, performed numerical simulations for Figures Figures11166,,77,,8,8, and edited the manuscript. All authors read and approved the final manuscript.
Acknowledgements
We thank Shoshana Magnet, Penelope Ironstone-Catterall and Tim Reluga for technical discussions. JMT was supported through a postdoctoral appointment funded by a grant from the Ontario Ministry of Research and Innovation awarded to CTB. RJS is supported by an NSERC Discovery Grant, an Early Researcher Award and funding from MITACS. Handling editor for this manuscript was Jane Heffernan.
This article has been published as part of BMC Public Health Volume 11 Supplement 1, 2011: Mathematical Modelling of Influenza. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2458/11?issue=S1.
• Laxminarayan R, Mills AJ, Breman JG, Measham AR, Alleyne G, Claeson M, Jha P, Musgrove P, Chow J, Shahid-Salles S, Jamison DT. Advancement of global health: key messages from the Disease Control Priorities Project. Lancet. 2006;367:1193–1208. doi: 10.1016/S0140-6736(06)68440-7. [PubMed]
• Liu R, Wu J, Zhu H. Media/psychological impact on multiple outbreaks of emerging infectious disease. Comput. 2007;8(3):153–164. doi: 10.1080/17486700701425870.
• Liu Y, Cui J. The impact of media coverage on the dynamics of infectious disease. Int. 2008;1:65–74. doi: 10.1142/S1793524508000023.
• Simpson CR. Nature as News: Science Reporting in The New York Times 1898 to 1983. The International Journal of Politics, Culture and Society. 1987;1(2):218–241.
• Slater MD, Rasinski KA. Media Exposure and Attention as Mediating Variables Influencing Social Risk Judgments. Journal of Communication. 2005;55(4):810–827. doi: 10.1111/j.1460-2466.2005.tb03024.x.
• Ma R. Media, Crisis, and SARS: An Introduction. Asian Journal of Communication. 2005;15(3):241–246. doi: 10.1080/01292980500260656.
• Roche JP, Muskavitch MAT. Limited Precision in Print Media Communication of West Nile Virus Risks. Science Communication. 2003;24:353–365. doi: 10.1177/1075547002250300.
• Lasswell H. Propaganda Technique in the World War. Cambridge: MIT Press; 1971.
• Lazarsfeld P. Radio and the Printed Page: An Introduction to the Study of Radio and Its Role in the Communication of Ideas. New York: Duell Sloan, and Pearce; 1940.
• Radway J. Reading the Romance. Chapel Hill: The University of North Carolina Press; 1984.
• Hall S. Cultural Representations and Signifying Practices. London: Sage Publications; 1997.
• Williams R. Television: Technology and Cultural Form. London: Fontana; 1974.
• Hebdige D. Subculture: The meaning of style. London: Metkuen & Co; 1979.
• Treichler P. How to Have Theory in an Epidemic. Durham: Duke University Press; 1999.
• Klein E, Laxminarayan R, Smith DL, Gilligan CA. Economic incentives and mathematical models of disease. Envir. Dev. Economics. 2007;12:707–732.
• Ramalingaswami V. Psychosocial effects of the 1994 plague outbreak in Surat, India. Military Med. 2001;166:29–30. [PubMed]
• Li Y, Cui J. The effect of constant and pulse vaccination on SIS epidemic models incorporating media coverage. Commun. 2009;14:2353–2365. doi: 10.1016/j.cnsns.2008.06.024.
• Cui J, Sun Y, Zhu H. The impact of media on the spreading and control of infectious disease. J. 2003;20:31–53. doi: 10.1007/s10884-007-9075-0.
• Reluga TC, Bauch CT, Galvani AP. Evolving public perceptions and stability in vaccine uptake. Math. 2006;204:185–198. doi: 10.1016/j.mbs.2006.08.015. [PubMed]
• Samba E, Nkrumah F, Leke R. Getting polio eradication back on track in Nigeria. N. 2004;350:645–646. doi: 10.1056/NEJMp038210. [PubMed]
• Gross L. A Broken Trust: Lessons from the Vaccine-Autism Wars. PloS Biology. 2009;7(5):e1000114. doi: 10.1371/journal.pbio.1000114. [PubMed]
• Barry M. The tail end of Guinea Worm – global eradication without a drug or a vaccine. New England Journal of Medicine. 2007;356(25):2561–2563. doi: 10.1056/NEJMp078089. [PubMed]
• Smith? RJ, Li J, Gordon R, Heffernan JM. Can we spend our way out of the AIDS epidemic? A world halting AIDS model. BMC Public Health. 2009;9(Suppl 1):S15. [PubMed]
• Capasso V, Serio G. A generalization of the Kermack-McKendrick deterministic epidemic model. Math. 1978;42:43–62. doi: 10.1016/0025-5564(78)90006-8.
• Chierichetti F, Lattanzi S, Panconesi A. Rumor Spreading in social networks. Lecture Notes in Computer Science. 2009;5556:375–386. full_text.
• Pitell B. On spreading a rumour. SIAM J. 1987;47:213–223. doi: 10.1137/0147013.
• Alexander ME, Bowman C, Moghadas SM, Summers R, Gumel AB, Sahai BM. A Vaccination model for transmission dynamics of influenza. SIAM J. Appl. Dyn. Syst. 2004;3:503–524. doi: 10.1137/030600370.
• Rwezaura H, Mtisi E, Tchuenche JM. In: Recent Progress in Infectious Disease Modelling. Tchuenche JM, Chiyaka C, editor. New York: Nova Science Publishers; 2009. A Mathematical Model of Influenza with Vaccination and Treatment; pp. 31–83.
• Diekmann O, Heesterbeek JAP, Metz JAJ. On the definition and computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations. J. 1990;28(4):365–382. doi: 10.1007/BF00178324. [PubMed]
• Heffernan JM, Smith RJ, Wahl LM. Perspectives on the basic reproductive ratio. J. R. Soc. Interface. 2005;2(4):281–293. doi: 10.1098/rsif.2005.0042. [PubMed]
• van den Driessche P, Watmough J. Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Math. 2002;180:29–48. doi: 10.1016/S0025-5564(02)00108-6. [PubMed]
• Cui J, Tao X, Zhu H. An SIS infection model incorporating media coverage. Rocky Mountain J. Math. 2008;38(5):1323–1334. doi: 10.1216/RMJ-2008-38-5-1323.
• Castillo-Chavez C, Feng Z, Huang W. In: Mathematical Approaches for Emerging and Reemerging Infectious Diseases: An Introduction. Castillo-Chavez C, van den Driessche P, Kirschner D, Yakubu AA, editor. Berlin-Heidelberg-New York: Springer-Verlag; 2002. On the computation of R0 and its role on global stability; pp. 229–250.
• Miller Neilan RLM, Schaefer E, Gaff H, Fister KR, Lenhart S. Modeling optimal intervention strategies for cholera. Bull Math Biol. 2010;72(8):2004–2018. doi: 10.1007/s11538-010-9521-8. [PubMed]
• Tchuenche JM, Khamis SA, Agusto FB, Mpeshe CS. Optimal control and sensitivity analysis of an influenza model with treatment and vaccination. Acta Biotheor. 2011;59(1):1–28. doi: 10.1007/s10441-010-9095-8. [PubMed]
• Joshi HR. Optimal Control of an HIV Immunology Model. Optim. 2002;23:199–213. doi: 10.1002/oca.710.
• Fister KR, Lenhart S, McNally JS. Optimizing chemotherapy in an HIV model. Electronic J. Diff. Equations. 1998;1998(32):1–12.
• Fleming WH, Rishel RW. Deterministic and Stochastic Optimal Control. New York: Springer Verlag; 1975.
• Lukes DL. Differential Equations: Classical to Controlled. New York: Academic Press; 1982.
• Pontryagin LS, Boltyanskii VG, Gamkrelize RV, Mishchenko EF. The Mathematical Theory of Optimal Processes. New York: Wiley; 1967.
• Nuño M, Chowell G, Gumel AB. Assessing the role of basic control measures, antivirals and vaccine in curtailing pandemic influenza: Scenarios for the US, UK and the Netherlands. J. R. Soc. Interface. 2006;4(14):505–521. doi: 10.1098/rsif.2006.0186. [PubMed]
Articles from BMC Public Health are provided here courtesy of
BioMed Central

 PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's Canada Institute for Scientific and Technical Information in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers.