Spatial early warning signals of social and epidemiological tipping points in a coupled behaviour-disease network

The resurgence of infectious diseases due to vaccine refusal has highlighted the role of interactions between disease dynamics and the spread of vaccine opinion on social networks. Shifts between disease elimination and outbreak regimes often occur through tipping points. It is known that tipping points can be predicted by early warning signals (EWS) based on characteristic dynamics near the critical transition, but the study of EWS in coupled behaviour-disease networks has received little attention. Here, we test several EWS indicators measuring spatial coherence and autocorrelation for their ability to predict a critical transition corresponding to disease outbreaks and vaccine refusal in a multiplex network model. The model couples paediatric infectious disease spread through a contact network to binary opinion dynamics of vaccine opinion on a social network. Through change point detection, we find that mutual information and join count indicators provided the best EWS. We also show the paediatric infectious disease natural history generates a discrepancy between population-level vaccine opinions and vaccine immunity status, such that transitions in the social network may occur before epidemiological transitions. These results suggest that monitoring social media for EWS of paediatric infectious disease outbreaks using these spatial indicators could be successful.

Resurgences of vaccine-preventable diseases severely stress public health systems, interrupt tourism and public services, and disrupt economies through the huge costs of large-scale interventions 1 . These impacts motivate the study of factors that support vaccine uptake. Undervaccination may be attributed to vaccine refusal 2 , the spread of anti-vaccine opinion facilitated by media coverage and its sensationalisation of true adverse vaccine effects 3 , the expectation of adverse effects 4 , misstatement of the cause of illnesses 5 , the spread of other rumours and false information 6 , and the effect of social norms 7 .
These phenomena show how the social diffusion of information is heavily responsible for the trajectory of disease spread through its ability to alter individual behaviour. Much work has modelled opinion dynamics for different applications through the use of voter models 8 and majority opinion models 9 , among other frameworks, and their combination with network structure has revealed much about the occurrence of opinion cascades 10 and forecasting 11 . For instance, models coupling behavioural dynamics and spreading processes change the predicted dynamics of influenza transmission 12 and climate change 13 alike.
Opinion propagation can be represented by information diffusion through social networks 14 . Similarly, infection spread is often conceptualized as spreading through a physical contact network 15 . Many models explore the dynamics of n-layer multiplex networks, where each layer represents a different aspect of the dynamics of a single coupled system [16][17][18][19] . In these cases, the theory of phase transitions in spatially structured systems is important. For instance, epidemic regimes have previously been modelled as the outcome of phase transitions in physical systems and are featured widely in the epidemiology literature [20][21][22][23][24][25][26] . Generally, phase transitions occur when a physical system moves from one state to another after going through some critical point. Second-order transitions occur when a macroscopic variable varies continuously 27 and are called critical transitions 28 .
Systems approaching these critical transitions sometimes display characteristic spatial or temporal behaviours called early warning signals (EWS) that can predict coming epidemic disease outbreaks and other events 29 . EWS are perhaps better defined as statistically significant, recognisable and characteristic behaviours known to precede critical transition in dynamical systems 30,31 . For instance, critical slowing down can precede both first-and second-order transitions 31,32 and is accompanied by the divergence of correlation length in a physical system 33 . Many statistics have been used to study EWS in spatially extended systems; temporal 34,35 and spatial correlation 36,37 have been found to precede transitions in spreading processes. Other measurements have been applied to spin systems, where each site in a lattice may be in one of two possible states, possibly partially dependent on the state of neighbouring sites. The spin model has also been applied to opinion dynamics; a simple voter model with binary opinion dynamics is analogous to a physical spin system, where particles represent agents and spins represent different opinions 38 . Consensus formation can be seen as a second-order phase transition to an ordered state (where all spins are aligned). In this regime, knowledge of the opinion of a single agent predicts the opinion of all other agents in the system 39 . Since the transition in finite networks is smooth 40 , the distance across the network over which the opinions of connected agents are strongly correlated increases smoothly; this is analogous to divergence of the correlation length of a physical system 41 .
Above some critical temperature, disordered systems take the form of a spin glass. In a spatial opinion model, this describes a state where opinions between neighbours are generally uncorrelated 42 . On a static network, this state induces a larger number of edges between dissimilar neighbours as compared to that of consensus regimes. This is related to join count statistics, where the numbers of edges between like neighbours are compared to the number between dislike neighbours as a test of geographical distribution. This is arguably the most natural and well-defined measure for graphs presenting binary data and is used for spatial analysis 43 .
The necessity of disease surveillance and early warning signals for outbreaks has been discussed in multiple contexts, from epidemic mitigation to bioterrorism prevention [44][45][46] . Potential mitigation of unnecessary expense motivates us to find reliable EWS that remain easily computable on large high-resolution data sets. Furthermore, the study of EWS in coupled disease-behaviour multiplex networks has received relatively little attention, suggesting a significant gap in the literature. Our objective is to evaluate and compare the relative merits of the mutual information, Moran's I, Geary's C and join count statistics as EWS of the occurrence of epidemics and changes in aggregate opinion on a coupled disease-behaviour network model. We use three differently parametrised models (V1, V2 and V3) coupling a binary vaccination opinion dynamic to an SIRV epidemic process. The resulting trends in the EWS for model V2 will be explored in the Results and Discussion sections, with V1 and V3 presented in Supplementary Information S5.
The outline of this paper is as follows: the Methods section will present the EWS and their derivations and give the details of the model used. The Results section will analyse the trends in the warning signals and the Discussion section will present a review of the study and any shortcomings of our approach, with further results pertinent to the study presented in the Supplementary Information.

Methods
We assume an acute, self-limiting infection that confers lifelong natural immunity upon recovery, and for which a vaccine is readily available. Similar premises have been used to represent the natural history of many paediatric infectious diseases such as measles 47 . In particular, we assume an SIRV p natural history consisting of four mutually exclusive disease states. Agents are initially susceptible to infection (S). Upon infection the agent enters the infected state ( → S I), which we treat as a combination of both the latent, ill and infectious periods 48 . Upon clearing the infection, agents enter the recovered state of lifelong immunity ( → I R); additionally, susceptible agents may be vaccinated and so enter the vaccinated state ( → S V p ) 49 . We also include injunctive social norms (i.e. peer pressure) as well as a perceived cost of vaccination that captures both economic costs and the fear of perceived adverse vaccine effects 50 . As in some models 36 , we include a noise parameter ξ to account for environmental and demographic stochasticity 51 with the simplifying assumption of perfect vaccination 52 (reversion from the recovered state to the susceptible only through agent death). During simulation, each time step represents a single week. the model. We model a multilayer network where each layer is given an identical undirected Erdös-Rényi random graph with size N and mean node degree Q n . Each agent n can be described by a pair of states; for instance, each agent is assigned the joint state V V ( , ) s p at the start of the simulation with probability α (they are a pro-vaccine vaccinated agent), else they are initialised with joint state N S ( , ) (an anti-vaccine susceptible agent) with probability α − 1 . The social process follows an NV s dynamic (Fig. 1b), representing pro-(V s ) and anti-vaccine (N ) opinion for each agent n. ξ represents the probability of any agent switching opinion randomly in each week and  → N V ( ) n s represents the probability of switching from anti-vaccine opinion to pro-vaccine opinion ( → N V s ) upon interaction with a disagreeing neighbour. We introduce an imitation dynamic by having each agent n compare its opinion with a single randomly chosen social contact (a neighbouring agent on the social layer) each week; n then changes its vaccination opinion only if there is disagreement (the agent and the neighbour have different vaccine opinions). This change of opinion depends on the perceived risk of vaccine adverse effects κ ("vaccine risk") and I n (the number of infected physical neighbours of n) according to the rules  Here, σ represents the strength of an injunctive social norm to maintain the current opinion (referred to as "social norm"), while ⁎ Q n represents the number of neighbours of n with opinion ⁎ and Q n represents the total number of social contacts of n.
The epidemiological dynamics (which we term 'physical dynamics') follow an SIRV p process (Fig. 1a), in which an agent n can progress through each of four disease compartments: S (susceptible), I (infected), R (recovered) and V p (vaccinated). Every week (i.e. in each time step), each susceptible agent n interacts with all its physical neighbours; each effective interaction carries the probability p of infection ( → S I), so that each susceptible agent faces the total probability − − p 1 (1 ) I n of infection in a single week. The duration of the illness is  weeks (with no impact on mortality), after which n gains lifelong natural immunity ( → I R). Alternately, if a susceptible agent n adopts a pro-vaccine opinion, they are immediately vaccinated → ( ) S V p and gain lifelong vaccine-derived immunity. We also assume that only susceptible agents are vaccinated. Thus, individual agents may change their opinion about vaccination multiple times in their life ( → → N V N s ), but once they are vaccinated they may not become unvaccinated ( → S V p ). This in turn creates an asymmetry between disease dynamics and social dynamics that will have implications for the model predictions. We will discuss this in the Results section.
Each agent n has probability µ of dying each week, upon which they are replaced by a new susceptible individual who is a pro-vaccine vaccinated agent V V ( , ) s p with probability α, or an anti-vaccine susceptible agent N S ( , ) with probability α − 1 , keeping the same physical and social contacts as the agent they replaced (that is, the network is static). Case importation is accounted for by infecting a randomly selected proportion ι of susceptible agents at the start of each week, and noise is introduced to the model by changing the vaccine opinions of a randomly selected proportion ξ of the entire population weekly.
At the start of the simulation, some susceptible agent is randomly selected as an index patient and infected; subsequent disease spread is governed solely by environment and inter-agent interaction. Models V1, V2 and V3 all use these model dynamics, and a flowchart and detailed written description of the model transitions are given in Supplementary Information S2. A complete list of the variables used is given in Supplementary Table S1. early warning signals. Mutual information M is defined as respectively upon interaction with a dissenting neighbour. α gives the probability of being birthed with provaccine opinion V s .
x y 2 where X and Y are discrete random variables; x takes value on the set X x x x { , , , } 1 2 3 = … and y on set y y y { , , , } 1 2 3 Y = …, with  a joint probability mass function of X and Y 53 . Mutual information is an entropy-based quantification of the "shared information" of two random variables quantifying how knowledge of one decreases the uncertainty of the other and vice versa 54 . Mutual information peaks at the critical temperature of spin systems during second-order transitions and has been widely used in detecting phase transitions 55,56 ; an advantage of this statistic is its ability to quantify non-linear dependence, unlike Moran's I and covariance which only account for linear dependence.
Join counts quantify the degree of clustering by giving the number of adjacencies between agents of different types. We divide the population into two attributive classes, with V s the compartment of pro-vaccine agents and N the compartment of anti-vaccine agents. Let [Ψ, Ω] be the number of social interactions between agents with vaccine opinions Ψ and Ω; then N N [ , ] represents the number of nearest-neighbour interactions between anti-vaccine agents, V V [ , ] s s the number of interactions between pro-vaccine agents and N V [ , ] s the number of interactions between pro-and anti-vaccine agents. These can be written as if agents j and k are social neighbours, with ω = 0 j k , otherwise (ω is the adjacency matrix of the social network); x n represents the opinion score of agent n, defined as In an opinion model, clustering manifests as agents consistently having a higher number of like-minded neighbours than expected based on the global prevalence of the opinion; join counts are then used to test the null hypothesis of positive correlation 57 . Join counts are used in many fields as a categorical test of spatial autocorrelation, including ecology 58 and geographical information systems 59 . In all parameter realisations here, the number of joins are counted naïvely rather than calculated. Joins between like-minded agents (e.g. N N [ , ] and V V [ , ] s s joins) will be called similar joins, and edges between disagreeing neighbours (e.g. N V [ , ] s ) will be called dissimilar joins.
The Moran's I coefficient J quantifies spatial correlation and is defined represents the mean opinion score of the population 43 . Used as a global statistic, Moran's I gives the degree of correlation between the values of neighbouring patches (agents and their social contacts); here, the numerical value of the vaccine opinion is the same as described in Eq. (5). Algebraic manipulation of Eq. (6) using Eq. (5) gives s s s s s (full derivation given in Supplementary Information S4); we can then consider Moran's I as a measure derived from the linear combination of join counts. Positive values signify spatial correlation, with negative values signifying anticorrelation. The Geary's C coefficient C is yet another measure of spatial correlation based on the cross-product (like Moran's I) 43 , but unlike Moran's I it accounts for the difference in opinion between two neighbours 60 . It is given as Lower values show spatial correlation, and large values represent anticorrelation. Like Moran's I (Eq. (7)), Geary's C can also be expressed as a linear combination of join counts www.nature.com/scientificreports www.nature.com/scientificreports/ proportion ratio ι = . × − 2 5 10 5 was added to provide periodic impulses of infection as a test of resilience in endemic disease regimes. Here, an ensemble of 100 simulations using parameters κ = 0, σ = 0 and α = .
0 05 returned the values < . Supplementary Information S3), where Ψ represents the mean number of agents with (social or physical) state Ψ, averaged over all realisations of that combination of parameter values.
p 0 2 was chosen to reflect the reproduction rate of a measles infection commonly estimated from empirical data 61 ; effective contacts occur in the simulation once per week during the period of infection. The probability of randomly switching opinion ξ = × − 1 10 1 4 was included as a source of noise. We found that the parameter ranges for vaccine risk κ ∈ − [ 1,1] and social norm σ ∈ [0, 3] were sufficiently broad to capture transitions in both social and physical dynamics (Fig. 2a,b), as well as the corresponding trends in the mutual information (Fig. 2c) and dissimilar join count (Fig. 2d).
The contours in each panel of Fig. 2 show the obvious correspondence between transitions in the social and physical dynamics of the model, and substantial changes in M and N N , ; here, the dissimilar join count N V , s (Fig. 2d) and mutual information M (Fig. 2c) increase while the vaccine risk κ and social norm σ parameters increase towards their respective (pre-transition) threshold values. These trends are generally asymmetric about both transitions; this can be seen in Fig. 4, where post-transition trends do not exhibit similarly detectable warnings (if any). This Parametrisation applies to model V2. The corresponding parametrisations and contour plots of models V1 and V3, as well as their post-transition trends are presented through the Supplementary Information.

Results
Due to the low initial vaccine coverage α = .
0 05, all realisations demonstrated an initial epidemic spread (defined in Supplementary Information S3) over the first 7 weeks. After this period, the dynamics settled down to a quasi-equilibrium state characterized by fluctuations around a mean value that is the focus of our study -the following subsections are grouped by major findings of the model. The term model variables refer to the outputs S , I , R , V p , N and V s . Population vaccine immunity status can differ from aggregate vaccine opinion. Because only susceptible individuals are vaccinated and individuals cannot become 'unvaccinated' (but may change their opinion about the vaccine over their lifetime), the population-averaged vaccine opinion is not equal to the population-averaged vaccine immunity status, even at the quasi-equilibrium state. With no social pressure (σ = 0), a small increase in vaccine risk κ → . 0 03125 pushes the system to endemic infection and anti-vaccine consensus (Fig. 3a) despite a high vaccination rate (Fig. 3d). Towards an explanation, if an agent n is newly birthed into this regime, the probability of having an infected neighbour vanishes 〈 〉 →  under the assumptions that the average agent with anti-vaccine opinion is almost certain to interact with a disagreeing contact. A similar calculation explains the phenomenon of high vaccination rate (Fig. 3e) despite mixed consensus (Fig. 3b) when vaccine risk becomes neutral (κ → 0). A vaccine then perceived as beneficial (κ → − . 0 03125) intuitively results in a high vaccination rate (Fig. 3f) and pro-vaccine consensus (Fig. 3c); in both these regimes, the disease survives only through case importation. In the absence of social norms and vaccine risk, the population's aggregate vaccine opinion may not be a good indicator of its vaccine immunity profile (and vice versa); Eq. (2) shows that the probability of changing opinion depends only on I n when social norm σ = 0. In this region, the pattern of disease spread will be determined by the initial conditions of the physical dynamics; slight changes in vaccine risk κ will push the network towards either of the consensuses, with minimal www.nature.com/scientificreports www.nature.com/scientificreports/ effect on the high vaccination rate (Fig. 3d). This phenomenon is shared by models V1 and V3, as shown in Supplementary Information S5.1.

eWS trends identify approaching transitions in both social and physical layers. The trends in
both model dynamics and the proposed EWS are shown in Fig. 4, with social norm σ = 0 (column A) and σ = .
0 25 (column B). The first vertical black line in all panels of Fig. 4 represents the transition in the social dynamics K s , defined as the smallest κ value at which V N s 〈 〉 ≈ 〈 〉 (the mean number of pro-vaccine agents equals the number of anti-vaccine agents); the second vertical black line represents the transition in the physical dynamics K p , similarly defined as the earliest κ where R V p 〈 〉 ≈ 〈 〉. Multiple physical and social transitions were found for some parameter combinations; these trends and the attendant behaviour of the EWS for models V1 and V3 can be seen in Supplementary Information S5.3. We also note that the equilibrium values of EWS and model variables were averaged over − 15 20 realisations of all parameter combinations. All EWS show recognisable trends preceding both transitions for both social norm values σ = 0 (column A) and σ = .
0 25 (column B); for instance N N , (Fig. 4c,d), N V , s (Fig. 4e,f), J (Fig. 4i,j) and C (Fig. 4k,l) increase sharply preceding K s with all but N V , s (Fig. 4e,f) approaching some maximum value preceding K p , while V V , s s (Fig. 4c,d) and M (Fig. 4g,h) sharply decrease and approach some minimum value before K s and K p respectively. Mutual information M (Fig. 4g,h) in particular shows clear changes in trend well before the social transition K s occurs. Though N V , s shows a similar rising-falling pattern for both σ = . 0, 0 25, its maximum value with σ = .
0 25 (Fig. 4f) is much lower than that for σ = 0 (Fig. 4e). For σ = . 0 25, the mean of the Geary's C C (Fig. 4l) shows almost no change, though its envelope broadens post-transition; we see this as a failure of the EWS (no forewarning given). Similar observations hold for model V1 (Supplementary Information S5.2), with the failure of the Geary's C coefficient C shown clearly in Supplementary Figs. S6 and S7. As stated in the Methods section, the pre-and post-transition trends of the EWS do not generally resemble each other; asymmetry of the EWS about K s can be seen Fig. 4 and Supplementary Figs. S6 and S7, showing that (in general) less of a warning is given (if any) when the κ-series is reversed. This is explicitly demonstrated in Supplementary Figs. S15 and S16, where skewness γ 1 is used to quantify asymmetry of the trend of each EWS.
We can then say that all proposed EWS other than C give appreciable signals approaching K s and K p when σ = .
0 25 (Fig. 4B). K s precedes K p (Fig. 4a,b), showing that a shift in consensus will always precede a crisis in vaccination coverage in this model. Also shown is a marked decrease in − K K p s (the gap between the two transitions K s and K p , which we call the intertransition distance) as the social norm strengthens (for example, σ → . 0 25 in column (B) of Fig. 4). The generalisation of these trends to all tested values of σ is confirmed in Fig. 5a, where − K K p s is everywhere positive, though the distance between K s and K p vanishes with increasing σ ; the inset of Fig. 5a shows the location of K s (blue) and K p (red) with respect to σ, so that − K K p s (purple) gives the width of the area between the two curves in the inset graph at each σ. Other disparate models of the disease display largely similar concave decreases in the intertransition distance, suggesting that this behaviour arises generally from the model dynamics rather than in some specific subspace of the parameter space (see Supplementary Information S5.2).

Stronger social norms result in decreased lead time for all eWS.
The findings of the preceding subsection are intuitive, as vaccination depends more heavily on individual vaccination opinion than the number of infected neighbouring agents in Eq. (2), so that the opinion dynamics exert more influence than any feedback effect occurring in the physical dynamics. However, the vanishing intertransition distance − K K p s presents a problem if we depend on predictions of K s to enact interventions avoiding the collapse of the system to a non-vaccinated regime (i.e. avoiding K p ). For social norms of increasing strength, we can therefore look at the  K s , where Π Ψ represents some κ value at which we can assert that a signal occurs in some sequence Ψ of κ values; since we've established that K s precedes K p everywhere, then necessarily any warning of a social transition also warns of the following physical transition, so the quantity − Π Ψ K p is not discussed here.
One way for us to quantify this lead time is to use a change point detection tool to find κ values at which the two classes EWS (subscript WS) and model variables (subscript S/P) give signals (i.e. display statistically significant changes in trend/behaviour). Here, the Lanzante test 62 from the trend 63 (Fig. 5b), the positivity of some curves shows that some of the proposed EWS do indeed give early warnings of coming transitions (largely for σ ≤ . 1 875). Failures of N V , s and M occur in the range σ .
≤ ≤ . 1 875 2 125, while all other tests give valid warning signals everywhere 2 5 σ < . ; model variables R and 〈 〉 V p fail in the range σ . ≤ ≤ . 1 65 2 125 (Fig. 5c). The failure of all the tests after 2 5 σ = . likely results from insufficient length of the EWS' κ-series; the inset of Fig. 5a shows that → − K 1 s as σ increases. 2 5 (as in Fig. 5b). In line with our focus on social dynamics as a predictor, V s appears to be the best performing signal of all the model variables; as was reasonably expected, I appears to perform badly, since its role as a transitory compartment in the disease dynamics means that it never "gathers" sufficient agents over the course of each realisation to give a true indication of the state of the system (other than indicating the presence or absence of endemic infection).
Since the perceived cost of vaccination κ ≥ −1 in this study, our method of detecting the change point will not accurately predict a change point Π Ψ close to −1. Since no one warning signal gives the highest lead time for any large contiguous range of σ values, there is unfortunately no single objective way to choose a "strongest" signal; they are all suitable tools to predict coming crises in aggregate opinion and vaccination dynamics. However, it is worth noting that mutual information M and the dissimilar join count N V , s perform better than the other indicators; Fig. 5b shows that mutual information M gave the largest lead time measured for 45% of tested social norm σ values, and V V ,

eWS can provide better forewarning than trends in model variables. One final question is whether
the proposed EWS (mutual information, Moran's I, Geary's C, join counts) give earlier warnings than simply monitoring trends in model variables (such as the number of infections or pro-vaccine agents, using a change point test for prediction in both cases). There are many ways to quantify this, including maximin comparison (finding the larger of minimum values of classes WS and S/P) and maximax comparison (finding the larger of the maxima of each class) at each value of the social norm σ.
To compare the minima of the EWS and model variable lead times, we define χ min Lan as   36 4% of tested σ values. The second part of the comparison (Fig. 6, red curve) is between the maxima of the lead times; as above, we define the comparison variable χ min Lan and tolerance ε max Lan as   Fig. 6, there is no apparent pattern to the positivity/negativity of χ min Lan and χ max Lan . These comparisons are given for other tests and models in Supplementary Information S4.

Discussion
Here we studied a range of early warning signals for critical transitions in a two-layer coupled behaviour-disease model for paediatric infectious diseases. We compared the indicators to one another and the approach of simply monitoring trends in model variables. We found that the performance of the indicators was variable depending on model parameters, but the mutual information statistic and the dissimilar join count showed consistently high pre-transition lead times over various strengths of the social norm, many times giving the highest lead times of all the EWS. Through maximin and maximax comparisons, we found that using EWS provide more advance warning than simply monitoring trends in model variables in a clear majority of cases.
We note that join counts have the additional advantage of easy computability, since they require only counting pairs of a given type. This contrasts with other more computationally intensive indicators such as autocorrelation which require making decisions about whether to study lag-1 or higher order lags, as well as choosing parameter values governing computation of residuals. Moran's I was also shown to predict the approach of transitions, although perhaps this finding is trivial considering that it is a linear combination of similar and dissimilar join counts. Its predictive power was not as strong as many of the other indicators such as join count and mutual information, hence the added complexity of its calculation may not justify its use. Potential downfalls of the mutual information statistic include its computational complexity and the availability of a suitable data set pairing the personal health of each agent with their social activity.
We also showed that a population may have relatively high vaccine coverage despite a low pro-vaccine opinion. This discrepancy between social and physical dynamics is due to the paediatric infectious disease natural history we assumed in the model. Unlike influenza, where revaccination must occur seasonally, an individual who receives a sufficient number of measles or chickenpox vaccine doses generally has lifelong immunity and therefore the opinion towards the vaccine can decline well before the level of vaccine immunity does. (Individuals can change their opinion but never become 'unvaccinated' .) The implication of this asymmetry is that monitoring social media networks for changes in opinion using early warning signals like mutual information might provide advance warning of outbreak hot-spots.
The distance between the change point in the EWS indicators and the critical transition in the social dynamics decreases as the social norm grows stronger, as does the distance between the transitions in social and physical dynamics of the model. Given the relative scale of the social norm and vaccine risk parameter values used, stronger social norms decrease the time interval between birth and vaccination decision (the vaccination rate converges to its equilibrium value in fewer time steps than in other regimes); feedback between this and the disease incidence in the network (which affects the number of infected neighbours in each agent's neighbourhood) alters the probability function controlling the vaccination decisions, effecting faster alignment of majority opinion and vaccination coverage.
This study only lays the foundation for the investigation of spatial EWS for such coupled-behaviour systems. There is much work to be done before they can be meaningfully applied to empirical data. For instance, we assumed that the network was static. This simplifying assumption could be relaxed in future work by using an evolving social dynamic in which agents are allowed to form or break links with new agents based on their node degree 64 or vaccine opinion and associated social pressures 65 . The rate of interaction between agents was also assumed fixed in our model (relative to the speed of other dynamics, such as the birth/death interval and the length of illness). A valid extension of the model would be a variable rate of communication between agents, since the rate of communication has been shown to influence the rate and efficiency of opinion formation 66 . A further avenue of research would explore interventions to turn populations away from critical transitions. This could answer research questions such as: How far in advance must we act to prevent a social or physical critical transition, and does this change our interpretation of the EWS? If we assume that any of the EWS can be used for monitoring, how would this change in trend alter the reliability of the EWS?
Many researchers and public health bodies are drawing attention to global resurgences of vaccine-preventable illness, and speak to the vast efforts and multiple approaches taken to mitigating outbreaks. A few of these approaches have focused on human behaviour and opinion dynamics, either by directly tracking aggregate vaccine opinion, or monitoring alerts and media reports. Our work demonstrates the potential uses of early warning systems of critical transitions in preventative epidemiology. In particular, our work provides proof-of-concept for the idea of monitoring social networks for early warning signals of both social and epidemiological shifts, and also suggests several EWS indicators that might work well for this purpose.

Data availability
The datasets generated during and analysed during the current study are available from the corresponding author upon reasonable request.