Spatial model for risk prediction and sub-national prioritization to aid poliovirus eradication in Pakistan

BackgroundPakistan is one of only three countries where poliovirus circulation remains endemic. For the Pakistan Polio Eradication Program, identifying high risk districts is essential to target interventions and allocate limited resources.MethodsUsing a hierarchical Bayesian framework we developed a spatial Poisson hurdle model to jointly model the probability of one or more paralytic polio cases, and the number of cases that would be detected in the event of an outbreak. Rates of underimmunization, routine immunization, and population immunity, as well as seasonality and a history of cases were used to project future risk of cases.ResultsThe expected number of cases in each district in a 6-month period was predicted using indicators from the previous 6-months and the estimated coefficients from the model. The model achieves an average of 90% predictive accuracy as measured by area under the receiver operating characteristic (ROC) curve, for the past 3 years of cases.ConclusionsThe risk of poliovirus has decreased dramatically in many of the key reservoir areas in Pakistan. The results of this model have been used to prioritize sub-national areas in Pakistan to receive additional immunization activities, additional monitoring, or other special interventions.

[1]  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.

[2]  Analysis of vaccination campaign effectiveness and population immunity to support and sustain polio elimination in Nigeria , 2016, BMC Medicine.

[3]  Leonhard Held,et al.  Posterior and Cross-validatory Predictive Checks: A Comparison of MCMC and INLA , 2010 .

[4]  O. Diop,et al.  Surveillance Systems to Track Progress Toward Polio Eradication — Worldwide, 2015–2016 , 2017, MMWR. Morbidity and mortality weekly report.

[5]  Ahmed S Bux,et al.  Polio in Pakistan's North Waziristan. , 2015, The Lancet. Global health.

[6]  J. Hanley,et al.  The meaning and use of the area under a receiver operating characteristic (ROC) curve. , 1982, Radiology.

[7]  J. Alexander,et al.  Progress and peril: poliomyelitis eradication efforts in Pakistan, 1994-2013. , 2014, The Journal of infectious diseases.

[8]  Pakistan. Pakistan Demographic and Health Survey 2012-13 , 2013 .

[9]  J. Glass,et al.  Differential diagnosis of acute flaccid paralysis and its role in poliomyelitis surveillance. , 2000, Epidemiologic reviews.

[10]  M. Greiner,et al.  Principles and practical application of the receiver-operating characteristic analysis for diagnostic tests. , 2000, Preventive veterinary medicine.

[11]  Pulak Ghosh,et al.  A spatial Poisson hurdle model for exploring geographic variation in emergency department visits , 2013, Journal of the Royal Statistical Society. Series A,.

[12]  Jiqiang Guo,et al.  Stan: A Probabilistic Programming Language. , 2017, Journal of statistical software.

[13]  Leonhard Held,et al.  Spatio‐temporal disease mapping using INLA , 2011 .

[14]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[15]  E. Ionides,et al.  Dynamics affecting the risk of silent circulation when oral polio vaccination is stopped , 2017, Epidemics.

[16]  Ali Arab,et al.  Spatial and Spatio-Temporal Models for Modeling Epidemiological Data with Excess Zeros , 2015, International journal of environmental research and public health.

[17]  M. Pallansch,et al.  Duration of poliovirus excretion and its implications for acute flaccid paralysis surveillance: a review of the literature. , 1997, The Journal of infectious diseases.

[18]  Y. Manor,et al.  Molecular epidemiology of silent introduction and sustained transmission of wild poliovirus type 1, Israel, 2013. , 2014, Euro surveillance : bulletin Europeen sur les maladies transmissibles = European communicable disease bulletin.

[19]  Finn Lindgren,et al.  Bayesian Spatial Modelling with R-INLA , 2015 .

[20]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[21]  J. Corander,et al.  Implications of a circulating vaccine-derived poliovirus in Nigeria. , 2010, The New England journal of medicine.

[22]  Kaushik Banerjee,et al.  Development of a District-Level Programmatic Assessment Tool for Risk of Measles Virus Transmission. , 2017, Risk analysis : an official publication of the Society for Risk Analysis.

[23]  Guillaume Chabot-Couture,et al.  Predictive spatial risk model of poliovirus to aid prioritization and hasten eradication in Nigeria , 2014, BMC Medicine.

[24]  A. Khurshid,et al.  Detection of multiple cocirculating wild poliovirus type 1 lineages through environmental surveillance: impact and progress during 2011-2013 in Pakistan. , 2014, The Journal of infectious diseases.

[25]  K. Hajian‐Tilaki,et al.  Receiver Operating Characteristic (ROC) Curve Analysis for Medical Diagnostic Test Evaluation. , 2013, Caspian journal of internal medicine.

[26]  Progress toward poliomyelitis eradication--Pakistan, 1994-1998. , 1999, MMWR. Morbidity and mortality weekly report.

[27]  N. Grassly,et al.  Asymptomatic wild-type poliovirus infection in India among children with previous oral poliovirus vaccination. , 2010, The Journal of infectious diseases.

[28]  Rebecca Martin,et al.  World Health Organization Regional Assessments of the Risks of Poliovirus Outbreaks , 2013, Risk analysis : an official publication of the Society for Risk Analysis.

[29]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[30]  J. Wakefield,et al.  Bayesian inference for generalized linear mixed models. , 2010, Biostatistics.

[31]  H. Rue,et al.  An explicit link between Gaussian fields and Gaussian Markov random fields: the stochastic partial differential equation approach , 2011 .

[32]  J. Stockman Effectiveness of Immunization against Paralytic Poliomyelitis in Nigeria , 2010 .

[33]  H. Rue,et al.  Improving the INLA approach for approximate Bayesian inference for latent Gaussian models , 2015, 1503.07307.

[34]  N. Nathanson,et al.  From Emergence to Eradication: The Epidemiology of Poliomyelitis Deconstructed , 2010, American journal of epidemiology.

[35]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .