Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Local risk perception enhances epidemic control

  • José L. Herrera-Diestra ,

    Contributed equally to this work with: José L. Herrera-Diestra, Lauren Ancel Meyers

    Roles Conceptualization, Formal analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing

    jdiestra@ictp-saifr.org

    Affiliations ICTP South American Institute for Fundamental Research, São Paulo, Brazil, IFT-UNESP, São Paulo, Brazil, CeSiMo, Facultad de Ingeniería, Universidad de Los Andes, Mérida, Venezuela

  • Lauren Ancel Meyers

    Contributed equally to this work with: José L. Herrera-Diestra, Lauren Ancel Meyers

    Roles Conceptualization, Investigation, Methodology, Writing – original draft, Writing – review & editing

    Affiliation Department of Integrative Biology, The University of Texas at Austin, Austin, Texas, United States of America

Abstract

As infectious disease outbreaks emerge, public health agencies often enact vaccination and social distancing measures to slow transmission. Their success depends on not only strategies and resources, but also public adherence. Individual willingness to take precautions may be influenced by global factors, such as news media, or local factors, such as infected family members or friends. Here, we compare three modes of epidemiological decision-making in the midst of a growing outbreak using network-based mathematical models that capture plausible heterogeneity in human contact patterns. Individuals decide whether to adopt a recommended intervention based on overall disease prevalence, the proportion of social contacts infected, or the number of social contacts infected. While all strategies can substantially mitigate transmission, vaccinating (or self isolating) based on the number of infected acquaintances is expected to prevent the most infections while requiring the fewest intervention resources. Unlike the other strategies, it has a substantial herd effect, providing indirect protection to a large fraction of the population.

Introduction

As outbreaks emerge, public health agencies often implement a variety of pharmaceutical and non-pharmaceutical interventions to prevent epidemic expansion, including vaccination and medical prophylaxis, school closures and other social distancing measures, and information campaigns to promote awareness, hygienic precautions and voluntary isolation [14]. However, such measures require population adherence and are often hindered by failure to take recommended actions [5]. Around the globe, for example, seasonal influenza vaccine coverage falls significantly below the 75% baseline recommended by the World Health Organization, but varies widely between countries and across age groups [6]. In the USA, 2015-2016 uptake was only 59.3% in children and 41.7% in adults [7]. For measles, routine childhood vaccination is declining in Texas and other areas of the United States where personal belief and other non-medical vaccination exemptions are allowed [810]. Parental decision-making regarding childhood vaccines is complex and context dependent [11], but likely influenced by false claims regarding vaccine safety, low perceived risks of infectious diseases, and other forms of misinformation from the “anti-vaxxer” movement [8, 1113]. Recently, there have been calls for a special government commission on vaccine safety, despite overwhelming scientific consensus that vaccines are both safe and effective [12, 1416].

As outbreaks unfold, people can take a variety precautionary measures to avoid infection, including immunization and social distancing [1, 1719]. They often judge personal risk based on their impressions of overall disease prevalence and severity [2, 2023]. When infection risk appears low, small risks of adverse effects from the vaccine can seem relatively important and cause vaccine coverage to drop below levels required to control transmission. A variety of other factors can influence the perceived utility of disease prevention, including epidemiological news and rumors, costs of vaccination and other control measures, trust in health professionals, government agencies, media and non-official information sources, as well as societal pressure to ensure the health of one’s children [2429].

Studies have shown that media reports about outbreaks that specify numbers of cases, hospitalizations or deaths, can influence avoidance behavior and contact patterns at both individual and community levels. In some cases, oversimplified or erroneous media reports regarding flawed vaccines can trigger panic and increases in vaccine hesitancy [3033]. For both seasonal and pandemic influenza, such interactions between vaccination decision-making and transmission dynamics can profoundly shape the course of epidemics [20, 3436].

Individual intervention decisions can have far-reaching effects. For example, vaccination protects not only the immunized individual, but also social contacts who might have been infected by the individual. Social distancing decisions may break chains of transmission by protecting the decision-maker and more generally disrupting social dynamics. Following [37, 38], we refer to this indirect protection as a herd effect. Although previously equated to the reduction in incidence in the unimmunized population [38], we quantify the herd effect of an intervention effort by estimating the number of infections averted per vaccine administered (or per individual social distancing action). The general phenomenon in which individual intervention actions reduce the risk of infection to others has also been called vaccination efficiency, vaccination effectiveness, herd immunity, and indirect protection [3740].

The magnitude of the herd effect will critically depend on contact patterns [41, 42]. Measures taken by gregarious individuals may have greater immediate benefits than those taken by solitary individuals, with downstream epidemiological consequences modulated by the full social network [43, 44]. Contact patterns may also influence the decision-making process itself, by constraining epidemiological perspectives. When gauging infection risk, individuals may consider global information (e.g., from news media) or local first-hand encounters with disease (e.g., infected acquaintances, friends or family members) [30, 31]. While traditional compartmental models assume homogeneity in both epidemiological risks and intervention benefits, network-based models provide a tractable framework for studying the complex interplay between contact networks, intervention decision making and disease transmission [34, 4551].

Here, we investigate the epidemiological impacts of different decision paradigms using a network-based SIR epidemic model, in which individuals also make vaccination or social distancing choices based on their perceived epidemiological risks. Depending on the decision model, they estimate either overall disease prevalence, their number of infected social contacts, or their fraction of infected social contacts. When the perceived threat is sufficiently high, they take a measure that immediately affords full protection for the duration of the epidemic. We compare the efficacy of these three different paradigms across a range of diseases in a realistically heterogeneous network, and show that the most naive model–simply counting one’s infected contacts–affords the most epidemiological protection using the least amount of resources (e.g., vaccinations or economic costs associated with social distancing).

Materials and methods

We simulate the spread of an infectious disease in a network (population) with an exponential degree distribution–as estimated for typical urban populations [46, 52]–using a susceptible-infected-recovered (SIR) chain-binomial model [41, 53] with an additional immunized state. At each time step, individuals decide whether or not to take an instantaneously protective action to avoid infection, based on their perceived risk of infection, as defined by the given local or global decision model. We assume that there are sufficient resources to immediately protect any willing individuals.

Contact network

We model contact patterns in the population using an exponential network with N = 10000 nodes and mean number of contacts μ = 10 [54], generated according to the configuration model [55], unless otherwise is specified. We assume that this network constrains both disease transmission and local risk perceptions, when individuals monitor infected social contacts. To evaluate the impact of network topology, we also analyze the SIR-intervention dynamics on a homogeneous random graph (all nodes have same degree) and Barabási-Albert scale-free network [56], with degree distributions constrained to achieve the same epidemic threshold as the focal exponential network. All three networks share where 〈k〉 and 〈k2〉 are the average degree and squared degree in the network, respectively.

Epidemic dynamics

We model SIR transmission dynamics of a flu-like disease using chain-binomial stochastic simulations [41, 53]. Epidemics begin at time t = 0 by infecting a single randomly chosen node in an otherwise completely susceptible population and terminate when there are no remaining infected individuals. Individuals remain infectious for l = 7 days before recovering with full immunity to future infection [57, 58]. Infected individuals transmit disease to each susceptible contact at a rate β. Immunized and recovered individuals are assumed to be fully resistant to infection. Results are averages over 500 simulations.

The basic reproduction number (R0) is the expected number of secondary cases when a single case of disease is introduced into a naive population, and is related to the epidemic growth rate. To study the impact of transmission rate on the vaccination-epidemiological dynamics, we consider R0 values ranging between one and ten, which spans the range for many common human pathogens, including influenza, Ebola, SARS, Pertussis, HIV/AIDS, Mumps, Rubella, Polio, Smallpox, Diphtheria, etc [52, 5966]. For each value of R0, we determine the corresponding β using [45] (1) where T = 1 − (1 − β)l is the transmissibility over the entire infectious period, also known as the secondary attack rate [67].

Immunization models

We assume that individuals make daily decisions regarding whether or not to take precautionary measures based on their perceived risk (Fig 1). In other words, they have a general sense of the contagiousness of a given disease based on past experience or conventional wisdom. We use T as a rough proxy for perceived risk. If and when they choose to take action, they instantly gain full resistance to infection for the remainder of the simulation. Although these models apply to any transmission preventing measure, we henceforth refer to the interventions as vaccinations.

thumbnail
Fig 1. Three decision strategies.

Individuals decide to vaccinate based on one of three risk measures: The number of infected contacts (local count), the fraction of infected contacts (local prevalence), or the overall fraction infected (global prevalence). In this example, six of the 24 nodes are infected, yielding a global prevalence of 0.25. The white node towards the top has a single contact that happens to be infected; the white node towards the bottom has two of its five contacts infected.

https://doi.org/10.1371/journal.pone.0225576.g001

We model three different individual decision strategies in which individuals consider either the disease states of their direct social contacts or the global situation, perhaps gleaned through news or social media. Let vX(t) denote the willingness of a individual to vaccinate under strategy X at time t.

Local decision strategies.

In the first model, local prevalence, individuals assess infection risk by tracking the fraction of their social contacts that are currently infected. The probability that a susceptible individual i vaccinates at time t is given by (2) where ηi(t) is the number of neighbors of i that are infected at time t, ki is the total number of neighbors (degree) of i, and 〈k〉 is the average degree of the network.

In the second model, local count, individuals track their number rather than proportion of infected neighbors, and decide to vaccinate according to (3)

Local prevalence is arguably a less plausible strategy than local count, given that the decisions require the additional knowledge of total number of contacts (degree) of each individual.

Global decision strategy.

The global prevalence model assumes that individuals base their vaccination decisions on the epidemiological state of the entire population, as given by (4) where I(t) is the total number of infected individuals in the population at time t and N is the size of the population. This assumes general knowledge of the evolving dynamics of the epidemic, perhaps through news, social media or public health messaging.

The mean degree (〈k〉) appears in the global prevalence and local prevalence as a normalizer. If node i has the average degree (ki = 〈k〉) and its local prevalence mirrors global prevalence (), then it will have the same probability of vaccinating across all three models.

In all three models, we assume that individuals will vaccinate with a probability equal to their perceived real-time probability of being infected. For example, if an individual’s perception of infection in the immediate future is 25%, then a precautionary measure will be taken with probability 0.25. The local count model comes closest to estimating actual risk of infection. Specifically, vlc is the probability that any currently infected contact will transmit disease to the focal node at some point during his or her infections period. This exactly estimates risk if all infected contacts were just infected and at the beginning of their infectious period, but overestimates risk if some are nearing recovery. (S1 Fig illustrates how vaccination decisions change under each of these models as disease prevalence increases.).

Herd effect

To assess the indirect and direct protection afforded by a given decision strategy D at a given R0, we calculate a quantity we call the herd effect, given by (5) where 〈C0(R0)〉 and 〈CD(R0)〉 are the expected total number infections in epidemics without vaccination and with vaccination decision strategy D, respectively, and 〈VD(R0)〉 is the expected total number of individuals vaccinated under D. We estimate these expected values by averaging over 500 simulations with the specified R0 and decision model. Barring extreme stochasticity, we expect H > 0 for any reasonably protective vaccine strategy. When H is between zero and one, more vaccines are given than infections averted, suggesting that vaccines may be mistimed or misplaced. This could happen, for example, if risk is underestimated early in the epidemic and overestimated late in the epidemic. An H near one indicates that approximately one infection is averted for every vaccine given. Note that this is an average, and does not necessarily mean that every vaccination prevents infection of the recipient. If each vaccine averts, on average, greater than one infection (H > 1), then the value of H corresponds to the level of indirect protection achieved by the decision strategy.

Results

The decision models yield distinct vaccine adoption and disease transmission dynamics (Fig 2). As disease begins to spread, individuals perceive increasing risks and vaccinate according to the decision model, thereby protecting themselves and interrupting potential chains of transmission to others. While all three strategies reduce the total number of infections, the local count strategy affords the greatest and most efficient protection of the three. Under the global prevalence strategy, perceived risk is homogeneous. As cases mount, the vaccination rate rises synchronously throughout the population, arguably resulting in too much too late vaccine coverage.

thumbnail
Fig 2. Disease and vaccination dynamics under different decision models.

Shading indicates the fraction of nodes in each state: susceptible (black), vaccinated (light gray), recovered (dark gray), infected (white). Columns corresponds to different strategies, as labeled above; rows correspond to R0 = 2 (top) and R0 = 5 (bottom).

https://doi.org/10.1371/journal.pone.0225576.g002

The local strategies avert more infections with fewer vaccinations than the global strategy. As epidemics unfold, risk is both heterogeneous and dynamic, with some portions of the network experiencing greater forces of infection than others. Local decision-making allows earlier detection and response to increasing personal risk, and prevents unnecessary vaccination in lower risk settings, both prior to and following epidemic waves. The local count strategy is more protective than the local prevalence strategy. By tracking the number rather than proportion of infected contacts, individuals more accurately assess the local force of infection. For example, compare a solitary individual with just two social contacts and a gregarious individual with 20. If they both have two infected contacts, then their risk of infection will be similar (assuming that time spent with each contact is sufficient for transmission). Under local count, their perceived risk and consequent vaccination probability will be identical; under local prevalence, the solitary individual will perceive higher risk (i.e., 100% of contacts infected) than the gregarious individual. Under all models, overall vaccine coverage increases as R0 increases, with the global prevalence achieving near universal coverage by R0 = 5.

The relative and absolute impacts of each strategy are remarkably robust to the transmissibility of the pathogen (Fig 3). Without vaccines, the expected epidemic size increases non-linearly with R0, reaching almost 100% by R0 = 10 (Fig 3A). All of the vaccine strategies avert a large and increasing fraction of cases, as R0 increases. In fact, the total epidemic size is non-monotonic, with slightly more expected infections around R0 = 4 than around R0 = 10. The local count strategy consistently yields the greatest protection, followed by local prevalence.

thumbnail
Fig 3. Epidemiological impacts of decision strategies change with R0.

(A) The proportion of individuals that are infected increases and then declines slightly as R0 increases, across all three decision strategies. (B) The proportion of individuals that choose to vaccinate initially increases sharply with R0 across all three strategies, but then declines for only the local count strategy. (C) Herd effect is the proportion of infections averted per vaccinated individual, and is highest for the local count strategy across all R0. All values are averages across 500 stochastic simulations.

https://doi.org/10.1371/journal.pone.0225576.g003

The decision models result in dramatically different vaccination rates, with the global prevalence strategy leading to near universal vaccination, consistently more than double the coverage produced by the local count strategy (Fig 3B). The population-level protection afforded by the local count strategy exhibits a non-trivial trend with R0 (Fig 3C). Between R0 = 1 and R0 = 2, its impact grows logarithmically from less than one infection averted per vaccinator to nearly two infections averted per vaccinator. Thereafter, the indirect benefits continue to grow slowly, reaching three infections averted per vaccinator when R0 = 9.

To explore the dynamic interactions between behavior and epidemiology in the three models, we consider individual nodes based on their degree (number of contacts). In general, the higher the degree of a node, the higher their risk for becoming infected and infecting others, and the greater the number of local infections they could potentially perceive when making vaccination decisions. Indeed, across all decision models, higher degree individuals vaccinate earlier, in terms of the fraction of the population infected at time of vaccination (Fig 4A, bottom). However, the fraction of individuals vaccinated in each degree class does not necessarily increase with degree (Fig 4A, top). Local count is the only strategy under which vaccination coverage monotonically increases with degree. Under global decision-making, coverage is inversely related to degree, and under local prevalence, coverage peaks for moderately connected individuals. By the time individuals choose to vaccinate under either of the two suboptimal strategies, their local risk of infection is already quite high (Fig 4A, middle), particularly for more gregarious individuals. Although the vaccinating individuals are immediately protected, comparable individuals (of the same degree class and local risk) who stochastically fail to make the same low probability vaccination decision are likely to become infected. Consequently, the risk of infection increases steeply with degree under all models except local count (Fig 4B). In a sensitivity analysis, we find that these qualitative results are robust to our assumptions about the efficacy of the vaccine, the time lag between an individual deciding to vaccinate and becoming protected against infection, the size of the network and the duration of the infectious period (S3 Fig).

thumbnail
Fig 4. Vaccination decisions vary with degree.

(A) For each degree class, we graph the proportion of individuals that vaccinate (top) along with the epidemiological situation at the time a node chooses to vaccinate in terms of the number of its infected neighbors (middle) and the overall disease prevalence in the population (bottom). Epidemiological risks—the chances of both becoming infected and infecting others—generally increase with degree. The local count strategy (light blue) is the only strategy for which the probability of vaccinating consistently increases with risk. Compared to the two other strategies, high degree individuals vaccinate earlier in terms of both local and global disease prevalence and, consequently, are less likely to become infected. (B) The number of individuals infected in each degree class under each decision strategy. The top of the curve indicates the number of individuals in each degree class (i.e., the underlying degree distribution) on a log scale. For each degree class, the stacked values indicate the expected number of individuals infected under the various strategies. The top of the gray area indicates the expected number of infections in the baseline scenario without vaccines; red indicates individuals expected to remain uninfected. Generally, local count has the lowest expected attack rate (light blue), followed by local prevalence (dark blue), and finally global prevalence. This ranking does not hold for the lowest degree individuals; instead, local prevalence has a lower expected attack rate than local count, as indicated by the dark blue line cutting through the light blue region. For all graphs, values are averages across 500 simulations, assuming R0 = 5.

https://doi.org/10.1371/journal.pone.0225576.g004

Finally, we consider the impact of the underlying contact network on the interplay between transmission and vaccination dynamics (Fig 5). We compare our focal exponential network to a uniform random network in which all nodes have the same degree and a Barabási-Albert scale-free network. The local count strategy robustly affords the most efficient population-level protection, averting the maximum number of infections (or nearly maximum in the case of the scale-free network and low R0) with the fewest vaccines.

thumbnail
Fig 5. Decision dynamics across social networks.

We compare three different network structures–homogeneous (top), exponential (middle), and scale-free (bottom)–and show the susceptible (black), infected (dark gray) and vaccinated (light gray) for two different values of R0. Results are averages across 500 stochastic simulations.

https://doi.org/10.1371/journal.pone.0225576.g005

Discussion and conclusions

Public health interventions and individual-level adherence decisions can profoundly influence the fate of unfolding epidemics. In this study, we assume that individuals have access to a fully protective measure, such as self-isolation, medical prophylaxis, or an immediately and completely efficacious vaccine. They continuously make real-time risk assessments and thereby decide whether or not to adopt the intervention, based on either direct knowledge of infected friends and family (number or fraction of infected social contacts) or indirect information about population-level prevalence, perhaps gleaned through news media.

Of the three decision models, global risk assessments prove least effective across a large range of disease scenarios (R0 ranging from one to ten). Nearly all non-infected individuals eventually vaccinate, yet the total cases more than double those occurring under the alternative strategies. There is a mismatch between risk and action. Risk is highly variable in time and space, given the heterogeneity of the underlying contact network and branching nature of transmission. Yet, the global model assumes that perceived risk and the consequent likelihood of adherence is homogeneous throughout the network, though variable in time. By the time global prevalence triggers wide-spread action, the highest risk individuals have already been exposed and the lowest risk individuals may still not, and might never, require protection.

In contrast, when individuals make decisions based on local risk assessments, the intervention efforts more closely track the epidemiological dynamics. Tallying infected contacts rather than estimating the fraction of infected contacts provides a more accurate indication of real-time risk and results in more efficient intervention. Assuming that all social contacts are equally likely to transmit disease, two out of three infected contacts carries the same immediate risk as two out of ten. The advantage of local risk assessment stems from two sources of variation in risk. First, disease transmission is an inherently local process in which risk aggregates around currently infected individuals. Second, this is magnified in realistically heterogeneous networks, by the concentration of risk around the center (most connected individuals) of the network.

Although several prior studies have also explored the epidemiological impacts of local and globally-informed vaccination decisions [4851], ours is the first to consider a decision-model based on counts rather than fractions of infected contacts and to systematically compare three different decision paradigms across a range of network structures. Massaro et al. model two networks–the epidemiological contact network through which disease spreads and the social network through which risk information spreads. The more similar the two networks, the greater the individual and population-level protection achieved by vaccination [49]. This is consistent with our finding that locally-sourced decisions provide greater protection than globally-sourced decisions. Bagnoli et al. compare the local fraction strategy across contact networks with different degree distributions, and likewise found that the herd effect is magnified by heterogeneity in degree. [48].

Given infinite resources, all three of the decision paradigms would markedly diminish an emerging outbreak. However, interventions may be constrained by limited supplies or lack of population access to medical countermeasures, such as vaccines or antimicrobials. Even social distancing measures, such as self-isolation, may be limited by economic necessity—the need to go to work, school or daycare—or care-giving obligations for extended family. While such limitations should be formally analyzed, our simple analysis suggests that the best paradigm for averting infections also requires the fewest resources. For example, for a flu-like R0 of two, compare the local count strategy, where individuals protect themselves as their number of infected friends and family climb, to the global strategy, where decisions are based on population prevalence. For every individual that takes action, almost two infections are averted under the local strategy whereas less than one infection is averted under the global strategy. Local counting results in far fewer total infections (3% versus 13%) while requiring far less intervention resources (23% versus 60% of individuals taking protective action).

Several studies suggest that immunizing or isolating interventions should target the most connected individuals in a population [4244, 52]. However, we rarely know the full contact network of a population. As proxies, we can target populations subgroups that tend to have high numbers of potentially disease-spreading contacts, such as young and school-aged children or health-care workers. We can also use biased sampling to identify highly connected individuals, such as the random acquaintance strategy in which random individuals are asked to name one of their social contact; individuals with more contacts are more likely to be named [6871]. In a sense, the winning paradigm of our study—counting infected contacts—similarly biases interventions efforts towards more connected parts of the network. The more connected one is, the more likely one is to have several infected contacts.

The model is intentionally simplistic, providing a best case scenario for each of the three strategies. We assume that resources are unlimited, protection is immediate and complete, and adherence probabilities perfectly mirror perceived risks. Furthermore, depending on the decision paradigm, individuals fairly accurately estimate the infectiousness of the disease, their number or fraction of infected social contacts, or the population average risk of infection. The model also assumes that individuals are short-sighted and make reactive decisions to avert immediate threat. We conjecture that the qualitative results of our analysis—the optimality of assessing risk based on the numbers of infected friends and family—are robust for a large class of ‘on-the-fly’ interventions that afford relatively rapid protection in the heat of an epidemic, but may not apply to preventative measures taken early in an outbreak or those with long efficacy lags. (For example, see alternative models presented in S1 and S2 Figs).

As a final caveat, we highlight our assumption that all edges (contacts) in our networks are equally likely to transmit disease. In reality, contacts can be highly heterogeneous, with household and health care contacts far more likely to transmit disease than casual social acquaintances. Our results should be robust when such heterogeneity is distributed randomly throughout the network. However, if individuals with more contacts tend to spend less time with each one, then epidemiological risk may be more homogeneous throughout the network and the advantage of the local decision strategies reduced. Although we do not model this scenario directly, we considered a homogeneous network where all individuals have the same number of contacts. This is roughly equivalent to mass action models that assume homogeneous contact rates and complete mixing [54]. The local strategies still prevail, but their relative efficiency is reduced, with far more vaccines required to achieve the same benefit (Fig 5). Conversely, in a network with greater heterogeneity (scale-free), the advantages of the local strategies are magnified.

This study prompts two practical questions. First, how do people actually make intervention decisions? Perhaps individuals fall nicely into one of these three decision-making camps. More likely, individual risk assessments are constrained by historical inertia [21, 24, 34, 46, 72], influenced by decisions of friends and family [1, 3, 10, 26, 28, 46], and integrate information from a combination of local and global data sources of variable reliability. Realistic decision models, driven by sociological survey data, can elucidate vaccine campaign failures and identify key pressure points for increasing uptake. Second, how can we streamline intervention campaigns to achieve efficient, rather than universal, adherence? This study reminds us that more intervention is not necessarily better intervention. The decision paradigm that most reduced transmission also required the least resources. Given the simplicity of our model, we do not suggest that public health agencies should promote ‘infection-counting’. Rather, we conclude that public health agencies should prioritize local disease surveillance and risk communications efforts and believe that data-driven models can be instrumental in designing effective outbreak information campaigns.

Supporting information

S1 Fig. Probability of vaccination as a function of the proportion of individuals infected.

Assuming R0 = 5, an individual with k = 65 in our exponential network, we plot the probability of vaccination versus the proportion infected, which indicates either the fraction of neighbors infected (local prevalence and local count) or the overall fraction of the population infected (global prevalence). The behavior of the strategies that use prevalence is similar (blue); however, during the curse of a simulation, they have different values.

https://doi.org/10.1371/journal.pone.0225576.s001

(TIF)

S2 Fig. Epidemiological outcomes for threshold-based vaccine decisions.

Individuals determine their infection risk according to the original local count, local prevalence or global prevalence equations. However, they vaccinate if and when their perceived risk crosses a specified vaccination threshold, rather than vaccinating probabilistically according to risk. This decision threshold (x-axes) impacts the proportion of the population that (A) vaccinates (and are not infected), (B) becomes infected (with or without vaccination), and (C) remains susceptible, under the three different decision models. Y-axes values are means across 500 stochastic simulations, assuming R0 = 5.

https://doi.org/10.1371/journal.pone.0225576.s002

(TIF)

S3 Fig. Sensitivity analysis with respect to (A) the delay between vaccination and protection and (B) both the network size and duration of the infectious period.

In both sets of graphs, the columns correspond to the three different strategies. Each point indicates a mean over 500 stochastic simulations that assume R0 = 5. (A) The x-axes indicate the delay between an individual deciding to vaccinate and becoming immune to infection. Each simulation assumes an exponential network with N = 10000 individuals and that 80% of vaccinated individuals become protected against infection (the remaining 20% remain fully susceptible). (B) The rows correspond to different network sizes (1000, 5000, 15000 and 20000).

https://doi.org/10.1371/journal.pone.0225576.s003

(TIF)

Acknowledgments

These project was supported by NIGMS MIDAS grant U01-GM087719. José Luis Herrera Diestra is supported by the São Paulo Research Foundation (FAPESP) under grants 2016/01343-7 and 2017/00344-2.

References

  1. 1. Serpell L, Green J. Parental decision-making in childhood vaccination. Vaccine 2006; 24: 4041–4046. pmid:16530892
  2. 2. Koh D, Takahashi K, Lim M-K, Imai T, Chia S-E, Qian F et al. SARS risk perception and preventive measures, Singapore and Japan. Emerg. Infect. Dis., 11 (4) (2005), pp. 641–642. pmid:15834989
  3. 3. Bhattacharyya S, Bauch CT (2015) Parental Decisions Unfold in Layers during a Vaccine Scare: Insights from Measles Vaccine Uptake Data. JSM Math Stat 2(1): 1007.
  4. 4. National Public Health Information Coalition: Promote the importance of immunizations with this communications toolkit. From: https://www.nphic.org/niam. Accessed 8 October 2018.
  5. 5. Larson HJ, de Figueiredo A, Xiahong Z, Schulz WS, Verger P, Johnston IG, Jones NS (2016). The State of Vaccine Confidence 2016: Global Insights Through a 67-Country Survey. EBioMedicine, 12, 295–301. http://doi.org/10.1016/j.ebiom.2016.08.042 pmid:27658738
  6. 6. Fifty-Sixth World Health Assembly. Agenda item 14.14. Available from: http://www.who.int/immunization/sage/1_WHA56_19_Prevention_and_control_of_influenza_pandemics.pdf. Accessed 8 September 2017.
  7. 7. Available from: https://www.cdc.gov/flu/fluvaxview/coverage-1516estimates.htm. Accessed 8 October 2018.
  8. 8. Hotez PJ. Texas and its measles epidemics. PLoS Med. 2016;13(10):e1002153. pmid:27780206
  9. 9. Omer SB, Enger KS, Moulton LH, Halsey NA, Stokley S, Salmon DA. Geographic clustering of nonmedical exemptions to school immunization requirements and associations with geographic clustering of pertussis. Am J Epidemiol. 2008;168(12):1389–1396. pmid:18922998
  10. 10. National Conference of State Legislatures. States with religious and philosophical exemptions from school immunization requirements. http://www.ncsl.org/research/health/school-immunization-exemption-state-laws.aspx. Published August 23, 2016. Accessed November 16, 2017.
  11. 11. Larson HJ, Jarrett C, Eckersberger E, Smith DM, Paterson P. Understanding vaccine hesitancy around vaccines and vaccination from a global perspective: a systematic review of published literature, 2007-2012. Vaccine. 2014;32(19):2150–2159. pmid:24598724
  12. 12. Report of the SAGE Working Group on Vaccine Hesitancy. World Health Organization. http://www.who.int/immunization/sage/meetings/2014/october/1_Report_WORKING_GROUP_vaccine_hesitancy_final.pdf. Published October 1, 2014. Accessed November 16, 2017.
  13. 13. Chapman GB, Coups EJ. Predictors of influenza vaccine acceptance among healthy adults. Prev. Med., 29 (1999), pp. 249–262. pmid:10547050
  14. 14. American Academy of Pediatrics. Vaccine safety: examine the evidence. https://www.aap.org/en-us/Documents/immunization_vaccine_studies.pdf. Updated April 2013. Accessed November, 2017.
  15. 15. Jain A, Marshall J, Buikema A, Bancroft T, Kelly JP, Newschaffer CJ. Autism occurrence by MMR vaccine status among US children with older siblings with and without autism. JAMA. 2015;313(15):1534–1540. pmid:25898051
  16. 16. Taylor LE, Swerdfeger AL, Eslick GD. Vaccines are not associated with autism: an evidence-based meta-analysis of case-control and cohort studies. Vaccine. 2014;32(29):3623–3629. pmid:24814559
  17. 17. Lau JT, Yang X, Pang E, Tsui HY, Wong E, Wing YK, SARS-related perceptions in Hong Kong. Emerg. Infect. Dis. 11 (3).
  18. 18. Philipson T, Private vaccination and public health: an empirical examination for U.S. measles. J. Hum. Resour. 31 (3), 1996.
  19. 19. Ahituv A, Hotz VJ, Philipson T, The responsiveness of the demand for condoms to the local prevalence of AIDS. J. Hum. Resour. 31 (4), 1996.
  20. 20. Funk S, Salathé M, Jansen VAA. Modeling the influence of human behaviour on the spread of infectious diseases: a review. J. R. Soc. Interface, 7 (50) (2010), pp. 1247–1256. pmid:20504800
  21. 21. Durham DP, Casman EA. Incorporating individual health-protective decisions into disease transmission models: a mathematical framework. J. R. Soc. Interface, 9 (68) (2012), pp. 562–570. pmid:21775324
  22. 22. De Zwart O, Veldhuijzen I, Elam G. Perceived threat, risk perception, and efficacy beliefs related to SARS and other (emerging) infectious diseases: results of an international survey. Int. J. Behav. Med., 16 (2009), pp. 30–40. pmid:19125335
  23. 23. Perra N, Balcan D, Gonçalves B, Vespignani A (2011) Towards a Characterization of Behavior-Disease Models. PLoS ONE 6(8): e23084. pmid:21826228
  24. 24. Thomas EA, Goldstone SE (2011) Should I or shouldn’t I: Decision making, knowledge and behavioral effects of quadrivalent HPV vaccination in men who have sex with men. Vaccine 29: 570–6. pmid:20950728
  25. 25. Hackett AJ (2008) Risk, its perception and the media: The MMR controversy. Community Pract 81: 22–5. pmid:18655642
  26. 26. Brown KF, Kroll JS, Hudson MJ, Ramsay M, Green J, Long SJ et al. (2010) Factors underlying parental decisions about combination childhood vaccinations including MMR: A systematic review. Vaccine 28: 4235–48. pmid:20438879
  27. 27. Samba E, Nkrumah F, Leke R (2004) Getting polio eradication back on track in Nigeria. N Engl J Med 350: 645–6. pmid:14960740
  28. 28. Myers LB, Goodwin R (2011) Determinants of adults’ intention to vaccinate against pandemic swine flu. BMC Public Health 11: 15. pmid:21211000
  29. 29. Fabry P, Gagneur A, Pasquier JC (2011) Determinants of A (H1N1) vaccination: Cross-sectional study in a population of pregnant women in Quebec. Vaccine 29: 1824–9. pmid:21219988
  30. 30. Liu R, Wu J, Zhu H (2007) Media/psychological impact on multiple outbreaks of emerging infectious diseases. Comput Math Method M 8: 153–64.
  31. 31. Cui J, Sun Y, Zhu H (2008) The impact of media on the control of infectious diseases. J Dyn Differ Equ 20: 31–53.J.
  32. 32. Li Y, Cui J (2009) The effect of constant and pulse vaccination on SIS epidemic models incorporating media coverage. Commun Nonlinear Sci 14: 2353–65.
  33. 33. Tchuenche J, Dube N, Bhunu C (2011) The impact of media coverage on the transmission dynamics of human influenza. BMC Public Health 11: S5. pmid:21356134
  34. 34. Cornforth DM, Reluga TC, Shim E, Bauch CT, Galvani AP, Meyers LA (2011) Erratic Flu Vaccination Emerges from Short-Sighted Behavior in Contact Networks. PLoS Comput Biol 7(1): e1001062. https://doi.org/10.1371/journal.pcbi.1001062 pmid:21298083
  35. 35. Fu F, Christakis NA and Fowler JH. (2017) Dueling biological and social contagions. Scientific Reports volume 7, Article number: 43634
  36. 36. Andrews MA, Bauch CT. The impacts of simultaneous disease intervention decisions on epidemic outcomes. J Theor Biol. 2016 Apr 21;395:1–10. pmid:26829313
  37. 37. Fine P, Eames K, Heymann DL. “Herd Immunity”: A Rough Guide, Clinical Infectious Diseases, Volume 52, Issue 7, 1 April 2011, Pages 911–916, https://doi.org/10.1093/cid/cir007
  38. 38. John TJ and Samuel R. Herd Immunity and Herd Effect: New Insights and Definitions. European Journal of Epidemiology, Vol. 16, No. 7 (2000), pp. 601–606.
  39. 39. Yamin D, Jones FK, DeVincenzo JP, Gertler S, Kobiler O, Townsend JP, and Galvani AP (2016). Vaccination strategies against respiratory syncytial virus. Proceedings of the National Academy of Sciences of the United States of America, 113(46), 13239–13244. pmid:27799521
  40. 40. Chao DL, and Dimitrov DT (2016). Seasonality and the effectiveness of mass vaccination. Mathematical biosciences and engineering: MBE, 13(2), 249–59. pmid:27105983
  41. 41. Ferrari MJ, Bansal S, Meyers LA, Bjornstad ON; Network frailty and the geometry of herd immunity; Proc R Soc B Biol Sci 2006, 273: 2743–2748.
  42. 42. Pastor-Satorras R, Vespignani A. (2002) Immunization of complex networks. Phys. Rev. E 65, 036104.
  43. 43. Pastor-Satorras R, Castellano C, Van Mieghem P, and Vespignani A. (2015) Epidemic processes in complex networks. Rev. Mod. Phys. 87, 925.
  44. 44. Goltsev AV, Dorogovtsev SN, Oliveira JG, and Mendes JFF. (2012) Localization and Spreading of Diseases in Complex Networks. Phys. Rev. Lett. 109, 128702. pmid:23006000
  45. 45. Meyers LA. (2007) Contact network epidemiology: Bond percolation applied to infectious disease prediction and control. Bulletin of the American Mathematical Society 44: 63–86.
  46. 46. Ndeffo Mbah ML, Liu J, Bauch CT, Tekel YI, Medlock J, Meyers LA, et al. (2012) The Impact of Imitation on Vaccination Behavior in Social Contact Networks. PLoS Comput Biol 8(4): e1002469. https://doi.org/10.1371/journal.pcbi.1002469 pmid:22511859
  47. 47. Herrera JL, Srinivasan R, Brownstein JS, Galvani AP, Meyers LA (2016) Disease Surveillance on Complex Social Networks. PLoS Comput Biol 12(7): e1004928. https://doi.org/10.1371/journal.pcbi.1004928 pmid:27415615
  48. 48. Bagnoli F Lió P, Sguanci L (2007) Risk perception in epidemic modeling. Phys. Rev. E 76, 061904.
  49. 49. Massaro E, Bagnoli F (2014) Epidemic spreading and risk perception in multiplex networks: A self-organized percolation method. Phys. Rev. E 90, 052817.
  50. 50. Rizzo A, Frasca M, and Porfiri M (2014) Effect of individual behavior on epidemic spreading in activity-driven networks. Phys. Rev. E 90, 042801.
  51. 51. Moinet A, Pastor-Satorras R, and Barrat A (2018) Effect of risk perception on epidemic spreading in temporal networks. Phys. Rev. E 97, 012313. pmid:29448478
  52. 52. Meyers LA, Pourbohloul B, Newman MEJ., Skowronski DM, Brunham RC. (2005). Network theory and SARS: Predicting outbreak diversity. Journal of Theoretical Biology 232: 71–81. pmid:15498594
  53. 53. Abbey H; An examination of the Reed-Frost theory of epidemics. Hum Biol 1952, 24:201–233. pmid:12990130
  54. 54. Bansal S, Grenfell BT, Meyers LA. (2007). When individual behavior matters: homogeneous and network models in epidemiology. Journal of the Royal Society Interface 4: 879–891.
  55. 55. Newman MEJ, Strogatz SH, and Watts DJ (2001) Random graphs with arbitrary degree distributions and their applications. Phys. Rev. E 64.
  56. 56. Barabási AL, Albert R (1999). Emergence of scaling in random networks. Science. 286 (5439): 509–512. pmid:10521342
  57. 57. Huai Y, Xiang N, Zhou L, Feng L, Peng Z, Chapman RS, et al. Incubation period for human cases of avian influenza A (H5N1) infection, China, Emerging Infectious Diseases • www.cdc.gov/eid • Vol. 14, No. 11, November 2008.
  58. 58. De Serres G, Rouleau I, Hamelin M-E, Quach C, Skowronski D, Flamand L, et al. Contagious period for pandemic (H1N1) 2009, Emerging Infectious Diseases • www.cdc.gov/eid • Vol. 16, No. 5, May 2010.
  59. 59. Chowell G, Hengartner NW, Castillo-Chavez C, Fenimore PW, Hyman JM. The basic reproductive number of Ebola and the effects of public health measures: the cases of Congo and Uganda. J Theor Biol. 2004 Jul 7;229(1):119–26. pmid:15178190
  60. 60. Althaus CL. Estimating the Reproduction Number of Ebola Virus (EBOV) During the 2014 Outbreak in West Africa. PLoS Currents. 2014;6. pmid:25642364
  61. 61. Kretzschmar M, Teunis PFM, Pebody RG (2010) Incidence and Reproduction Numbers of Pertussis: Estimates from Serological and Social Contact Data in Five European Countries. PLoS Med 7(6): e1000291. pmid:20585374
  62. 62. Wallinga J, Teunis P. Different Epidemic Curves for Severe Acute Respiratory Syndrome Reveal Similar Impacts of Control Measures, American Journal of Epidemiology, Volume 160, Issue 6, 15 September 2004, Pages 509–516, https://doi.org/10.1093/aje/kwh255
  63. 63. Keeling MJ, Rohani P. Modeling Infectious Diseases in Humans and Animals (2008) Princeton University Press: 21.
  64. 64. Du Z, Zhang W, Zhang D, Yu S, and Hao Y. (2017). Estimating the basic reproduction rate of HFMD using the time series SIR model in Guangdong, China. PLoS ONE, 12(7), e0179623. http://doi.org/10.1371/journal.pone.0179623 pmid:28692654
  65. 65. Bacaer N, Abdurahman X, Ye J, Auger P. On the basic reproduction number R0 in sexual activity models for HIV/AIDS epidemics: example from Yunnan, China. Math Biosci Eng. 2007 Oct;4(4):595–607. pmid:17924713
  66. 66. Nsubuga RN, White RG, Mayanja BN, Shafer LA (2014) Estimation of the HIV Basic Reproduction Number in Rural South West Uganda: 1991–2008. PLoS ONE 9(1): e83778. https://doi.org/10.1371/journal.pone.0083778 pmid:24404138
  67. 67. Halloran ME. (2005). Secondary Attack Rate. In Encyclopedia of Biostatistics (eds P. Armitage and T. Colton).
  68. 68. Newman MEJ (2001), Ego-centered networks and the ripple effect. Social Networks 25, 83–95.
  69. 69. Cohen R, Havlin S, and ben-Avraham D. (2003) Efficient Immunization Strategies for Computer Networks and Populations. Phys. Rev. Lett. 91, 247901. pmid:14683159
  70. 70. Kitsak M, Gallos LK, Havlin S, Lijeros F, Muchnik L, Stanley HE, Makse HA. (2010) Identification of influential spreaders in complex networks, Nature 6, 888–893.
  71. 71. Christley RM, Pinchbeck GL, Bowers RG, Clancy D, French NP, Bennett R and Turner J. Infection in Social Networks: Using Network Analysis to Identify High-Risk Individuals, Am J Epidemiol 2005;162:1024–1031. pmid:16177140
  72. 72. Taylor E, Atkins KE, Medlock J, Li M, Chapman GB, Galvani AP. Cross-Cultural Household Influence on Vaccination Decisions. Medical Decision Making Vol 36, Issue 7, pp. 844—853. First published date: June-17-2015.