The complex interplay between risk tolerance and the spread of infectious diseases ================================================================================== * Maximilian Nguyen * Ari Freedman * Matthew Cheung * Chadi Saad-Roy * Baltazar Espinoza * Bryan Grenfell * Simon Levin ## Abstract Risk-driven behavior provides a feedback mechanism through which individuals both shape and are collectively affected by an epidemic. We introduce a general and flexible compartmental model to study the effect of heterogeneity in the population with regards to risk tolerance. The interplay between behavior and epidemiology leads to a rich set of possible epidemic dynamics. Depending on the behavioral composition of the population, we find that increasing heterogeneity in risk tolerance can either increase or decrease the epidemic size. We find that multiple waves of infection can arise due to the interplay between transmission and behavior, even without the replenishment of susceptibles. We find that increasing protective mechanisms such as the effectiveness of interventions, the number of risk-averse people in the population, and the duration of intervention usage reduces the epidemic overshoot. When the protection is pushed past a critical threshold, the epidemic dynamics enter an underdamped regime where the epidemic size exactly equals the herd immunity threshold and overshoot is eliminated. Lastly, we can find regimes where epidemic size does not monotonically decrease with a population that becomes increasingly risk-averse. ## 1 Introduction Recent outbreaks such as the COVID-19 pandemic, the 2014 Ebola outbreak, the 2009 influenza A (H1N1) pandemic, and the 2002 SARS epidemic brought to light many of the challenges of mounting an effective and unified epidemic response in a country as large and as diverse as the United States. Particularly during the COVID-19 pandemic, people were split in opinion on questions such as the origin of the virus [1], whether they would social distance or wear a mask [2–4], or whether the country should even have a pandemic response at all [5]. As time progressed, the situation became more dire and the death toll accumulated. People then had a new battery of questions to address, such as whether or not they would adhere to mandatory lock-downs [6–8] or whether they felt comfortable using the new mRNA vaccines [9, 10]. Compounding the issue were the multiple streams of information and potential misinformation spread through social media and other channels [11–14]. People’s stances on the questions and issues were diverse, arising from the milieu of differences in culture, geography, scientific education, sources of information, political leanings, and individual identity [15–17]. Taken altogether, these differences within the population reflect a spectrum of people’s risk tolerances to a circulating infectious disease. For any given intervention, such as social distancing, wearing a mask, or taking a vaccine, each person in the population falls somewhere on a spectrum of willingness to adopt the intervention. Given a threat level of an infectious disease in the population, some people will readily wear masks, whereas other people will refuse to. In this study, we aim to analyze the impact of heterogeneity in risk tolerance and the resulting behavioral response on the dynamics of epidemics. We seek to add to a burgeoning literature on the impact of human behavior in epidemic response [18–32], which the recent pandemic highlighted as an area for further exploration in preparation for the next large scale global health crisis [33–35]. To study the impact of heterogeneity in risk tolerance on epidemic dynamics, we introduce a simple and flexible modeling framework based on ordinary differential equations that can be used for different interventions and an arbitrary partitioning of the population with regard to risk tolerance and behavioral responses. We will examine and discuss potential interesting outcomes that can arise from coupling individual-level preferences and population-level epidemiology. ## 2 Results ### Model of Adaptive Intervention Usage under Heterogeneous Risk Tolerance Here we assume people’s risk aversion manifests as the rate at which they adopt individual interventions in response to an infectious disease outbreak. The intuition underlying this paradigm is that more risk-averse individuals are more sensitive to becoming sick and thus will adopt interventions at a faster rate than more risk-tolerant people. We consider the following SPIR compartmental model of a population with *n* differing levels of risk tolerance (1-4). This model features four types of classes: unprotected susceptible (*S*), protected susceptible (*P*), infectious (*I*), and recovered with permanent immunity (*R*). Since there are *n* differing levels of risk tolerance, we subdivide the susceptible population into *n* discrete groups indexed by *i*, where *i* ∈{1, 2, …, *n*}. Each tolerance level is characterized by an intervention adoption rate parameter (*λ**i*) and an intervention relaxation rate parameter (*δ**i*). Transitions of susceptibles between their unprotected class (*S**i*) and their corresponding protected class (*P**i*) are governed by the corresponding parameters of the same index (*λ**i*, *δ**i*). Overall, the system is governed by 3 + 2*n* parameters: a transmission rate parameter (*β*), a recovery rate parameter (*γ*), an intervention effectiveness parameter (*ϵ*), and an intervention adoption rate (*λ**i*) and intervention relaxation rate (*δ**i*) for each tolerance level. Since we always assume there are no protected individuals initially, the dynamics of the infected compartment is initially only a function of the transmission and recovery parameters. Thus we can use the typical definition for the basic reproduction number ![Graphic][1]. ![Formula][2] ![Formula][3] ![Formula][4] ![Formula][5] The transition from the unprotected susceptible state to the protected susceptible state represents individuals implementing an intervention that confers them protection against disease transmission from an infected individual. The rate at which intervention adoption occurs may be driven by individuals considering information such as the epidemic incidence rate (e.g. cases per day), the total number of infected individuals in the population (e.g. total number of active cases), and mortality rate (e.g. deaths per day) [22]. Here we assume that individuals have knowledge about the total number of infected individuals (*I*) and respond accordingly. Parameterizing each person’s individual risk tolerance by *λ**i*, we assume each individual person adopts an intervention at a rate *λ**i**I*. Then, if there are *S**i* number of people that behave exactly the same (i.e. have the same level of risk-aversion), then at the population scale there is a collective adoption rate of *λ**i**S**i**I*. The same reasoning holds for each of the *n* tolerance levels. We also consider a model where the adoption rate is driven by individuals reacting to the incidence rate (Supplemental Materials); while this produces a more complex mathematical model, the results are qualitatively similar. The effectiveness of the intervention being used is captured by the parameter *ϵ*, which linearly scales down the transmission rate between infected and protected susceptibles. In the limit of *ϵ* = 1, the intervention is perfectly effective and protected individuals cannot become infected. In the limit of *ϵ* = 0 then the intervention is completely ineffective, which reduces the model to an SIR model without interventions. For simplicity, we assume each epidemic features only a single type of intervention (whether that be masking, social distancing, vaccines, etc.) and that the effectiveness of an intervention is identical across the population. In reality, multiple interventions may be available concurrently, which would drive additional variation in behavior due to differences in risk sensitivity across the population. This model allows for protected individuals to relax their usage of interventions, becoming unprotected in the process. Here, individuals in the protected class can relax back to the unprotected class through an infection-independent rate (*δ**i**P**i*) that is governed by the intervention relaxation rate parameter (*δ**i*) for each tolerance level. In the limit of *δ**i* = 0, an intervention is irreversible, which would represent an intervention such as vaccines with permanent immunity. When *δ**i* is non-zero, individuals are using interventions such as masking or social distancing. This relaxation rate is motivated by factors such as psychological fatigue of social distancing [36, 37] and physical discomfort with wearing masks [38]. In general, we will consider the regime where the relaxation rate *δ* is of comparable scale or smaller than the transmission scale (i.e. *δ* ≤ *β*). This reflects intuition that people are likely to continue to protect themselves with interventions even beyond an initial outbreak [39]. For simplicity, we consider the model for the case when *n* = 1 and *n* = 2. A schematic for these two cases is shown in Figure 1. However, the framework is general and can be extended to any discrete number of groups. ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F1) Figure 1: Flow diagram for an SIR model with adaptive interventions for either (a) a population with homogeneous risk tolerance or (b) a heterogeneous population with two different levels of risk tolerance. Susceptible individuals can access a more protected susceptible state through usage of interventions. The transition rate to the protected state depends on the incidence level. The protected state offers a 1 − *ϵ* reduction in transmission rate over the normal susceptible state. For convention, when there are two susceptible classes, we assume the first susceptible class (*S*1) has a lower risk tolerance for becoming infected (i.e. more risk-averse). As a result, these individuals more readily adopt the intervention (i.e. *λ*1 > *λ*2), making individuals in this class transition more rapidly to the protected susceptible state (*P*1). The second susceptible class (*S*2) is more risk tolerant (i.e. more risk-taking), and thus is less eager to use the intervention, making individuals in this class transition more slowly to their protected susceptible state (*P*2). The model introduced here has several similarities and distinctions compared to existing behavior-disease models in the literature[19, 22, 28, 40–42], particularly to the global, prevalence-based spread of fear of disease of Perra et al. [43]. Our model shares the spirit of a disease contagion occuring at the same time as a social contagion. While the Perra model couches the social contagion in the language of individuals and fear, here we take the viewpoint that individuals have a personal characteristic called risk-tolerance, that encapsulates fear and other underlying biopsychosocial factors. The main feature of our model is the consideration of heterogeneity in the population, allowing for the population to be split into arbitrary factions of risk-tolerance. In addition, compared to the Perra model, the model here allows for people to remove protection at a rate independent of the disease. ### Adaptive Adoption of Interventions Can Produce Damped Oscillations The coupling of intervention usage to the incidence rate and the resulting adaptive changes enables the epidemic dynamics to display a much richer set of behavior over the simple SIR model. From Figure 2, we see this particular set of conditions can deterministically produce multiple waves of infection, even when vital dynamics (i.e. birth and death processes) are not considered. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F2) Figure 2: Time series for population with homogeneous risk tolerance and adaptive intervention usage and the corresponding phase space trajectory indicate the possibility for damped oscillations. Parameter values: *β* = 3, *γ* = 1, *ϵ* = 0.8, *λ* = 10, *δ* = 0.01. Initial conditions: *I*(0) = 10*−*6, *S*(0) = 1 − *I*(0), *P* (0) = *R*(0) = 0. Evidence for cycling of individuals between using interventions and not using interventions during the COVID pandemic can be seen in longitudinal usage [44–47] data. The possibility for these oscillations highlight the intimate connection between individual human behavior and intervention usage in shaping the dynamics of epidemics, while also be affected by the collective decision of everyone in the population. The coupling of behavior and epidemiology here provides a feedback mechanism where an increasing incidence rate prompts more individuals to adopt an intervention, which lowers the overall incidence rate; however, as the epidemic wanes and factors such as fatigue or discomfort set in, people begin dropping their usage of interventions, which may eventually lead to another wave of outbreaks if enough people become unprotected while infected individuals still remain, and then the cycle can be repeated. The oscillatory phenomenology here is reminiscent of several other behavioral models in the literature [48, 49]. ### Protective Mechanisms Saturate in Underdamped Regime that Eliminates Epidemic Overshoot One might have the intuition that having more people that will more readily adopt an intervention (i.e. mask, social distance, or vaccinate) or increasing the effectiveness of the intervention in reducing transmission will further decrease the size of the epidemic. Here we define epidemic size to be the cumulative number of infections over the course of the epidemic, or equivalently, the total depletion of susceptible people over the course of the epidemic. While we find this intuition to be mostly correct, we unexpectedly find that the protection conferred by either of these mechanisms can saturate once a critical parameter threshold has been passed. In Figure 3, *left*, we see that increasing the effectiveness of the intervention or increasing the fraction of the population that are risk-averse monotonically decreases the epidemic size. However, in the dark blue region (which we will refer to as the underdamped regime) where both protection mechanisms are at their highest, we see no further reduction in the epidemic size. This regime corresponds to an epidemic where the epidemic size exactly equals the herd immunity threshold, which is the minimum number of susceptibles that must be eliminated from the pool to prevent further outbreak from occurring. Since the epidemic size is the sum of the herd immunity threshold and the epidemic overshoot, this further implies that epidemic overshoot is eliminated in the underdamped regime. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F3) Figure 3: Left. Epidemic size as a function of varying the fraction of the population that are low-risk tolerance (i.e. those with higher *λ*). Right. Corresponding orbits in the I versus S+P plane for the sampled points in parameter space. Parameter values: *β* = 2, *γ* = 1, *λ*1 = 100, *δ*1 = 0.1, *λ*2 = 1, *δ*2 = 0.1,. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. The orbits of the dynamics from different areas of this parameter space are shown in Figure 3, *right*. We notice that orbits in the overdamped region (Orbits A, B, C) are parabolas and that each of those orbits ends at a different final number of susceptibles. In contrast, orbits in the underdamped regime (Orbits D, E, F) are qualitatively different in that they display long tails that converge to the same final number of susceptibles. We see that even though the final epidemic size is the same throughout the underdamped region, the trajectories to reach the same final epidemic size can look qualitatively different. Figures S4-S5 show a larger sampling of trajectories if one fixes either the fraction of the population with low risk tolerance or the intervention effectiveness respectively. It becomes clear that at the border of the underdamped region, we can see a clear change in the qualitative behavior of the trajectories as the threshold is crossed. Under some assumptions, one can prove that that the epidemic overshoot is eliminated in the underdamped regime (Supplemental Information). However, we should make the point that this is not evidence that highly effective interventions are a waste or that the overall population should tolerate risky behavior. As this is a model with a large parameter space, we cannot visualize all of it. If we could, we would find many parameter regimes where the protection mechanisms never reach a critical threshold, which implies the conventional intuition of increasing intervention effectiveness and having more risk-averse people always being beneficial applies. ### The Threshold to the Underdamped Regime is Reduced when Intervention Usage is Prolonged The transition to the underdamped regime is more easily accessed when the usage of interventions is prolonged (or equivalently when the rate at which protected individuals relax back into the regular susceptibility classes decreases). Consider the following scenario which is identical to the previous setup, except now the intervention reversion rate (*δP**i*) has been reduced by an order of magnitude (Figure 4). This corresponds to a scenario where people continue to use the intervention (i.e. such as wearing masks or social distancing) on a timescale significantly longer than the transmission timescale ![Graphic][6]. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F4) Figure 4: Final epidemic size versus fraction of population that are risk-averse (*S*1). Simulations in the left column have a higher *δ* than simulations in the right column. Parameter values: *R* = *β, γ* = 1, *λ*1 = 100, *λ*2 = 1. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. This suggests that increasing the timescale at which individuals continue to use interventions decreases the number of risk-averse individuals needed to achieve the same epidemic size. This is reflected in the horizontal shift of the transition region to the left when comparing Figures in the left column and the right column (Figure 4). We also briefly investigate what happens if we allow for heterogeneous relaxation rates between the difference groups, which is more reflective of what occurs in reality, finding the epidemic size to be more attune to the slower relaxation time scale (Supplemental Information). ### Increasing Heterogeneity in Risk Tolerance can Either Increase or Decrease the Epidemic Size There are many types of heterogeneity that may be present in the population, arising from differences in characteristics such as biological, cultural factors, and socioeconomic factors. The effect of each type of heterogeneity on disease dynamics is an area of active study. Some previous literature has suggested that increasing heterogeneity in the population through increasing the variation in their contact patterns, age, or general susceptibility results in a either a reduction or no change in the epidemic size [50–53]. We find in this model of heterogeneous behavior that it is possible to switch from a regime where increasing the heterogeneity in risk tolerance results in a decrease in epidemic size to a regime where increasing the heterogeneity in risk tolerance results in an increase in epidemic size. This result also does not have to necessarily be due to a dramatic shift in parameters. From Figure 5, we see this shift can arise from solely varying the fraction of the population with low risk tolerance by a small amount. ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F5) Figure 5: Epidemic size under differing levels of heterogeneity in the adoption rate for interventions. The mean adoption rate of the two groups (i.e. geometric average of *λ*1, *λ*2) is compared to the difference in the two adoption rates as parameterized by a homogeneity index (see Methods for definition). *Left* is when the fraction of the population with low risk tolerance (*x*1) is 0.2, *center* is when *x*1 = 0.35, *right* is when *x*1 = 0.5. Parameter values: *β* = 2, *γ* = 1, *ϵ* = 0.7, *δ*1 = *δ*2 = 0.5. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. The intuition underlying this result is that when the average adoption rate (*λ**average*), which is expressed as a (geometric or arithmetic) weighted mean of the adoption rates of the two groups, is fixed at a constant level, then the epidemic size can be suppressed through either varying the fraction of the population in each group or through varying each group’s adoption rate. When risk-averse people make up a smaller fraction of the population than risk-taking people, then it would be more beneficial in reducing epidemic size to have the adoption rates of the two groups be more similar (i.e. more homogeneous) as that would imply risk-taking people (which are then the majority of the population) would have a similar adoption rate to risk-averse people. As an increasingly larger fraction of the population is composed of risk-averse people, then it becomes increasingly beneficial in reducing epidemic size to have the adoption rates of the two groups be more different (i.e. more heterogeneous) as the deleterious effects of highly risk-taking people (which are then the minority of the population) can be mitigated by the large presence of risk-averse people. ### Epidemic Size Does Not Necessarily Decrease with an Increasing Number of Risk-Averse People General intuition would suggest that as one increases the number of risk-averse people in the population that the overall epidemic size would go down. However, the introduction of the adaptive behavior mechanism allows for regimes where this is no longer strictly the case. Thus, it is no longer a guarantee that decreasing the population’s overall risk tolerance will always improve epidemic outcomes. Consider Figure 6a, where we find a small region after the transition to the underdamped regime where there is an increase in the epidemic size when increasing the fraction of the population that are risk-averse. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F6) Figure 6: (a) Final epidemic size as a function of the proportion of the population that are risk averse (*S*1). Parameter values: *λ*1 = 10, *λ*2 = 0.5, *ϵ* = 1, *δ*1 = *δ*2 = 0.1. Initial conditions: *I*(0) = 10*−*7, *P*1(0) = *P*2(0) = 0, *R*(0) = 0. (b) Same as in (a) except now the effectiveness of the intervention (*ϵ*) is allowed to vary. If we also vary the effectiveness of the intervention as an additional axis (Figure 6b), we observe that there is a small trench in the threshold region surrounding the plateau area. This double descent suggests that the landscape can potentially be quite complicated when risk tolerance in the population is partitioned into even more groups. While the mechanism and intuition underlying this double descent remains to be determined, it is possible that it is caused by peculiar transient dynamics that are occurring at the beginning of outbreaks in that parameter regime [54]. ## 3 Discussion and Conclusions In this paper, we have proposed a simple model to model heterogeneity in risk tolerance levels in the population. We find that including a behavioral mechanism for adopting interventions that adapts with the level of infections greatly expands the variety in epidemic dynamics and outcomes that can occur. The general picture from the findings suggest that epidemic dynamics under adaptive intervention adoption fall into either an underdamped regime or an overdamped regime. The underdamped regime has a special property in which the epidemic size equals the herd immunity threshold exactly, which means no epidemic overshoot occurs. The system can be driven into this regime when protection mechanisms (such as numbers of risk-averse people, intervention effectiveness, and duration of intervention usage) are increased to a sufficiently high level. This regime is also marked by damped oscillations in the phase space of infecteds and susceptibles. In direct contrast, the overdamped regime closely resembles the dynamics of a simple *SIR* model without behavior, in which there are no oscillations and a non-zero overshoot, which makes the epidemic size greater than the herd immunity threshold. While we have limited our numerical exploration of the model to two groups of different risk-tolerance levels, in principle our framework can accommodate an arbitrary number of risk-tolerance groups. While the combinatorial increase in parameter space make it difficult to thoroughly explore larger models, we might naively expect there to be a complex landscape of outcomes, which might reveal other counter-intuitive phenomenology akin to the non-monotonic dependence of epidemic size with increasing number of risk-averse people for two groups. We have looked for some evidence in the historical data on outbreaks for these damped oscillations due to cycles in adoption and relaxation of interventions. While such data is in very limited supply, previous analysis suggested that relaxation of social distancing measures may have led to multiple waves of infection in the Spanish flu of the early 20th century [55, 56]. Dating back to the time of the bubonic plague, there is data from an outbreak in 1636 in the parish of St. Martin in the Fields, which showed how relaxation of quarantining and isolation measures lead to a smaller secondary wave of infections [57]. More recently, premature relaxation of interventions during the COVID pandemic might have led to secondary outbreaks [58]. The complex milieu of disease, biological, behavioral, political, economic, and cultural factors however make it is difficult to isolate cause-and-effect in real-world epidemics [59]. The results here are reminiscent of feedback control systems commonly studied in control theory. Here the set point is the herd immunity threshold, which is determined by the basic reproduction number (*R*). The ability for the population to reach this set point for epidemic size without additional overshoot depends on the effectiveness of the feedback mechanism from coupling intervention usage to the number of infected people. In the model presented, the adoption of interventions is a continuous process, in which the different groups are constantly reacting to the level of infections without requiring any notion of time or thresholds. In contrast, existing research on mitigation have considered more active control where activation and intervention timing play a key role [60–62]. Future work may explore how to synergistically utilize both active and continuous mechanisms for control. The inclusion of heterogeneity in risk tolerance and adaptive adoption of interventions leads to several unexpected conclusions. We find that increasing heterogeneity in risk tolerance levels in the population can lead to either an increase or decrease in the epidemic size. The direction of the trend depends nonlinearly on the composition of the population in terms of the ratio of risk-averse to risk-taking individuals and their respective intervention adoption rates. This adds to a small literature that demonstrates how heterogeneity can actually lead to a larger epidemic [27, 63]. Interestingly, these results on heterogeneity also can be used to address the question of whether distributed or centralized control of mitigation results in a smaller outbreak. Control of factors such as mobility may be more practically achieved in a more centralized and unified fashion [64, 65], whereas a distributed approach may be more appropriate when considering factors such as speed and individual agency. In a centralized scenario, a single entity controls the dynamics. That situation has an exact correspondence to the homogeneous population considered here, where all individuals respond in unison. Distributed control allows for more localized control, such as individuals or small groups deciding if they want to mask or social distance. This corresponds to the heterogeneous scenarios considered here, where there are multiple groups each with differing risk tolerance level. The results suggest that in some scenarios a single, coordinated response would be better for mitigation, whereas in other parameter regimes, a more decentralized strategy would be more optimal. We also find that increasing overall protection mechanisms does not always result in a monotonic decrease in epidemic size. In scenarios when the adoption rate begins to approach the transmission rate, near the critically damped boundary a nonmonotonicity can arise. This suggests that when intervention usage and effectiveness are tenuous, the dynamics become more complex and predicting what epidemic outcomes will result becomes significantly more difficult. Understanding how these nonlinear effects combine with other biological and behavioral heterogeneities will be important to explore in future work. ## 4 Methods ### Defining the *λ* Homogeneity Index The *λ* homogeneity index is defined as follows. We will assume the initial condition that at the beginning of the dynamics, the total population is composed of the fraction of the population in the low risk tolerance group *x*1, the high risk tolerance group *x*2, or the infected class. We will assume the fraction of the population initially infected is sufficiently small so that size of the two susceptible compartments is given by *x*1 and 1 − *x*1 respectively. Define a homogeneity index parameter (*c*) that captures the degree of homogeneity between the two adoption rates, *λ*1 and *λ*2. We define the parameter domain to be the unit interval, *c* ∈ (0, 1]. An index value of *c* = 1 indicates *λ*1 = *λ*2, while decreasing the index value towards 0 increases the difference between *λ*1 and *λ*2. We can define the average adoption rate as being either a geometric or arithmetic mean of the two adoption rates. The choice one makes is arbitrary, so we present prescriptions for both routes. In both cases, we will map the level of homogeneity to the unit interval. ### Geometric Average Let the geometric average of the two adoption rates be given by *λ*Geometric Average. ![Formula][7] Define *λ*2 as a fraction between 0 and 1 of the average *λ*. ![Formula][8] These two equations combine to define *λ*1. ![Formula][9] ### Arithmetic Average Let the arithmetic average of the two adoption rates be given by *λ*Arithmetic Average. ![Formula][10] Define *λ*2 as a fraction between 0 and 1 of the average *λ*. ![Formula][11] These two equations combine to define *λ*1. ![Formula][12] Again, an index value of 1 indicates *λ*1 = *λ*2, while decreasing the index value towards 0 increases the difference between *λ*1 and *λ*2. ### Homogeneity Index in Larger Groups The above parameterization of homogeneity is useful for in the setting of a two-group model because the entirety of the possible difference in heterogeneity between *λ*1 and *λ*2 can been mapped to the unit interval. This allows us to more easily see the qualitative shift in trend in phenomenology that prompted this section of text in the first place. An alternative parameterization might have used the variance between the two adoption rates. However, the domain of such a parameterization is unbounded, making it impossible to explore fully numerically. Moving beyond two groups would necessitate a different parameterization (such as through using a variance or coefficient-of-variation approach). ## Author Contributions M.M.N., B.E.C., C.M.S.-R., B.T.G., and S.A.L. designed the study. M.M.N., A.S.F., M.A.C. performed the simulations and calculations. M.M.N. wrote the manuscript. All authors contributed to interpreting the results and editing the manuscript. ## Data Availability The code used to generate all Figures is provided in the Supplemental Material. ## Funding M.M.N., A.S.F., M.A.C., and S.A.L. would like to acknowledge funding from NSF (CCF1917819, CNS-2041952, DMS-2327711), Army Research Office (W911NF-18-1-0325), and a gift from William H. Miller III. C.M.S.-R. acknowledges funding from the Miller Institute for Basic Research in Science of UC Berkeley via a Miller Research Fellowship. B.E.C. would like to acknowledge funding from NSF (IHBEM grant 2327710 and Expeditions NSF 1918656). B.T.G. would like to acknowledge the Princeton Catalysis Initiative and Princeton Precision Medicine. ## Competing Interests The authors declare no competing interests. ## Supplemental Material ### Analysis of the Model When Individuals Respond to the Incidence Rate The rate at which intervention adoption occurs may be driven by individuals considering information such as the epidemic incidence rate (e.g. cases per day), the total number of infected individuals in the population (e.g. total number of active cases), and mortality rate (e.g. deaths per day) [22]. Here we will consider the first case, where individuals adopt interventions based on the incidence rate for infections. Recall that the incidence rate is given by ![Graphic][13]. Parameterizing each person’s individual risk tolerance by *λ**i*, let us assume each individual adopts an intervention at rate ![Graphic][14]. Then, if there are *S**i* number of people that behave exactly the same (i.e. have the same level of risk-aversion), then at the population scale there is a collective adoption rate of ![Graphic][15]. The same reasoning holds for each of the *n* tolerance levels. The corresponding equations for this model are given by(11-14). ![Formula][16] ![Formula][17] ![Formula][18] ![Formula][19] When comparing the final epidemic size and epidemic trajectories between this model and the model in the main text where individuals adopt interventions at a rate that is based on the total number of infected people (Figures S4-S5), the results are indistinguishable (Figures S1-S2). ![Figure S1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F7.medium.gif) [Figure S1:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F7) Figure S1: Left. Epidemic size as a function of varying the fraction of the population that are low-risk tolerance (i.e. those with higher *λ*) when individuals react to incidence rate. Right. Corresponding orbits in the I versus S+P plane for the sampled points in parameter space when the fraction of the population that are low-risk tolerance has been fixed. Parameter values: *β* = 1.5, *γ* = 1, *λ*1 = 100, *δ*1 = 0.1, *λ*2 = 1, *δ*2 = 0.1,. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. Thus, given the equivalence in results, we have gone with the mathematically simpler and cleaner model based on total number of infected people in the main text. ![Figure S2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F8.medium.gif) [Figure S2:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F8) Figure S2: Left. Epidemic size as a function of varying the fraction of the population that are low-risk tolerance (i.e. those with higher *λ*) when individuals react to incidence rate. Right. Corresponding orbits in the I versus S+P plane for the sampled points in parameter space when the intervention effectiveness has been fixed. Parameter values: *β* = 1.5, *γ* = 1, *λ*1 = 100, *δ*1 = 0.1, *λ*2 = 1, *δ*2 = 0.1,. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. ### The Herd Immunity Threshold is Set by a Complex Interplay Between Transmission (*R*), Behavior, and Intervention Effectiveness (*ϵ*) While we could make some analytical calculations for the epidemic size in the homogeneous model, the heterogeneous two group case requires a numerical approach to find the plateau region. In general, it is set by a highly nonlinear interaction between the transmission (*R*), behavior as determined by the fraction of the population that are risk-averse and risk-taking, and the effectiveness of the intervention (*ϵ*). Consider the following progression of Figures (Figure S3), where in each subsequent Figure, the effectiveness of the intervention in blocking transmission (*ϵ*) is increasing. ![Figure S3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F9.medium.gif) [Figure S3:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F9) Figure S3: Final epidemic size versus fraction of population that are risk-averse (*S*1) with a progressive increase in intervention effectiveness (*ϵ*). Parameter values: *λ*1 = 100, *δ*1 = 1, *λ*2 = 1, *δ*2 = 1. Initial conditions: *I*(0) = 10*−*7, *P*1(0) = *P*2(0) = 0, *R*(0) = 0. ### Comparing Orbits Inside and Outside of the Plateau Region of Herd Immunity The following Figures sample more orbits for the Figure considered in the main text. ![Figure S4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F10.medium.gif) [Figure S4:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F10) Figure S4: Left. Epidemic size as a function of varying the fraction of the population that are low-risk tolerance (i.e. those with higher *λ*). Right. Corresponding orbits in the I versus S+P plane for the sampled points in parameter space when the fraction of the population that are low-risk tolerance has been fixed. Parameter values: *β* = 1.5, *γ* = 1, *λ*1 = 100, *δ*1 = 0.1, *λ*2 = 1, *δ*2 = 0.1,. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. ![Figure S5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F11.medium.gif) [Figure S5:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F11) Figure S5: Left. Epidemic size as a function of varying the fraction of the population that are low-risk tolerance (i.e. those with higher *λ*). Right. Corresponding orbits in the I versus S+P plane for the sampled points in parameter space when the intervention effectiveness has been fixed. Parameter values: *β* = 1.5, *γ* = 1, *λ*1 = 100, *δ*1 = 0.1, *λ*2 = 1, *δ*2 = 0.1,. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. ### Proving Underdamped Regime Eliminates Epidemic Overshoot In general, the nonlinear feedback between the protected classes and infected individuals make it difficult to make analytical calculations in the full model. Under some simplifications however, we can make some progress. In this section, we derive the epidemic size at which the protection saturates in a restricted case of the homogeneous model (*n* = 1). Let us consider the homogeneous model in the limit of an intervention with perfect effectiveness (i.e. *ϵ* = 1). We will consider the case where the recovery rate from infection and relaxation rate for interventions are comparable (i.e. *γ* = *δ*). Since the equation for recovered individuals can be ignored since the population is closed (*S* + *P* + *I* + *R* = 1), this reduces the dynamics to the following system: ![Formula][20] ![Formula][21] ![Formula][22] ![Formula][23] The initial conditions sets the number of individuals initially susceptible to be *f*, the number of initially infected is assumed to be small *I*(0) *<<* 1, and the remainder of the population is already immune to infection (i.e. recovered). To ensure that the epidemic initially grows in size, we assume that ![Graphic][24]. This follows from comparing the incidence term (*βSI*) to the recovery term (*γI*) in17. To find the asymptotic behavior for *S*, we will attempt to eliminate *P* from(15). We start first by seeking an equation that relates the *P* and *I* compartments. Consider the following ansatz that considers the difference between the two compartments: ![Formula][25] Differentiating this equation with respect to time and using(16)-(17) yields: ![Formula][26] Using(19), this simplifies to ![Graphic][27], which has a critical point at *x* = 0. At *x* = 0, we obtain that ![Graphic][28], indicating a regime where the behavior of *P* and *I* scale linearly with each other. Combining this with(15) yields: ![Formula][29] Now that we have an equation that is linear in *I*, we are in a good position to find a final size relationship for the number of susceptibles. To start we take the ratio of(17) and(21). ![Formula][30] Using the partial fractions ![Graphic][31] and ![Graphic][32], we get upon indefinite integration of(22) that ![Graphic][33], where *k* is a constant that holds throughout the trajectory of the dynamics. Thus, considering the values of *S* and *I* at the beginning of the epidemic (*t* = 0) and the end of the epidemic (*t* =) and using the conditions that *I*(∞) = 0 and *I*(0) ≈ 0 yields the following transcendental equation for the final epidemic size. ![Formula][34] Since ![Graphic][35], then the argument of the logarithm on the left hand side must be positive, and subsequently the left hand side evaluates to a real number. Due to the equality, the right hand side must also evaluate to a real number, implying the argument of the logarithm on the right hand side must also be positive. Positivity implies the following inequality for the lower bound for the final number of susceptibles: ![Formula][36] An upper bound can be given by simply noting that in the long time limit, the recovery term (*γI*) must be at least as large as the incidence term (*βS*(∞)*I*) in(17), otherwise the epidemic would still be growing. This implies ![Graphic][37]. To summarize, when the interventions are perfectly effective, the rate of relaxation from the protected class is equal to the rate of recovery from infection, and the number in the protected class scales linearly with the number of infected, then the final fraction of susceptibles is bounded as follows: ![Formula][38] We see in the parameter limit of when the adoption rate of interventions is very fast compared to the transmission rate (i.e. *λ* >> *β*), that the lower bound reduces to ![Graphic][39]. Since both bounds now coincide, then *S**∞* must equal that value. Interestingly this corresponds to the herd immunity threshold of the standard SIR model. As the overshoot is the excess number of cases beyond the herd immunity threshold, we see that in this parameter limit there is no overshoot. This analysis for the homogeneous case also carries over to the heterogeneous case for two groups when the adoption rate between the two groups is significantly different (i.e. *λ*1 >> *λ*2). This results in a separation of time scales in which the faster adopters quickly transition to the protected state and can essentially be treated as immune over the course of the remaining epidemic over the slow adopters. This amounts to effectively reducing the dynamics to the homogeneous model considered here where *f* and 1− *f* fractions of the population in the susceptible and recovered respectively correspond to the fraction of the population in the slow adopter (*λ*2) and fast adopter groups (*λ*1). ### Heterogeneity in Risk Tolerance through Arithmetic Averaging There are two ways for calculating the difference (heterogeneity) in adoption rates, either through geometric or arithmetic averaging. Both can be justified, and we presented the geometric formulation in the main text. We find that the arithmetic formulation gives qualitative similar results under a suitable parameter shift (Figure S6). ![Figure S6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F12.medium.gif) [Figure S6:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F12) Figure S6: Epidemic size under differing levels of heterogeneity in the adoption rate for interventions. The mean adoption rate of the two groups (i.e. arithmetic average of *λ*1, *λ*2) is compared to the difference in the two adoption rates as parameterized by a homogeneity index (see Methods for definition). *Left* is when the fraction of the population with low risk tolerance (*x*1) is 0.2, *center* is when *x*1 = 0.5, *right* is when *x*1 = 0.8. Parameter values: *β* = 2, *γ* = 1, *ϵ* = 0.7, *δ*1 = *δ*2 = 0.5. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. ### Heterogeneity in Relaxation Rates of Protection In the main text we primarily assumed that the relaxation rate in both groups were homogeneous for simplicity (i.e. *δ*1 = *δ*2). However, in the real world, one might actually expect for them to be inversely correlated with the intervention adoption rates, such that if *λ*1 > *λ*2, then *δ*1 *< δ*2. The intuition is that a more risk-averse person will adopt protection faster and relax the protection at a slower rate. We numerically study if we allow for heterogeneous relaxation rates. As one might expect, our numerical results that the behavior of the epidemic size is intermediate between the two homogeneous extremes. (Figure S7). ![Figure S7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F13.medium.gif) [Figure S7:](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F13) Figure S7: Final epidemic size versus fraction of population that are risk-averse (*S*1). Simulations in the left column have a higher *δ* than simulations in the right column. Parameter values: *R* = *β, γ* = 1, *λ*1 = 100, *λ*2 = 1. Initial conditions: *I*(0) = 10*−*8, *P*1(0) = *P*2(0) = *R*(0) = 0. ## Code to Generate Figures Code executed in MATLAB R2023a. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-39.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-39) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-40.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-40) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-41.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-41) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-42.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-42) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-43.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-43) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-44.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-44) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-45.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-45) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-46.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-46) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-47.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-47) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-48.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-48) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-49.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-49) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-50.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-50) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/09/18/2024.07.01.24309771/F14/graphic-51.medium.gif) [](http://medrxiv.org/content/early/2024/09/18/2024.07.01.24309771/F14/graphic-51) ## Acknowledgments The ideas and conversations that sparked this work originated during a special session on the Mathematics of Infectious Diseases: A Session in Memory of Dr. Abdul-Aziz Yakubu, III at the 2024 Spring Sectional Meeting of the American Mathematical Society. The authors would like to thank the organizers of that session: Abba Gumel, Daniel Cooney, and Chadi Saad-Roy. ## Footnotes * mmnguyen{at}princeton.edu * arisf{at}princeton.edu * matthew.cheung{at}princeton.edu * csaadroy{at}berkeley.edu * baltazar.espinoza{at}virginia.edu * grenfell{at}princeton.edu * slevin{at}princeton.edu * Updated figures and expanded discussion throughout the manuscript. * Received July 1, 2024. * Revision received September 18, 2024. * Accepted September 18, 2024. * © 2024, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Maxmen, A. & Mallapaty, S. The COVID lab-leak hypothesis: what scientists do and don’t know. Nature 594, 313–315 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/d41586-021-01529-3&link_type=DOI) 2. 2.Betsch, C. et al. Social and behavioral consequences of mask policies during the COVID-19 pandemic. Proceedings of the National Academy of Sciences 117, 21851–21853 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzM2LzIxODUxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTgvMjAyNC4wNy4wMS4yNDMwOTc3MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 3. 3.Fischer, C. B. et al. Mask adherence and rate of COVID-19 across the United States. PLOS ONE 16, e0249891 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0249891&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33852626&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 4. 4.Yang, L. et al. Sociocultural determinants of global mask-wearing behavior. Proceedings of the National Academy of Sciences 119, e2213525119 (2022). 5. 5.Morens, D. M., Folkers, G. K. & Fauci, A. S. The Concept of Classical Herd Immunity May Not Apply to COVID-19. The Journal of Infectious Diseases 226, 195–198 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiac109&link_type=DOI) 6. 6.Wong, C. M. L. & Jensen, O. The paradox of trust: perceived risk and public compliance during the COVID-19 pandemic in Singapore. Journal of Risk Research 23, 1021–1030 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/13669877.2020.1756386&link_type=DOI) 7. 7.Brzezinski, A., Kecht, V., Van Dijcke, D. & Wright, A. L. Science skepticism reduced compliance with COVID-19 shelter-in-place policies in the United States. Nature Human Behaviour 5, 1519–1527 (2021). 8. 8.Kleitman, S. et al. To comply or not comply? A latent profile analysis of behaviours and attitudes during the COVID-19 pandemic. PLOS ONE 16, e0255268 (2021). 9. 9.Machingaidze, S. & Wiysonge, C. S. Understanding COVID-19 vaccine hesitancy. Nature Medicine 27, 1338– 1339 (2021). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 10. 10.Fedele, F. et al. COVID-19 vaccine hesitancy: a survey in a population highly compliant to common vaccinations. Human Vaccines & Immunotherapeutics 17, 3348–3354 (2021). 11. 11.Gallotti, R., Valle, F., Castaldo, N., Sacco, P. & De Domenico, M. Assessing the risks of ‘infodemics’ in response to COVID-19 epidemics. Nature Human Behaviour 4, 1285–1293 (2020). 12. 12.Loomba, S., de Figueiredo, A., Piatek, S. J., de Graaf, K. & Larson, H. J. Measuring the impact of COVID-19 vaccine misinformation on vaccination intent in the UK and USA. Nature Human Behaviour 5, 337–348 (2021). 13. 13.Offer-Westort, M., Rosenzweig, L. R. & Athey, S. Battling the coronavirus ‘infodemic’ among social media users in Kenya and Nigeria. Nature Human Behaviour, 1–12 (2024). 14. 14.Towers, S. et al. Mass Media and the Contagion of Fear: The Case of Ebola in America. PLOS ONE 10, e0129179 (2015). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0129179&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26067433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 15. 15.1. Olson, J. M. & 2. Zanna, M. P. Murray, D. R. & Schaller, M. in Advances in Experimental Social Psychology (eds Olson, J. M. & Zanna, M. P. 75–129 (Academic Press, 2016). 16. 16.Kramer, P. & Bressan, P. Infection threat shapes our social instincts. Behavioral Ecology and Sociobiology 75, 47 (2021). 17. 17.Lu, J. G., Jin, P. & English, A. S. Collectivism predicts mask use during COVID-19. Proceedings of the National Academy of Sciences 118, e2021793118 (2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxODoiMTE4LzIzL2UyMDIxNzkzMTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTgvMjAyNC4wNy4wMS4yNDMwOTc3MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 18. 18.Bansal, S., Grenfell, B. T. & Meyers, L. A. When individual behaviour matters: homogeneous and network models in epidemiology. Journal of The Royal Society Interface 4, 879–891 (2007). 19. 19.Funk, S., Salathé, M. & Jansen, V. A. A. Modelling the influence of human behaviour on the spread of infectious diseases: a review. Journal of The Royal Society Interface 7, 1247–1256 (2010). 20. 20.Fenichel, E. P. et al. Adaptive human behavior in epidemiological models. Proceedings of the National Academy of Sciences 108, 6306–6311 (2011). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA4LzE1LzYzMDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8xOC8yMDI0LjA3LjAxLjI0MzA5NzcxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 21. 21.Edmunds, W. J., Medley, G. F. & Nokes, D. J. Evaluating the cost-effectiveness of vaccination programmes: a dynamic perspective. Statistics in Medicine 18, 3263–3282 (1999). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/(SICI)1097-0258(19991215)18:23<3263::AID-SIM315>3.0.CO;2-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10602150&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000084077400010&link_type=ISI) 22. 22.Weitz, J. S., Park, S. W., Eksin, C. & Dushoff, J. Awareness-driven behavior changes can shift the shape of epidemics away from peaks and toward plateaus, shoulders, and oscillations. Proceedings of the National Academy of Sciences 117, 32764–32771 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzUxLzMyNzY0IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTgvMjAyNC4wNy4wMS4yNDMwOTc3MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 23. 23.Wagner, C. E. et al. Economic and Behavioral Influencers of Vaccination and Antimicrobial Use. Frontiers in Public Health 8 (2020). 24. 24.Tyson, R. C., Hamilton, S. D., Lo, A. S., Baumgaertner, B. O. & Krone, S. M. The Timing and Nature of Behavioural Responses Affect the Course of an Epidemic. Bulletin of Mathematical Biology 82, 14 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11538-019-00684-z&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31932981&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 25. 25.Espinoza, B., Marathe, M., Swarup, S. & Thakur, M. Asymptomatic individuals can increase the final epidemic size under adaptive human behavior. Scientific Reports 11, 19744 (2021). 26. 26.Wagner, C. E., Saad-Roy, C. M. & Grenfell, B. T. Modelling vaccination strategies for COVID-19. Nature Reviews Immunology 22, 139–141 (2022). 27. 27.Espinoza, B., Swarup, S., Barrett, C. L. & Marathe, M. Heterogeneous adaptive behavioral responses may increase epidemic burden. Scientific Reports 12, 11276 (2022). 28. 28.Qiu, Z. et al. Understanding the coevolution of mask wearing and epidemics: A network perspective. Proceedings of the National Academy of Sciences 119, e2123355119 (2022). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 29. 29.Traulsen, A., Levin, S. A. & Saad-Roy, C. M. Individual costs and societal benefits of interventions during the COVID-19 pandemic. Proceedings of the National Academy of Sciences 120, e2303546120 (2023). 30. 30.Saad-Roy, C. M. & Traulsen, A. Dynamics in a behavioral–epidemiological model for individual adherence to a nonpharmaceutical intervention. Proceedings of the National Academy of Sciences 120, e2311584120 (2023). 31. 31.Smith, R. A. et al. COVID-19 Mitigation Among College Students: Social Influences, Behavioral Spillover, and Antibody Results. Health Communication 38, 2002–2011 (2023). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/10410236.2022.2049047&link_type=DOI) 32. 32.Kiss, I. Z., Miller, J. C. & Simon, P. L. Mathematics of Epidemics on Networks: From Exact to Approximate Models (Springer International Publishing, Cham, 2017). 33. 33.Morse, S. S. et al. Prediction and prevention of the next pandemic zoonosis. The Lancet 380, 1956–1965 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0140-6736(12)61684-5&link_type=DOI) 34. 34.Osterholm, M. T. in The Covid-19 Reader (Routledge, 2020). 35. 35.Bergstrom, C. T. & Hanage, W. P. Human behavior and disease dynamics. Proceedings of the National Academy of Sciences 121, e2317211120 (2024). 36. 36.Franzen, A. & Wöhner, F. Fatigue during the COVID-19 pandemic: Evidence of social distancing adherence from a panel study of young adults in Switzerland. PLOS ONE 16, e0261276 (2021). 37. 37.Jørgensen, F., Bor, A., Rasmussen, M. S., Lindholt, M. F. & Petersen, M. B. Pandemic fatigue fueled political discontent during the COVID-19 pandemic. Proceedings of the National Academy of Sciences 119, e2201266119 (2022). 38. 38.Cheok, G. J. W. et al. Appropriate attitude promotes mask wearing in spite of a significant experience of varying discomfort. Infection, Disease & Health 26, 145–151 (2021). 39. 39.Barak, D., Gallo, E., Rong, K., Tang, K. & Du, W. Experience of the COVID-19 pandemic in Wuhan leads to a lasting increase in social distancing. Scientific Reports 12, 18457 (2022). 40. 40.Funk, S. et al. Nine challenges in incorporating the dynamics of behaviour in infectious diseases models. Epidemics. Challenges in Modelling Infectious DIsease Dynamics 10, 21–25 (2015). 41. 41.Wang, Z., Andrews, M. A., Wu, Z.-X., Wang, L. & Bauch, C. T. Coupled disease–behavior dynamics on complex networks: A review. Physics of Life Reviews 15, 1–29 (2015). 42. 42.Verelst, F., Willem, L. & Beutels, P. Behavioural change models for infectious disease transmission: a systematic review (2010–2015). Journal of The Royal Society Interface 13, 20160820 (2016). 43. 43.Perra, N., Balcan, D., Gonçalves, B. & Vespignani, A. Towards a Characterization of Behavior-Disease Models. PLOS ONE 6, e23084 (2011). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0023084&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21826228&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 44. 44.Mathieu, E. et al. Coronavirus Pandemic (COVID-19). Our World in Data (2020). 45. 45.Hale, T. et al. A global panel database of pandemic policies (Oxford COVID-19 Government Response Tracker). Nature Human Behaviour 5, 529–538 (2021). 46. 46.Rader, B. et al. Mask-wearing and control of SARS-CoV-2 transmission in the USA: a cross-sectional study. The Lancet Digital Health 3, e148–e157 (2021). 47. 47.Salomon, J. A. et al. The US COVID-19 Trends and Impact Survey: Continuous real-time measurement of COVID-19 symptoms, risks, protective behaviors, testing, and vaccination. Proceedings of the National Academy of Sciences 118, e2111454118 (2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxODoiMTE4LzUxL2UyMTExNDU0MTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTgvMjAyNC4wNy4wMS4yNDMwOTc3MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 48. 48.Liu, M., Liz, E. & Röst, G. Endemic Bubbles Generated by Delayed Behavioral Response: Global Stability and Bifurcation Switches in an SIS Model. SIAM Journal on Applied Mathematics 75, 75–91 (2015). 49. 49.Manrubia, S. & Zanette, D. H. Individual risk-aversion responses tune epidemics to critical transmissibility (R = 1). Royal Society Open Science 9, 211667 (2022). 50. 50.Miller, J. C. Epidemic size and probability in populations with heterogeneous infectivity and susceptibility. Physical Review E 76, 010101 (2007). 51. 51.Britton, T., Ball, F. & Trapman, P. A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2. Science 369, 846–849 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNjkvNjUwNS84NDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8xOC8yMDI0LjA3LjAxLjI0MzA5NzcxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 52. 52.Gomes, M. G. M. et al. Individual variation in susceptibility or exposure to SARS-CoV-2 lowers the herd immunity threshold. Journal of Theoretical Biology 540, 111063 (2022). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2022.111063&link_type=DOI) 53. 53.Allard, A., Moore, C., Scarpino, S. V., Althouse, B. M. & Hébert-Dufresne, L. The Role of Directionality, Heterogeneity, and Correlations in Epidemic Risk and Spread. SIAM Review 65, 471–492 (2023). 54. 54.Hastings, A. et al. Transient phenomena in ecology. Science 361, eaat6412 (2018). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjEvNjQwNi9lYWF0NjQxMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA5LzE4LzIwMjQuMDcuMDEuMjQzMDk3NzEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 55. 55.Hatchett, R. J., Mecher, C. E. & Lipsitch, M. Public health interventions and epidemic intensity during the 1918 influenza pandemic. Proceedings of the National Academy of Sciences 104, 7582–7587 (2007). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTA0LzE4Lzc1ODIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOS8xOC8yMDI0LjA3LjAxLjI0MzA5NzcxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 56. 56.Caley, P., Philp, D. J. & McCracken, K. Quantifying social distancing arising from pandemic influenza. Journal of The Royal Society Interface 5, 631–639 (2007). 57. 57.Newman, K. L. S. Shutt Up: Bubonic Plague and Quarantine in Early Modern England. Journal of Social History 45, 809–834 (2012). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jsh/shr114&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22611587&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 58. 58.Tsai, A. C., Harling, G., Reynolds, Z., Gilbert, R. F. & Siedner, M. J. Coronavirus Disease 2019 (COVID-19) Transmission in the United States Before Versus After Relaxation of Statewide Social Distancing Measures. Clinical Infectious Diseases 73, S120–S126 (Supplement_2 2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa1502&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33009800&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) 59. 59.Caulkins, J. P. et al. The optimal lockdown intensity for COVID-19. Journal of Mathematical Economics. The economics of epidemics and emerging diseases 93, 102489 (2021). 60. 60.Bussell, E. H., Dangerfield, C. E., Gilligan, C. A. & Cunniffe, N. J. Applying optimal control theory to complex epidemiological models to inform real-world disease management. Philosophical Transactions of the Royal Society B: Biological Sciences 374, 20180284 (2019). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2018.0284&link_type=DOI) 61. 61.Lauro, F. D., Kiss, I. Z. & Miller, J. C. Optimal timing of one-shot interventions for epidemic control. PLOS Computational Biology 17, e1008763 (2021). 62. 62.Morris, D. H., Rossine, F. W., Plotkin, J. B. & Levin, S. A. Optimal, near-optimal, and robust epidemic control. Communications Physics 4, 1–8 (2021). 63. 63.Volz, E. M., Miller, J. C., Galvani, A. & Meyers, L. A. Effects of Heterogeneous and Clustered Contact Patterns on Infectious Disease Dynamics. PLOS Computational Biology 7, e1002042 (2011). 64. 64.Bonaccorsi, G. et al. Economic and social consequences of human mobility restrictions under COVID-19. Proceedings of the National Academy of Sciences 117, 15530–15535 (2020). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzI3LzE1NTMwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDkvMTgvMjAyNC4wNy4wMS4yNDMwOTc3MS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 65. 65.Espinoza, B., Castillo-Chavez, C. & Perrings, C. Mobility restrictions for the control of epidemics: When do they work? PLOS ONE 15, e0235731 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0235731&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F09%2F18%2F2024.07.01.24309771.atom) [1]: /embed/inline-graphic-1.gif [2]: /embed/graphic-1.gif [3]: /embed/graphic-2.gif [4]: /embed/graphic-3.gif [5]: /embed/graphic-4.gif [6]: /embed/inline-graphic-2.gif [7]: /embed/graphic-11.gif [8]: /embed/graphic-12.gif [9]: /embed/graphic-13.gif [10]: /embed/graphic-14.gif [11]: /embed/graphic-15.gif [12]: /embed/graphic-16.gif [13]: /embed/inline-graphic-3.gif [14]: /embed/inline-graphic-4.gif [15]: /embed/inline-graphic-5.gif [16]: /embed/graphic-17.gif [17]: /embed/graphic-18.gif [18]: /embed/graphic-19.gif [19]: /embed/graphic-20.gif [20]: /embed/graphic-26.gif [21]: /embed/graphic-27.gif [22]: /embed/graphic-28.gif [23]: /embed/graphic-29.gif [24]: /embed/inline-graphic-6.gif [25]: /embed/graphic-30.gif [26]: /embed/graphic-31.gif [27]: /embed/inline-graphic-7.gif [28]: /embed/inline-graphic-8.gif [29]: /embed/graphic-32.gif [30]: /embed/graphic-33.gif [31]: /embed/inline-graphic-9.gif [32]: /embed/inline-graphic-10.gif [33]: /embed/inline-graphic-11.gif [34]: /embed/graphic-34.gif [35]: /embed/inline-graphic-12.gif [36]: /embed/graphic-35.gif [37]: /embed/inline-graphic-13.gif [38]: /embed/graphic-36.gif [39]: /embed/inline-graphic-14.gif