BaySTDetect: detecting unusual temporal patterns in small area data via Bayesian model choice.

Space-time modeling of small area data is often used in epidemiology for mapping chronic disease rates and by government statistical agencies for producing local estimates of, for example, unemployment or crime rates. Although there is typically a general temporal trend, which affects all areas similarly, abrupt changes may occur in a particular area, e.g. due to emergence of localized predictors/risk factor(s) or impact of a new policy. Detection of areas with "unusual" temporal patterns is therefore important as a screening tool for further investigations. In this paper, we propose BaySTDetect, a novel detection method for short-time series of small area data using Bayesian model choice between two competing space-time models. The first model is a multiplicative decomposition of the area effect and the temporal effect, assuming one common temporal pattern across the whole study region. The second model estimates the time trends independently for each area. For each area, the posterior probability of belonging to the common trend model is calculated, which is then used to classify the local time trend as unusual or not. Crucial to any detection method, we provide a Bayesian estimate of the false discovery rate (FDR). A comprehensive simulation study has demonstrated the consistent good performance of BaySTDetect in detecting various realistic departure patterns in addition to estimating well the FDR. The proposed method is applied retrospectively to mortality data on chronic obstructive pulmonary disease (COPD) in England and Wales between 1990 and 1997 (a) to test a hypothesis that a government policy increased the diagnosis of COPD and (b) to perform surveillance. While results showed no evidence supporting the hypothesis regarding the policy, an identified unusual district (Tower Hamlets in inner London) was later recognized to have higher than national rates of hospital readmission and mortality due to COPD by the National Health Service, which initiated various local enhanced services to tackle the problem. Our method would have led to an early detection of this local health issue.

[1]  P. Müller,et al.  Optimal Sample Size for Multiple Testing , 2004 .

[2]  M. Kulldorff,et al.  A Space–Time Permutation Scan Statistic for Disease Outbreak Detection , 2005, PLoS medicine.

[3]  L Knorr-Held,et al.  Bayesian modelling of inseparable space-time variation in disease risk. , 2000, Statistics in medicine.

[4]  Jon Wakefield,et al.  A Bayesian measure of the probability of false discovery in genetic epidemiology studies. , 2007, American journal of human genetics.

[5]  Sylvia Richardson,et al.  Use of Space–Time Models to Investigate the Stability of Patterns of Disease , 2008, Environmental health perspectives.

[6]  Interrupted Time Series Quasi-Experiments , 2022 .

[7]  A. Lawson,et al.  Review of methods for space–time disease surveillance , 2010, Spatial and Spatio-temporal Epidemiology.

[8]  B. Miller,et al.  Cause-specific mortality in British coal workers and exposure to respirable dust and quartz , 2009, Occupational and Environmental Medicine.

[9]  A. Hansell,et al.  Geographic Variations in Risk: Adjusting for Unmeasured Confounders Through Joint Modeling of Multiple Diseases , 2009, Epidemiology.

[10]  P. Green,et al.  Modelling spatially correlated data via mixtures: a Bayesian approach , 2002 .

[11]  E G Knox,et al.  The Detection of Space‐Time Interactions , 1964 .

[12]  L Knorr-Held,et al.  Modelling risk from a disease in time and space. , 1998, Statistics in medicine.

[13]  Martin Kulldorff,et al.  Prospective time periodic geographical disease surveillance using a scan statistic , 2001 .

[14]  Y. MacNab,et al.  Autoregressive Spatial Smoothing and Temporal Spline Smoothing for Mapping Rates , 2001, Biometrics.

[15]  A. Gelman Prior distributions for variance parameters in hierarchical models (comment on article by Browne and Draper) , 2004 .

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

[17]  C. D. Mathers,et al.  Chronic obstructive pulmonary disease: current burden and future projections , 2006, European Respiratory Journal.

[18]  Sylvia Richardson,et al.  A comparison of Bayesian spatial models for disease mapping , 2005, Statistical methods in medical research.

[19]  Kenneth Rice,et al.  FDR and Bayesian Multiple Comparisons Rules , 2006 .

[20]  D. Rubin,et al.  Inference from Iterative Simulation Using Multiple Sequences , 1992 .

[21]  Deepayan Sarkar,et al.  Detecting differential gene expression with a semiparametric hierarchical mixture method. , 2004, Biostatistics.

[22]  D. Cocchi,et al.  Multiple testing on standardized mortality ratios: a Bayesian hierarchical model for FDR estimation. , 2011, Biostatistics.

[23]  Andrew Gelman,et al.  General methods for monitoring convergence of iterative simulations , 1998 .

[24]  L. Fahrmeir,et al.  Bayesian inference for generalized additive mixed models based on Markov random field priors , 2001 .

[25]  Andrew Thomas,et al.  WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility , 2000, Stat. Comput..

[26]  D. Strachan,et al.  Validity and interpretation of mortality, health service and survey data on COPD and asthma in England , 2003, European Respiratory Journal.

[27]  Robert Haining,et al.  Spatial Data Analysis: Theory and Practice , 2003 .

[28]  Ying C. MacNab,et al.  Spline smoothing in Bayesian disease mapping , 2007 .

[29]  N. Mantel The detection of disease clustering and a generalized regression approach. , 1967, Cancer research.

[30]  A. Hansell,et al.  What do chronic obstructive pulmonary disease patients die from? Amultiple cause coding analysis , 2003, European Respiratory Journal.

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

[32]  R. Rudd Coal miners’ respiratory disease litigation , 1998, Thorax.

[33]  Bradley P. Carlin,et al.  Hierarchical Spatio-Temporal Mapping of Disease Rates , 1997 .

[34]  D. Coggon,et al.  Coal mining and chronic obstructive pulmonary disease: a review of the evidence , 1998, Thorax.

[35]  A. Seaton The new prescription: industrial injuries benefits for smokers? , 1998, Thorax.

[36]  Leonhard Held,et al.  Gaussian Markov Random Fields: Theory and Applications , 2005 .

[37]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .