Home | About | Journals | Submit | Contact Us | Français |

**|**BMC Public Health**|**v.11(Suppl 1); 2011**|**PMC3317585

Formats

Article sections

- Abstract
- Introduction
- Model framework
- The optimal control model
- Numerical simulations
- Conclusion
- Competing interests
- Authors contributions
- References

Authors

Related links

BMC Public Health. 2011; 11(Suppl 1): S5.

Published online 2011 February 25. doi: 10.1186/1471-2458-11-S1-S5

PMCID: PMC3317585

Jean M Tchuenche,^{1} Nothabo Dube,^{2} Claver P Bhunu,^{2} Robert J Smith?,^{}^{3} and Chris T Bauch^{1}

Jean M Tchuenche: ku.oc.oohay@shtamoib_tmj; Nothabo Dube: moc.liamg@57ehtmam; Claver P Bhunu: moc.liamg@unuhbpc; Robert J Smith?: ac.awattou@34htimsr; Chris T Bauch: ac.hpleugou@hcuabc

Mathematical Modelling of Influenza

Jean M Tchuenche and Robert Smith?

1471-2458-11-S1-info.pdfCopyright ©2011 Tchuenche et al; licensee BioMed Central Ltd.

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article has been cited by other articles in PMC.

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.

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.

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.

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.

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 *m _{I}* > 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

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.

The disease-free equilibrium of the system is given by

The endemic equilibrium of the system is given by

It satisfies and

where *h*_{1}(*I*) = *m*_{1} + *I*, *h*_{2}(*I*) = (*θ* + *µ* + *β*_{1}*I*)*h*_{1}(*I*) – *β*_{2}*I*^{2}, *h*_{3}(*I*) = (*θ + ω* + (1 – *γ*)*β*_{1}*I*)*h*_{1}(*I*) – *β*_{3}(1 – *γ*)*I*^{2}. 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, *R _{v}*, 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

and

The basic reproductive ratio is the spectral radius *ρ*(*FV*^{–1}) which is

(5)

**Lemma 1***The disease-free equilibrium E _{v}*

**Proof.** The Jacobian of the system evaluated at *E _{v0}* is given by

The eigenvalues of *J*_{EV0} 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 *R _{v}* < 1 and, consequently, all the eigenvalues of the Jacobian matrix above have negative real part. Thus, the disease-free equilibrium is locally asymptotically stable.

We adopt the method of Castillo-Chavez *et al,*[33] and we rewrite the set of model equations in the form

with *G*(*X _{G}*,0) = 0.

(H1) For is globally asymptotically stable.

(H2) *G*(*X _{G}*,

If the above two conditions are satisfied, then the following theorem holds.

**Theorem 2 ***(Castillo-Chavez et al,*[33]*): The fixed point**is a globally stable equilibrium of (2.28) provided that R _{v}* < 1

Therefore, *E _{v}*

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 *u _{v}* and

(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 *B*_{1} and *B*_{2} 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 *u _{v}*(

where *U* = {(*u _{v}*,

The existence of an optimal control can be obtained by using a result by Joshi [36] and Fister *et al.*[37].

**Theorem 3***Consider the control problem with the system of Equations (4.1)-(4.4). There exists an optimal control**such 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 *c*_{2} – *c*_{1}(*|u _{v}|^{k}* +

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, *c*_{1}, *c*_{2} > 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 *w*_{11}(*t*) ≥ 0, *w*_{12}(*t*) ≥ 0 are penalty multipliers satisfying *w*_{11}(*t*)(*a*_{11} – *u _{v}*(

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}*[

To determine an explicit expression for our controls , (without *w*_{11}*, w*_{12}, *w*_{21}, *w*_{22}), 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 *w*_{12} ≥ 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 *w*_{22} ≥ 0.

3. On the set . Hence

or

Combining all the three sub-cases in a compact form gives

(13)

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) = *S*_{0}, *I*(0) = *I*_{0}, *V*(0) = *V*_{0}, *R*(0) = *R*_{0} and *λ _{i}*[

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.

The general model with pulse vaccination is given as

for *t* ≠ *t _{k}*, where

Δ*S* = –*θS*

Δ*V* = *θS*

when *t* = *t _{k}.* 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, *I*_{crit}*.*

For *I* <*I*_{crit}, this model would look like

For *I* >*I*_{crit}, 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* <*I*_{crit}, we assume that there is no mixing, but rather that new infectives arrive impulsively into the system at fixed times *t _{k}* and in numbers

For *I* >*I*_{crit}, fear of the disease keeps susceptibles from mixing with the infected, but the vaccinated will.

Thus *β*_{4} = *β*_{6} = 0. Since *I ^{i}*

The model then becomes

(14)

(15)

(16)

(17)

for *I* <*I*_{crit} and

(18)

(19)

(20)

for *I* >*I*_{crit}.

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 *t _{k}*

where is the value immediately after the *k*th impulse. Then, since the period is constant, we have

This is a recursion relation with solution

Consequently,

Thus, if *m*^{+} >*I*_{crit}, 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 *Î* >*I*_{crit}. 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*^{+} >*I*_{crit} will, as a result of the media, instead equilibrate at a much larger value *I* >*m*^{+} >*I*_{crit}. 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).

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 *B*2 (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 *B*1 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 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.

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.

The authors declare that they have no competing interests.

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.

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] [Cross Ref]
- 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. [Cross Ref]
- Liu Y, Cui J. The impact of media coverage on the dynamics of infectious disease. Int. 2008;1:65–74. doi: 10.1142/S1793524508000023. [Cross Ref]
- 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. [Cross Ref]
- Ma R. Media, Crisis, and SARS: An Introduction. Asian Journal of Communication. 2005;15(3):241–246. doi: 10.1080/01292980500260656. [Cross Ref]
- 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. [Cross Ref]
- 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. [Cross Ref]
- 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. [Cross Ref]
- 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] [Cross Ref]
- Samba E, Nkrumah F, Leke R. Getting polio eradication back on track in Nigeria. N. 2004;350:645–646. doi: 10.1056/NEJMp038210. [PubMed] [Cross Ref]
- Gross L. A Broken Trust: Lessons from the Vaccine-Autism Wars. PloS Biology. 2009;7(5):e1000114. doi: 10.1371/journal.pbio.1000114. [PMC free article] [PubMed] [Cross Ref]
- 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] [Cross Ref]
- 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. [PMC free article] [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. [Cross Ref]
- 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. [Cross Ref]
- 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. [Cross Ref]
- 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] [Cross Ref]
- 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. [PMC free article] [PubMed] [Cross Ref]
- 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] [Cross Ref]
- 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. [Cross Ref]
- 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] [Cross Ref]
- 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] [Cross Ref]
- Joshi HR. Optimal Control of an HIV Immunology Model. Optim. 2002;23:199–213. doi: 10.1002/oca.710. [Cross Ref]
- 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. [PMC free article] [PubMed] [Cross Ref]

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 national science library 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. |