On fitting spatio-temporal disease mapping models using approximate Bayesian inference

Spatio-temporal disease mapping comprises a wide range of models used to describe the distribution of a disease in space and its evolution in time. These models have been commonly formulated within a hierarchical Bayesian framework with two main approaches: an empirical Bayes (EB) and a fully Bayes (FB) approach. The EB approach provides point estimates of the parameters relying on the well-known penalized quasi-likelihood (PQL) technique. The FB approach provides the posterior distribution of the target parameters. These marginal distributions are not usually available in closed form and common estimation procedures are based on Markov chain Monte Carlo (MCMC) methods. However, the spatio-temporal models used in disease mapping are often very complex and MCMC methods may lead to large Monte Carlo errors and a huge computation time if the dimension of the data at hand is large. To circumvent these potential inconveniences, a new technique called integrated nested Laplace approximations (INLA), based on nested Laplace approximations, has been proposed for Bayesian inference in latent Gaussian models. In this paper, we show how to fit different spatio-temporal models for disease mapping with INLA using the Leroux CAR prior for the spatial component, and we compare it with PQL via a simulation study. The spatio-temporal distribution of male brain cancer mortality in Spain during the period 1986–2010 is also analysed.

[1]  Milton C. Chew Distributions in Statistics: Continuous Univariate Distributions-1 and 2 , 1971 .

[2]  Valerie Isham,et al.  Non‐Negative Matrices and Markov Chains , 1983 .

[3]  L. Tierney,et al.  Accurate Approximations for Posterior Moments and Marginal Densities , 1986 .

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

[5]  N. Breslow,et al.  Approximate inference in generalized linear mixed models , 1993 .

[6]  D. Clayton,et al.  Bayesian analysis of space-time variation in disease risk. , 1995, Statistics in medicine.

[7]  David G. Clayton,et al.  Generalized linear mixed models , 1996 .

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

[9]  Norman E. Breslow,et al.  Estimation of Disease Rates in Small Areas: A new Mixed Model for Spatial Dependence , 2000 .

[10]  R. Assunção,et al.  Diffusion and prediction of Leishmaniasis in a large metropolitan area in Brazil with a Bayesian space–time model , 2001, Statistics in medicine.

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

[12]  C B Dean,et al.  Detecting Interaction Between Random Region and Fixed Age Effects in Disease Mapping , 2001, Biometrics.

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

[14]  M. Pollán,et al.  Geographical pattern of brain cancer incidence in the Navarre and Basque Country regions of Spain , 2003, Occupational and environmental medicine.

[15]  Volker J Schmid,et al.  Bayesian Extrapolation of Space–Time Trends in Cancer Registry Data , 2004, Biometrics.

[16]  Charmaine B. Dean,et al.  Penalized quasi-likelihood with spatially correlated data , 2004, Comput. Stat. Data Anal..

[17]  S. Richardson,et al.  Interpreting Posterior Relative Risk Estimates in Disease-Mapping Studies , 2004, Environmental health perspectives.

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

[19]  Charmaine B. Dean,et al.  Approximate inference for disease mapping , 2006, Comput. Stat. Data Anal..

[20]  Jon Wakefield,et al.  Disease mapping and spatial regression with count data. , 2007, Biostatistics.

[21]  H. Rue,et al.  Approximate Bayesian inference for hierarchical Gaussian Markov random field models , 2007 .

[22]  Ying C MacNab,et al.  Regression B‐spline smoothing in Bayesian disease mapping: with an application to patient safety surveillance , 2007, Statistics in medicine.

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

[24]  M A Martínez-Beneito,et al.  An autoregressive approach to spatio-temporal disease mapping. , 2008, Statistics in medicine.

[25]  E. Seneta Non-negative Matrices and Markov Chains , 2008 .

[26]  Håvard Rue,et al.  Norges Teknisk-naturvitenskapelige Universitet Implementing Approximate Bayesian Inference Using Integrated Nested Laplace Approximation: a Manual for the Inla Program Implementing Approximate Bayesian Inference Using Integrated Nested Laplace Approximation: a Manual for the Inla Program , 2022 .

[27]  María Dolores Ugarte,et al.  Prediction error estimators in Empirical Bayes disease mapping , 2008 .

[28]  María Dolores Ugarte,et al.  Empirical Bayes and Fully Bayes procedures to detect high-risk areas in disease mapping , 2009, Comput. Stat. Data Anal..

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

[30]  María Dolores Ugarte,et al.  Spatio‐temporal modeling of mortality risks using penalized splines , 2009 .

[31]  María Dolores Ugarte,et al.  Evaluating the performance of spatio‐temporal Bayesian models in disease mapping , 2009 .

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

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

[34]  Ying C MacNab,et al.  On Gaussian Markov random fields and Bayesian disease mapping , 2011, Statistical methods in medical research.

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

[36]  Duncan Lee,et al.  A comparison of conditional autoregressive models used in Bayesian disease mapping. , 2011, Spatial and spatio-temporal epidemiology.

[37]  Leonhard Held,et al.  Using integrated nested Laplace approximations for the evaluation of veterinary surveillance data from Switzerland: a case‐study , 2011 .

[38]  Leonhard Held,et al.  A primer on disease mapping and ecological regression using $${\texttt{INLA}}$$ , 2011, Comput. Stat..

[39]  María Dolores Ugarte,et al.  A P-spline ANOVA type model in space-time disease mapping , 2012, Stochastic Environmental Research and Risk Assessment.

[40]  M. D. Ugarte,et al.  Comparing CAR and P-spline models in spatial disease mapping , 2012, Environmental and Ecological Statistics.

[41]  Duncan Lee,et al.  Locally adaptive spatial smoothing using conditional auto‐regressive models , 2012, 1205.3641.

[42]  M. D. Ugarte,et al.  Projections of cancer mortality risks using spatio-temporal P-spline models , 2012, Statistical methods in medical research.

[43]  P Botella-Rocamora,et al.  Spatial moving average risk smoothing , 2013, Statistics in medicine.

[44]  Andrew B. Lawson,et al.  Bayesian disease mapping , 2013 .

[45]  Gianluca Baio,et al.  Spatial and spatio-temporal models with R-INLA. , 2013, Spatial and spatio-temporal epidemiology.

[46]  Finn Lindgren,et al.  Bayesian computing with INLA: New features , 2012, Comput. Stat. Data Anal..

[47]  Ying C MacNab MacNab spatial regression models − On identification in Bayesian disease mapping and ecological , 2014 .