Abstract
The spread of disease can be slowed by certain aspects of real-world social networks, such as clustering1,2 and community structure3, and of human behaviour, including social distancing4 and increased hygiene5, many of which have already been studied. Here, we consider a model in which individuals with essential societal roles—be they teachers, first responders or health-care workers—fall ill, and are replaced with healthy individuals. We refer to this process as relational exchange, and incorporate it into a dynamic network model to demonstrate that replacing individuals can accelerate disease transmission. We find that the effects of this process are trivial in the context of a standard mass-action model, but dramatic when considering network structure, featuring accelerating spread, discontinuous transitions and hysteresis loops. This result highlights the inability of mass-action models to account for many behavioural processes. Using empirical data, we find that this mechanism parsimoniously explains observed patterns across 17 influenza outbreaks in the USA at a national level, 25 years of influenza data at the state level, and 19 years of dengue virus data from Puerto Rico. We anticipate that our findings will advance the emerging field of disease forecasting and better inform public health decision making during outbreaks.
Similar content being viewed by others
Main
Consider the school teacher who is infected with influenza by a student. At some point, they may stay home from work due to the illness and a replacement instructor will fill their role. For the ill teacher, this is probably a benefit, in that they have time to recover; however, the replacement teacher is now in a social situation where infection may be much more likely. A similar narrative relates both to other essential members of society and to common situations such as asking a co-worker to attend a meeting because you have fallen ill. There is overwhelming evidence that this behavioural process occurs consistently and regularly during outbreaks. For example, numerous studies report that ∼4% of teachers are absent and have to be replaced with a substitute instructor on any given day during the influenza season in the United States of America6,7,8,9. Furthermore, ref. 7 found that teacher absenteeism peaks on the peak week of the influenza season and that vaccinating students can lower the rate of teacher absenteeism by up to 10%. For nurses tracked during the 2009 H1N1 pandemic, 24% had influenza-like-illness and 74% of these individuals stayed home from work and had to be replaced10. Clearly such relational exchange happens; nevertheless, it remains to be quantified how this behavioural process affects transmission and whether evidence exists for relational exchange in real-world outbreaks.
In its most basic form, relational exchange is defined as a node replacement process where some individuals (for example, teachers, custodians, health workers, and even children on a hockey team) will be replaced by susceptible individuals if they are ever infected. This replacement process occurs at some rate, termed γ in our equations, to account for a potential delay between when individuals become infectious and when they are diagnosed. Once replaced, a new susceptible individual is given some of the connections of the essential individual (such as students or patients). This relational exchange is important because: the new susceptible node is introduced into what is most probably a more dangerous situation with respect to disease risk; and bringing susceptible nodes from a different region of the contact network reduces the diameter of the population.
To begin, we investigated a standard mass-action model where nodes are distinguished by their state in a susceptible–infectious–recovered (SIR) model. If we do not explicitly tag nodes as being essential or non-essential, we can assume that each infectious node is replaced by a susceptible node at an effective rate γ. This simple model is described in the Methods, and always leads to smaller epidemic peaks and final epidemic sizes than the equivalent model with no relational exchange (that is, γ = 0). The logic behind this result is straightforward, essentially the system behaves like a hybrid between SIR and susceptible–infectious–susceptible (SIS) models, with an effective recovery rate equal to the sum of the recovery and relational exchange rates.
What becomes apparent is that the critical feature of relational exchange, namely that the replacement individual is put into a more dangerous situation than they were in before, is simply not captured by mass-action models. Replacement individuals are thus not equivalent to a random susceptible individual. To properly account for this effect, we introduce network structure into the population. In the network model, individuals in contact with infectious nodes will see their links rewired (the relation is exchanged) from an infectious to a susceptible node. We test a more realistic rewiring scheme through simulations in the Supplementary Information. In both the analytical analysis and the simulation models, the effect of relational exchange is simple: for every exchange the degree of the infectious node (for example, a nurse) is reduced by one, the degree of the node on the other end (for example, a patient) remains the same, and the degree of the replacement susceptible node is increased by one. Again, the events occur at a rate γ for each link with infectious nodes, and the dynamics otherwise follow the standard epidemiological model. To solve this model analytically, we can use pair approximations to describe the dynamics of both the disease and the network. To further facilitate analytic treatment, we consider the SIS model in the main text; however, qualitatively equivalent results apply to SIR models, see Fig. 1b and the Supplementary Information. The pair approximation approach consists of following the fraction of individuals that are susceptible [S] or infectious [I]
![](http://media.springernature.com/lw306/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ1_HTML.gif)
![](http://media.springernature.com/lw306/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ2_HTML.gif)
where r is the recovery rate, β is the transmission rate, and γ is the replacement rate; as well as the fraction of pairs (that is, links)
![](http://media.springernature.com/lw351/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ3_HTML.gif)
![](http://media.springernature.com/lw424/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ4_HTML.gif)
![](http://media.springernature.com/lw371/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ5_HTML.gif)
a, Expected final epidemic size, which, unlike the classic model, can undergo either a continuous or discontinuous transition at the epidemic threshold. b, Time evolution at β/r = 0.175 for the SIS (solid lines) and SIR (dashed lines) models for both the classic model (green) and the acceleration caused by relational exchange (blue). The solid black line represents exponential growth for the SIS model. c, Illustration of the hysteresis loop that can be encountered when varying the rates at which infectious individuals are replaced.
For relational exchange, we find both that the final epidemic size and rate of spread just before the peak can be higher, see Fig. 1. Put simply, relational exchange can parsimoniously account for faster transmission near an outbreak’s peak than would be predicted given early data on transmission. We should stress that replacement does not occur only near the peak. In the model, replacement occurs throughout the time evolution of the disease, yet we find an interesting trade-off between its early effect and its effect near the peak. This trade-off arises from the state of nodes involved in the exchanged relation; for instance, whether or not the students of the original teacher are themselves infectious. Thus, accelerating transmission is an emergent feature of relational exchange and is not due to a change in model rules. As we demonstrate with a series of agent-based models in the Supplementary Methods, accelerating transmission also occurs in considerably more complex models, including those with complex social structure, heterogeneity in contact between age groups, and empirically derived social networks.
With the form used in equations (1)–(5), our model is fully solvable. We use conservation equations for nodes, [S] + [I] = 1, and links, [SS] + [SI] + [II] = 〈k〉/2, where 〈k〉 is the average degree. This leads to a simplified system whose solution is presented in the Methods. We find that, unlike the standard SIS model on a static network, whose endemic state can only emerge via a continuous transcritical bifurcation, the endemic state in the relational exchange model can also appear through a discontinuous saddle-node bifurcation. As shown in the Methods, the transcritical bifurcation occurs when
![](http://media.springernature.com/lw342/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ6_HTML.gif)
whereas the saddle-node bifurcation, which may imply the presence of a hysteresis loop, occurs at
![](http://media.springernature.com/lw387/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ7_HTML.gif)
This analysis helps us highlight a second important feature of the model. Namely that, although the rules of the process are local, they model a global policy: the replacement rate is a function of how closely we are surveying the state of essential nodes and of how quickly we wish to replace them. However, a fast replacement rate is effective only at low prevalence, such that secondary infections because of relational exchange are rare. Now, let us assume that the replacement rate is high enough to keep an outbreak under control, but that after some time the rate is slightly reduced. This might occur, for example, after the initial fear wears off. Such a mechanism can push the system over a discontinuous transition, such that a microscopic change in γ can lead to a macroscopic change in disease prevalence. We then wish to bring the system back to its previous state, so we increase γ back to its previous value. Unfortunately, as shown in Fig. 1c, the system exhibits a hysteresis loop. The replacement rate must be increased well beyond its previous value for the system to return to the initial state.
We now turn our attention to empirical data. The question we wish to answer is whether evidence exists for relational exchange in real-world diseases. We select as our metric for relational exchange the presence of accelerating exponential transmission near the outbreak peak. This phenomenon is a ubiquitous feature of relational exchange, but very uncommon in all other general models of disease spread published to date. For this investigation, we selected two pathogens: influenza and dengue. The rationale for selecting these diseases is that they both have strong seasonality—which can drive accelerating spread in some models11,12—and have rich historical data sets. However, because dengue is a vectored pathogen, which induces a separation of timescales between transmission and behaviour, and influenza is not, we expect the importance of relational exchange to be far greater for influenza13,14,15. All models were fitted to individual seasons because both influenza and dengue lead to sterilizing immunity during a single outbreak, which is true even if there are multiple circulating strains of dengue16,17. To test for the presence of relational exchange, we fit exponential growth models to the first eight weeks of each season and then compared the predicted number of cases to the model prediction from week nine until the peak (see Methods).
In short, we find the predicted pattern: evidence for relational exchange for influenza, but not for dengue. Figure 2 shows sixteen years of national-level influenza data from the United States of America’s Influenza-Like-Illness surveillance network. In all sixteen years, and during the 2009 pandemic, we find strong evidence that the exponential rate of transmission accelerates before the peak. For dengue, with nineteen years of data, we see only a few scattered instances of accelerating transmission–aside from the 1995 outbreak, which was the most dramatic outbreak observed since 1990.
a, We examined seventeen outbreaks, sixteen years and the 2009 pandemic, of national Influenza-Like-Illness (ILI) data in the USA (grey dashed lines) by fitting exponential growth models to the first eight weeks of the influenza season, defined by the USA Centers for Disease Control and Prevention (CDC) as weeks 40–20, and then comparing the observed number of cases to the exponential prediction up until the peak. For these data, ILI was defined as a fever ≥100 °F (approx. 37.778 °C) and one additional upper-respiratory symptom, for example, sore throat, during the CDC-defined influenza season without another known non-influenza cause. If more cases are reported than predicted by this exponential fit, we considered this evidence for accelerating spread and relational exchange (solid red lines). We find evidence for relational exchange in all sixteen years and during the fall wave of the 2009 pandemic. b, Example of the fitting process for a, with training data in blue and the exponential fit in solid black. c, Same analysis as in a, but for the dengue virus using nineteen years of data reported to the USA CDC in Puerto Rico, in which we found only a few isolated instances of accelerating spread. d, Example of the fitting process for c, with training data in blue and the exponential fit in solid black.
We chose to contrast influenza and dengue precisely because both viruses are affected by seasonal forcing. The observation that influenza shows evidence for accelerating spread, and dengue largely does not, indicates that seasonal forcing alone cannot account for the observed patterns. This observed difference between dengue and influenza can parsimoniously be accounted for by relational exchange. Furthermore, existing models of influenza, even those with seasonal forcing, fail to capture the accelerating growth seen in real data13,18. For example, both papers reported evidence that absolute humidity affected the start of influenza season, more so than accelerating the rate of spread leading up to the peak, and consistently under-predicted the peak magnitude and predicted a later peak than observed. The similarity to these failed predictions and the comparison of the standard SIR model to the SIR model with relational exchange in our Fig. 2b is striking. To date, no alternative model can explain this difference, without substantially increased complexity or unrealistic assumptions. Therefore, we conclude that evidence exists for relational exchange in real-world outbreaks.
Because the Influenza-Like-Illness surveillance network is at the national level and the dengue data were the state level, we also fit the relational exchange model to state-level influenza data. Our analysis suggests that across 25 influenza seasons, and in all of the USA (we were unable to obtain data from Vermont and Alaska/Hawaii were not states during the period with historical data), greater than 70% of each season in each state (that is, >900 seasons) support relational exchange. The heatmap of the USA in Fig. 3 plots the proportion of influenza seasons between 1921 and 1951, showing evidence for relational exchange, with darker blues indicating more support than lighter blues and yellows. The data used for this analysis were not obtained from the Influenza-Like-Illness surveillance network, and were instead obtained from the USA National Notifiable Diseases Surveillance System as digitized by Project Tycho19. As a result, we now have more reliable support derived from multiple data sources and at multiple geopolitical scales, indicating that our results on relational exchange are robust to sampling errors and/or synchronization effects. We furthered evaluated the support for relational exchange by fitting models with and without relational exchange to the state-level influenza data and found strong support for the presence of relational exchange in all seasons and statistical support for the relational exchange model across more than 400 influenza outbreaks, see Supplementary Methods.
A heatmap of the USA illustrating the proportion of influenza seasons between 1921 and 1951 showing evidence for relational exchange, with darker blues indicating more support than lighter blues and yellows. On average, over 70% of seasons across all states provide support for relational exchange. The data were obtained from the USA National Notifiable Diseases Surveillance System as digitized by Project Tycho19.
It is worth highlighting that we restricted our analysis to the USA Centers for Disease Control and Prevention influenza season and that we do find evidence for exponential growth during the first eight weeks of every influenza season, which was our training period. However, we find an even higher rate of exponential growth in the weeks leading up to the peak. Lastly, and perhaps most importantly, we see evidence for accelerating spread during the 2009 H1N1 pandemic. This pandemic occurred in the USA outside of the typical influenza season—thus probably altering any effect of seasonal forcing—and generated substantial changes in health-seeking behaviour. The fact that we still see evidence for accelerating spread during the 2009 H1N1 pandemic suggests that the mechanism must transcend seasonality and behaviours associated with health-care seeking. Relational exchange parsimoniously accounts for both seasonal and pandemic influenza. It is also worth noting that dengue and influenza do not have the same functional type of seasonal forcing. However, because dengue is also seasonally forced, this difference between the two diseases demonstrates that accelerating spread can occur with, but is not an expected consequence of, seasonal forcing. Accelerating spread is an expected consequence for relational exchange. Seasonal forcing and behaviour should thus be included and perhaps coupled for more precise predictions.
This work is not without caveats. First, in our analytical treatment, we assumed the population was infinitely large. While this is a standard assumption, clearly the model cannot be directly applied to real-world populations. Second, relational exchange modifies the degree of nodes; however, the degree distribution will never be broader than an exponential distribution, such that using a heterogeneous mean-field approach would be marginally more precise, without qualitatively changing the observed phenomena20. Third, the compartmental approach distinguishes nodes only on the basis of their state, such that while in the real world an individual might return to their past connections after recovering, our model cannot capture this behaviour, because nodes are undistinguishable. Similarly, while random rewiring increases the accuracy of compartmental models, it could also dampen the effects of relational exchange in epidemics occurring over longer timescales. Fourth, although we demonstrate that accelerating spread still occurs in age-structured populations, see Supplementary Methods, we did not use such a model to analyse the influenza and dengue time series. Future work should focus on evaluating age-structured relational exchange models with empirical data. Finally, we considered only two diseases in our empirical investigation. Future work should consider classifying a variety of different diseases, on the basis of whether they exhibit evidence for relational exchange.
There are four additional implications of this study, which will affect the broader scientific, medical, and public health communities. First, as recently demonstrated, the evolution of increased pathogen transmissibility decreases as a pathogen spreads through a heterogeneous population21. This occurs because highly connected individuals are typically infected early in an outbreak and, as these individuals recover, the epidemic potential of the population decreases22. Therefore, any process, such as relational exchange, which maintains highly connected, susceptible individuals can increase the chance of more transmissible strains evolving and persisting. Second, if individuals entering high-risk societal roles can be vaccinated or selected from resistant individuals, then relational exchange will no longer have negative effects on outbreak progression. Third, one clear prediction from this model is that replacement workers, for example, substitute teachers, may have higher rates of illness. Future work should focus on empirical studies to evaluate this prediction. Finally, methods for forecasting disease spread must include both realistic population structure and salient aspects of human behaviour. Without these key features, we cannot hope for robust, actionable models for predicting epidemics.
Methods
Mass-action compartmental approach.
We use a standard mass-action model where nodes are distinguished by their state in a SIR model (variables in brackets correspond to the fraction of the population in each state). The epidemic is described by the following system of equations,
![](http://media.springernature.com/lw298/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ8_HTML.gif)
![](http://media.springernature.com/lw315/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ9_HTML.gif)
![](http://media.springernature.com/lw268/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ10_HTML.gif)
where r is the recovery rate, β is the transmission rate, and γ is the replacement rate. All variables in brackets are dynamical, and we use the dot notation to identify time derivatives.
Network model with pair approximations.
We consider the network structure by introducing pair approximations in the SIR ordinary differential equation system
![](http://media.springernature.com/lw278/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ11_HTML.gif)
![](http://media.springernature.com/lw292/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ12_HTML.gif)
with an implicit via the conservation condition
. The pairs are followed by
![](http://media.springernature.com/lw348/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ13_HTML.gif)
![](http://media.springernature.com/lw438/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ14_HTML.gif)
![](http://media.springernature.com/lw358/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ15_HTML.gif)
where all missing terms go to links involving recovered nodes (which do not need to be explicitly followed). To allow analytical treatment, we consider SIS dynamics as per equations (1)–(5) for the remainder of the analysis.
Analytical solution of the network model.
The evolution of the dynamical variables in the model is constrained by the conservation of nodes and links, which allows us to effectively reduce equations (1)–(5) to the following system of three equations:
![](http://media.springernature.com/lw292/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ16_HTML.gif)
![](http://media.springernature.com/lw337/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ17_HTML.gif)
![](http://media.springernature.com/lw358/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ18_HTML.gif)
We can readily see that a disease-free state, in which every node is susceptible,
![](http://media.springernature.com/lw414/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ19_HTML.gif)
is a steady state of the network model. To see whether there exists an endemic state in which a non-zero fraction of individuals are infectious, we set in equations (16)–(18). This yields
![](http://media.springernature.com/lw278/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ20_HTML.gif)
![](http://media.springernature.com/lw311/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ21_HTML.gif)
![](http://media.springernature.com/lw347/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ22_HTML.gif)
where we have defined the dimensionless parameters a ≡ β/r and b ≡ γ/r. Substituting these last equations in the link conservation equation, we obtain
![](http://media.springernature.com/lw432/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ23_HTML.gif)
whose solutions,
![](http://media.springernature.com/lw462/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ24_HTML.gif)
correspond to the possible values of the endemic state. The stability of the three possible steady states are obtained via a standard linear stability analysis of equations (16)–(18).
The disease-free steady state undergoes a transcritical bifurcation whenever either [I]+∗ = [I]df∗ = 0 or [I]−∗ = [I]df∗ = 0, which happens when the constant term in equation (23) equals zero.
![](http://media.springernature.com/lw326/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ25_HTML.gif)
Fixing either a or b, we find the following threshold values
![](http://media.springernature.com/lw337/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ26_HTML.gif)
Similarly, we see that the endemic steady state appears through a saddle-node bifurcation that occurs when [I]+∗ = [I]−∗ which, from equation (24), happens when
![](http://media.springernature.com/lw387/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ27_HTML.gif)
Fixing either a or b, we find the following threshold values
![](http://media.springernature.com/lw402/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ28_HTML.gif)
Notice that whenever a〈k〉 = 4, the threshold b−sn diverges, but b+sn equals −1 − 1/〈k〉(〈k〉 − 1) is always negative and therefore can be discarded. From equations (26) and (28), we see that a bistable region appears or disappears when atr = asn or btr = b±sn, which respectively yield the criteria
![](http://media.springernature.com/lw337/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ29_HTML.gif)
for the existence of a bistable region.
Empirical analysis.
For exponential growth, the corresponding equation for the number of infected individuals as a function of elapsed time is:
![](http://media.springernature.com/lw102/springer-static/image/art%3A10.1038%2Fnphys3832/MediaObjects/41567_2016_Article_BFnphys3832_Equ30_HTML.gif)
We estimated the parameters of both models from data using a nonlinear least squares algorithm coded in the R programming language. For both dengue and influenza we used data from the first eight weeks of each season. We further explored the sensitivity of our results to the size of the training data and found them to be robust.
Data availability.
Empirical data obtained from the USA National Notifiable Diseases Surveillance System as digitized by Project Tycho19. All other data that support the plots within this paper and other findings of this study are available from the corresponding author upon reasonable request.
References
Volz, E. M., Miller, J. C., Galvani, A. & Meyers, L. A. Effects of heterogeneous and clustered contact patterns on infectious disease dynamics. PLoS Comput. Biol. 7, e1002042 (2011).
Scarpino, S. V. et al. Epidemiological and viral genomic sequence analysis of the 2014 ebola outbreak reveals clustered transmission. Clin. Infect. Dis. 60, 1079–1082 (2015).
Salathé, M. & Jones, J. H. Dynamics and control of diseases in networks with community structure. PLoS Comput. Biol. 6, e1000736 (2010).
Glass, R. J. et al. Targeted social distancing design for pandemic influenza. Emerg. Infect. Dis. 12, 1671–1681 (2006).
Fewtrell, L. et al. Water, sanitation, and hygiene interventions to reduce diarrhoea in less developed countries: a systematic review and meta-analysis. Lancet Infect. Dis. 5, 42–52 (2005).
Ervasti, J. et al. Association of pupil vandalism, bullying and truancy with teachers’ absence due to illness: a multilevel analysis. J. Sch. Psychol. 50, 347–361 (2012).
Graitcer, S. B. et al. Effects of immunizing school children with 2009 influenza A (H1N1) monovalent vaccine on absenteeism among students and teachers in Maine. Vaccine 30, 4835–4841 (2012).
Alker, H. J., Wang, M. L., Pbert, L., Thorsen, N. & Lemon, S. C. Impact of school staff health on work productivity in secondary schools in Massachusetts. J. Sch. Health 85, 398–404 (2015).
Goldhaber, D. & Hansen, M. Is it just a bad class? Assessing the long-term stability of estimated teacher performance. Economica 80, 589–612 (2012).
Mitchell, R. et al. Impact of the 2009 influenza A (H1N1) pandemic on Canadian health care workers: a survey on vaccination, illness, absenteeism, and personal protective equipment. Am. J. Infect. Control 40, 611–616 (2012).
Colón-González, F. J., Lake, I. R. & Bentham, G. Climate variability and dengue fever in warm and humid Mexico. Am. J. Trop. Med. Hyg. 84, 757–763 (2011).
Shaman, J., Karspeck, A., Yang, W., Tamerius, J. & Lipsitch, M. Real-time influenza forecasts during the 2012–2013 season. Nature Commun. 4, 2837 (2013).
Shaman, J., Pitzer, V. E., Viboud, C., Grenfell, B. T. & Lipsitch, M. Absolute humidity and the seasonal onset of influenza in the continental United States. PLoS Biol. 8, e1000316 (2010).
Johansson, M. A. et al. Local and global effects of climate on dengue transmission in Puerto Rico. PLoS Negl. Trop. Dis. 3, e382 (2009).
Wearing, H. J. & Rohani, P. Ecological and immunological determinants of dengue epidemics. Proc. Natl Acad. Sci. USA 103, 11802–11807 (2006).
Reich, N. G. et al. Interactions between serotypes of dengue highlight epidemiological impact of cross-immunity. J. R. Soc. Interface 10, 20130414 (2013).
Rothman, A. Immunity to dengue virus: a tale of original antigenic sin and tropical cytokine storms. Nature Rev. Immunol. 11, 532–543 (2011).
Shaman, J., Goldstein, E. & Lipsitch, M. Absolute humidity and pandemic versus epidemic influenza. Am. J. Epidemiol. 173, 127–135 (2011).
van Panhuis, W. G. et al. Contagious diseases in the united states from 1888 to the present. N. Engl. J. Med. 369, 2152–2158 (2013).
Marceau, V., Noël, P.-A., Hébert-Dufresne, L., Allard, A. & Dubé, L. J. Adaptive networks: coevolution of disease and topology. Phys. Rev. E 82, 036116 (2010).
Leventhal, G. E., Hill, A. L., Nowak, M. A. & Bonhoeffer, S. Evolution and emergence of infectious diseases in theoretical and real-world networks. Nature Commun. 6, 6101 (2015).
Ferrari, M. J., Bansal, S., Meyers, L. A. & Bjørnstad, O. N. Network frailty and the geometry of herd immunity. Proc. R. Soc. Lond. B 273, 2743–2748 (2006).
Acknowledgements
S.V.S. received funding support from the Santa Fe Institute and the Omidyar Group. A.A. received funding support from the Fonds de recherche du Québec—Nature et technologies. L.H.-D. received funding support from the James S. McDonnell Foundation Postdoctoral Fellowship and the Santa Fe Institute.
Author information
Authors and Affiliations
Contributions
S.V.S. and L.H.-D. conceived the project; L.H.-D. and A.A. performed the simulations and calculations; S.V.S. analysed the empirical data; all authors interpreted the results and produced the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary information
Supplementary information (PDF 118026 kb)
Rights and permissions
About this article
Cite this article
Scarpino, S., Allard, A. & Hébert-Dufresne, L. The effect of a prudent adaptive behaviour on disease transmission. Nature Phys 12, 1042–1046 (2016). https://doi.org/10.1038/nphys3832
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/nphys3832
This article is cited by
-
The unintended consequences of inconsistent closure policies and mobility restrictions during epidemics
BMC Global and Public Health (2023)
-
Exploring the interaction of influenza A subtypes H1N1 and H3N2 based on an evolution-driven transmission model
Nonlinear Dynamics (2022)
-
Effect of Human Behavior on the Evolution of Viral Strains During an Epidemic
Bulletin of Mathematical Biology (2022)
-
Interventions to control nosocomial transmission of SARS-CoV-2: a modelling study
BMC Medicine (2021)
-
A review and agenda for integrated disease models including social and behavioural factors
Nature Human Behaviour (2021)