Anticipating epidemic transitions with imperfect data
There are numerous causative factors linked with disease emergence , including pathogen evolution , ecological change and variation in host demography and behavior 1–5 ., Combined , they can make each pathogen’s emergence seem idiosyncratic ., In spite of this apparent particularity , there is a recent literature on the possibility of anticipating epidemic transitions using model-independent metrics 6–14 ., Referred to as early-warning signals ( EWS ) , these metrics are summary statistics ( e . g . the variance and autocorrelation ) which undergo characteristic changes as the transition is approached ., In addition to infectious disease transmission , EWS have been investigated for transitions in a broad range of dynamical systems , including ecosystem collapse and climate change 15–21 ., The motivation for EWS comes from the theories of dynamical systems and stochastic processes , in particular the slowing down that universally occurs in the vicinity of dynamical critical points 22–24 ., Theoretical results for disease emergence are promising , and suggest that the transition from limited stuttering chains of transmission ( R0 < 1 ) to sustained transmission and outbreaks ( R0 > 1 ) is preceded by detectable EWS 8 , 13 , 14 ., A major obstacle to deploying early-warning systems is the type of data available to calculate the EWS ., Theoretical predictions assume the data will be sequential recordings ( or “snapshots” ) of the true number of infected in the population through time 8–13 ., In this paper we refer to this as snapshots data ., However , epidemiological data originate instead from notifications by public health practitioners whenever a new case is identified ., Public health bodies aggregate individual cases into regular case reports ( e . g . the US Centers for Disease Control and Prevention’s Morbidity and Mortality Weekly Report ) , as shown in Fig 1 ., Different combinations of serial interval ( difference in time of symptom onset between primary and secondary cases ) and aggregation period lead to time series which have very different appearances ., Even assuming perfect reporting , variability in both the incubation period and onset of clinical symptoms mean that snapshots data cannot be reconstructed from case report data ., In addition to aggregation , case reports are subject to reporting error ( see Fig 2 ) ., Underreporting may occur due to asymptomatic infection , poorly implemented notification protocols , or socio-political factors 25–29 ., Misdiagnoses and clerical errors in the compilation of reports can result in both under- and over-reporting 30–32 ., Due to self reporting and contact tracing , once an index case has been positively identified secondary cases are more likely to be diagnosed , which may lead to clustering in case reports ., The combination of case aggregation and reporting error results in a mismatch between snapshots and imperfect epidemiological data ., EWS , such as the variance ( Fig 3 , top panel ) , are affected by imperfect data ( Fig 3 , bottom panel ) and may not display the characteristic trends that form the basis for detecting disease emergence ., This provides reason to question the direct application of EWS to observed data ., In this paper we report on a simulation study aimed at investigating the robustness of a range of EWS to case report data ., We simulated a stochastic SIR model of a pathogen emerging via increasing R0 , and corrupted the simulated case reports by applying a negative binomial reporting error ., The area under the curve ( AUC ) statistic was computed to quantify how well trends in an EWS identify emergence ., We find that performance depends on both the EWS and the reporting model ., Broadly , the mean , variance , index of dispersion and first differenced variance perform well ., The autocorrelation , autocovariance and decay time perform well unless either, i ) the data are highly overdispersed or, ii ) the aggregation period is less than the infectious period ., The coefficient of variation , kurtosis , and skewness have a more subtle dependence on the reporting model , and are not reliable ., We conclude that seven of ten EWS perform well for most realistic reporting scenarios ., The dynamics of disease spread in a host population are modeled as a stochastic process using an SIR model with birth and death 37 ., The model compartments and parameters are listed in Table 1 . Transition rates and effects are listed in Table 2 . The basic reproductive number for the SIR model is R0 ( t ) = β ( t ) / ( γ+ α ) , where β ( t ) varies due to nondescript secular trends in the transmissibility ., Simulated data are generated using the Gillespie algorithm 33 , which simulates a sequence of transition events ( infection , recovery , birth and death ) , and returns the number of individuals in each model compartment through time ., The SIR simulations are of a population with average size N0 = 106 ., The parameter ζ gives the rate at which new cases arise due to external sources , and is set to ζ = 1 per week ., The death rate , α , is the reciprocal of the life expectancy , set to 70 years ., Case counts , Ct , are given by the number of recovery events ( at rate γIt ) within each aggregation period , and are included in the model as an additional variable ( see Table 1 ) ., Reporting error is applied to the case count at the end of each aggregation period by sampling a negative binomial distribution ,, P ( K t = k | C t ) = Γ ( ϕ + k ) k !, Γ ( ϕ ) ( ρ C t ρ C t + ϕ ) k ( ϕ ρ C t + ϕ ) ϕ , ( 1 ), with reporting probability ρ and dispersion parameter ϕ 38 ., Given Ct cases , the mean number reported is μt = ρCt ., The variance is specified by the dispersion parameter via the relation σ t 2 = μ t + μ t 2 / ϕ ., Increasing ϕ reduces the overdispersion of the data , so that for large ϕ the distribution of reports is approximately Poisson ., Previous work has proposed a range of different EWS to anticipate dynamical transitions 8 , 12–15 , 17 , 18 ., The ten candidate EWS considered in this paper are listed in Table 3 . We consider additional indicators to those most frequently studied in the EWS literature ( the variance , autocorrelation and coefficient of variation ) ., As R0 approaches 1 , the mean number of cases caused by introductions rises , making it a potential EWS ., The index of dispersion is a similar measure to the coefficient of variation , and is defined as the variance to mean ratio ., The decay time ( or correlation time ) is a log-transform of the autocorrelation , which diverges as R0 approaches 1 ( the definition of critical slowing down ) ., In addition to the autocorrelation , which is normalized by the variance , we consider the unnormalized autocovariance ., As both the autocorrelation and variance increase , the autocovariance may outperform these two measures ., Theoretical results show the increase in variance accelerates as R0 approaches 1 , suggesting the first differenced variance as a complementary EWS ., Additionally we investigate the performance of two higher-order moments , the skewness and kurtosis ., Functional expressions for the dependence of each EWS on R0 can be found using the Birth-Death-Immigration ( BDI ) process , a variation of the SIR model which neglects susceptible depletion ( i . e . St = N0 ) ., The BDI process is a one-dimensional stochastic process , depending only on the number of infected It , and possesses an exact mathematical solution ( for full details see 13 ) ., This allows expressions for the moments and correlation functions of It to be found ( Table 3 , fourth column ) ., BDI theory predicts that most EWS ( the mean , variance , index of dispersion , autocovariance , decay time and first differenced variance ) are expected to grow hyperbolically as R0 approaches one ., The autocorrelation is expected to grow exponentially , the kurtosis quadratically and the skewness linearly ., The coefficient of variation is the only EWS which does not grow , instead remaining constant ., We propose observing these trends in data as a basis for anticipating disease emergence ., The numerical estimators used in this paper are listed in Table 3 , third column , discussed in more depth below ., Theoretical predictions from the BDI process are based on It and do not take into account effects of reporting error and aggregation ., The focus of this paper is to examine the robustness of each EWS to reporting process parameters , using simulated case report data , Kt ., BDI theory predicts that 9 out of 10 EWS increase as the transition is approached ., We quantify the association of each EWS with time using Kendall’s rank correlation coefficient 19 ., A coefficient close to ( +/− ) 1 implies consistent increases/decreases of the EWS in time ., As the underlying dynamics of the case reports are stochastic , the value of the rank correlation coefficient is itself a random variable ., Multiple simulations of the test ( emerging ) and null ( stationary/not emerging ) scenarios result in two distributions of correlation coefficients for each EWS ., We measure performance using the AUC statistic , defined as the overlap of the two distributions , and may be interpreted as the probability that a randomly chosen test coefficient is higher than a randomly chosen null coefficient , AUC = P ( τtest > τnull ) 39 , 40 ., The name comes from one method of calculating it , the area under the receiver operating characteristic ( ROC ) curve , a parametric plot of the false positive rate against true positive rate as the decision threshold is varied 41 ., Instead of explicitly calculating the ROC curve , the AUC can be efficiently calculated after ranking the combined set of test and null correlation coefficients by value 40 ,, AUC =  r test - n test ( n test + 1 ) / 2  / ( n test n null ) , ( 2 ), where rtest is the sum of the ranks of test coefficients and ntest and nnull are the number of realizations of the test and null models respectively ., In this paper the AUC statistic quantifies how successfully an EWS distinguishes whether or not a disease is approaching an epidemic transition ., An AUC = 0 . 5 implies that an observed rank coefficient value conveys no information about whether or not the disease is emerging , i . e . the EWS is ineffective ., If the AUC < 0 . 5 then a decreasing trend in the EWS indicates emergence , whereas if AUC > 0 . 5 an increasing trend indicates emergence ., A larger |AUC − 0 . 5| implies better performance , if |AUC − 0 . 5| = 0 . 5 the rank coefficient value classifies the two scenarios perfectly ., The mathematical definitions of the EWS depend on expectations of the stochastic process , Ef ( X )  ( Table 3 , second column ) ., To calculate EWS from non-stationary time series data we use centered moving window averages with bandwidth b as estimators for expectation values ., For example , the mean at time t is estimated using, μ ^ t = ∑ s = t - ( b - 1 ) δ t + ( b - 1 ) δ X s 2 b - 1 , ( 3 ), where δ is the size of one time step ., Near the ends of the time series ( t < bδ and t > T − bδ ) , the normalization factor 2b − 1 is reduced to ensure it remains equal to the number of data points within the window ., Applying Eq 3 to the time series for X results in a time series for μ ^ ., Certain EWS depend on others , for example the variance depends on the mean ., EWS are therefore calculated iteratively , for example μ ^ is first calculated using Eq 3 , and then σ ^ 2 is found using, σ ^ t 2 = ∑ s = t - ( b - 1 ) δ t + ( b - 1 ) δ ( X s - μ ^ s ) 2 2 b - 1 ., ( 4 ), Estimators for each EWS are in Table 3 . For snapshots data Xt = It , and for case report data Xt = Kt ., Throughout this paper we use a bandwidth of b = 35 time steps ( weeks or months depending on aggregation period ) ., Results have been found to be similar for a bandwidth of b = 100 time steps ., To quantify the sensitivity of each EWS to reporting process , we calculate the AUC from simulated data for a range of different model parameter combinations ., The experimental design is fully factorial ( i . e . considers all parameter value combinations ) ., The following four parameters are varied:, ( i ) the infectious period , 1/γ , which can be either 7 or 30 days ,, ( ii ) the reporting probability , ρ = 2−8x for x in {0 , 0 . 05 , 0 . 1 , … , 1} ,, ( iii ) the dispersion parameter , ϕ , which is one of {0 . 01 , 1 , 100} ,, ( iv ) the aggregation period , δ , which is either weekly or monthly ., For the test model , the disease emerges over T = 20 years , via an increase R0 ., For the null model , R0 is constant ., One epidemiological interpretation for the test scenario is it models transmission in a population with high vaccine coverage , where gradual pathogen evolution results in increasing evasion of host immunity ., An alternative interpretation is it models zoonotic spillover , where pathogen evolution within an animal reservoir results in gradually increasing human transmissibility 42 ., In both interpretations , the null model assumes no change in host-to-host transmissibility ., The transmission dynamics were simulated using the Gillespie algorithm 33 ., The Gillespie algorithm assumes all model parameters ( including the transmissibility ) are constant ., To simulate disease emergence we modify the Gillespie algorithm , discretely increasing β at the end of each day and after each reaction to ensure an approximately linear increase in R0 over T = 20 years , from R0 ( 0 ) = 0 to R0 ( T ) = 1 . For the null model , transmission is simulated for 20 years at a constant rate , R0 = 0 ., Our choice of null has no secondary transmission , making the classification problem easy under perfect reporting ., This enables clearer identification of responses to reporting process effects as results span the full range of the AUC statistic ., We repeated the experiment with null model R0 = 0 . 5 , and found no qualitative differences ., For both scenarios transmission is subcritical , with disease presence maintained by reintroduction from an external reservoir ., For each parameter combination 1000 replicates of both scenarios are generated ., We perform these computational experiments in R using the pomp package 43 to simulate the SIR model and the spaero package 44 to calculate the EWS ., Code was written to simulate aggregation and reporting error ., All code to reproduce the results is archived online at doi:10 . 5281/zenodo . 1185284 ., Provided the data are aggregated monthly , with high reporting probability and low overdispersion , the coefficient of variation , skewness and kurtosis have similar AUC values when calculated from snapshots data ( Fig, 4 ) and case report data ( Fig 5 , right column ) ., Unlike the other seven EWS , this it is not the case for weekly data ., If calculated from weekly snapshots data with 1/γ = 1 week , the coefficient of variation has an AUC = 0 . 18 ( Fig 4 , bottom right ) ., With reporting , if ρ = 1 , ϕ = 100 the AUC = 0 . 005 ( Fig 6 , top right ) ., By switching to case report data the performance of the coefficient of variation has improved dramatically ., Similar improvements are seen for the skewness and kurtosis ., In addition , and perhaps counterintuitively , these three EWS’s performances are further enhanced at lower reporting probabilities ( compare the top right and bottom right panels of Fig 6 ) ., At low overdispersion and low reporting probability , the coefficient of variation ( |AUC − 0 . 5| = 0 . 5 ) is joint with the mean and variance as the best performing statistic , closely followed by the skewness ( |AUC − 0 . 5| = 0 . 497 ) and kurtosis ( |AUC − 0 . 5| = 0 . 491 ) ., The improvement in performance at low reporting probability is acutely sensitive to other model parameters ., Both overdispersion in the reporting ( for example Fig 6 , left column ) and larger aggregation period ( Fig 5 , right column ) severely dampen the sensitivity to ρ ., All three EWS perform poorly if ϕ = 0 . 01 , regardless of the other model parameters ., This group of EWS are all measures of the correlation between neighboring data points ., At high reporting probability ( ρ > 0 . 33 ) and low overdispersion ( ϕ = 100 ) , all three perform well ( AUC > 0 . 77 ) , regardless of infectious and aggregation periods ( see Fig 5 ) ., Performance is comparable with snapshots data ( Fig 4 ) ., Overall , they perform best if 1/γ = 1 week ( Fig 5 , top row ) and worst if 1/γ > δ ( Fig 5 , bottom left ) ., At low overdispersion , decreasing the reporting probability reduces the AUC ( compare the top right and bottom right panels of Fig 6 , AUC = 1 . 000 vs 0 . 831 ) ., The performance drop is largest if 1/γ = 1 month and δ = 1 week ., The performance of all three EWS is negatively affected by overdispersion ., Sensitivity to overdispersion is least for 1/γ = δ = 1 week , performance is only poor if ϕ = 0 . 01 and/or ρ ≲ 0 . 036 ( Fig 5 , top left ) ., These three EWS are reliable indicators of emergence provided δ ≥ 1/γ and ϕ = 100 ., Unless reporting error is highly overdispersed ( ϕ = 0 . 01 ) , the mean , variance and first differenced variance perform extremely well ( AUC ≈ 1 , see Fig 5 ) ., If case reports are aggregated weekly and have high overdispersion ( ϕ = 0 . 01 ) , they are among the best performing EWS ., The mean and variance have AUC > 0 . 85 , and the first differenced variance has AUC ≈ 0 . 66 , but is largely unaffected by reporting probability and infectious period ., However , if case reports are aggregated monthly and ϕ = 0 . 01 , then all three perform poorly ., This holds regardless of reporting probability and infectious period , and is in line with the results for other EWS ., The index of dispersion ( unrelated to the dispersion parameter ) has a similar performance to the previous group of EWS , however with certain differences ., We first consider low overdispersion ( ϕ = 100 ) ., At low reporting probabilities the index of dispersion performs best if 1/γ = 1 week and δ = 1 month ( Fig 5 , top right ) ., For other combinations of infectious period and aggregation period , performance suffers a sharp drop as reporting probability decreases ., This drop occurs at a reporting probability dependent on the infectious period and aggregation period , around ρ = 0 . 047 for δ = 1 week , and around ρ = 0 . 027 for δ = 1/γ = 1 month ., Unique among the EWS , the index of dispersion performs best at intermediate overdispersion ( ϕ = 1 ) , in particular at small reporting probability ., This is true for all infectious and aggregation periods , although most pronounced if 1/γ = 1 month and δ = 1 week ., For ϕ = 0 . 01 the index of dispersion performs better if the data are aggregated weekly , and best if the infectious period is also one week , with AUC ≈ 0 . 71 for ρ = 0 . 047 ( Fig 6 , bottom left ) ., Provided ρ ≳ 0 . 05 and ϕ > 0 . 01 , performance is good for all aggregation and infectious periods ., Overall performance is best if 1/γ = 1 week and δ = 1 month ., Taken in isolation , the mean and variance are the EWS least impacted by reporting ., Unless the overdispersion in the observation process is high ( ϕ = 0 . 01 ) , their performance is largely unaffected by reporting process parameters ., At low reporting probabilities they outperform the autocorrelation , autocovariance , decay time and index of dispersion , and are independent of aggregation period and infectious period ., EWS sensitive to correlation between neighboring data points perform well unless, i ) ϕ = 0 . 01 and/or, ii ) 1/γ > δ and ρ ≲ 0 . 06 ., While it is clear how high overdispersion in reporting reduces correlation in the data , an explanation for, ii ) is less clear ., If calculated from snapshots data , the coefficient of variation , kurtosis and skewness are the worst performing statistics ( |AUC − 0 . 5| ≈ 0 ) ., Using case report data improves performance under certain conditions ., If cases are aggregated weekly with low reporting probability and low overdispersion then they are among the best performing EWS , with |AUC − 0 . 5| ≈ 0 . 5 ., In addition the trends of the skewness and kurtosis ( both decreasing ) are opposite those given by the BDI process ( both increasing ) ., Overall , we conclude that these three EWS are unreliable indicators of disease emergence as their performance is conditional on a limited range of reporting process parameters ., For mathematical reasons , proposed EWS for disease emergence have assumed access to regular recordings ( “snapshots” ) of the entire infectious population 8–13 ., However , epidemiological data are typically aggregated into periodic case reports subject to reporting error ., To examine the practical consequences of this mismatch between theory and data , in this paper we calculated EWS from case report data ., We performed extensive numerical simulations to determine the sensitivity of each candidate EWS to imperfect data ., Case aggregation and reporting error change the statistical properties of the data , and can have subtle effects on an EWS’s performance ., We identified four groups of EWS based on their sensitivity to the various reporting process parameters ., The performance of one group , consisting of the EWS with either polynomial or no growth with R0 , has a nuanced relationship with the reporting process parameters ., We therefore conclude that the coefficient of variation , kurtosis and skewness perform poorly as EWS ., In general , the other EWS ( the mean , variance , first differenced variance , index of dispersion , autocorrelation , autocovariance and decay time ) all performed well and are strong candidates for incorporation in monitoring systems intended to provide early warning of disease emergence ., Surprisingly , the combination of reporting error and aggregation of data does not always have a detrimental effect on EWS performance ., The coefficient of variation , kurtosis and skewness perform best when both reporting probability and overdispersion are low ., At first glance this result appears counter-intuitive: as an increasingly large fraction of cases are missed , performance improves ., The point to stress here is that by changing the parameters of the reporting process we are systematically changing the statistical properties of the time series ., For instance , the BDI process predicts no trend in the coefficient of variation , due to the standard deviation and mean increasing with R0 at an identical rate 13 ., With aggregation and reporting error this identity does not necessarily hold , introducing a trend in the coefficient of variation and improving its performance ., To fully explain this phenomenon requires an analytical solution for the statistics of Kt , which requires solving the stochastic process including aggregation and reporting error ., However , it can be seen to be plausible if we focus only on stochasticity resulting from reporting error ., For low overdispersion ( e . g . ϕ = 100 ) , the reporting probability distribution can be approximated by a Poisson distribution with parameter λ = ρCt ., Ignoring demographic stochasticity , we replace Ct with ECt = ηδ ( 1 − R0 ) −1 ., Both the coefficient of variation and skewness for this distribution are λ−1/2 = { ( 1−R0 ) /ρηδ}1/2 and the kurtosis is λ−1 = ( ρηδ ) −1 ( 1 − R0 ) ., These two expressions both decrease as R0 increases from 0 to 1 , consistent with the experimentally observed AUC < 0 . 5 ., The improved performance at low ρ is a consequence of the increased stochasticity in reporting outweighing demographic stochasticity ., Can these EWS be used to anticipate disease emergence ?, If overdispersion and reporting probability are known to be low , then yes ., However , it is unlikely that the reporting process is sufficiently understood for an emerging disease ., We conclude that these three EWS are unreliable and therefore not good indicators of emergence ., There is a similar reason for why the index of dispersion has a peak in performance at intermediate reporting overdispersion ., The negative binomial reporting distribution , conditioned on ECt as above , has index of dispersion given by σ2/μ = 1 + ρηδ{ϕ ( 1−R0 ) }−1 ., Therefore increasing reporting overdispersion ( i . e increasing ϕ−1 ) amplifies the response of σ2/μ to changes in R0 ., This leads to a greater differential , improving performance of the index of dispersion as an EWS ., However , increased reporting overdispersion also implies increased volatility of data within a finite sized window , which reduces reliability ., These two countervailing factors provide an explanation for the optimal performance at intermediate overdispersion values ., In our analysis we considered an SIR model with epidemiologically plausible parameters ., The negative binomial distribution is meant to provide a stringent test on EWS performance , and the parameter ranges are conservative ( especially for overdispersion ) ., For instance , if there are 10 actual cases in a week , and reporting error is negative binomially distributed with ρ = 0 . 1 and ϕ = 0 . 01 , then the mean number of reported cases is 1 ., However , the probability of no cases being reported is P ( K = 0 ) = 0 . 955 whereas the variance in reported cases is σ2 = 101 ., The resulting time series is highly volatile , with little similarity in appearance to the underlying time series of actual cases ., It is unlikely that case reports for an emerging disease will have such high overdispersion ., In addition , for a highly pathogenic emerging disease , such as Middle East respiratory syndrome ( MERS ) or H7N9 avian influenza , the reporting probability is likely much higher than ρ = 1/256 ( the smallest value we studied ) ., Nonetheless , one of the encouraging findings of this study is that high reporting is not essential for reliable early warning ., Clear trends in the EWS can still be identified , provided there are sufficiently many introductions for cases to be sporadically detected prior to emergence ., These dynamics are typical for a reemerging vaccine controllable disease , such as measles , where cases are continually introduced into disease-free regions from endemic regions 45 , 46 ., Performance of EWS which depend on correlation between neighboring case reports was found to be contingent on the aggregation period being larger than the serial interval ( equal to the infectious period for the SIR model ) ., If this is not the case , there is a smaller probability that successive links in a chain of transmission fall into neighboring case reports ., We speculate that this reduces the impact of fluctuations in a particular report on the subsequent report , diminishing their correlation ., This effect is exacerbated if the reporting probability is low ., A more rigorous explanation requires a full solution to the stochastic process with aggregation and reporting error ., For many known pathogens the serial interval is larger than one week , for example measles virus and Bordetella pertussis 47 ., For other pathogens it is less than one week , such as SARS coronavirus 48 and influenza virus 47 , 49 ., In order for the autocorrelation , autocovariance and decay time to be reliable EWS , our results suggest the data need to be aggregated by periods larger than the serial interval ., The performance boost outweighs costs associated with having fewer data points ., We expect that these three EWS will work best for pathogens with short serial intervals; for pathogens with extremely long serial intervals ( such as HIV ) reliable use of these EWS is unlikely ., The purpose of this study was not to identify the best EWS , but to investigate the robustness of this approach to the reporting process ., In order to isolate the effects of incomplete reporting and aggregation error , we ignored parameter uncertainty by fixing epidemiological parameters ( e . g . the infectious period and the introduction rate ) , rather than drawing them from a distribution ., As shown in Table 3 , both the mean and the variance scale with the introduction rate , which is a product of the per capita introduction rate and the susceptible population size ., On the other hand , the index of dispersion , autocorrelation and decay time are all independent of introduction rate ., Uncertainty of important factors , such as the susceptible population size , is a key challenge to anticipating emergence , and these three EWS may outperform the mean and the variance if uncertainty is included ., Thus , while the mean and the variance are most robust to imperfect data , they are not necessarily the best EWS ., Instead , our results suggest that imperfect data is not a barrier to the use of EWS ., One challenge to early-warning stems from the potential suddenness of novel pathogen emergence , for example SARS was unknown prior to the global outbreak in 2002-2003 ., For known pathogens , intermittent data availability presents a separate challenge ., Mumps was excluded from the US CDC’s MMWR in 2002 following a period of low incidence ., Subsequently , there was a series of large outbreaks , notably in 2006 in the Midwest , and mumps was reincluded ., Methods such as EWS are contingent on surveillance efforts being maintained ., In addition to underlining the importance of disease surveillance , our work suggests ways it can be improved ., Case reports sometimes include additional metadata , for example whether all suspected cases are counted or only clinically confirmed cases ., The reporting error of case reports with differing case identification criteria is expected to be very different , as has been seen for instance with MERS 50 ., This paper shows that EWS depend on the reporting process , and cross-validating EWS calculated from each data stream could improve performance ., Provided it is available , how to leverage metadata is a promising avenue for future research into enhancing EWS ., These results provide an essential stepping stone from previous theoretically focused works to implementable early-warning systems ., Our findings further reinforce the hypothesis that disease emergence is preceded by detectable EWS ., While epidemiological factors preclude early-warning for certain pathogens , for example Ebola virus ( estimates of R0 have consistently been greater than one 51 ) and HIV ( see above ) , they do not rule out many others , including reemerging childhood diseases 52 , H7N9 avian influenza virus 53 , and MERS coronavirus 54 ., These pathogens all present public health risk , and EWS may be able to play an important role in monitoring for their emergence .
Epidemic transitions are an important feature of infectious disease systems ., As the transmissibility of a pathogen increases , the dynamics of disease spread shifts from limited stuttering chains of transmission to potentially large scale outbreaks ., One proposed method to anticipate this transition are early-warning signals ( EWS ) , summary statistics which undergo characteristic changes as the transition is approached ., Although theoretically predicted , their mathematical basis does not take into account the nature of epidemiological data , which are typically aggregated into periodic case reports and subject to reporting error ., The viability of EWS for epidemic transitions therefore remains uncertain ., Here we demonstrate that most EWS can predict emergence even when calculated from imperfect data ., We quantify performance using the area under the curve ( AUC ) statistic , a measure of how well an EWS distinguishes between numerical simulations of an emerging disease and one which is stationary ., Values of the AUC statistic are compared across a range of different reporting scenarios ., We find that different EWS respond to imperfect data differently ., The mean , variance and first differenced variance all perform well unless reporting error is highly overdispersed ., The autocorrelation , autocovariance and decay time perform well provided that the aggregation period of the data is larger than the serial interval and reporting error is not highly overdispersed ., The coefficient of variation , skewness and kurtosis are found to be unreliable indicators of emergence ., Overall , we find that seven of ten EWS considered perform well for most realistic reporting scenarios ., We conclude that imperfect epidemiological data is not a barrier to using EWS for many potentially emerging diseases .
Anticipating disease emergence is a challenging problem , however the public health ramifications are clear ., A proposed tool to help meet this challenge are early-warning signals ( EWS ) , summary statistics which undergo characteristic changes before dynamical transitions ., While previous theoretical studies are promising , and find that epidemic transitions are preceded by detectable trends in EWS , they do not consider the effects of imperfect data ., To address this , we developed a simulation study which assesses how case aggregation and reporting error impact on 10 different EWS’s performance ., Case report data were simulated by combining a stochastic SIR transmission model with a model of reporting error ., Temporal trends in an EWS were used as a method of distinguishing between an emerging disease ( R0 approaching 1 ) and a stationary disease ( constant R0 ) ., We investigated the robustness of EWS to reporting process parameters , namely the aggregation period , reporting probability and overdispersion of reporting error ., Seven of ten EWS perform well for realistic reporting scenarios , and are strong candidates for incorporation in disease emergence monitoring systems .
medicine and health sciences, pathology and laboratory medicine, infectious disease epidemiology, pathogens, mathematical models, epidemiological methods and statistics, probability distribution, mathematics, stochastic processes, skewness, research and analysis methods, statistical distributions, infectious diseases, epidemiology, mathematical and statistical techniques, probability theory, statistical dispersion, epidemiological statistics, physical sciences