key: cord-0496089-onxh6dri authors: Epstein, Joshua M.; Hatna, Erez; Crodelle, Jennifer title: Coupled Contagion: A Two-Fears Epidemic Model date: 2021-02-22 journal: nan DOI: nan sha: 8398f2bc9354f5606034cc4aa74db1b7c91b520b doc_id: 496089 cord_uid: onxh6dri We present a differential equations model in which contagious disease transmission is affected by contagious fear of the disease and contagious fear of the control, in this case vaccine. The three contagions are coupled. The two fears evolve and interact in ways that shape distancing behavior, vaccine uptake, and their relaxation. These behavioral dynamics in turn can amplify or suppress disease transmission, which feeds back to affect behavior. The model reveals several coupled contagion mechanisms for multiple epidemic waves. Methodologically, the paper advances infectious disease modeling by including human behavioral adaptation, drawing on the neuroscience of fear learning, extinction, and transmission. In classical mathematical epidemiology -the venerable tradition of the 1927 Kermack-McKendrick modelindividuals do not adapt their contact behavior during epidemics [1, 2] . Specifically, they do not endogenously engage in social distancing based on fear. Yet, such behavior is well-documented in true epidemics. In 2008, Epstein et al. published Coupled Contagion Dynamics of Fear and Disease [3] , a model that introduced the idea of two interacting contagions: one of the disease proper and one of fear of the disease. The model's core narrative is that epidemic growth induces fear. Contagious fear among healthy susceptible people, in turn, induces self-isolation. By depriving the epidemic of fuel, in the form of susceptibles, this self-isolation suppresses the disease. When disease prevalence becomes low, however, so does the fear. Thus, susceptible people (no longer fearful) come out of hiding. But, because there are still infectious individuals in circulation, this pours gasoline (susceptible individuals) on the remaining embers (the infectives), igniting a second wave. Precisely this occurred historically, in the 1918 Pandemic Flu (see [4] ), and history repeated itself in the multi-wave COVID-19 pandemic [5] . Recent work on the neuroscience of fear lends scientific support to the postulate of fear contagion, and a recent agent-based model explicitly includes fear modules grounded in that neuroscience; see [6, 7] . The literature on behavioral adaptation in epidemics has grown in several directions; see [8] [9] [10] [11] . In the present work, we modify and extend the original coupled contagion model [3] in light of recent advances, subsuming it in a more general framework that-while including contagious fear of disease-adds contagious fear of vaccine. The World Health Organization recently included vaccine refusal in the top ten threats to global health [12] . It is responsible for the resurgence of several deadly vaccine-preventable diseases, including measles and pertussis in the US and even polio in several countries [13, 14] . During the swine flu pandemic of 2009, roughly 40 percent of Americans refused the vaccine [15] . And, writing as COVID-19 vaccination is beginning, there is concern that refusal will undermine the attainment and maintenance of herd immunity to the SARS-CoV-2 virus and its variants. In our model, as discussed in [16] , "Everything turns on the relationship between the two fears, one of disease, the other of vaccine." If fear of the disease exceeds fear of the vaccine in the population, the rate of vaccine acceptance rises, and the disease may be suppressed. However, if the prevalence of the disease is suppressed enough, fear of the disease may fall below fear of the vaccine (as might happen when a disease recedes from our collective memory). Now the vaccine is scarier than the disease, people eschew the vaccine, and a new disease cycle can explode. This narrative also rings true historically. Smallpox, one of the great scourges of human history, kills roughly 30 percent of those infected [17] . Yet, even when inoculation (with cowpox) was discovered, cycles of vigilance and complacency kept smallpox alive. In her social history of smallpox, the Speckled Monster, Jennifer Lee Carrell [18] recounts, "In London, inoculation's popularity waxed and waned through the 1730s, with the force of the disease: in bad years, people flocked to be inoculated; in lighter years, the practice shrank. Inoculation was a security-the only security-to cling to within the terror of an epidemic; in times of good health, however, it looked like a foolish flirtation with danger." Our two-fear model generates such cycles and related dynamics. As discussed in [3] , our approach differs from the so-called "rational epidemics" tradition (stemming from [19, 20] ) which does not model contagious fears. Rather, as in microeconomics and game theory, agents maximize utility conditional on the disease's prevalence. While illuminating in several important settings, prevalence-elastic optimal adaptation in the rational choice tradition is not well-suited to capture prevalence independent fear contagions such as Morgellan's disease, an internet-disseminated delusional parasitosis [21] , or the mass panics that occurred in Surat India (1994), or during Ebola [22] , or in recent episodes of vaccine refusal [23] . Indeed, cognitive neuroscience demonstrates that the fear response and fear learning generally are not fundamentally choice-like, or even necessarily conscious, none of which means they can't be modeled or counteracted [6, 7, [24] [25] [26] [27] [28] ]. After the model's technical exposition (Section 2), we offer four base scenarios with discussions of their dynamics and importance (Section 3). Analytical results and extensive sensitivity analyses are provided and discussed in Section 4. We begin with the pure compartmental Susceptible-Infected-Recovered (SIR) version of a contagious disease alone. Every subsequent scenario subsumes the preceding one, as follows: Scenario 1: contagious disease Scenario 2: contagious disease + fear of the disease Scenario 3: contagious disease + fear of the disease + vaccination Scenario 4: contagious disease + fear of the disease + vaccinations + fear of the vaccinations These four scenarios are of central concern to public health. All numerical assumptions (parameter settings and initial conditions) are given in the Supplementary Information (SI), ensuring replicability. Several mathematical conditions for growth are derived there as well. The SI also includes a pure fear "Salem Witches" scenario, where fear propagates in the absence of any disease, further distinguishing the approach from prevalence elastic rational adaptation. On emotional contagion and its mechanisms, see [6, 29] . Proceeding with the exposition, we first define all state variables and parameters of the model in Tables 1 and 2 respectively. It may be of interest that we use an average infectious period of seven days ( 1 γ = 7) and a basic reproduction number (R 0 ) of two ( β γ = 2) for the scenarios discussed. The proportion of susceptible individuals with no fear The proportion of susceptible individuals who fear the disease The proportion of susceptible individuals who fear the vaccine The proportion of (pathogen) infectious individuals The proportion of recovered individuals (persons who had the disease and gained immunity) The proportion vaccinated individuals The proportion of recently-vaccinated individuals who fear the vaccine because of an adverse reaction v(t) The rate of vaccination (1/days) The mathematical model relating these variables and parameters consists of the eight coupled nonlinear ordinary differential equations shown below. For expository efficiency, we use a well-mixed model. Natural extensions would include social networks [30] [31] [32] [33] and agent-based formulations [34] . Capital letters indicate infection states, while subscripts indicate fear states. For example, the S f d compartment is the fraction of the population that is susceptible to the disease and fears the disease (subscript f d). The S f v compartment is the fraction of the population that is susceptible to the disease and fears the vaccine (subscript f v). While this dynamical system is rich, as discussed below, Eqs (1), (4), and (5) reduce to the familiar SIR model when there are no fears and the terms subscripted by f d or f v are set to zero. Moreover, each of the fear contagions also propagates in classical fashion, as can be seen in Eqs (2) and (3). A flow diagram of the model is shown in Figure 1 . The equations include six population compartments, each representing the proportion of individuals in the given state at any time. The sum of these six compartments is always 1. A susceptible individual may acquire the disease by being exposed to an infectious person. The effective contact parameter represents the rate of transmission. When effective contact occurs, a susceptible person becomes infectious for an average period of 1 γ days. Once the infectious period ends, the individual recovers and gains permanent immunity to the pathogen (compartment R nat ). We simplified the model by not including a pre-infectious (latent) period, although this is an obvious extension. We consider three types of susceptible individuals: persons without fear (S), persons who fear the disease (S f d ), and persons who fear the vaccine (S f v ). A susceptible person can retain only one fear at a given time, and all fearful persons share the same fear intensity. A non-fearful (susceptible) individual may acquire fear of the disease by interacting with infectious or disease-fearful persons. These processes represent scenarios in which a susceptible person observes or communicates with an infectious (sick) individual or with a disease fearful person. Unlike the transmission of the pathogen, such interactions could occur at a distance (as on social media) and thus require a dedicated effective contact rate parameter (β f d ). Note that an infectious individual can infect a susceptible person with either the pathogen or fear of the disease, but again, not both. Fear of the disease affects the behavior of susceptible individuals. These persons may take protective actions, such as self-isolation, mask-wearing, social distancing, avoidance of travel and mass gatherings, and improved personal hygiene. In the interest of simplicity, such actions are modeled using the relative risk parameter p, which is used to scale down β. This is a fundamentally different representation than in [3] , where distanced individuals were a separate compartment. Here, they are not. In the present model, a value of p = 0.25 represents a 75% decrease in the likelihood of a disease-fearful individual becoming infected with the disease compared to a susceptible individual with no fear. Disease-fearful individuals may also choose to gain permanent immunity through vaccination. We assume that a small proportion (σ) of vaccinated individuals experience adverse effects or associate an unrelated discomfort with the vaccine. These individuals (A) acquire a (transmissible) fear of the vaccine while gaining full immunity. The rest of the vaccinated individuals, a proportion of 1 − σ, gain immunity without acquiring the fear (R nat ). Susceptible individuals acquire fear of the vaccine by interacting with vaccinated persons who had an adverse experience (A) or with vaccine fearful susceptible persons (S f v ). The effective contact rate of such interactions is β f v . We know from neuroscience that, post-traumatic stress notwithstanding, fear is not permanent but decays in the absence of an aversive stimulus. In this model, susceptible people may naturally overcome both fears (of disease and vaccine) and join the compartment of non-fearful susceptible individuals (S). Our model contains two paths for such fear decay, or "extinction" as it is called in behavioral neuroscience [35] . Specfically, we think of exposures (direct or indirect) to disease-infected people as classical associative fear-conditioning trials. A classic example of a fear conditioning trial is as follows. If a person is simply shown a benign blue light, no manifestations of fear (e.g., freezing, pupil dilation, adrenaline spikes, increased heart rate, electrodermal activity) or neural correlates of fear, such as activation of (e.g., oxygenation and recruitment of blood to) the amygdala, as seen in fMRI [36] are observed. By contrast, if the subject is unexpectedly given an aversive electric shock, the amygdala is immediately stimulated, triggering a suite of fear responses. Importantly, if the two stimuli are repeatedly paired -blue light followed shortly by shock -the subject will come to associate (not necessarily consciously) the light with the shock, to the point where the blue light alone elicits the same amygdala response as the shock. By a process of associative learning, the subject has been "conditioned" to fear the blue light. If these light-shock pairings are discontinued, the fear of the blue light will decay. Both the fear acquisition phase and fear extinction phases can be modeled mathematically [37] . Consider a person whose fear of the disease has prompted self-isolation. This person's fear may decay in two ways. The first is by eliminating direct (aversive) exposures to disease-infected individuals; conditioning trials are thereby suspended, and "fear extinction" commences. In the model, this natural decay is exponential, consistent with the simple seminal Rescorla-Wagner model [38] . We assume that, in the absence of a fear stimulus, a person will retain fear for an average duration of 1 γ f days. On the widespread use of the Rescorla-Wagner model, see [39] . For other learning models, see for example [40] and [41] . The second path to overcoming fear is social and distinct from extinction through stimulus deprivation. Individuals may lose fear by communicating with persons who have recovered from the fearful event. These reassuring exposures (think of repeated blue light and candy pairings) can damp the conditioned fear. This would be called counter-conditioning, over-writing a negative response with a positive one. (On the relative effectiveness of extinction and counter-conditioning in diminishing children's fear, see [42] ). By interacting with a recovered person (R nat ), a disease fearful person (S f d ) may lose their fear. Similarly, a vaccine-fearful person (S f v ) may lose their fear by interacting with a protected vaccinated person (R vac ). In our model, vaccinated persons who had gained fear due to a negative vaccine experience (A) abandon the fear only via the first path: natural exponential decay. We assume that their first-hand experience with the vaccine makes them resistant to social influence. Analogous to disease fear, they retain their fear of vaccine for an average duration of 1 γ f days and then join the compartment of vaccinated individuals (R vac ). Widespread distancing and vaccination also cut the disease's growth rate and can even make it negative-the herd immunity condition-which amplifies their suppressive effects. The daily rate at which fearful disease-susceptible persons vaccinate, v(t), may change over time due to a mechanism of social influence; see Eq (8) . Specifically, we assume that the growth rate of v(t) increases dv dt > 0 when fear of the disease exceeds the fear of the vaccine. It decreases dv dt < 0 when the reverse obtains -when vaccine-fear is more prevalent than disease-fear. We represent this effect using the difference between the two prevalences (S f d − S f v ): Clearly, S f d − S f v = 0, is a tipping point of the dynamics. Several mechanisms can affect the fear ordering. If the model begins with disease fear exceeding vaccine fear (S f d − S f v > 0), vaccination expands. However, this itself can endogenously suppress the disease to the point where fear of disease falls below fear of vaccine. At this point, the fear ordering switches, reversing the sign of dv dt , opening the door for disease resurgence through vaccine refusal. Of course, two other mechanisms can drive fear of vaccine to exceed fear of disease. One is an accumulation of adverse vaccine events represented by the A compartment. Another mechanism (not included here) would be exogenous suppression of disease fear (S f d ) through statements by officials underestimating the threat. We turn now to the core scenarios of the model. Again, all numerical assumptions are provided in the text or the SI. Here, we "dock" the model to the classic case, an SIR epidemic with no fears, with a disease transmission rate β, and a single recovery (and subsequently immune) rate, γ. In this case, Eqs (1), (4), and (5) reduce to the Kermack-McKendrick model. A reference plot of the main dynamics is given in Figure 2 , which illustrates our graphical strategy. To reduce clutter, it will prove useful to have four plots focused on different aspects of the coupled contagions: susceptibles, vaccine uptake, infection, and removals, as shown in Figure 2 . With all fears and all vaccinations clamped at zero (see purple and green curves in Figure 2a and black curve in Figure 2b) , we see the classical one peak curve of infectives in Figure 2c , the falling susceptible curve in 2a, and the familiar recovered curve in Figure 2d . Now we add contagious fear of the pathogen, so there are two contagions, as depicted in Figure 3 . The core narrative here is that the initial spike of infections (the blue curve) stimulates a fear spike (the purple curve). People reduce their contacts out of fear (this is modeled through p), which suppresses disease spread. As the disease wanes, however, so does the fear of it. Now, susceptibles go back into circulation, which pours fuel on the infective embers, and a second wave ensues. The second wave is larger than the first. Why? Because in our model, there are two mechanisms of fear decay, and they amplify one another. One mechanism is the "natural decay" governed by the parameter γ f . The second is the "contagious," fear-reversal mechanism. People who have recovered from the disease are in contact with those who are still fearful. The recovereds' low fear is also transmitted, emboldening the fearful people in hiding to come 'out of the basement' when it is still unsafe. This 'complacency contagion,' if you will, amplifies the natural fear decay rate to produce a very sharp fear reduction. This pours a larger number of susceptibles onto the circulating infectives than would either fear decay mechanism alone. The result is that the second wave of the disease can be larger than the first, as occurred in 1918 [43] . All of this is chronicled in Figure 3 . Its robustness is explored in Section 4. As the data science of social media shows [44] , fear can spread much faster and much farther than the disease itself (a good thing when it induces preventive measures). For the earlier 2008 model, an analytic expression for the R 0 of fear, and conditions for fear of disease to spread faster than the disease itself, are given in [3] . The mathematics are different here and several analytical growth conditions for the present model are given in the SI text. An obvious reason for fear to outpace disease is that disease transmission requires direct physical contact while fear transmission does not. Indeed, there are two channels to acquire disease fear in our model-through contact with an infectious person (in the I compartment) or contact with a frightened susceptible person, in the S f d compartment. Scared individuals-whether sick or not-remove themselves from circulation, social distancing with an effectiveness governed by the parameter p. This endogenously affects the contact dynamic, and thus the disease epidemic itself. Sometimes, the self-isolation is sufficient to produce herd immunity and epidemic fade-out. (See Section 4). In other cases, because disease prevalence is low, individuals recover from fear at a rate α f despite the presence of disease. This releases fresh susceptibles onto the still-circulating infectives, generating a second wave, as shown in Figure 3 . We now extend the model further, adding vaccination, but not yet the fear of it. Vaccination can mitigate the second wave generated in Scenario 2, a beneficial result from a public health perspective. In Figure 4a , fear of disease (purple) exceeds fear of vaccine (green), which is clamped at zero. Vaccine uptake thus increases, as shown in Figure 4b . Now we see both "natural" and vaccine-induced removals (Figure 4d) . The combined effect is to suppress the second wave, as evident from the lower-left infection curve. Notice that, while the second wave of Figure 4 is clearly suppressed, there is still a small second wave. The mechanism for the two peaks here lies in the assumed effectiveness of social distancing. We have assumed that the probability p of contracting the disease while fearful is low. If we increase p, the peaks get closer together until (in the limit) they converge to a single peak again. The sensitivity of the phenomenon to variations in p is given in Section 4 below. In Scenario 4, the fear of vaccination "wins," and the outbreak is again unmitigated. People do vaccinate at the beginning of the outbreak but stop too soon because the fear ordering reverses. In Figure 5a , we see that the fear ordering changes at roughly 100 days, at which point fear of vaccine (green) rises above the fear of disease (purple). This reverses the sign of the rate of change of v(t) [(Eq 8)], and a second wave of infections ensues. The base scenarios exhibit several mechanisms for the emergence, timing, size, and decay of multiple waves. We now explore their sensitivity to various parameters. To begin, we return to the case of the disease and fear of the disease only (Scenario 2, Section 3.1.2) and study the effect of changing the relative risk p of acquiring the disease for pathogen-fearful individuals. Figure 6 shows that if p is decreased to 0, meaning that those who are fearful of the disease go into hiding and have a 0% chance of contracting the disease, the epidemic will be prevented (we see that only about 10% of the population gets the disease). By contrast, as p increases, fearful individuals become more risk-neutral (increasing their likelihood of contracting the pathogen), and the epidemic worsens; at its worst, with the fearful individuals not altering their behavior at all (p = 1), we see that about 80% of the population becomes infected. We see also that there are single-wave and two-wave regimes, depending on p. If p is small (less than about 0.4) but positive, we see a second wave emerge as fearful individuals hide away and then return to circulation. If p is larger than about 0.4, then the fearful individuals don't lower their risk enough to preserve a susceptible population sufficient to produce a second wave. As noted earlier, an important extension (among several) of the original coupled contagion model [3] is our inclusion of a second mechanism of fear loss. In addition to spontaneous loss of fear, disease fearful persons (S f d ) may lose their fear by interacting with recovered persons (R nat ). The effective contact rate for this interaction is α f . We can see that the infectives curve of different α f values overlap until the first peak is reached; see Figure 7a . The curves do not differ because the number of recovered in the beginning is too low to reduce fear significantly. The differences become apparent once the infectives curve drops. Higher values of α f cause people to lose their fear and abandon their protective measures. This process increases the number of persons that are infected and in turn, the number of recovered. The larger number of recovered causes a larger fraction of persons to lose their fear of the disease and so on. The result is a second wave when the contact rate is sufficiently high. As α f increases, the second peak is higher and occurs sooner. This process increases the fraction of infected persons; see Figure 7b . Table S1 , column 2. We now focus on the two-fears scenario (Scenario 4) and explore how the contact rate of the fear of the disease (β f d ) and the contact rate of the fear of the vaccine (β f v ) influence the model behavior. The two effective contact rates determine how fast the fears are transmitted in the population. The fraction of vaccinated and infected persons as a function of (β f d , β f v ) is shown in Figure 8 . When fear of the vaccine is transmitted sufficiently faster than the fear of the disease (dark blue region in Figure 8a ), the population eschews vaccine, and a large portion of the population becomes infected with the disease (over 75%, as shown in bright orange region of Figure 8b ). When the fear the disease is transmitted fast with a sufficiently low transmission rate of the fear of the vaccine, a small proportion of the population is infected (the blue region in 8b). The fear of the disease spreads too slowly, and the number of persons who adopt protective behavior is too low to create two peaks. We now explore how the fears' contact rates (β f d , β f v ) affect the number of peaks that occur in the infectives curve (I). We define a peak as a local maximum with a proportion of infectives above 0.01. Figure 9a shows that the model produces zero, one, or two peaks. Point A is located in a region of fast transmission of the vaccine fear, which leads to few vaccinated individuals. Only one peak is generated (see Figure 9b ) because the disease's level of fear is too low for temporarily decreasing the disease spread, and a majority of the population is infected. Two peaks are generated in the purple area of Figure 9a . The two peaks region in (β f d , β f v ) space can be well-represented by the equation β f v = 4.6(β f d − 0.45) + 1.3 (see black arrow). Notice that this area includes the interface between the red and blue regions of Figure 8 which indicate the transition between the high and low fractions of vaccinated and infected persons. As we move along this arrow toward larger values of β f d and β f v , we see that the second peak in the infection increases while the first peak decreases (Figure 9d ). As more fear enters the population, we see individuals fearful of the disease hiding out (decreasing the first peak of infection), driving the vaccination rate up, and then a disproportionate fear of the vaccine (increasing the second peak). Along this arrow, the total fraction of the population that becomes infected remains relatively constant (Figure 9e ) as the arrow is parallel to the interface of the two regions. If β f d is high enough (above ∼0.45 per day), there will be no epidemic unless β f v is sufficiently high. Even then, there will be only one peak because there are too few susceptibles to cause an initial peak high enough to reverse the order of fears (that is, too many are fearful of the disease, reducing their transmission rate); see Figure 9c . A proportion σ of individuals experience adverse reactions to the vaccine and develop a temporary fear of the vaccine, which may be transmitted to susceptible individuals. The more people that have such reactions, the easier it is for the fear to spread. When we increase the proportion of adverse reactions, the second peak in the proportion of infectives rises while the first peak remains unaffected (Figure 10a ). The first peak remains the same because too few persons fear the vaccine (Figure 10b) , and they have a negligible effect on the vaccination rate (Figure 10c ). After the first peak, the fear of the vaccine spreads, the vaccination rate drops, and more people become infected, leading to the second peak. Further sensitivity analyses can, of course, be conducted. But these demonstrate how the scenario dynamics respond to variations in several key parameters. The SI gives parameter values and initial conditions for the scenarios and derives analytical expressions for the growth (the R n values) of the disease and fear epidemics. The last of these specialize to give the condition for a fear epidemic in the absence of disease, which is also shown. We have extended earlier work on coupled contagion dynamics [3] . In addition to a contagious disease and contagious fear of it, we have added a second contagion: fear of the control, in this case, vaccine. In addition, unlike [3] , we include both the classical extinction of fear and its contagious evaporation. The interaction of these entangled contagions reveals several novel behavioral mechanisms for multiple waves of infection and for their timing, size, and form. To adequately capture such mechanisms, infectious disease modeling must begin to incorporate behavioral neuroscience. Human behavior is complex and involves interacting affective, deliberative, and social components. To be sure, some health decisions qualify as canonically rational. But often, as Hume noted, "Reason is a slave to the passions," the more so in settings of extreme stress like pandemics. Simple models grounded in the neuroscience of fear and its transmission can deepen epidemic modeling. The parameters employed in each scenario are shown in Table S1 below, illustrating again the cumulative nature of the exercise. Table S2 : Initial values used in the scenarios in Section 3.1. Growth conditions for the no-fear model As shown in the text, in the absence of fear, our model reduces to the classical SIR model. The resulting R n values are no exception. Here, we calculate the condition for continued growth of an epidemic. If we consider the case of no fear by setting S f d = S f v = 0, we arrive at the typical growth condition: R n = βS γ > 1. If we plot this value for the simulation shown in Section 3.1.1, Figure 2 , we see that when the R n value is above 1, the disease is increasing; see Figure S1 In a similar way, we can calculate an R n value for the spread of the fear of the disease by finding a condition on which the derivative of S f d is positive Finally, the extreme case is fear contagion in the absence of actual disease. Baseless fear contagions are common outside public health. The Salem Witch Trials of 1692 come immediately to mind. But in public health there are purely psychogenic contagions like Morgellan's disease, noted earlier. If we assume that there is no disease (I = 0), the growth condition above simplifies to Finally, if we begin with one baselessly fearful individual, we see that for values of disease-free R n greater than 1, we have an epidemic of fear without disease; see Figure S2 . A contribution to the mathematical theory of epidemics The kermack-mckendrick epidemic model revisited Coupled contagion dynamics of fear and disease: mathematical and computational explorations The effect of public health measures on the 1918 influenza pandemic in u.s. cities CDC covid data tracker. trends in number of covid-19 cases and deaths in the us reported to cdc, by state/territory Agent Zero: Toward neurocognitive foundations for generative social science Advancing agent zero Adaptive human behavior in epidemiological models Modelling the influence of human behaviour on the spread of infectious diseases: a review Towards a characterization of behavior-disease models Infection prevention behaviour and infectious disease modelling: a review of the literature and recommendations for the future WHO releases list of 10 threats to global health Increase in measles cases -united states Vaccine hesitancy, vaccine refusal and the anti-vaccine movement: influence, impact and implications Compliance with anti-h1n1 vaccine among healthcare workers and general population Are we already missing the next epidemic? Smallpox and its eradication (world health organization The Speckled Monster: A Historical Tale of Battling the Smallpox Epidemic Rational epidemics and their public control Integrating behavioral choice into epidemiological models of aids Delusional parasitosis facilitated by web-based dissemination Mass media and the contagion of fear: the case of ebola in america Weaponized health communication: Twitter bots and russian trolls amplify the vaccine debate Emotion circuits in the brain The emotional brain, fear, and the amygdala Synaptic self: How our brains become who we are. Penguin Emotion circuits in the brain Rethinking the emotional brain Emotional contagion Coupled disease-behavior dynamics on complex networks: A review Epidemics on networks Modelling disease outbreaks in realistic urban social networks Mathematics of epidemics on networks A distributed platform for global-scale agent-based models of disease transmission Timing of extinction relative to acquisition: a parametric analysis of fear extinction in humans Neural signatures of human fear conditioning: an updated and extended meta-analysis of fmri studies Context, learning, and extinction A theory of pavlovian conditioning: Variations in the effectiveness of reinforcement and nonreinforcement The widespread influence of the rescorla-wagner model A model for pavlovian learning: variations in the effectiveness of conditioned but not of unconditioned stimuli The relative effectiveness of extinction and counter-conditioning in diminishing children's fear Mechanistic modelling of the three waves of the 1918 influenza pandemic The pandemic of social media panic travels faster than the covid-19 outbreak