Skip to main content
  • Other Publications
    • Philosophical Transactions B
    • Proceedings B
    • Biology Letters
    • Open Biology
    • Philosophical Transactions A
    • Proceedings A
    • Royal Society Open Science
    • Interface
    • Interface Focus
    • Notes and Records
    • Biographical Memoirs

Advanced

  • Home
  • Content
    • Latest issue
    • All content
    • Subject collections
    • Videos
  • Information for
    • Authors
    • Reviewers
    • Readers
    • Institutions
  • About us
    • About the journal
    • Editorial board
    • Author benefits
    • Policies
    • Citation metrics
    • Publication times
    • Open access
  • Sign up
    • Subscribe
    • eTOC alerts
    • Keyword alerts
    • RSS feeds
    • Newsletters
    • Request a free trial
  • Submit
Open Access

Anticipating the emergence of infectious diseases

Tobias S. Brett, John M. Drake, Pejman Rohani
Published 5 July 2017.DOI: 10.1098/rsif.2017.0115
Tobias S. Brett
Odum School of Ecology, University of Georgia, Athens, GA, USACenter for the Ecology of Infectious Diseases, University of Georgia, Athens, GA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tobias S. Brett
  • For correspondence: tsbrett@uga.edu
John M. Drake
Odum School of Ecology, University of Georgia, Athens, GA, USACenter for the Ecology of Infectious Diseases, University of Georgia, Athens, GA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Pejman Rohani
Odum School of Ecology, University of Georgia, Athens, GA, USACenter for the Ecology of Infectious Diseases, University of Georgia, Athens, GA, USADepartment of Infectious Diseases, University of Georgia, Athens, GA, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

In spite of medical breakthroughs, the emergence of pathogens continues to pose threats to both human and animal populations. We present candidate approaches for anticipating disease emergence prior to large-scale outbreaks. Through use of ideas from the theories of dynamical systems and stochastic processes we develop approaches which are not specific to a particular disease system or model, but instead have general applicability. The indicators of disease emergence detailed in this paper can be classified into two parallel approaches: a set of early-warning signals based around the theory of critical slowing down and a likelihood-based approach. To test the reliability of these two approaches we contrast theoretical predictions with simulated data. We find good support for our methods across a range of different model structures and parameter values.

1. Introduction

Tipping points—where small changes in circumstances precipitate dramatic shifts in state—are a feature of many natural systems. These sudden transitions can have devastating consequences, for example, irreversible climate change [1] or ecological collapse [2]. Forewarning, with sufficient time to act, is of clear importance.

In principle, careful scientific research and detailed mechanistic understanding could lead to the formulation of predictive models capable of anticipating such transitions. However the combination of nonlinearity, non-stationarity, noise and data availability is a barrier to mechanistic modelling for many natural systems. There has been a trend towards looking for alternative methods of anticipating transitions, which do not rely on an empirically validated model [3]. One promising avenue is the development of early-warning signals (EWS) [4], which are summary statistics derived from dynamical systems theory and are calculable directly from observed data.

Dynamical systems theory states that as a stable system approaches a tipping point, the time taken to recover after a perturbation increases; ultimately diverging at the transition [4]. This phenomenon, known as critical slowing down, is expected to be observable in systems which are continually subject to shocks—whose origin can be either extrinsic (for instance, environmental fluctuations) or intrinsic (such as demographic noise). The effects of critical slowing down are manifest in the time series, leading to systematic changes in the summary statistics. For example, due to the persistence of perturbations the time series becomes increasingly correlated, which can be measured by the autocorrelation and correlation time. Because EWS are informed by generic properties derived from dynamical systems theory they have a key advantage over model fitting approaches: detailed calibration is not required. EWS have been applied to anticipate transitions in a range of complex systems, including ecosystem collapse in lakes and microbial communities [5,6], desertification in arid ecosystems [7] and changes in climate [8].

Identifying reliable warning signals is especially timely for emerging infectious diseases. Examples of pathogens with pandemic potential include influenza virus [9] and SARS-like coronaviruses [10]. Similarly, pathogens which have been successfully controlled can re-emerge, for example, pertussis [11].

Constructing validated models for emerging diseases is complicated in part by complexity, and also by virtue of the pathogen's novelty in a new host. Key epidemiological determinants, such as the infectious period, the mode of transmission, infectiousness and rate of spillover, are typically unknown or poorly quantified. This makes them a prime candidate for model-free approaches, such as EWS. To date, work on anticipating disease emergence has predominantly focused on specific systems and models, for instance, in studies on the effects of climate change on disease emergence [12], with a previous attempt at applying EWS to simulated data finding the approach to emergence difficult to detect [13].

The scenario considered in this paper is based on the slow emergence of a pathogen. This can occur for a variety of reasons: as a result of changing immunological landscape (increasing the susceptible pool—cf. monkeypox [14] or measles [15]), pathogen adaptation to new host [16], immune evasion [17,18] or long-term transient dynamics (the ‘honeymoon effect’ [19]). These mechanisms all lead to increasing inter-individual transmission of the pathogen—characterized by an increase in R0, a threshold quantity. The primary concern for disease containment is averting the tipping point which occurs when R0 = 1.

Formally, R0 is defined as the average number of secondary infections caused by an infectious individual [20]. From theory we know that if R0 < 1 transmission is subcritical and the pathogen is incapable of long-term persistence. Instead, continued circulation can only be sustained through repeated reintroduction from an external reservoir. In this regime, the transmission dynamics are characterized by what are known as ‘stuttering chains’ [21]. On the other hand, if R0 > 1 the disease will typically invade successfully and grow in a fully susceptible population, leading to a large-scale outbreak and possible endemicity [22]. Such dynamics are referred to as supercritical. It is this transition, from subcritical to supercritical, that we seek to anticipate.

The birth–death–immigration (BDI) process can be viewed as the simplest possible model of emerging diseases with direct transmission, modelling spillover of a pathogen from an external source, together with secondary chains of transmission fed by a large susceptible pool. A major benefit in using the BDI process is that it has a known analytical solution [23], which we use as the basis for the two approaches presented. For the first approach, we derive key metrics which may be used as EWS. The second approach uses a likelihood-based method to formally test emergence as a hypothesis.

As they are derived from a very generic model of disease emergence, the approaches are expected to have applicability across directly transmitted emerging disease systems. We test the robustness of measures to changes in model structure, simulating disease emergence using the SIS and SIRS compartmental models subject to varying immunology and demographics. Additionally, we comment on how the speed at which the transition is crossed affects the reliability of the two approaches, something which has been considered poorly understood in previous work on EWS [24].

2. Early-warning signals using the birth–death–immigration process

Previous work exploring EWS for infectious disease dynamics has shown that performance depends on the direction the critical transition is traversed, with the two contrasting cases shown in figure 1 [13]. In studies employing a combination of analytical results and simulated data, EWS have successfully been identified for diseases approaching elimination—for instance, due to increasing vaccine uptake—where R0 approaches the critical transition from above [13,27]. The case of disease emergence has thus far proven more challenging [13].

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

(a) Simulation of disease elimination through vaccination. At time T = 8 years the vaccine threshold for elimination is reached. Weekly data were generated using the SIR model with birth and death [20]. The mean infectious period 1/γ = 14 days, the death rate μ = 0.02 yr−1, R0 = β/(γ + μ) = 5 and the population size N = 4 × 105. The birth rate of susceptible individuals is affected by vaccination via ν = μ(1 − v), where v is the vaccine uptake, which linearly increases from 0 to 1 over the 10 year period. The elimination threshold is passed when v = 1 − R−10 = 0.8. (b) Simulation of disease emergence through increasing infectiousness. Simulations were performed using an SIS model with 1/γ = 7 days, N = 103 and with R0 = β/γ increasing linearly from 0 to 1.25 over the 10 year period. Susceptible individuals acquire infection from external sources with per capita rate ζ/N = 7 × 10−4 d−1. All simulations in this paper are performed using the NRM algorithm [25,26]. (Online version in colour.)

This challenge can be understood from a theoretical perspective. The conventional theory of EWS, as detailed, for example, by Scheffer et al. [4] and used by O'Regan & Drake [13], relies on the presence of critical slowing down, presupposing that the system is at a stable fixed point and subject to exponentially decaying perturbations. When the number of infectious individuals is sufficiently large the transmission dynamics can be well approximated as a linearized system subject to Gaussian noise [28,29].

For endemic diseases prior to elimination the number of infected individuals is large enough for this approximation to be valid, and EWS perform well [13]. Conversely, for subcritical emerging diseases there are typically very few infected individuals present in the population. This means that the effects of demographic stochasticity feature strongly in the dynamics [29]. The stochastic dynamics are highly non-Gaussian and the conventional theory of critical slowing down is invalid. In the light of these findings, we instead develop EWS starting from an alternative model, the BDI process, which explicitly incorporates transmission between individuals. By starting from a model which captures the conditions preceding disease emergence, we show how EWS can be expected to behave for this class of transition, and assess their prospects as generic indicators.

The BDI process is one of the simplest models of a subcritical disease [23]. It neglects the effects of susceptible depletion, assuming that only a small fraction of the population are infected at any one time, and that the pool of susceptibles is replenished with sufficient speed. The BDI process focuses on three simple transitions: (i) an infected individual is introduced into the population, (ii) the infection can spread to other individuals with a rate proportional to the number of infected individuals present, and (iii) infected individuals may recover from the disease.

We assume that at a constant rate ζ individuals import the infection due to contact with external sources, for example a zoonotic reservoir. We focus on diseases which are emerging due to changes in the transmission rate of the disease, β(t). If there are n infected individuals present in the population, then the total force of infection is λ(t) = β(t)n + ζ, and the rate of recovery is γn. Throughout the paper we assume that changes in β(t) occur on a much slower time scale than the typical duration of an infection, 1/γ.

The transition rate Tm,n is defined as the probability per unit time of transitioning from a state with n infected individuals to a state with m infected individuals. The BDI process is a one-step stochastic process which consists of two possible transitions; for infectious diseases these are infection and recovery. Respectively, the transition rates for infection and recovery areEmbedded Image 2.1In the BDI process, the chain of transmission stemming from a particular introduced index case is given by a branching process [30]. A particular outbreak can be considered as a superposition of the separate chains of transmission caused by each introduced case during the outbreak [31]. The basic reproductive number is defined as the average number of secondary cases, R0 = β/γ, found using the offspring distribution of the associated branching process [30,32].

The probability of n individuals being infected at time t is Pn(t). The change in probability with time is found by solving the master equation, a set of coupled linear differential equations built from the transition rates. For the BDI process,Embedded Image 2.2The advantage in considering such a simple model is that the master equation, which determines how the probability distribution of the number of infectious individuals changes in time, can be solved exactly without the need for any approximations [23]. This can be achieved through use of the moment generating function, Embedded Image. The variable ψ is used to find the moments and correlation functions of the stochastic process via differentiation, for example, the ith moment, Embedded Image. We find from the master equation that Z(ψ, t) solves the partial differential equationEmbedded Image 2.3Further details of the master equation calculation are presented in appendix A.

Using that R0 = β/γ, the mean-field behaviour of the system solves dμ1/dt = (R0 − 1)γμ1 + ζ. If R0 > 1 then μ1 grows exponentially, whereas if R0 < 1 the disease persists at a low level, μ1 = (ζ/γ)/(1 − R0), sustained by immigration.

Table 1 shows the list of candidate EWS calculated from the moment generating function of the stationary BDI process which will be considered in this paper. The mean, variance, index of dispersion (also known as the variance to mean ratio), and correlation time all diverge as R0 → 1. The autocorrelation approaches one, whereas the coefficient of variation remains constant below the transition.

View this table:
  • View inline
  • View popup
Table 1.

Measurable quantities calculated from the BDI process. All expressions are valid for 0 ≤ R0 < 1. See equations (A 3), (A 4) and (A 7) for details.

The correlation time gives a measure of the persistence of correlations. Its divergence as R0 → 1 implies that correlations persist for extremely long times, and is the mathematical definition of critical slowing down. Intuitively, critical slowing down can be understood by considering the lengths of transmission chains. Although on average each imported case will infect fewer than one individual (as R0 < 1), there is a possibility that a large chain of secondary cases is sparked, which then takes a long time to die out. As R0 increases, the probability of larger chains of transmission increases, in turn leading to an increasingly autocorrelated dynamics. The increased probability of larger chains of transmission also causes an increase in the mean and variance of the number of infected.

In addition to the EWS derived from the BDI process presented in table 1 we also study two signals from information theory: the Kolmogorov complexity, and the Shannon entropy [33]. The Kolmogorov complexity is a measure whose origin lies in the length of the minimal computer program needed to reproduce the sequence of data, with longer computer programs implying more complex data [33–35]. The increased correlation near the transition may therefore result in a decrease in complexity, which can serve as an EWS. The Shannon entropy is a measure of the information content of the time series. As R0 approaches one, the probability of rare large outbreaks increases, which suggests the entropy will also increase near the transition. For details of how these two measures are computed, see appendix B.

Two simplifying assumptions are implicit in the model. Firstly, the BDI process neglects susceptible depletion. Secondly, the speed with which β changes is ignored, assuming that it varies infinitely slowly. The effects of relaxing the first of these assumptions will be investigated in the next section. Disease emergence over faster time scales will be investigated in §5.

3. Comparison with models featuring susceptible depletion

For EWS to be useful they must be robust to changes in model structure, and not just specific to the BDI process. Using results from simulations of common disease models, in this section we explore the sensitivity of the EWS to a range of demographic and immunological properties.

The SIS model is a model of infectious disease spread in which individuals move from the susceptible class to the infected class upon infection, and then return to the susceptible class when they cease being infectious. The total population size is assumed constant. The SIRS model is similar to the SIS model, but in addition the pathogen confers temporary immunity to reinfection, which wanes at rate ρ [20]. The dynamics of both models can be represented using a master equation, as with the BDI process. Exact numerical realizations of these models can be obtained using the next reaction method (NRM) algorithm [25,26]. The candidate EWS for the SIS and SIRS models are shown in figure 2. The results for the BDI process listed in table 1 are also plotted.

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

EWS for the stationary SIS and SIRS models. Red lines correspond to the theoretical results found from the BDI model (table 1). Symbols depict the results of simulations with the duration of immunity as indicated in the figure. Independent simulations are performed for each value of R0. EWS are calculated by time averaging over 4 × 103 weeks of data (see appendix D). The force of infection in the SIS and SIRS models is λ = (βn + ζ)/N, where n is the number of infected and N the total population size. The remaining model parameters are γ = 1 and ζ = 1 week−1. The autocorrelation is calculated at lag one week, AC(1). (Online version in colour.)

The only additional parameter in the SIS model is the population size, N, allowing the effects of susceptible depletion on EWS to be studied in isolation (figure 2a). We find that typically, the effect of reducing the population size is to decrease the responsiveness of a signal to increases in R0. The mean and variance always increase prior to the transition; however, the magnitude of the increase is reduced. The value of R0 for which the variance is maximized is seen to increase. A similar effect is also observed for the index of dispersion, autocorrelation and correlation time. These three EWS all peak around R0 = 1 for large population sizes, but for N = 102 there is a noticeable shift in the peak position and height to the right. For the coefficient of variation, decreasing the population size smooths the abrupt change in value at R0 = 1. The Kolmogorov complexity falls as the transition is approached, in line with the increased correlation: if there is a fluctuation above (or below) zero within the detrended data it is more likely to be sustained in further data points. Regular sequences allow for increased use of the copy operation (see appendix B), reducing the complexity of the time series. Perhaps counterintuitively, the complexity of the time series is lower near the transition for larger populations. The entropy has a very similar behaviour to the variance, peaking in the vicinity of the transition. Overall, the results from the BDI process appear to provide an upper bound on the values of the EWS in the SIS model.

Figure 2b,c shows results from the SIRS where the mean duration of immunity, 1/ρ, is 5 and 25 years. For the smallest population size, N = 102, the mean number of infected individuals is less than one, irrespective of R0. As might be expected, EWS vary little with R0 and do not display the trends observed for the BDI process. For larger population sizes, the behaviour of some EWS above the transition is markedly different from the SIS model, for example, the autocorrelation remains close to one and the complexity decreases with R0. Below the transition, the opposite is the case for N > 102, with the EWS having similar behaviour to the BDI process. For a given population size, the responsiveness of the EWS decreases with immune duration, cf. the similarity in the index of dispersion and Kolmogorov complexity for N = 104 in the SIS model and N = 106 in the SIRS model. Apart from the coefficient of variation, the BDI process results again provide upper bounds on the EWS, as with the SIS model.

Below the transition, decreasing the population size and increasing the immune duration reduce the responsiveness of the EWS. This can be understood by returning to the theory of critical slowing down. As the susceptible pool is depleted the probability of further disease transmission is diminished, reducing the probability of extremely long chains of transmission. The impact of susceptible depletion is greater for smaller population sizes. Increasing the duration of immunity means individuals who acquire the infection are absent from the susceptible pool for longer, magnifying the impact of susceptible depletion. The longer the duration of immunity the larger the susceptible pool required for agreement with the BDI results.

Below the transition each EWS behaves similarly regardless of the model, provided that the population size is large. This independence of model structure supports the general applicability of BDI results to more complex dynamics. We note that there is a large distinction in behaviour between models above the transition; for instance the coefficient of variation strictly decreases for the SIS model, whereas for the SIRS model it may decrease, increase or remain approximately flat depending on the model parameters. The generality of the EWS is strictly for disease emergence.

4. Cox's δ: a likelihood-ratio test

An alternative model-specific approach to detecting critical transitions has been proposed by Boettiger and Hastings, focusing on anticipating the loss of stability of fixed points [24]. They propose likelihood-ratio testing, using a statistic they refer to as Cox's δ. The test determines whether the dynamics are approaching a transition or are instead stationary. Each hypothesis is represented by a model, and a maximum-likelihood estimate (MLE) for the data is found. The value of Cox's δ provides a measure of how much disease emergence is favoured. The procedure by which Cox's δ is calculated also allows for an estimation of the time remaining until the transition is reached.

Boettiger and Hastings assume that the Ornstein–Uhlenbeck (OU) process is an appropriate model for the dynamics. The OU process is a continuous-time stochastic process in which there is deterministic reversion to the mean, with Gaussian white noise perturbing the system [28]. It has a known mathematical expression for the likelihood of a time series, which Boetigger and Hastings use to efficiently calculate Cox's δ and thereby determine whether the system is approaching a transition [24].

A barrier to using Cox's δ as a method for anticipating disease emergence is this reliance on the OU process to calculate the likelihood. Dynamically, the OU process does not constrain the number of infected individuals to be a non-negative integer. When there are a small number of infected individuals present in the data, this means biologically impossible paths are included in the likelihood calculation. The problem can be overcome by instead using the BDI process, which presents a more biologically plausible model of disease emergence.

An exact solution to the transition probability exists for the BDI process from which the likelihood of a time series can be constructed (see appendix C). We assume R0(t) can be decomposed into two parts, a baseline which is set as R0(0) and a linear trend, ΔR0, i.e.Embedded Image 4.1For the test hypothesis (the disease is emerging), the MLE is calculated from a likelihood surface allowing both of these parameters to vary. In calculating the MLE for the null hypothesis (no emergence), we fix ΔR0 = 0, still allowing R0(0) to vary. For an observed time series {n}Tt=0, the log-likelihood of the MLE for the test model isEmbedded Image 4.2and for the null model isEmbedded Image 4.3The expressions for PBDI are presented in appendix C.

Cox's δ statistic is defined as twice the difference in the log-likelihoods of the MLEs for the two nested models,Embedded Image 4.4Because the models are nested δ≥0. Increasingly positive δ implies stronger support for the test model over the null model. A δ > 2 is taken as the criterion for concluding that the disease is emerging, on the basis of the AIC score [24].

Figure 3 shows the log-likelihood surface calculated using the BDI process for a time series generated using the SIS model. The MLE of the test model is R0(0) = 0.00 and ΔR0 = 1.83 × 10−3 week−1. The critical transition is predicted to occur after 10.52 years, very close to the true value of 10 years. Cox's δ statistic is 10.78, implying that the hypothesis of disease emergence is statistically supported [24]. One caveat is that this strong result has been achieved under idealized model conditions, assuming knowledge of the rates of immigration and recovery, and calculated using a long sequence of data (312 weeks of data). In the following section, we explore further the reliability of Cox's δ at anticipating transitions and compare with the EWS.

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Estimating the rate of emergence using maximum-likelihood methods. Simulation generated using an SIS model with R0(t) = R0(0) + ΔR0t increasing linearly from 0 to 1.2 over 12 years. Data are output weekly. The likelihood surface is calculated using only data from the first 6 years for parameters on a 100 × 100 grid of uniformly spaced data points with ΔR0(0) ∈ [0, 0.008] and R0(0) ∈ [0, 0.8]. The gold star indicates the location of the true parameter set and the red dot indicates the parameters corresponding to the maximum-likelihood estimate (MLE). Simulation performed using the modified next reaction method [26]. The remaining model parameters are fixed to their true values to reduce computation time; 1/γ = 1 week, ζ = 1 week−1 and N = 104 individuals. The MLE is that R0(0) = 0.00 and ΔR0 = 1.83 × 10−3 week−1. The true values are R0(0) = 0.00 and ΔR0 = 1.92 × 10−3 week−1. Side panels show the likelihood profiles. (Online version in colour.)

5. Impact of time scale on reliability of predictors

In §3, we presented results for stationary dynamics, appropriate for a disease emerging over an extremely long time scale. In practice, we are interested in anticipating disease emergence over faster time scales, as shown in figure 3. In this section, we investigate how the EWS and Cox's δ perform as indicators of disease emergence. Results for the EWS are shown in figures 4 and 5, and results for Cox's δ are shown in figures 6, 7 and 8.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Performance of candidate EWS for various window sizes and transition speeds. Simulations are performed using an SIS model. Window sizes are 100, 200 and 300 weeks. The population consists of 1000 individuals, all other parameters are the same as in figure 2. Red lines correspond to the stationary results for the BDI model, green lines to simulated time series. EWS are calculated from the number of infected individuals in the population, data output weekly. The autocorrelation is calculated at lag one week, AC(1). (Online version in colour.)

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

(a) 95% CIs of EWS for the SIS model with N = 1000. R0 increases linearly from 0 to 2. The total simulation time is T = 20 years and T = 100 years. Data are output weekly, EWS are calculated using a window size W = 52 weeks for both columns. The average of the replicates is shown in black and a sample replicate is shown in dark blue. (b) Numerical estimate of the probability that R0 > 1, given an observed value of an EWS. The noisiness in the estimate for more extreme values of the EWS is due to the infrequency of these events. (Online version in colour.)

Figure 6.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 6.

Cox's δ calculated using the BDI process solution. The time-series data are the same as used in figure 5, with the additional stationary time series generated by running the SIS model for 40 years with R0 = 0.5. Model parameters are N = 1000, γ = 1 and ζ = 1. Time to emergence is given as a fraction of the initial time to emergence. In the bottom row, the log-likelihood is shown for the test (green) and null (red) models for window size W = 100 weeks. (Online version in colour.)

Figure 7.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 7.

Probability that Cox's δ > 2, calculated from 200 replicates. Cox's δ computed using the likelihood-ratio test and BDI process solution, see §4. Fixed model parameters are N = 1000, γ = 1 and ζ = 1. (Online version in colour.)

Figure 8.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 8.

Mean and 95% CIs of the remaining time to emergence. Time to emergence calculated using the MLE, Te = W + (1 − R0(0))/ΔR0, where W is the window size. True value of the time to emergence shown in black. Calculated using the same data as in figure 7. (Online version in colour.)

Figure 4 shows the EWS calculated for the SIS model, with the stationary BDI process results also shown for comparison. The EWS are calculated from the time-series data; for simplicity, we use unweighted moving averages (see appendix D). The time series were generated by an SIS model with linearly increasing R0; from 0 to 2 over T = 10, 20, 30 and 40 years. The moving window average is calculated using three window sizes: W = 100, 200 and 300 weeks.

As shown in figure 4, all EWS bar one have a clear response to the approaching transition as long as the time scale is longer than 10 years (T > 10). For instance, the variance increases for all window sizes and speeds prior to R0 = 1, in line with predictions. The exception is the coefficient of variation, which performs poorly as an indicator of transitions, remaining largely unchanged until R0 > 1. It is also highly sensitive to the choice of window size, the sharp peak following the transition observed for T = 10 and 20 years increases with window size. The autocorrelation performs particularly well, reaching AC(1) > 0.85 for T = 20, 30 and 40 years. It does not appear particularly sensitive to window size below the transition.

As the speed of emergence is increased, the lag behind the BDI process results increases. Overall, we see that the transition speed has three notable effects which adversely impact all EWS. The first effect is dynamical: the incidence level is continually responding to the changing equilibrium resulting from increasing R0, leading to a lag between the values of the EWS calculated at equilibrium and those observed in the time series. Combined with demographic stochasticity, this leads to an apparent ‘bifurcation delay’ as the transition is traversed [36,37].

The second effect arises in computing the EWS. The theoretical results for the EWS are calculated by finding the moments and autocorrelation of a stochastic process. Owing to the absence of multiple replications, these are calculated over a moving window [4]. However, when the disease emerges over a finite time scale, the process is non-stationary, and therefore also non-ergodic, due to changing R0(t). In practice, this means that averages are calculated using data points which were generated at different values of R0; for example, if T = 10 years and W = 100 weeks, the difference in R0 between the start and end of the window is 0.38. This reduces the responsiveness of the EWS, also contributing to the lag behind the BDI process results. Furthermore, there can be other artifacts arising from non-ergodicity, for instance, increased variance of the data. Increasing both the window size and transition speed exacerbates these effects due to the larger range in values of R0.

Thirdly, there is Monte Carlo error due to the limited number of data points used. This is manifest in the stochastic variation of the observed EWS as R0 changes. This is particularly apparent in the first column of figure 4, where a chance drop in prevalence immediately prior to R0 = 1 leads to a failure of most EWS to behave as theory predicts. Increasing the size of the window has the effect of reducing this error (evidenced by the smoother curves for W = 300), however, at the expense of the signal's responsiveness, as detailed above.

Given the existence of Monte Carlo error, to quantify the performance of the EWS further, we compute 90% CIs for the EWS for disease emergence over 20 years and over 100 years using a window size of 52 weeks (figure 5a). A measure of how effective an EWS is at detecting emergence is the range in R0 over which a specific reading of an EWS lies within the 90% CI. A wide range means that the EWS gives little information about the value of R0. For the variance, a reading of σ2 = 40 falls within the 90% CI for 0.8 ≲ R0 ≲ 1.1, implying a strong chance that the system is near the transition. A reading of K = 0.8 for the Kolmogorov complexity falls within the 90% CI for all values of R0 shown, implying a poor predictive power.

To quantify this further, figure 5b shows an estimate of probability that R0 > 1 given a particular reading of an EWS, y (appendix E provides details on how this estimate is calculated). The sensitivity of P(R0 > 1 | y) to changes in y depends heavily on the EWS. Large values for the mean, variance and entropy strongly indicate that transmission of the pathogen is supercritical. On the other hand, for EWS which peak near R0 = 1, such as the index of dispersion, autocorrelation and correlation time, there is low indication whether the disease is supercritical or not, unless the observed value is especially large.

Because of their narrow confidence intervals, the variance and entropy are exquisitely sensitive, with P(R0 > 1 | y) increasing from near 0 to 1 over a small range of y. Others, such as the index of dispersion and autocorrelation, have an intermediary range of values of y where P(R0 > 1 | y) ≈ 0.5.

We see a trade-off between how reliable an EWS is at classifying whether R0 > 1 and the range of values of R0 over which P(R0 > 1 | y) notably increases. A practical conclusion from this is that sequential readings of an EWS are needed to ascertain whether a disease is emerging.

Using a moving window, in figure 6 we show Cox's δ calculated for the same set of time series studied in figure 4. As the transition is approached Cox's δ typically increases in significance, with large responses for window sizes W = 200 and 300 weeks. There are spurious increases in Cox's δ for the stationary time series, due to the intrinsic variability of the dynamics, but these are dampened by using a larger window.

To test its reliability as an indicator of disease emergence, figure 7 shows the probability that δ > 2 as the transition is approached.

Large window sizes lead to an improved performance, as does a faster speed of emergence. For faster speeds of emergence (T ≤ 20 years) and larger windows (W ≥ 200 weeks), Cox's δ performs well, with P(δ > 2) > 0.75 prior to the transition. From figure 7, the recommendation is to use the largest window size feasible given the data.

The effect of the speed of emergence and the window size on the performance of Cox's δ is more subtle than for the EWS shown in figure 4. For EWS, the performance improves when the disease emerges at a slower rate due to the smaller change in R0 between the start and end of the window. The opposite is the case for Cox's δ. Instead, slower emergence reduces the difference in likelihood between the emerging and non-emerging models, reducing the power of Cox's δ.

From the estimates of the baseline R0 and ΔR0 the time to emergence can be estimated (figure 8). The mean estimated time to emergence generally does agree well with the true value, with typically marginally better agreement for smaller window sizes. However, the 95% CIs are large, in line with the variation in the MLE of ΔR0 seen in figure 6. Increasing the window size reduces to the confidence intervals, with a more pronounced reduction in the probability of underestimating the time to emergence, cf. the difference between W = 100 and 300 weeks for T = 20 years. With sufficiently large window size given the speed of emergence, the estimated time to emergence can be taken as an overestimate.

6. Discussion/conclusion

This work is intended to provide a theoretical base for the development of methods to anticipate disease emergence. We have presented two parallel approaches by which this may be accomplished. The first approach is to extend the burgeoning literature on EWS, reformulating them in a manner appropriate for emerging diseases. The second approach makes use of Cox's δ, a likelihood-ratio test for emergence. For both approaches, the application to emerging diseases is achieved through use of the BDI process, a well-understood model from the theory of stochastic processes.

Both the EWS and Cox's δ perform best when calculated using a large number of data points as this reduces statistical uncertainty. The time scales imposed by real-world disease emergence likely mean that available data are more limited. We therefore apply both approaches to weekly time series data of diseases emerging over a range of different intervals. We find that—apart from the coefficient of variation—all EWS undergo similar behaviour regardless of the time scale of emergence. More advanced time averaging methods may be needed for very fast emergence to detect a strong signal. Although larger window sizes do reduce uncertainty they are typically not necessary, with W = 100 weeks being sufficient. Cox's δ performs best with the largest window size possible so that there is the largest possible difference in R0 between the start and end of the window.

Overall we conclude that the two methods are complementary, performing best under differing conditions. Cox's δ reliably detects emergence over fast time scales, where the null hypothesis is strongly disfavoured. EWS perform best at detecting long term trends in the time series, where there is a more gradual approach to the epidemic transition. Additionally, our results favour the use of some EWS over others. The behaviour of the coefficient of variation and Kolmogorov complexity prior to the transition means they both perform poorly as indicators of emergence. The remaining EWS (the mean, variance, index of dispersion, autocorrelation, correlation time and entropy) have a behaviour which is resiliant to parameter changes and reliable in the face of stochasticity, making them strong candidate EWS.

By grounding our results in theory, rather than a detailed model for a particular disease we expect our results will have applicability to a wide range of emerging infectious diseases. This generality is of particular importance for emerging infectious diseases, where there are likely many unknowns regarding transmission of the pathogen.

The key simplifying assumption of the BDI process which makes it mathematically solvable is that it neglects the impact of disease transmission on the availability of susceptible individuals. In reality, infection diminishes the size of the susceptible pool and can also confer immunity to reinfection, delaying the individual's return to susceptibility. Our results show how these demographic and immunological considerations impact on the performance of EWS. Although the signals' strengths are reduced, for all parameters considered the performance is still informative for a population of 106 individuals.

Our work is not the first to investigate the signature of criticality prior to the epidemic transition; in particular, we note the work of Jansen and collaborators [15]. These previous works seek to identify the approach to the critical transition through changes in the tail of the outbreak size distribution [15,38]. In contrast, both approaches in this paper are based around summary statistics. This has advantages: the sizes of individual distinct outbreaks can be difficult to identify, a large number of outbreaks are required to accurately reconstruct the tail of the outbreak size distribution, and summary statistics calculated using moving windows clearly display temporal variation in the dynamics.

The SIS and SIRS models considered in this paper have only one susceptible and one infected class, assuming all individuals are equally infectious and interact homogeneously. Risk structure has been shown to significantly change the necessary conditions for epidemics, with R0 = 1 being an underestimate for the epidemic threshold [39]. For diseases emerging in populations structured as many loosely connected smaller community units, such as Ebola in sub-Saharan Africa, this may significantly impact on predictability. R0 is a measure based on the average context; however heterogeneities in the contact structure of individuals can lead to instances of superspreading [22], particularly relevant for sexually transmitted diseases and something not captured in the BDI process. The large variations in transmission rates reduce the worth of R0 as a measure, and a universal increase in R0 is unlikely to be the driver of emergence. Superspreading events have been identified for SARS and more recently for MERS [21,40]. In addition, there is the challenge of incomplete and unreliable data, modelled through explicit inclusion of a reporting process [41]. Further work is ongoing studying EWS and Cox's δ for these dynamics.

This work lays theoretical foundations, but more research is needed to develop actionable technologies applicable to actual disease data. For instance, improved reporting would lead to an increase in the mean number of case reports, but not an increased autocorrelation. A simultaneous increase in both presents stronger evidence of disease emergence. The development of ways in which readings of multiple EWS and also Cox's δ can be combined and leveraged to improve the quality of predictions is therefore desirable. One motivation in studying EWS is to ultimately develop comprehensive EWS software packages. These packages would be plug and play ready: public health practitioners could use them to give quantifiable information about the risk of a critical transition and disease emergence in a particular context.

As a final note, while this paper has focused on anticipating infectious disease emergence this does not limit the applicability of our methods. Emerging diseases are an example of an invasion process, where a new class of individual attempts to establish itself in a population. Other examples are found across the biological sciences. The BDI process is an appropriate model for a large range of invasion processes. We propose that the methods in this paper will have a similarly general status for anticipating successful invasion. Our work, therefore, opens up a new class of transitions which it may be possible to anticipate.

Data accessibility

All code and data are available at Dryad: http://dx.doi.org/10.5061/dryad.7kt3h [42].

Author's contributions

T.S.B., J.M.D. and P.R. designed the study. T.S.B. designed and wrote the software. T.S.B. and P.R. analysed and interpreted the data. T.S.B. drafted the manuscript. T.S.B., J.M.D. and P.R. reviewed and edited versions of the manuscript. All authors agree with the conclusions of the paper.

Competing interests

We declare we have no competing interests.

Funding

Research was funded by the National Institute of General Medical Sciences of the National Institutes of Health (award no. U01GM110744).

Appendix A. Critical slowing down for the stationary birth–death–immigration process

The master equation for the BDI process, equation (2.3), can be solved both when R0, γ and ζ are constant or when they are functions of time [23]. For mathematical simplicity, we will focus on the case in which the disease is emerging over a long time scale relative to the infectious period. Under these conditions, we can find exact theoretical values for the EWS. If there are y infectious individuals initially present then Z(ψ, 0) = eψy and equation (2.3) has the solution [23]Embedded Image A 1with A(ψ) = (1 − R0 eψ) and B(ψ, t) = e(R0−1)γt(eψ − 1).

We are interested in dynamics prior to the emergence of disease so restrict ourselves to considering R0 < 1. As τ → ∞, we see that B(ψ, t) → 0. The generating function converges onto the stationary solutionEmbedded Image A 2for any number of initial infected. This is the moment generating function of a negative binomial distribution with mean and varianceEmbedded Image A 3andEmbedded Image A 4Both equation (A 3) and equation (A 4) diverge as R0 → 1. Using equations (A 1) and (A 2), we can calculate the correlation function in the stationary state. The conditional mean, Embedded Image, isEmbedded Image A 5and therefore the correlation function Embedded Image isEmbedded Image A 6The autocorrelation, defined for a stationary process as AC(τ) = (C(τ) − μ21)/σ2, isEmbedded Image A 7As the critical transition is approached, the correlation time Embedded Image diverges. This is the hallmark of critical slowing down.

Appendix B. Shannon entropy and Kolmogorov complexity

The Shannon entropy is a measure of the information content of a random variable, and is defined as Embedded Image [33]. For time series data, an approximation to the probability distribution over n is calculated numerically via Embedded Image, where n(s) is the number of infected individuals present at time s and Embedded Imagen=n(s) is the indicator function, equal to one if n = n(s) and zero otherwise. As mentioned above, when the process is non-stationary this method is inexact.

Another information theoretic measure is the Kolmogorov complexity [33]. The Kolmogorov complexity of a sequence of data is the shortest computer algorithm required to reconstruct it. A quantification of the Kolmogorov complexity was proposed by Kaspar & Schuster [34]. Before applying their algorithm the data are converted to sequence of bits, 0 or 1. Their algorithm proceeds by considering two operations: ‘copy’ and ‘insert’. Each operation adds one to the complexity measure. The Kolmogorov complexity is the minimum number of operations required to reconstruct the data sequence. For further details, see [34]. The specific algorithm used to calculate the Kolmogorov complexity can be found in [35]. As with the other EWS, we calculate the Kolmogorov complexity over a moving window. To convert the time series to a sequence of bits, we perform a linear detrending over the data within the window. The residual time series is discretized into 1 whenever the residual is larger than zero, and 0 otherwise.

Appendix C. Transition probability distribution and likelihood

The generating function given in equation (A 1) can also be used to find the the probability of there being nt+δ infected at time t + δ given nt infected at time t [28]. Assuming that R0 does not vary significantly between time t and t + δ, the transition probability for the BDI process isEmbedded Image C 1whereEmbedded Image C 2andEmbedded Image C 3with ρ = [1 − AC(δ)]/[1 − R0AC(δ)] and p = μ1/(μ1 + ζ/γ) = 1/(2 − R0). Given a time series {nt}Tt=0, the likelihood of the data for a given rate of emergence, R0(t) = R0(0) + ΔR0t, isEmbedded Image C 4The MLE for (R0(0), ΔR0) is given by Embedded Image Embedded Image.

Appendix D. Calculating early-warning signals from single time-series data

If the BDI process is stationary then the process is ergodic, allowing expectation values to be calculated from a time series by time averaging. For example, the correlation function isEmbedded Image D 1The sum is performed over integer multiples of the timestep δ, with the approximation becoming exact as W → ∞. Although not formally equivalent, and therefore potentially introducing errors (see the discussion in §5), we will also use this approach to calculate expectation values for non-stationary processes, in line with the early-warnings literature.

Appendix E. Estimating the probability of supercriticality

Suppose we observe that EWS Θ has a value y. We would like to know the probability that R0 > 1 given this observation, P(R0 > 1 | Θ = y).

Monte Carlo sampling allows for the estimation of P(Θ = y | R0 = r). Through use of Bayes' theorem,Embedded Image E 1Assuming a prior P(R0 = r) = 1/Λ if r ∈ [0, Λ] and 0 otherwise,Embedded Image E 2Finally, P(R0 > 1 | Θ = y) can be found numerically by binning recorded values of Θ and approximating the integrals as summations over all bins.

  • Received February 16, 2017.
  • Accepted June 9, 2017.
  • © 2017 The Author(s)

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.

References

  1. ↵
    1. Lenton TM,
    2. Held H,
    3. Kriegler E,
    4. Hall JW,
    5. Lucht W,
    6. Rahmstorf S,
    7. Schellnhuber HJ
    . 2008 Tipping elements in the Earth's climate system. Proc. Natl Acad. Sci. USA 105, 1786–1793. (doi:10.1073/pnas.0705414105)
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Scheffer M,
    2. Carpenter SR
    . 2003 Catastrophic regime shifts in ecosystems: linking theory to observation. Trends Ecol. Evol. 18, 648–656. (doi:10.1016/j.tree.2003.09.002)
    OpenUrlCrossRefWeb of Science
  3. ↵
    1. Drake JM,
    2. Griffen BD
    . 2010 Early warning signals of extinction in deteriorating environments. Nature 467, 456–459. (doi:10.1038/nature09389)
    OpenUrlCrossRefPubMedWeb of Science
  4. ↵
    1. Scheffer M et al.
    2009 Early-warning signals for critical transitions. Nature 461, 53–59. (doi:10.1038/nature08227)
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. Carpenter SR et al.
    2011 Early warnings of regime shifts: a whole-ecosystem experiment. Science 332, 1079–1082. (doi:10.1126/science.1203672)
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Dai L,
    2. Vorselen D,
    3. Korolev KS,
    4. Gore J
    . 2012 Generic indicators for loss of resilience before a tipping point leading to population collapse. Science 336, 1175–1177. (doi:10.1126/science.1219805)
    OpenUrlAbstract/FREE Full Text
  7. ↵
    1. Kéfi S,
    2. Rietkerk M,
    3. Alados CL,
    4. Pueyo Y,
    5. Papanastasis VP,
    6. ElAich A,
    7. De Ruiter PC
    . 2007 Spatial vegetation patterns and imminent desertification in Mediterranean arid ecosystems. Nature 449, 213–217. (doi:10.1038/nature06111)
    OpenUrlCrossRefPubMedWeb of Science
  8. ↵
    1. Dakos V,
    2. Scheffer M,
    3. van Nes EH,
    4. Brovkin V,
    5. Petoukhov V,
    6. Held H
    . 2008 Slowing down as an early warning signal for abrupt climate change. Proc. Natl Acad. Sci. USA 105, 14 308–14 312. (doi:10.1073/pnas.0802430105)
    OpenUrlAbstract/FREE Full Text
  9. ↵
    1. Lam TT-Y et al.
    2013 The genesis and source of the H7N9 influenza viruses causing human infections in China. Nature 502, 241–244. (doi:10.1038/nature12515)
    OpenUrlCrossRefPubMedWeb of Science
  10. ↵
    1. Menachery VD et al.
    2015 A SARS-like cluster of circulating bat coronaviruses shows potential for human emergence. Nat. Med. 21, 1508–1513. (doi:10.1038/nm.3985)
    OpenUrlCrossRefPubMed
  11. ↵
    1. Rohani P,
    2. Drake JM
    . 2011 The decline and resurgence of pertussis in the US. Epidemics 3, 183–188. (doi:10.1016/j.epidem.2011.10.001)
    OpenUrlCrossRefPubMedWeb of Science
  12. ↵
    1. Altizer S,
    2. Ostfeld RS,
    3. Johnson PT,
    4. Kutz S,
    5. Harvell CD
    . 2013 Climate change and infectious diseases: from evidence to a predictive framework. Science 341, 514–519. (doi:10.1126/science.1239401)
    OpenUrlAbstract/FREE Full Text
  13. ↵
    1. O'Regan SM,
    2. Drake JM
    . 2013 Theory of early warning signals of disease emergence and leading indicators of elimination. Theor. Ecol. 6, 333–357. (doi:10.1007/s12080-013-0185-5)
    OpenUrlCrossRef
  14. ↵
    1. Rimoin AW et al.
    2010 Major increase in human monkeypox incidence 30 years after smallpox vaccination campaigns cease in the Democratic Republic of Congo. Proc. Natl Acad. Sci. USA 107, 16 262–16 267. (doi:10.1073/pnas.1005769107)
    OpenUrlAbstract/FREE Full Text
  15. ↵
    1. Jansen VA,
    2. Stollenwerk N,
    3. Jensen HJ,
    4. Ramsay M,
    5. Edmunds W,
    6. Rhodes C
    . 2003 Measles outbreaks in a population with declining vaccine uptake. Science 301, 804–804. (doi:10.1126/science.1086726)
    OpenUrlFREE Full Text
  16. ↵
    1. Grenfell BT,
    2. Pybus OG,
    3. Gog JR,
    4. Wood JL,
    5. Daly JM,
    6. Mumford JA,
    7. Holmes EC
    . 2004 Unifying the epidemiological and evolutionary dynamics of pathogens. Science 303, 327–332. (doi:10.1126/science.1090727)
    OpenUrlAbstract/FREE Full Text
  17. ↵
    1. Mooi FR,
    2. van der Maas NAT,
    3. De Melker HE
    . 2013 Pertussis resurgence: waning immunity and pathogen adaptation—two sides of the same coin. Epidemiol. Infect. 142, 685–694. (doi:10.1017/S0950268813000071)
    OpenUrlWeb of Science
  18. ↵
    1. Bart MJ et al.
    2014 Global population structure and evolution of Bordetella pertussis and their relationship with vaccination. mBio 5, e01074. (doi:10.1128/mBio.01074-14)
    OpenUrlCrossRefPubMed
  19. ↵
    1. Riolo MA,
    2. King AA,
    3. Rohani P
    . 2013 Can vaccine legacy explain the British pertussis resurgence? Vaccine 31, 5903–5908. (doi:10.1016/j.vaccine.2013.09.020)
    OpenUrlCrossRefPubMedWeb of Science
  20. ↵
    1. Keeling MJ,
    2. Rohani P
    . 2008 Modeling infectious diseases in humans and animals. Princeton, NJ: Princeton University Press.
  21. ↵
    1. Lloyd-Smith JO,
    2. Schreiber SJ,
    3. Kopp PE,
    4. Getz W
    . 2005 Superspreading and the effect of individual variation on disease emergence. Nature 438, 355–359. (doi:10.1038/nature04153)
    OpenUrlCrossRefPubMedWeb of Science
  22. ↵
    1. Anderson RM,
    2. May RM
    . 1991 Infectious diseases of humans: dynamics and control. Oxford, UK: Oxford University Press.
  23. ↵
    1. Bailey NTJ
    . 1990 The elements of stochastic processes with applications to the natural sciences, vol. 25. New York, NY: John Wiley & Sons.
  24. ↵
    1. Boettiger C,
    2. Hastings A
    . 2012 Quantifying limits to detection of early warning for critical transitions. J. R. Soc. Interface 9, 2527–2539. (doi:10.1098/rsif.2012.0125)
    OpenUrlAbstract/FREE Full Text
  25. ↵
    1. Gibson MA,
    2. Bruck J
    . 2000 Efficient exact stochastic simulation of chemical systems with many species and many channels. J. Phys. Chem. A 104, 1876–1889. (doi:10.1021/jp993732q)
    OpenUrlCrossRefWeb of Science
  26. ↵
    1. Anderson DF
    . 2007 A modified next reaction method for simulating chemical systems with time dependent propensities and delays. J. Chem. Phys. 127, 214107. (doi:10.1063/1.2799998)
    OpenUrlCrossRefPubMed
  27. ↵
    1. O'Regan SM,
    2. Lillie JW,
    3. Drake JM
    . 2016 Leading indicators of mosquito-borne disease elimination. Theor. Ecol. 9, 269–286. (doi:10.1007/s12080-015-0285-5)
    OpenUrl
  28. ↵
    1. van Kampen NG
    . 1981 Stochastic processes in physics and chemistry. Amsterdam, The Netherlands: North Holland.
  29. ↵
    1. Black AJ,
    2. McKane AJ
    . 2012 Stochastic formulation of ecological models and their applications. Trends Ecol. Evol. 27, 337–345. (doi:10.1016/j.tree.2012.01.014)
    OpenUrlCrossRefPubMedWeb of Science
  30. ↵
    1. Harris TE
    . 1963 The theory of branching processes. New York, NY: Dover Publications.
  31. ↵
    1. Singh S,
    2. Myers CR
    . 2014 Outbreak statistics and scaling laws for externally driven epidemics. Phys. Rev. E 89, 042108. (doi:10.1103/PhysRevE.89.042108)
  32. ↵
    1. Heesterbeek J,
    2. Dietz K
    . 1996 The concept of R0 in epidemic theory. Stat. Neerl. 50, 89–110. (doi:10.1111/j.1467-9574.1996.tb01482.x)
    OpenUrlCrossRef
  33. ↵
    1. Cover TM,
    2. Thomas JA
    . 1991 Elements of information theory. New York, NY: Wiley-Interscience.
  34. ↵
    1. Kaspar F,
    2. Schuster H
    . 1987 Easily calculable measure for the complexity of spatiotemporal patterns. Phys. Rev. A 36, 842–848. (doi:10.1103/PhysRevA.36.842)
    OpenUrlCrossRefPubMedWeb of Science
  35. ↵
    1. Rohani P,
    2. Miramontes O
    . 1995 Host-parasitoid metapopulations: the consequences of parasitoid aggregation on spatial dynamics and searching efficiency. Proc. R. Soc. Lond. B 260, 335–342. (doi:10.1098/rspb.1995.0101)
    OpenUrlAbstract/FREE Full Text
  36. ↵
    1. Ashcroft P,
    2. Galla T
    . 2013 Pattern formation in individual-based systems with time-varying parameters. Phys. Rev. E 88, 062104. (doi:10.1103/PhysRevE.88.062104)
    OpenUrl
  37. ↵
    1. Dibble CJ,
    2. O'Dea EB,
    3. Park AW,
    4. Drake JM
    . 2016 Waiting time to infectious disease emergence. J. R. Soc. Interface 13, 20160540. (doi:10.1098/rsif.2016.0540)
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Blumberg S,
    2. Lloyd-Smith JO
    . 2013 Inference of R0 and transmission heterogeneity from the size distribution of stuttering chains. PLoS Comput. Biol. 9, e1002993. (doi:10.1371/journal.pcbi.1002993)
    OpenUrlCrossRefPubMed
  39. ↵
    1. Cross PC,
    2. Johnson PL,
    3. Lloyd-Smith JO,
    4. Getz WM
    . 2007 Utility of R0 as a predictor of disease invasion in structured populations. J. R. Soc. Interface 4, 315–324. (doi:10.1098/rsif.2006.0185)
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Cowling BJ,
    2. Park M,
    3. Fang VJ,
    4. Wu P,
    5. Leung GM,
    6. Wu JT
    . 2015 Preliminary epidemiologic assessment of MERS-CoV outbreak in South Korea, May–June 2015. Euro Surveill. 20, 7–13.
    OpenUrlPubMed
  41. ↵
    1. He D,
    2. Ionides EL,
    3. King AA
    . 2009 Plug-and-play inference for disease dynamics: measles in large and small populations as a case study. J. R. Soc. Interface 7, 271–283. (doi:10.1098/rsif.2009.0151)
    OpenUrl
  42. ↵
    1. Brett TS,
    2. Drake JM,
    3. Rohani P
    . 2017 Data from: Anticipating the emergence of infectious diseases. Dryad Digital Repository. (doi:10.5061/dryad.7kt3h)
View Abstract
PreviousNext
Back to top
PreviousNext
July 2017
Volume 14
, issue 132
Journal of The Royal Society Interface: 14 (132)
  • Table of Contents
  • About the Cover
  • Index by author
  • Ed Board (PDF)

Keywords

modelling of infectious disease
epidemiology
complex systems
early-warning signals
Share
Anticipating the emergence of infectious diseases
Tobias S. Brett, John M. Drake, Pejman Rohani
J. R. Soc. Interface 2017 14 20170115; DOI: 10.1098/rsif.2017.0115. Published 5 July 2017
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Connotea logo Facebook logo Google logo Mendeley logo
Email

Thank you for your interest in spreading the word on Journal of The Royal Society Interface.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Anticipating the emergence of infectious diseases
(Your Name) has sent you a message from Journal of The Royal Society Interface
(Your Name) thought you would like to see the Journal of The Royal Society Interface web site.
Print
Manage alerts

Please log in to add an alert for this article.

Sign In to Email Alerts with your Email Address
Citation tools
Research article:

Anticipating the emergence of infectious diseases

Tobias S. Brett, John M. Drake, Pejman Rohani
J. R. Soc. Interface 2017 14 20170115; DOI: 10.1098/rsif.2017.0115. Published 5 July 2017

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Download

Article reuse

  • Article
    • Abstract
    • 1. Introduction
    • 2. Early-warning signals using the birth–death–immigration process
    • 3. Comparison with models featuring susceptible depletion
    • 4. Cox's δ: a likelihood-ratio test
    • 5. Impact of time scale on reliability of predictors
    • 6. Discussion/conclusion
    • Data accessibility
    • Author's contributions
    • Competing interests
    • Funding
    • Appendix A. Critical slowing down for the stationary birth–death–immigration process
    • Appendix B. Shannon entropy and Kolmogorov complexity
    • Appendix C. Transition probability distribution and likelihood
    • Appendix D. Calculating early-warning signals from single time-series data
    • Appendix E. Estimating the probability of supercriticality
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

See related subject areas:

  • computational biology
  • biocomplexity

Related articles

Cited by

Interface Focus To read similar articles, check out our sister journal

Open biology

  • Interface
    • About this journal
    • Contact information
    • Purchasing information
    • Submit
    • Author benefits
    • Open access membership
    • Recommend to your library
    • FAQ
    • Help

Royal society publishing

  • ROYAL SOCIETY PUBLISHING
    • Our journals
    • Open access
    • Publishing policies
    • Conferences
    • Podcasts
    • News
    • Blog
    • Manage your account
    • Terms & conditions
    • Cookies

The royal society

  • THE ROYAL SOCIETY
    • About us
    • Contact us
    • Fellows
    • Events
    • Grants, schemes & awards
    • Topics & policy
    • Collections
    • Venue hire
1742-5662

Copyright © 2018 The Royal Society