Regression B‐spline smoothing in Bayesian disease mapping: with an application to patient safety surveillance
- 14 March 2007
- journal article
- research article
- Published by Wiley in Statistics in Medicine
- Vol. 26 (24) , 4455-4474
- https://doi.org/10.1002/sim.2868
Abstract
In the context of Bayesian disease mapping, recent literature presents generalized linear mixed models that engender spatial smoothing. The methods assume spatially varying random effects as a route to partially pooling data and ‘borrowing strength’ in small‐area estimation. When spatiotemporal disease rates are available for sequential risk mapping of several time periods, the ‘smoothing’ issue may be explored by considering spatial smoothing, temporal smoothing and spatiotemporal interaction. In this paper, these considerations are motivated and explored through development of a Bayesian semiparametric disease mapping model framework which facilitates temporal smoothing of rates and relative risksviaregression B‐splines with mixed‐effect representation of coefficients. Specifically, we develop spatial priors such as multivariate Gaussian Markov random fields and non‐spatial priors such as unstructured multivariate Gaussian distributions and illustrate how time trends in small‐area relative risks may be explored by splines which vary in either a spatially structured or unstructured manner. In particular, we show that with suitable prior specifications for the random effects ensemble, small‐area relative risk trends may be fit by ‘spatially varying’ or randomly varying B‐splines. A recently developed Bayesian hierarchical model selection criterion, the deviance information criterion, is used to assess the trade‐off between goodness‐of‐fit and smoothness and to select the number of knots. The methodological development aims to provide reliable information about the patterns (both over space and time) of disease risks and to quantify uncertainty. The study offers a disease and health outcome surveillance methodology for flexible and efficient exploration and assessment of emerging risk trends and clustering. The methods are motivated and illustrated through a Bayesian analysis of adverse medical events (also known asiatrogenic injuries) among hospitalized elderly patients in British Columbia, Canada. Copyright © 2007 John Wiley & Sons, Ltd.Keywords
This publication has 29 references indexed in Scilit:
- An innovative application of Bayesian disease mapping methods to patient safety research: a Canadian adverse medical event studyStatistics in Medicine, 2006
- Bayesian P-SplinesJournal of Computational and Graphical Statistics, 2004
- Hierarchical Bayesian spatial modelling of small‐area rates of non‐rare diseaseStatistics in Medicine, 2003
- Bayesian Measures of Model Complexity and FitJournal of the Royal Statistical Society Series B: Statistical Methodology, 2002
- The Analysis of Designed Experiments and Longitudinal Data by Using Smoothing SplinesJournal of the Royal Statistical Society Series C: Applied Statistics, 1999
- Inference in Generalized Additive Mixed Models by Using Smoothing SplinesJournal of the Royal Statistical Society Series B: Statistical Methodology, 1999
- Smoothing Spline Models for the Analysis of Nested and Crossed Samples of CurvesJournal of the American Statistical Association, 1998
- Hierarchical Spatio-Temporal Mapping of Disease RatesJournal of the American Statistical Association, 1997
- Approximate Inference in Generalized Linear Mixed ModelsJournal of the American Statistical Association, 1993
- On the distinction between the conditional probability and the joint probability approaches in the specification of nearest-neighbour systemsBiometrika, 1964