Space-time interactions in Bayesian disease mapping with recent tools: Making things easier for practitioners

Abstract

Spatio-temporal disease mapping studies the distribution of mortality or incidence risks in space and its evolution in time, and it usually relies on fitting hierarchical Poisson mixed models. These models are complex for practitioners as they generally require adding constraints to correctly identify and interpret the different model terms. However, including constraints may not be straightforward in some recent software packages. This paper focuses on NIMBLE, a library of algorithms that contains among others a configurable system for Markov chain Monte Carlo (MCMC) algorithms. In particular, we show how to fit different spatio-temporal disease mapping models with NIMBLE making emphasis on how to include sum-to-zero constraints to solve identifiability issues when including spatio-temporal interactions. Breast cancer mortality data in Spain during the period 1990-2010 is used for illustration purposes. A simulation study is also conducted to compare NIMBLE with R-INLA in terms of parameter estimates and relative risk estimation. The results are very similar but differences are observed in terms of computing time.

Keywords

Disease mapping identifiability INLA NIMBLE spatio-temporal interactions sum-to-zero constraints

Get full access to this article

View all access options for this article.

References

Besag

. Spatial interaction and the statistical analysis of lattice systems (with discussion). Journal of the Royal Statistical Society: Series B (Statistical Methodology) 1974; 36: 192–236.

Goicoa

Ugarte

Etxeberria

et al. Comparing CAR and P-spline models in spatial disease mapping. Environ Ecol Stat 2012; 19: 573–599.

Bernardinelli

Clayton

Pascutto

et al. Bayesian analysis of space-time variation in disease risk. Stat Med 1995; 14: 2433–2443.

Assunçao

Reis

Oliveira

. Diffusion and prediction of leishmaniasis in a large metropolitan area in Brazil with a Bayesian space-time model. Stat Med 2001; 20: 2319–2335.

Knorr-Held

. Bayesian modelling of inseparable space-time variation in disease risk. Stat Med 2000; 19: 2555–2567.

Ugarte

Adin

Goicoa

et al. On fitting spatio-temporal disease mapping models using approximate bayesian inference. Stat Methods Med Res 2014; 23: 507–530.

Ugarte

Goicoa

Etxeberria

et al. A P-spline ANOVA type model in space-time disease mapping. Stochastic Environmental Research and Risk Assesment 2012; 26: 835–845.

Ugarte

Goicoa

Militino

. Spatio-temporal modelling of mortality risks using penalized splines. Environmetrics 2010; 21: 270–289.

Gelfand

Sahu

. Identifiability, improper priors, and gibbs sampling for generalized linear models. J Am Stat Assoc 1999; 94: 247–253.

10.

Goicoa

Adin

Ugarte

et al. In spatio-temporal disease mapping models, identifiability constraints affect PQL and INLA results. Stoch Environ Res Risk Assess 2018; 32: 749–770.

11.

deValpine

Turek

Paciorek

et al. Programming with models: writing statistical algorithms for general model structures with NIMBLE. J Comput Graph Stat 2017; 26: 403–417.

12.

Lindgren

Rue

. Bayesian spatial modelling with R-INLA. J Stat Softw 2015; 63: 1–25.

13.

Rue

Martino

Chopin

. Approximate Bayesian inference for latent Gaussian models by using integrated nested laplace approximations. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 2009; 71: 319–392.

14.

Martínez-Beneito

. A general modelling framework for multivariate disease mapping. Biometrika 2013; 100: 539–553.

15.

Besag

York

Mollié

. Bayesian image restoration, with two applications in spatial statistics. Annals of the institute of statistical mathematics. Ann Inst Stat Math 1991; 43: 1–20.

16.

Rue

Held

. Gaussian Markov Random Fields: Theory and Applications. Boca Raton: CRC press, 2005.

17.

Lunn

Thomas

Best

et al. WinBUGS - A Bayesian modelling framework: Concepts, structure, and extensibility. Stat Comput 2000; 10: 325–337.

18.

de Valpine

Paciorek

Ture

et al. NIMBLE User Manual, 2021; doi:10.5281/zenodo.1211190. https://r-nimble.org. R package manual version 0.11.1.

19.

Lawson

. NIMBLE for Bayesian disease mapping. Spat Spatiotemporal Epidemiol 2020; 33: 100323.

20.

Adin

Martínez-Beneito

Botella-Rocamora

et al. Smoothing and high risk areas detection in space-time disease mapping: a comparison of P-splines, autoregressive, and moving average models. Stoch Environ Res Risk Assess 2017; 31: 403–415.

21.

Harville

. Matrix Algebra From a Statistician’s Perspective. New York: Springer, 2008.

22.

Geman

. Stochastic relaxation, gibbs distributions, and the Bayesian restoration of images. IEEE Trans Pattern Anal Mach Intell 1984; PAMI-6: 721–741.

23.

Gelfand

Smith

AFM

. Sampling-based approaches to calculating marginal densities. J Am Stat Assoc 1990; 85: 398–409.

24.

Casella

George

. Explaining the gibbs sampler. Am Stat 1992; 46: 167–174.

25.

Roberts

Sahu

. Updating schemes, correlation structure, blocking and parameterization for the gibbs sampler. J R Stat Soc Ser B (Methodological) 1997; 59: 291–317.

26.

Neal

. Slice sampling. Ann Stat 2003; 31: 705–767.

27.

Brooks

Gelman

. General methods for monitoring convergence of iterative simulations. J Comput Graph Stat 1998; 7: 434–455.

28.

Gelman

Carlin

Stern

et al. Bayesian Data Analysis. Third Ed. Boca Raton: Chapman and Hall/CRC, 2013.

29.

Plummer

Best

Cowles

et al. CODA: Convergence diagnosis and output analysis for MCMC. R News 2006; 6: 7–11.

30.

Spiegelhalter

Thomas

Best

et al. WinBUGS User Manual version 1.4. 2003.

31.

Gelman

. Prior distributions for variance parameters in hierarchical models (comment on article by browne and draper). Bayesian Anal 2006; 1: 515–534.

32.

Spiegelhalter

Best

Carlin

et al. Bayesian measures of model complexity and fit. J R Stat Soc: Ser B (Statistical Methodology) 2002; 64: 583–639.

33.

Watanabe

. Asymptotic equivalence of Bayes cross validation and widely applicable information criterion in singular learning theory. J Mach Learn Res 2010; 11: 3571–3594.

34.

Vranckx M, Neyens T and Faes C. The (in)stability of bayesian model selection criteria in disease mapping. Spat Stat 2021; 43: 100502.

35.

MacNab

. On identification in Bayesian disease mapping and ecological-spatial regression models. Stat Methods Med Res 2014; 23: 134–155.

36.

Wood

Scheipl

Faraway

. Straightforward intermediate rank tensor product smoothing in mixed models. Stat Comput 2013; 23: 341–360.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

2.08 MB