Comparing methods for estimating R0 from the size distribution of subcritical transmission chains.

Many diseases exhibit subcritical transmission (i.e. 0<R0<1) so that infections occur as self-limited 'stuttering chains'. Given an ensemble of stuttering chains, information about the number of cases in each chain can be used to infer R0, which is of crucial importance for monitoring the risk that a disease will emerge to establish endemic circulation. However, the challenge of imperfect case detection has led authors to adopt a variety of work-around measures when inferring R0, such as discarding data on isolated cases or aggregating intermediate-sized chains together. Each of these methods has the potential to introduce bias, but a quantitative comparison of these approaches has not been reported. By adapting a model based on a negative binomial offspring distribution that permits a variable degree of transmission heterogeneity, we present a unified analysis of existing R0 estimation methods. Simulation studies show that the degree of transmission heterogeneity, when improperly modeled, can significantly impact the bias of R0 estimation methods designed for imperfect observation. These studies also highlight the importance of isolated cases in assessing whether an estimation technique is consistent with observed data. Analysis of data from measles outbreaks shows that likelihood scores are highest for models that allow a flexible degree of transmission heterogeneity. Aggregating intermediate sized chains often has similar performance to analyzing a complete chain size distribution. However, truncating isolated cases is beneficial only when surveillance systems clearly favor full observation of large chains but not small chains. Meanwhile, if data on the type and proportion of cases that are unobserved were known, we demonstrate that maximum likelihood inference of R0 could be adjusted accordingly. This motivates the need for future empirical and theoretical work to quantify observation error and incorporate relevant mechanisms into stuttering chain models used to estimate transmission parameters.

[1]  Barney S. Graham,et al.  Major increase in human monkeypox incidence 30 years after smallpox vaccination campaigns cease in the Democratic Republic of Congo , 2010, Proceedings of the National Academy of Sciences.

[2]  D. Earn,et al.  Vaccination and the theory of games. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[3]  N. Arinaminpathy,et al.  Evolution and emergence of novel human infections , 2009, Proceedings of the Royal Society B: Biological Sciences.

[4]  C P Farrington,et al.  Branching process models for surveillance of infectious diseases controlled by mass vaccination. , 2003, Biostatistics.

[5]  M E Woolhouse,et al.  Risk factors for human disease emergence. , 2001, Philosophical transactions of the Royal Society of London. Series B, Biological sciences.

[6]  C P Farrington,et al.  Epidemiology of transmissible diseases after elimination. , 2000, American journal of epidemiology.

[7]  Daniel B Di Giulio,et al.  Human monkeypox. , 2004, The Lancet. Infectious diseases.

[8]  M. McCauley,et al.  Has surveillance been adequate to detect endemic measles in the United States? , 2004, The Journal of infectious diseases.

[9]  Arlene King,et al.  Measles elimination in Canada. , 2004, The Journal of infectious diseases.

[10]  C. Fraser,et al.  Public Health Risk from the Avian H5N1 Influenza Epidemic , 2004, Science.

[11]  E. Cohen.  Vaccine Refusal, Mandatory Immunization, and the Risks of Vaccine-Preventable Diseases , 2010 .

[12]  P. E. Kopp,et al.  Superspreading and the effect of individual variation on disease emergence , 2005, Nature.

[13]  James O. Lloyd-Smith,et al.  Inference of R 0 and Transmission Heterogeneity from the Size Distribution of Stuttering Chains , 2013, PLoS Comput. Biol..

[14]  M. Dwass.  The total progeny in a branching process and a related random walk , 1969, Journal of Applied Probability.

[15]  J. Lloyd-Smith.  Maximum Likelihood Estimation of the Negative Binomial Dispersion Parameter for Highly Overdispersed Data, with Applications to Infectious Diseases , 2007, PloS one.

[16]  Daimin Shi,et al.  A new compounding life distribution: the Weibull–Poisson distribution , 2012 .

[17]  Mark E.J. Woolhouse,et al.  Host Range and Emerging and Reemerging Pathogens , 2005, Emerging infectious diseases.

[18]  Nathan D. Wolfe,et al.  Origins of major human infectious diseases , 2007, Nature.

[19]  N. Becker,et al.  Estimation for discrete time branching processes with application to epidemics. , 1977, Biometrics.

[20]  N. Becker,et al.  On parametric estimation for mortal branching processes , 1974 .

[21]  Frank Ball,et al.  Empty confidence sets for epidemics, branching processes and Brownian motion , 2002 .

[22]  Rustom Antia,et al.  The role of evolution in the emergence of infectious diseases , 2003, Nature.

[23]  N. Stollenwerk,et al.  Measles Outbreaks in a Population with Declining Vaccine Uptake , 2003, Science.

[24]  T. E. Harris,et al.  The Theory of Branching Processes. , 1963 .

[25]  R Bruce Aylward,et al.  The effect of mass immunisation campaigns and new oral poliovirus vaccines on the incidence of poliomyelitis in Pakistan and Afghanistan, 2001–11: a retrospective analysis , 2012, The Lancet.

[26]  T. E. Harris,et al.  The Theory of Branching Processes. , 1963 .

[27]  Kenneth Lange,et al.  Applied Probability , 2003 .

[28]  James O. Lloyd-Smith,et al.  Epidemic Dynamics at the Human-Animal Interface , 2009, Science.

[29]  C Paddy Farrington,et al.  Assessment of the status of measles elimination from reported outbreaks: United States, 1997-1999. , 2004, The Journal of infectious diseases.

[30]  L. Prior,et al.  Ecological Models and Data in R , 2011 .

[31]  C. J. Rhodes,et al.  The effect of superspreading on epidemic outbreak size distributions. , 2008, Journal of theoretical biology.

[32]  James O. Lloyd-Smith,et al.  Identifying genetic markers of adaptation for surveillance of viral host jumps , 2010, Nature Reviews Microbiology.

[33]  C. Mode,et al.  Estimating the transmission potential of supercritical processes based on the final size distribution of minor outbreaks , 2011, Journal of Theoretical Biology.