key: cord-0948503-mhw15krf authors: Khanjanianpak, Mozhgan; Azimi-Tafreshi, Nahid; Arenas, Alex; Gómez-Gardeñes, Jesús title: Emergence of protective behaviour under different risk perceptions to disease spreading date: 2022-07-11 journal: Philosophical transactions. Series A, Mathematical, physical, and engineering sciences DOI: 10.1098/rsta.2020.0412 sha: fc33bc3ef047553f7137d59eaab1daea5fa05dde doc_id: 948503 cord_uid: mhw15krf The behaviour of individuals is a main actor in the control of the spread of a communicable disease and, in turn, the spread of an infectious disease can trigger behavioural changes in a population. Here, we study the emergence of individuals’ protective behaviours in response to the spread of a disease by considering two different social attitudes within the same population: concerned and risky. Generally speaking, concerned individuals have a larger risk aversion than risky individuals. To study the emergence of protective behaviours, we couple, to the epidemic evolution of a susceptible-infected-susceptible model, a decision game based on the perceived risk of infection. Using this framework, we find the effect of the protection strategy on the epidemic threshold for each of the two subpopulations (concerned and risky), and study under which conditions risky individuals are persuaded to protect themselves or, on the contrary, can take advantage of a herd immunity by remaining healthy without protecting themselves, thanks to the shield provided by concerned individuals. This article is part of the theme issue ‘Emergent phenomena in complex physical and socio-technical systems: from cells to societies’. The behaviour of individuals is a main actor in the control of the spread of a communicable disease and, in turn, the spread of an infectious disease can trigger behavioural changes in a population. Here, we study the emergence of individuals' protective behaviours in response to the spread of a disease by considering two different social attitudes within the same population: concerned and risky. Generally speaking, concerned individuals have a larger risk aversion than risky individuals. To study the emergence of protective behaviours, we couple, to the epidemic evolution of a susceptible-infectedsusceptible model, a decision game based on the perceived risk of infection. Using this framework, we find the effect of the protection strategy on the epidemic threshold for each of the two subpopulations (concerned and risky), and study under which conditions risky individuals are persuaded to protect themselves or, on the contrary, can take advantage 2022 The Authors. Published by the Royal Society under the terms of the Creative Commons Attribution License http://creativecommons.org/licenses/ by/4.0/, which permits unrestricted use, provided the original author and source are credited. The prevention and control of the spread of communicable diseases have always constituted a fundamental challenge in human societies [1] [2] [3] [4] . Today, the COVID-19 pandemic has spread all over the world and, as of June 2021, SARS-CoV-2 has infected more than 175 million people, causing around 3.8 million deaths. Typically, vaccination is one of the most important preventive measures to prevent or reduce virus propagation and, when its availability is limited, one of the most important issues to study is the effectiveness of different kinds of vaccination strategies aimed at cutting off potential chains of transmission and avoiding as many potential deaths as possible [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] . However, when vaccines are unavailable or scarce, it is social behaviour combined with preventive measures that is the most effective way to reduce and control the spread of the disease [15] [16] [17] . Such preventive strategies include quarantine, self-isolation, social distancing and the use of prophylactic tools such as face masks. However, most of these behavioural changes and protective measures have associated social and economic costs. Therefore, when their use or application is not subject to law enforcement, regulatory policies or economic support, their application highly depends on an individual decision process. In this scenario, the dynamics of disease spread and the dynamics of individual decision making must be considered as two processes coevolving simultaneously. From a game theoretical point of view, it seems that behavioural adaptation could be coherently formulated as a game on a complex social network to better understand the evolution of human decisions in response to a given risk, here the spread of a pathogen [18] . Several studies have applied game theoretical frameworks coupled to epidemic spreading on populations wherein each rational individual tries to maximize his own payoff according to a self-evaluation of the cost-benefit ratio associated with protective measures [19] [20] [21] [22] [23] . Although the first studies along these lines constitute a combination of epidemic dynamics with a static game theory [24] , the evolutionary game benchmark provides a better description, since it allows agents to update their strategies by evaluating the most successful one according to the current epidemiological state [25] [26] [27] [28] [29] [30] [31] . This framework appears to be the most appropriate when individual decisions evolve in parallel with the development of an epidemic, allowing each individual to adopt the protection strategy that provides the greatest reward based on the perceived health risk [32, 33] . Most of the studies devoted to the adoption of prophylactic measures coupled to the unfolding of an epidemic implicitly assume a homogeneous perception of the risk of contagion throughout the population. However, people act differently, reacting heterogeneously to the same information, as evidenced during the current COVID-19 pandemic [34] . For example, some individuals are more concerned about the evolution of the incidence of cases, trying to be as up to date as possible to take action in an appropriate and responsible manner. On the contrary, others behave in a careless and unconcerned manner or, in the worst case, are critical and suspicious of the alarm messages issued by the health authorities. In [35, 36] , the authors have studied the effect of these two sub-populations on the spread of infectious diseases. They assume that the behavioural attitude of these two groups can be changed during the epidemic and find conditions under which the behavioural attitudes can mitigate the disease outbreak. In this paper, we consider the epidemic model based on risk perception introduced in [33] and generalize it to the case when individual reactions to the risk of infection and its associated cost are heterogeneous. In [33] , the interaction between individual decisions based on the perception of epidemic risk and the spread of the infectious disease was shown to lead to sustained oscillations over time in the degree of protection adopted by the community. In this study, we extend the model by including two different subpopulations: those who are concerned and those who are unconcerned/risky. Thus, faced with the same level of alarm, the first group perceives a much greater risk than the second and is therefore much more likely to take protective measures. To perform our study in a controlled manner, we assume that the partition between subpopulations is constant over time and introduce a parameter that controls the degree of diversity between the two attitudes to the same epidemiological situation. In this way, we show that risky agents start to be protected if the cost of contracting the disease, and the degree of diversity between the two attitudes, exceeds certain threshold values. We also show how the fraction of concerned agents can affect the protective and epidemic thresholds of each subpopulation. Finally, we show that herd immunity of risky agents is provided if the number of concerned agents is above a threshold. The article is organized as follows. In §2, we define our model. In the framework of the Microscopic Markov Chain Approach (MMCA), we derive the dynamic equations of the model. In §3 we present the main results of our work. In §3a, assuming that the fraction of concerned and risky individuals are equal, we find the phase diagrams of the model for the protected and infected fractions of each subpopulation. These results are generalized for different fractions of concerned and unconcerned individuals in §3b, and also the particular case in which risky agents never take protective measures (thus behaving like epidemic deniers) referred as the zealot's limit, is analysed in §3c. The article is concluded in §4. We consider a population of N individuals connected in pairs and forming a complex network. The infectious pathogen spreads from individual to individual following the connections defined by the social graph and according to the dynamical rules of the SIS (Susceptible-Infected-Susceptible) epidemic model, in which agents can be susceptible (S) or infected (I). Thus, a susceptible individual at time t becomes infectious from an infectious neighbour with probability λ while an infectious agent becomes susceptible with probability μ. The classical SIS epidemic model is modified to incorporate the possibility that some agents decide to adopt a prophylactic measure. In particular, we are interested in intermediate measures, such as the use of a mask, that offers partial protection against infection. To this end, we consider a parameter, γ ∈ [0, 1], that determines the descent in the probability of contagion, which is transformed into γ λ when in the encounter between an infectious agent and a susceptible one of the two adopts the prophylactic strategy. Thus, γ = 0 would imply that the preventive measure is perfect while γ = 1 means that the measure is completely useless. Similar to [33] , we assume that the preventive mechanism is not increased by the bilateral use of protective measures, and thus the reduction in infectivity is linear, γ λ, rather than quadratic, γ 2 λ, in the case where both individuals (infectious and susceptible) are protected. The use of the prophylactic measure depends on the individual free-will. To address the dynamics associated to this individual choice we couple, as in [33] , the SIS epidemic model to a two-strategy Protected-Unprotected (P-NP) decision game in which agents decide whether or not to adopt a prophylactic measure based on the associated costs and the disease incidence for each strategy. In particular, the choice of whether or not to take the prophylactic measure carries an associated cost and individuals must therefore assess the appropriateness of its use based on their perception of the risk of contagion. To account for the decision process, agents evaluate the strategic choice (P or NP) based on: (i) the cost of contracting the disease (T), (ii) the cost of the prophylactic measure (c), and (iii) the risk of contracting the disease, measured through its incidence (i.e. the fraction of infectious individuals I). The evaluation of the contagion risk at each time step is used to compute the expected benefit of each strategy, the so-called payoffs, P p and P np , that in their turn are used to construct the probabilities that drive strategic changes Γ p→np and Γ np→p . These probabilities change in time according to the epidemic incidence and thus govern the evolution of the strategic partition of the population. The strategic partition strongly influences the epidemic evolution, being the spreading of the pathogen favoured when NP strategy dominates or being its propagation mitigated when prophylaxis is widely adopted. A schematic plot of the coevolution of the spreading and decision processes is shown in figure 1 . Finally, the main novelty of this work lies in dividing the population into two groups: Concerned (C) and Risky/unconcerned (R). These groups have a different perception of risk and, therefore, under the same epidemic scenario show different behavioural responses. For example, concerned agents have a higher perception of the cost of contracting the disease than risky individuals (T C > T R ) while the cost associated with protection is perceived as lower by the former than by the latter (c R > c C ). Thus, on one hand, concerned individuals perceive a higher risk associated with the disease and a lower cost of protection so, consequently, they are inclined to adopt prophylactic measures. On the other hand, risk-takers are more relaxed in the face of epidemic alarm and are more likely to take risky decisions. In this paper, we assume that the populations of concerned and risky agents remain constant during the dynamics, i.e. neglecting possible changes in risk perception during the course of the epidemic dynamics. In particular, we fix a parameter f , the fraction of concerned individuals in the population so that the number of concerned and risky agents are N C = fN and N R = (1 − f )N. To perform the model analysis, we will translate the above basic evolution rules into discrete-time Markovian equations [37] . For each individual there are four accessible states (compartments): S p , S np , I p and I np . Therefore, the state at time t of an individual i (i = 1, . . . , N) belonging to group α (α ∈ {C, R}) is defined by the probabilities X i,α y (t), with X ∈ {S, I} and y ∈ {p, np} that agent i is in one of the former four states. Obviously, for a given agent i belonging to group α these probabilities satisfy the normalization condition: At each time step, each agent determines her state according to the values of the probabilities {X i,α y (t)} that are updated according to the possible transitions among the four compartments as plotted in figure 2. In particular, the Markovian equations governing the evolution of the four probabilities associated to a given individual i belonging to group α are Figure 2 . Transition probability trees for the states X i,α y (t) at each time step. The root shows initial state at t and the leaves represent the states at the next time t + 1. The first-row arrows show the P → NP and NP → P processes with probabilities Γ i,α p→np and Γ i,α np→p . The second-row arrows denote probabilities μ, q i,α np or q i,α p governing the changes in the epidemiological states of each node. (Online version in colour.) and indicates the probability that agent i of type α, being susceptible and protected (non-protected) at time t, is infected at this time step. Also, Γ α p→np (t) (Γ α np→p (t)) is the probability that an agent that is protected (non-protected) at time t decides to update the strategy to non-protected (protected) at this time step. Probabilities q i,α p (t), q i,α np (t), Γ α p→np (t) and Γ α np→p (t) are derived below. To derive the probabilities q i,α p (t) and q i,α np (t), let us interpret the interaction network as a twolayer graph, so that concerned and risky nodes belong to each of the two layers, respectively. This bilayer network is composed of N C = fN and N R = (1 − f )N nodes in each layer and contains interaction graph: intra-layer (R − R and C − C) and inter-layer (R − C) links as shown in figure 3 . Obviously, the number of links of each type is constant, since we assume that the fraction of concerned (f ) and risky (1 − f ) agents is constant. In this bilayer structure, a susceptible node in layer α can become infected through one of its infected neighbours in the same layer or through an inter-layer link connecting to an infected agent. Following the epidemic dynamics S I described above, a not-protected susceptible agent, i.e. having a state S α np (t), can contract the disease from each of those I α np (t) neighbours in both layers with probability λ, while infections coming from protected infected neighbours, I α p (t), in both layers occur with probability γ λ. However, a protected susceptible agent, S α p (t), will be infected with probability γ λ, regardless of whether the contact occurs through an intra-layer or an inter-layer link with agents in states I α p (t) and I α np (t). Consequently, we can write the following equations: and where α, β ∈ {C, R}, and A αβ denotes an N × N adjacency matrix whose element (i, j) is equal to 1 if there is a connection between agent i of type α and agent j of type β, while A αβ ij = 0 otherwise. Note that the overall adjacency matrix of the network A can be written as A = A CC + A RR + A CR . Thus, the former equations can be read as 1 minus the probability of not being infected by any infectious agent. In both expressions, this latter probability is split in two terms: the first product of brackets gives the probability of not being infected via intra-layer (R − R or C − C) links, while the second product indicates the probability of not contracting the virus from an agent of a different type, i.e. through R − C links. To round off the derivation of the Markovian equations, we now show the expression for the probabilities associated to the strategic updates, Γ α p→np (t) and Γ α np→p (t), for each type of agent α ∈ {C, R}. To do so, we first assign the payoff associated with each strategy at a given time t for an individual of class α: and are the expected (average) fraction of individuals in S p , S np , I p and I np compartments, respectively. As introduced above parameters c α and T α denote, respectively, the costs associated with protection and infection for an agent of class α. According to the behavioural responses of the concerned and risky agents, we assume that c R > c C and T C > T R . In particular, we consider the relation between parameters c and T for the concerned and risky people as follows: and where δ ∈ [0, 1] is a parameter that accounts for the awareness gap between concerned and risky agents. Finally, once agents have estimated the payoff associated with each strategy they can choose between acting as protected or not-protected in the next time step. In particular, for an agent of type α, Γ α p→np (t) and Γ α np→p (t) denote the update probabilities for the strategic choice P NP. Here we take the usual discrete-time and finite-population analogue of the replicator evolutionary rule [38] [39] [40] for the form of these update probabilities: and where Θ(x) is the Heaviside function: Once we have introduced the dynamical rules of the model and the associated Markovian equations, we proceed to analyse the behaviour arising from the interplay between the SIS spreading dynamics, the evolutionary dynamics for individual decisions and the partition of the population between concerned and risky agents. We do so by iterating numerically equations (2.2)-(2.3) from a given initial condition to construct the trajectory of the probabilities associated with each of the agents until a (static or dynamical) equilibrium is reached. All the results presented here consider a simplified network of social interactions represented by an Erdös-Rényi (ER) graph with size N = 2000 and mean degree k = 10. Although this is a very naive approximation to the structure of real social networks, it represents a starting point to understand the physical properties of the critical system. Our main interest is to study how contagion and the use of protective measures co-evolve in the groups of concerned and risky individuals so as to gain insights about the mutual influence between these antagonistic behaviours in a given population. To this end, we should define the following observables of interest: the relative fraction of risky agents that are infected (I R ), the relative fraction of risky agents that are protected (P R ), the relative fraction of concerned agents that are infected (I C ) and the relative fraction of concerned agents that are protected (P C ). Starting from the microscopic probabilities X i,α y (t) associated with each individual it is easy to compute the average (expected) occupation at a given time t of the four possible states accessible to each group α as: X α The corresponding four average quantities fulfil, at any time, the following conservation laws: Thus, it is possible to define the relative fractions of interest to our study as Let us first consider the case f = 1/2, so that the number of concerned and risky individuals in the population is equal (N C = N R = N/2), and analyse the behaviour of I α and P α as a function of the contagion probability, λ, and the failure probability of prophylaxis, γ . In figure 4 , we show the phase diagram for these quantities in the parameter plane (γ , λ) and fix the values of the remaining parameters as reported in the caption of figure 4 . For the concerned group, different regimes similar to those found in [33] are obtained. Namely: -Healthy state: wherein all concerned agents are susceptible but not protected. This happens for λ < λ c , i.e. for contagion probabilities below the epidemic threshold. In its turn, the behaviour of risky agents is totally different when, as in figure 4 , the value of δ is very small, so that the awareness gap between concerned and risky people becomes very large. In this case, the phase diagram for the fraction of protected risky individuals, P R (γ , λ), shows that risky agents refrain from protecting themselves regardless of the value of λ and γ . Even for high effective protection and low infection probabilities (γ 1 and λ λ c ), the chosen strategy is always NP. Thus, the phase diagrams correspond to the usual SIS with only two regimes: the disease-free (λ < λ c ) and the epidemic (λ > λ c ) ones. An interesting question is whether risky agents change their strategy by adopting protection for certain values of parameters. To unveil that, below we study the role of T C (or equivalently T R ), and δ, while keeping c R = 1 without loss of generality. In both cases, the reported curves P R (T C ) show a clear phase transition so that above a threshold value, T C c , risky agents start to take protection measures. The precise value of the cost threshold, T C c , depends on the values of γ and λ, increasing when the effectiveness of protection decreases and the contagion probability increases. Up to now, we have considered that the awareness gap between risky and concerned agents is extremely large, at δ = 0.01, thus being very far from an homogeneous response to epidemic risk in the entire population. As δ grows the awareness gap decreases, recovering the homogenous response limit studied in [33] for δ = 1. To monitor the transition between heterogeneous and homogeneous response to risk, we have explored the evolution of P R as a function of δ. Figure 6a , b shows the bifurcation diagram for P R as a function of δ for different values of the probability of protection failure γ and the infection probability λ. Again, a phase transition for the adoption of protection by risky agents shows up when δ equals some threshold value δ c . Interestingly, the value of δ c needed for risky agents to take protective measures notably increases as the measure becomes more inefficient, but it slightly changes with the probability of contagion λ. When protection is highly effective, i.e. for small values of γ , the increase of δ beyond δ c drives the system towards a supercritical Hopf bifurcation point, denoted as δ * c , in which the fixed point for the protection levels of risky (and also that of concerned) individuals loses its stability while a limit cycle in which the fraction P R (and P C ) oscillates in time in a sustained way. Solid lines in figure 6 show stable fixed-points for the range of 0 < δ < δ * c , while dashed-dotted and dashed lines identify lower and upper turning points of the stable limit cycles, respectively. From these diagrams, we observe that the increase of γ and λ shifts the Hopf bifurcation to larger values of δ, so that for large enough values of γ and λ, the oscillatory behaviour vanishes. Figure 6c , d shows the behaviour of both concerned and risky groups before (δ c < δ < δ * c ) and after (δ > δ * c ) the Hopf bifurcation. From these panels, we also observe that the amplitude of oscillations is smaller for concerned people than for the risky group. This implies that strategic changes are more likely for risky agents. For δ close enough to 1, behavioural responses of both groups are identical, as is expected (figure not shown). Our previous results have shown that when concerned and risky players equally populate the network, risky agents take advantage of the protective effort of concerned ones when the cost associated with the disease does not exceed a threshold, T C c or when δ is less than δ c . At this point, one would ask whether increasing the fraction of risky agents can lead to a change in their strategy, in particular of those threshold values corresponding to their onset of protection. To this end, we have studied how contagion and protection patterns change when the fraction f of concerned agents varies. In figure 7 , we show the phase diagrams for P R , I R , P C and I C as a function of f and T C . While for concerned agents protection is fully adopted, risky agents show their resistance to protect themselves unless T C is large enough, we consider the limit when δ = 0. In this limit, both risky and concerned agents act like zealots: on one hand, for risky agents the ratio between the cost of protection and the cost of the disease becomes infinite (since they neglect the cost of contracting the diseases T R = 0) and thus never accept protection, while for concerned agents the former ratio is 0 (since the cost associated with protection is c C = 0). This behaviour becomes evident by evaluating equations (2.8) and (2.9) for risky agents, for whom the pay-off associated to the NP strategy is always larger than that of the P strategy. Hence, no matter how severe the incidence of the disease, their decision is never changed. For the concerned agents c C vanishes in this limit, and only I p (t)/(I p (t) + S p (t)) and I np (t)/(I np (t) + S np (t)) fractions affect the protection dilemma for this group. We are particularly interested in the case that the two groups are complete zealots and thus the awareness gap is the largest possible one. With this aim, we assume the most effective protection, γ = 0. In this special case, all concerned agents are protected and healthy (P C = 1 and I C = 0) and the disease spreads only among risky agents for which P R = 0. The question thus is to what extent the protection provided by concerned agents extends to the risky population. This question is answered by the phase diagram for the fraction of risky agents that are infected (I R np ) as a function of f and λ, shown in figure 9 . Apart from the usual epidemic threshold λ c for small values of f the diagram shows a clear transition from healthy to endemic phase when f is varied. Again, we find a threshold value f c so that for f > f c risky individuals remain healthy for any value of λ (although not taking any protection) in a similar way to the results shown in figures 7 and 8 when, for f > f c , risky agents refrain from protecting themselves while being mostly susceptible sheltered by a kind of herd immunity provided by the protection taken by the concerned agents. Let us now derive the expression of f c in the case of a system composed of concerned and risky zealots. Since risky agents never change their NP strategy, we conclude Γ R np→p = 0 and Γ R p→np = 1. Inserting these conditions in equation (2.5), we obtain: where, according to equation (2.7) the infection probability for risky agents is given by In order to find the herd immunity threshold f c , let us assume that our system has reached an equilibrium state and that, in this equilibrium, the probability that a risky (not protected) agent is infected is small: I i,R np (t + 1) = I i,R np (t) ≡ R i . In this limit, we can write equation (3.5) as To satisfy this equation, the following equality must hold: where Λ max (A RR ) is the maximum eigenvalue of matrix A RR . Now let us recall that A RR is a N × N matrix that contains only those links connecting individuals of type R. Thus, considering that if two nodes i and j are connected then A ij = 1 (where A is the adjacency matrix of the network) and that the probability that these agents are both of type R is (1 − f ) 2 , we can approximate the elements of matrix A RR as A RR ij = (1 − f ) 2 A ij , so that A RR can be read as the probability that given a network a link connected two R agents. Thus, Λ max (A RR ) = (1 − f ) 2 Λ max (A) and equation (3.8) can be written as Now, by noting that Λ max (A) k 2 / k [41, 42] and, considering that for the ER networks we can approximate k 2 ∼ k 2 , we obtain the following expression for the epidemic threshold: The former equation defines a curve λ c (f ) (solid curve in figure 9 ) that pinpoints the border between the disease-free regime and the epidemic one. Finally, we consider the maximum value of the epidemic threshold λ c = 1, and find the critical fraction f c of concerned population that gives herd-immunity to the risky one as: ( Given the sharp behaviour of the curve λ c (f ) near f c a coarse-grained description of the phase diagram is as follows. For f < f c , the epidemic threshold of λ c is a function of the recovery probability and the fraction of the concerned agents. However, at critical point f c , the epidemic threshold jumps to 1. In other words, if the fraction of concerned people in the society exceeds a certain value, the disease-free phase is reachable for risky agents regardless of how large the probability of contagion is. The analytical curves are compatible with the numerical results obtained in the previous section, thus highlighting the generality of the herd-immunity threshold f c . In this work, we have studied the emergence of behavioural responses under the threat of the spread of a disease, how the adoption of protective measures impacts the course of epidemics, and how the interplay of different risk perceptions to the same disease shapes the collective response of the population. In this framework, people decide to protect or not depending on the perceived risk of infection and the costs associated with the protection measures and the disease. The main novelty of this work is the partition of the population into two groups with different social attitudes: concerned (C) and risky (R). In particular, both the cost of taking the protective measures, c, and that of contracting the disease, T, are assumed to be different for concerned and risky agents so that c C < c R and T C > T R . Within this framework, and considering a random distribution of concerned and risky agents in homogeneous Erdös-Rényi graphs, we have focused on the identification of the different equilibria that appear as a result of the coexistence and interplay of the two populations. We have shown that when the populations of concerned and risky players are identical, risky agents can capitalize on the protective effort of concerned ones, i.e. remain susceptible without being protected, provided the cost associated is less than a threshold or when the awareness gap is large enough. When the populations are not of equal size a herd-immunity effect over risky agents can be attained provided the fraction of concerned agents becomes larger than a threshold value. Other features such as the dynamical characterization or the coexistence equilibria or the effects of different mixing patterns between concerned and risky populations will be addressed elsewhere [43] . Our results, within the previous limitations, shed light on the difficulty of countering the existence of denialist minorities since they may find a false sense of security in the herd immunity created by the concerned majority. The danger of this effect is the potential spread of risky behaviour. This propagation has not been considered here, since the populations of the two groups were considered static, and is left for the future work. Data accessibility. All data used in this study are publicly available. Authors' contributions. J.G.G. and A.A. conceived the study. M.K. and N.A.-T. performed the numerical analysis and wrote the original draft. All authors edited the manuscript and approved its final version. All authors gave final approval for publication and agreed to be held accountable for the work performed therein. Infectious diseases of humans The mathematics of infectious diseases Modeling infectious diseases in humans and animals Spread of epidemic disease on networks Suppressing epidemics with a limited amount of immunization units Vaccination and epidemics in networked populations-an introduction Efficient immunization strategies for computer networks and populations Optimal vaccination program for two infectious diseases with cross immunity Imperfect vaccines and the evolution of pathogen virulence A simple vaccination model with multiple endemic states Vaccination intervention on epidemic dynamics in networks Prevention of infectious diseases by public vaccination and individual protection Imperfect vaccine and hysteresis 2020 Competition between vaccination and disease spreading Combined effects of prevention and quarantine on a breakout in SIR model Endogenous social distancing and its underappreciated impact on the epidemic curve Community lockdowns in social networks hardly mitigate epidemic spreading Behavioural change models for infectious disease transmission: a systematic review (2010-2015) Group interest versus self-interest in smallpox vaccination policy Group interest versus self-interest in smallpox vaccination policy Statistical physics of vaccination Game theoretic modelling of infectious disease dynamics and intervention methods: a review Evolutionary economics: analysis of traffic flow and epidemics The impact of imitation on vaccination behavior in social contact networks Mean-field analysis of an inductive reasoning game: application to influenza vaccination Hub nodes inhibit the outbreak of epidemic under voluntary vaccination Social contact networks and disease eradicability under voluntary vaccination Game theory of social distancing in response to an epidemic Policy resistance undermines superspreader vaccination strategies for influenza Evolutionary vaccination dilemma in complex networks Interplay between cost and benefits triggers nontrivial vaccination uptake An epidemiological model with voluntary quarantine strategies governed by evolutionary game dynamics Pulsating campaigns of human prophylaxis driven by risk perception palliate oscillations of direct contact transmitted diseases Using social and behavioural science to support COVID-19 pandemic response Coevolution of pathogens and cultural practices: a new look at behavioral heterogeneity in epidemics. Theor 2020 Co-evolution spreading of multiple information and epidemics on two-layered networks under the influence of mass media Nonperturbative heterogeneous mean-field approach to epidemic spreading in complex networks Spatial structure often inhibits the evolution of cooperation in the snowdrift game Evolutionary game theory: temporal and spatial effects beyond replicator dynamics Co-evolution of strategies and update rules in the prisoner's dilemma game on complex networks Spectra of random graphs with given expected degrees Approximating the largest eigenvalue of network adjacency matrices 2021Behavioural response to heterogeneous severity of COVID-19 explains temporal variation of cases among different age groups Competing interests. We declare we have no competing interests. Acknowledgements. We acknowledge useful discussions with A. Reyna-Lara, D. Soriano-Paños and B. Steinegger.