- Research article
- Open Access
- Open Peer Review
Modelling the effects of media during an influenza epidemic
BMC Public Healthvolume 14, Article number: 376 (2014)
Mass media is used to inform individuals regarding diseases within a population. The effects of mass media during disease outbreaks have been studied in the mathematical modelling literature, by including ‘media functions’ that affect transmission rates in mathematical epidemiological models. The choice of function to employ, however, varies, and thus, epidemic outcomes that are important to inform public health may be affected.
We present a survey of the disease modelling literature with the effects of mass media. We present a comparison of the functions employed and compare epidemic results parameterized for an influenza outbreak. An agent-based Monte Carlo simulation is created to access variability around key epidemic measurements, and a sensitivity analysis is completed in order to gain insight into which model parameters have the largest influence on epidemic outcomes.
Epidemic outcome depends on the media function chosen. Parameters that most influence key epidemic outcomes are different for each media function.
Different functions used to represent the effects of media during an epidemic will affect the outcomes of a disease model, including the variability in key epidemic measurements. Thus, media functions may not best represent the effects of media during an epidemic. A new method for modelling the effects of media needs to be considered.
Influenza causes annual epidemics and occasional pandemics, which have claimed millions of lives throughout history. In the past century four worldwide pandemic outbreaks of influenza have occurred: 1918, 1957, 1977 and 2009, [1, 2]. According to the Public Health Agency of Canada, inter-pandemic (or seasonal) influenza affects approximately 20,000 Canadians, with approximately 2,000 to 8,000 deaths annually . In the USA, it has been reported that flu-associated deaths can range from 3,000 to 49,000 individuals per year .
Mass media can affect disease transmission during an influenza epidemic or pandemic. Attention to health news has been increasing in importance during the last few decades, and thus, media reports can play an important role in defining health issues, since they serve as a major source of information and are able to incite changes in behaviour in the public . Individual response to a disease threat depends on risk perception that is gained largely through information reported by governments to mass media: number of infections, hospitalizations and deaths, as provided by the government [6, 7].
We have recently seen the use of mass media reports during two infectious disease outbreaks. The first novel infectious disease of the twenty-first century was SARS. It had distinct features such as rapid spatial spread and self-control, and vast media coverage [6, 7] that used to inform the public, provide a means of contract tracing, and advise isolation of exposed individuals.
Media coverage was substantial during the H1N1 epidemic in 2009 as well, which may have had an effect on the transmission of influenza by promoting social distancing and self-isolation . The media coverage of this influenza pandemic was widespread, with an increased sense of urgency since this influenza strain was related to the 1918 pandemic strain that caused approximately 50 million deaths worldwide .
Mathematical modelling has been used to study the effect of mass media on epidemics by employing the well-known Susceptible-Exposed-Infectious-Recovered (SEIR) model and various extensions [6–14]. In these studies, mass media has been incorporated using different, but qualitatively similar, functions that directly affect disease transmission and susceptibility. In general, the chosen functions are decreasing functions with respect to the current number of infected individuals in the population. However, the choice of function seems to be arbitrary. It is possible that the choice of function can change study results. For example, public health officials could be interested in epidemic measurements such as the peak number of infection, peak time, total number of infections and end of epidemic, which are all directly related to important public health resources (i.e. number of hospital bed, antiviral stockpile, vaccination doses). These key measurement may vary depending on the chosen media function. Therefore, a sensitivity analysis on these functions is required.
In addition to that mentioned above, there is a further drawback of the previous studies which include media functions, in that, deterministic systems of ordinary differential equations are employed. Deterministic models can describe the mean behaviour of an epidemic, but information surrounding any variability in key epidemic measurements cannot be made. A stochastic model is well suited to this task. Various methods of introducing stochasticity into disease models exist: [15–19]. Agent-Based Monte Carlo (ABMC) simulations, provide a way in which individuals with certain disease characteristics can be tracked in a population. This method also provides flexibility, as changes in biological assumptions can be easily incorporated, which are difficult to include in other methods.
In the sections that follow we give an overview of the functions used to describe media in the disease modelling literature. The functions are then incorporated into a standardized SEIR model, and model results are compared. A stochastic agent-based Monte Carlo (ABMC) simulation is then introduced, and is employed to study the variability within an epidemic depending on the media function chosen. A sensitivity analysis is also completed in order to determine the importance of certain model parameters on various epidemic outcomes for each media function.
where I is the number of infectious individuals in a population, γ is the recovery rate, and p i , i = 1, 2, 3 is a parameter used to represent media effects in these functions. In general, these functions are decreasing functions of I, which represents the fact that, as the number of infections increases in a population, and is reported by mass media, susceptible individuals will practice social distancing or control measures, which decreases susceptibility. Comparing the three functions, we see that, for a given I and γ, p 1 and p 2 can be written in terms of p 3:
These equations demonstrate that, if p 3 is increased, then p 1 and p 2 must decrease to have the functions remain the same value at a chosen I c .
In Figure 1, we plot Functions (1)-(3) (dotted, dashed, solid lines respectively) for media parameters p i , i = 1, 2, 3 (dotted, dashed, solid lines) with γ = 1/4 and I c = 300 (top) and I c = 1000 (bottom). Note that the functions are equal when I = 300 (top) and 1000 (bottom). Here, p 3 = 10 (left panel), 100 (middle panel), and 1000 (right panel), and p1 and p2 are determined using Equations (4) and (5).
To compare the effects of different media Functions (1)-(3), we must choose a standardized model. For the purposes of our study, we have chosen the basic Susceptible-Exposed-Infectious-Recovered (SEIR) model with a constant population size N:
Briefly, susceptible individuals (S) are infected by infectious individuals (I) with rate β and become exposed (E). Transmission can also be reduced as determined by f(I,p). Exposed individuals become infectious with rate σ, and infectious individuals recover with rate γ. The initial conditions and parameters values for this study can be found in Table 1.
Agent-based Monte Carlo simulation
For further comparison of the SEIR model with different media functions, we utilize an Agent-Based Monte Carlo (ABMC) simulation to capture the variability in the epidemic infection curve (that cannot be determined using a system of deterministic equations like that of System (6)). There are various ways of developing an ABMC simulation. We have employed a previous method as developed by Heffernan and Wahl, [16, 22]. Briefly, the ABMC simulation moves forward in time following event times: the next time that an individual changes state within the system. Agents in each of the susceptible, exposed, infectious and recovered compartments are assigned event times corresponding to Table 1 when they are introduced into the simulation, for each event that allows such an individual from that compartment to change state. For example, exposed individuals are assigned a time to become infectious, and infectious individuals are assigned a time to recovery and a time to infect a susceptible. The minimum event time determines the next event to occur in the simulation. The event is then carried out, and the next event is determined. To compare to the SEIR system as described above, we assume exponential distributions for all parameters. Table 1 lists the parameter values of the SEIR model and the mean value of the parameter distributions for the ABMC simulation. Note that, in general, the means of the lifetime distributions in the ABMC are simply the reciprocals of the rates used in System (6). However, so that infection event times always depend on the current size of the susceptible population S(t) (similar to what is assumed in System (6)), the infection time distribution mean is continuously updated. Figure 2, shows the progression of an individual through the epidemic.
Results and discussion
Basic reproductive ratio
The basic reproductive ratio R 0 is defined as the number of newly infectious individuals produced by one infectious individual in a totally susceptible population. For a description of different methods for calculating R 0 see . For System (6),
where N = S 0 is the total population size of susceptibles at time zero. Note that the calculation of R 0 is independent of f(I,p).
Comparison of media functions
Figure 3 plots System (6) using Media Functions (1)-(3) for p 3 = 10, 100, and 1000 (solid line), with and p 1 (dotted line) and p 2 (dashed line) determined by Equations (4) and (5), and γ = 1/4 and I c = 300. For comparison, the standard SEIR model with no media effect is also shown (dash-dotted line). This figure demonstrates that mass media will affect the epidemic curve. It also demonstrates that the epidemic curve varies depending on the media function used.
For a further comparison of the media functions, we determine the values of the basic reproductive ratio R 0 and the media constant p of f(I,p) based on epidemic data generated from our models. For an example, Table 2 lists the values of β and p determined through a fitting routine (in MATLAB - various routines were used and similar results were obtained) over the first 100 days of an epidemic, with data generated using Media Function (3) with σ and γ values listed in Table 1 and p 3 = 1000. Figure 4 plots the resulting epidemic curves of System (6) using Media Functions (1)-(3) (dotted, dashed and solid lines) and no media (dash-dotted line). Table 2 and Figure 4 demonstrate that even though all of the models have a very similar shape and basic reproductive ratio R 0 at the beginning of the epidemic, the resulting epidemic curves can still vary drastically depending on the media function chosen. Similar observations were made using different numbers of days of data used to fit β and p, and different Media Functions to generate the data (not shown). Therefore, it can be concluded that key epidemic measurements – peak number of infections, time of peak, end of epidemic, and total number of infections – will also vary depending on the Media Function chosen.
System (6) is useful in describing the mean behaviour of an epidemic, but it is unable to provide estimates of variation of important public health measures. A stochastic simulation lends itself well to demonstrating variation within an epidemic. Here, we employ an Agent-Based Monte Carlo (ABMC) simulation. Figure 5 shows 100 simulation runs of the ABMC simulation when no media is involved (top), and when Media Functions (1)-(3) are used (second to bottom rows). Each simulation run (gray line) is shown, as well the mean of the simulation runs (solid line) and solution of the System (6) with the corresponding media functions (dashed line). Note that the mean (solid line) and solution to System (6) are in agreement. This figure demonstrates that variability in the key epidemic measurement occurs. Table 3 lists the mean and standard error of each epidemic measurement when no media is considered (section a), and Media Functions (1)-(3) are included in the simulation (Sections b-d). In each section we list the simulation mean and standard error (top row), and the solution of System (6).
Within public health settings, a goal is to identify key characteristics of an epidemic that drive the infection with the hope of determining public health measures that can be implemented so that control or eradication of the pathogen can be achieved. By performing sensitivity analysis on System (6) with Equations (1)-(3), parameters that most affect epidemic outcomes for each media function can be identified, informing policy makers so that appropriate public health measures can be put into place. To conduct the sensitivity analysis, we use Latin Hypercube Sampling (LHS) and partial rank correlation coefficients (PRCC) .
We first conduct the sensitivity analysis with I c , p 3 and γ constant to directly compare the media functions. Figures 6 and 7 show the PRCC values determined for peak magnitude when I c = 300 (Figure 6) and 1000 (Figure 7) for Functions (1)-(3) (top to bottom) when p 3 = 10, 100 and 1000 (left, middle and right columns respectively). These figures demonstrate that for each set of parameters the PRCC value is similar. This is true for each key epidemic measurement - the PRCC values are similar between different values of p 3 and I c . However the PRCC values between key measurements are different i.e. the PRCC values are different if comparing peak magnitude to epidemic end time (not shown).
For further sensitivity analysis, all parameters are allowed to vary. Figure 8 shows the result of the sensitivity analysis for all four outcomes key to public health: peak number of infections, peak time, epidemic end time, and total number of infections. The results demonstrate that the SEIR model is more (or less) sensitive to certain parameters depending on what Mass Media Function (1)-(3) is used. For example, changes in β have a large effect on peak magnitude when Function (1) is used, but it has little to no effect when Functions (2-3) represent the effects of media. Variable p 3 also has a large effect on peak magnitude when media is represented by Function (1). It has a similar effect when Function (2) is used, but it has no effect on the system utilizing Function (3).
The sensitivity analysis indicates that models that include mass media influence will greatly depend on different parameters, depending on the media function chosen. This makes it very difficult for policy makers to determine an effective public health intervention strategy. This also explains the very different epidemic curves produced by System (6) when different media functions are employed, notwithstanding the similarities in the media functions when plotted at a specific level of media.
Technology and media play an increasing role in daily life. Mass media that is transmitted via technological media can therefore be used to inform the public during pandemics and epidemics. An understanding of the effects of media during an epidemic or pandemic threat can aid in the development of public health policy. Of particular interest to public health are the effects of media on key epidemic measurements - peak magnitude of infection, time of peak, end of epidemic, and total number of infections.
Mathematical modelling has been used to study the effect of media on epidemics by employing functions in the transmission terms of mathematical models [6–14]. A survey of the literature identified three functions that have been utilized to represent media in disease modelling [6–14]. We have conducted a comparison of these functions to determine the effects of media function on key epidemic measurement and variability within these measurements. We first demonstrated that by including mass media in System (6) the peak magnitude of the epidemic and the total number of infections would decrease. We also determined that the time to peak and the end of the epidemic would also occur earlier. However, we also demonstrated that, although the functions are similar in shape and magnitude, the resulting epidemic curve can be quite different (Figure 5). Therefore, the key epidemic measurements that are used to inform public health policy will be different. Furthermore, we demonstrated that variability in the key epidemic measurements also depends on the media function chosen (Figure 5 and Table 3).
A sensitivity analysis on System (6) with the different media functions was also conducted. Obtained from this analysis was the insight that some parameters are important for some outcomes and not for others. We can conclude that due to the different fixed functions resulting in very different epidemic behaviours, there is no clear control strategy present. Also, as a reult of the different behaviours from the different media functions, we are unsure as to which is the best function to use to model mass media. This suggests that a function representing media may not be the best course for modelling the effects of media during an epidemic. We suggest that perhaps, a separate model compartment representing media reports, such as those incorporated into surveillance data  could be used. Development of such a model is a course for future work.
Mass media reports can affect social behaviour, that ultimately, affects transmission of disease. However, an individual’s response to a media stimulus will wane over time [13, 25–29]. Models that employ a media function such as those studied here are difficult to modify to involve a waning response to a media stimulus over the age of an epidemic. This can be incorporated easily into a model whereby media is represented as a model compartment. A study of the effects of ‘waning media’ is a course for future work.
The current study has employed a simple SEIR model that includes three parameters. This model implicitly assumes that individuals mix at random in the population, the age and sex of individuals is unimportant, and that the population size remains constant over the epidemic. These assumptions are not true reflections of reality. An interesting direction for future work is to consider the effects of mass media reports on individual decision making strategies, and study how mass media can impact the contact structure of a population network .
The simple SEIR model employed in this study was used to study one season of influenza. Typically, seasonal forcing is used to model multiple seasons of influenza. Inclusion of mass media reports and media waning in a model of multiple seasons with seasonal forcing is an interesting direction for study.
Brauer F: Mathematical epidemiology is not an oxymoron. BMC Pub Health. 2009, 9 (Suppl1): S2-
Neumann G, Noda T, Kawaoka Y: Emergence and pandemic potential of swine-origin H1N1 influenza virus. Nature. 2009, 459: 931-939. 10.1038/nature08157.
Public Health Agency of Canada: Highlights from the Canadian Pandemic Influenza Plan for the Health Sector. 2006, Canada: Public Health Agency of Canada
Centers for Disease Control and Prevention: Estimating Seasonal Influenza-Associated Deaths in the United States: CDC Study Confirms Variability of Flu. 2011, Atlanta, GA, USA: Centers for Disease Control and Prevention
Viswanath K, Ramanadhan S, Kontos EZ: Mass media. Macrosocial Determinants of Population Health. Chap.13. Edited by: Galea S, Ann Arbor. 2007, MI, USA: Springer, 275-294.
Liu R, Wu J, Zhu H: Media/psychological impact on multiple outbreaks of emerging infectious diseases. Comput Math Methods Med. 2007, 8: 153-164. 10.1080/17486700701425870.
Tchuenche JM, Dube N, Bhunu CP, Smith RJ, Bauch CT: The impact of media coverage on the transmission dynamics of human influenza. BMC Pub Health. 2011, 11 (Suppl1): S5-
Pang J, Cui J-a J-a: An SIRS epidemiological model with nonlinear incidence rate incorporating media coverage. Second International Conference on Information and Computing Science. 2009, USA: IEEE, 116-119.
Cui J-a, Sun Y, Zhu H: The impact of media on the control of infectious diseases. J Dyn Differential Equations. 2008, 20: 31-53. 10.1007/s10884-007-9075-0.
Cui J-a, Tao X, Zhu H: An SIS infection model incorporating media coverage. Rocky Mt J Math. 2008, 38: 1323-1334. 10.1216/RMJ-2008-38-5-1323.
Li Y, Ma C, Cui J: The effect of constant and mixed impulsive vaccination on SIS epidemic models incorporating media coverage. Rocky Mot J Math. 2008, 38: 1437-1455. 10.1216/RMJ-2008-38-5-1437.
Schweitzer F, Mach R: The epidemics of donations: Logistic growth and power laws. PLoS ONE. 2008, 3: e1458-10.1371/journal.pone.0001458.
Sun C, Yang W, Arino J, Khan K: Effect of media-induced social distancing on disease transmission in a two patch setting. Math Biosci. 2011, 230: 87-95. 10.1016/j.mbs.2011.01.005.
Xiao D, Ruan S: Global analysis of an epidemic model with nonmonotone incidence rate. Math Biosci. 2007, 208: 419-429. 10.1016/j.mbs.2006.09.025.
Allen LJS: An Introduction to Stochastic Epidemic Models. Mathematical Epidemiology. Lecture Notes in Mathematics. Edited by: Brauer F, van den Driessche P, Wu J. 2008, Lubbock, TX, USA: Springer, 81-130.
Heffernan JM, Wahl LM: Monte Carlo estimates of natural variation in HIV infection. J Theor Biol. 2005, 236: 137-153. 10.1016/j.jtbi.2005.03.002.
Keeling MJ: Metapopulation moments: coupling, stochasticity and persistence. J Anim Ecol. 2000, 69: 725-736. 10.1046/j.1365-2656.2000.00430.x.
Keeling MJ: Multiplicative moments and measures of persistence in ecology. J Theor Biol. 2000, 205: 269-281. 10.1006/jtbi.2000.2066.
Krishnarajah I, Cook A, Marion G, Gibson G: Novel moment closure approximations in stochastic epidemics. Bull Math Biol. 2005, 67: 855-873. 10.1016/j.bulm.2004.11.002.
Pourbohloul B, Ahued A, Davoudi B, Meza R, Meyers LA, Skowronski DM, Villasenor I, Galvan F, Cravioto P, Earn DJD, Dushoff J, Fisman D, Edmunds WJ, Huper N, Scarpino SV, Trujillo J, Lutzow M, Morales J, Contreras A, Chavez C, Patrick DM, Brunham RC: Initial human transmission dynamics of the pandemic (H1N1) 2009 virus in North America. Influenza Other Respi Viruses. 2009, 3: 215-222. 10.1111/j.1750-2659.2009.00100.x.
Mostaço-Guidolin LC, Pizzi NJ, Moghadas SM: A classical approach for estimating the transmissibility of the 2009 H1N1 pandemic. Can Appl Math Q. 2011, 19: 185-194.
Heffernan JM, Smith RJ, Wahl LM: Perspectives on the basic reproductive ratio. J R Soc Interface. 2005, 2: 281-293. 10.1098/rsif.2005.0042.
Blower SM, Dowlatabadi H: Sensitivity and uncertainty analysis of complex models of disease transmission: An HIV model, as an example. Int Stat Rev. 1994, 62: 229-43. 10.2307/1403510.
Mykhalovskiy E, Weir L: The global public health intelligence network and early warning outbreak detection. Can J Public Health. 2006, 97 (1): 42-44.
Fanti KA, Vanman E, Henrich CC, Avraamides MN: Desensitization to media violence over a short period of time. Aggress Behav. 2009, 35: 179-187. 10.1002/ab.20295.
Griffiths MD, Shuckford GLJ: Desensitization to television violence: a new model. New Ideas Psychol. 1989, 7: 85-89. 10.1016/0732-118X(89)90039-1.
Kinnick KN, Krugman DM, Cameron GT: Compassion fatigue: Communication and burnout towards social problems. J Mass Commun Q. 1996, 73: 687-707.
McBurney DH, Balaban CD, Christopher DE, Harvey C: Adaptation to capsaicin within and across days. Physiol Behav. 1997, 61: 181-190. 10.1016/S0031-9384(96)00366-6.
Smith KC, Rimal RN, Sandberg H, Storey JD, Lagasse L, Maulsby C, Rhoades E, Barnett DJ, Omer SB, Links JM: Understanding newsworthiness of an emerging pandemic: International newspaper coverage of the H1N1 outbreak. Influenza Other Respi Viruses. 2012, 7: 847-853.
Wang Y, Cao J, Jin Z, Zhang H, Sun G-Q: Impact of media coverage on epidemic spreading in complex networks. Physica A: Stat Mech Appl. 2013, 392 (23): 5824-5835. 10.1016/j.physa.2013.07.067.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/14/376/prepub
The authors would like to thank Nicholas Geard, Amy Greer, James McCaw, Seyed Moghadas, Jianhong Wu and Huaiping Zhu for many helpful discussions. The authors would also like to thank the reviewers for their comments and suggestions that have added value to this study The work was funded by an Ontario Early Researcher Award, Mathematics of Information Technology and Complex Systems (MITACS), and the Natural Science and Engineering Research Council of Canada (NSERC). Simulations were conducted on Sharcnet.
The authors declare that they have no competing interests.
MSC performed the literature review. JMH and MSC developed the models, conducted the analysis and wrote the paper. Both authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.