Modelling and Bayesian analysis of the Abakaliki smallpox data.

The celebrated Abakaliki smallpox data have appeared numerous times in the epidemic modelling literature, but in almost all cases only a specific subset of the data is considered. The only previous analysis of the full data set relied on approximation methods to derive a likelihood and did not assess model adequacy. The data themselves continue to be of interest due to concerns about the possible re-emergence of smallpox as a bioterrorism weapon. We present the first full Bayesian statistical analysis using data-augmentation Markov chain Monte Carlo methods which avoid the need for likelihood approximations and which yield a wider range of results than previous analyses. We also carry out model assessment using simulation-based methods. Our findings suggest that the outbreak was largely driven by the interaction structure of the population, and that the introduction of control measures was not the sole reason for the end of the epidemic. We also obtain quantitative estimates of key quantities including reproduction numbers.

[1]  W. Foege,et al.  Faith Tabernacle smallpox epidemic, Abakaliki, Nigeria , 1968 .

[2]  N G Becker,et al.  Inference for an epidemic when susceptibility varies. , 2001, Biostatistics.

[3]  Youssef M. Marzouk,et al.  A Bayesian method for inferring transmission chains in a partially observed epidemic. , 2008 .

[4]  Xiao-Li Meng,et al.  POSTERIOR PREDICTIVE ASSESSMENT OF MODEL FITNESS VIA REALIZED DISCREPANCIES , 1996 .

[5]  Martin Eichner,et al.  Transmission potential of smallpox: estimates based on detailed data from an outbreak. , 2003, American journal of epidemiology.

[6]  I. Kurane,et al.  Recent advances in the study of live attenuated cell-cultured smallpox vaccine LC16m8 , 2015, Vaccine.

[7]  Fei Xiang,et al.  Efficient MCMC for temporal epidemics via parameter reduction , 2014, Comput. Stat. Data Anal..

[8]  A. Nizam,et al.  Containing Bioterrorist Smallpox , 2002, Science.

[9]  P. Yip,et al.  Estimating the initial relative infection rate for a stochastic epidemic model. , 1989, Theoretical population biology.

[10]  T. Kypraios.  A note on maximum likelihood estimation of the initial number of susceptibles in the general stochastic epidemic model , 2009 .

[11]  Andrew Golightly,et al.  Delayed acceptance particle MCMC for exact inference in stochastic kinetic models , 2014, Stat. Comput..

[12]  Rob Deardon,et al.  Computational Statistics and Data Analysis Simulation-based Bayesian Inference for Epidemic Models , 2022 .

[13]  Richard J Boys,et al.  Bayesian inference for stochastic epidemic models with time-inhomogeneous removal rates , 2007, Journal of mathematical biology.

[14]  P. O’Neill,et al.  Bayesian estimation of the basic reproduction number in stochastic epidemic models , 2008 .

[15]  J. Wallinga,et al.  Different Epidemic Curves for Severe Acute Respiratory Syndrome Reveal Similar Impacts of Control Measures , 2004, American journal of epidemiology.

[16]  G. Roberts,et al.  Bayesian inference for partially observed stochastic epidemics , 1999 .

[17]  I. Arita,et al.  The smallpox threat: a time to reconsider global policy. , 2014, Biosecurity and bioterrorism : biodefense strategy, practice, and science.

[18]  N. Becker,et al.  Estimation for an epidemic model. , 1976, Biometrics.

[19]  Changhyuck Oh.  Improved Estimation of the Initial Number of Susceptible Individuals in the General Stochastic Epidemic Model Using Penalized Likelihood , 2014, TheScientificWorldJournal.

[20]  R. Shanmugam.  Spinned Poisson distribution with health management application , 2011, Health care management science.

[21]  Gavin J. Gibson,et al.  Estimating parameters in stochastic compartmental models using Markov chain methods , 1998 .

[22]  Richard M. Huggins,et al.  A note on the estimation of the initial number of susceptible individuals in the general epidemic model , 2004 .

[23]  A. HendersonD.,et al.  The smallpox threat: a time to reconsider global policy. , 2014 .

[24]  J. D. Millar,et al.  Modeling potential responses to smallpox as a bioterrorist weapon. , 2001, Emerging infectious diseases.

[25]  N T Bailey,et al.  The estimation of parameters from population data on the general stochastic epidemic. , 1971, Theoretical population biology.

[26]  Theodore Kypraios,et al.  Bayesian non-parametric inference for stochastic epidemic models using Gaussian Processes , 2016, Biostatistics.

[27]  P. Yip,et al.  Estimating the Basic Reproductive Number in the General Epidemic Model with an Unknown Initial Number of Susceptible Individuals , 2008 .

[28]  Steve Leach,et al.  Transmission potential of smallpox in contemporary populations , 2001, Nature.