Big problems in spatio-temporal disease mapping: methods and software

Abstract

Fitting spatio-temporal models for areal data is crucial in many fields such as cancer epidemiology. However, when data sets are very large, many issues arise. The main objective of this paper is to propose a general procedure to analyze high-dimensional spatio-temporal count data, with special emphasis on mortality/incidence relative risk estimation. We present a pragmatic and simple idea that permits to fit hierarchical spatio-temporal models when the number of small areas is very large. Model fitting is carried out using integrated nested Laplace approximations over a partition of the spatial domain. We also use parallel and distributed strategies to speed up computations in a setting where Bayesian model fitting is generally prohibitively time-consuming and even unfeasible. Using simulated and real data, we show that our method outperforms classical global models. We implement the methods and algorithms that we develop in the open-source R package bigDM where specific vignettes have been included to facilitate the use of the methodology for non-expert users. Our scalable methodology proposal provides reliable risk estimates when fitting Bayesian hierarchical spatio-temporal models for high-dimensional data.

References in corpus (8)

• Asymptotic Equivalence of Bayes Cross Validation and Widely Applicable Information Criterion in Singular Learning Theory
• Programming with models: writing statistical algorithms for general model structures with NIMBLE
• A Unifying Framework for Parallel and Distributed Processing in R using Futures
• A general framework for Vecchia approximations of Gaussian processes
• A multi-resolution approximation for massive spatial datasets
• Highly Scalable Bayesian Geostatistical Modeling via Meshed Gaussian Processes on Partitioned Domains
• Alleviating confounding in spatio-temporal areal models with an application on crimes against women in India
• Scalable Bayesian modeling for smoothing disease risks in large spatial data sets

Cited by in corpus (4)

• High-dimensional order-free multivariate spatial disease mapping
• A scalable approach for short-term disease forecasting in high spatial resolution areal data
• A fast approach for analyzing spatio-temporal patterns in ischemic heart disease mortality across US counties (1999-2021)
• Improving Disease Risk Estimation in Small Areas by Accounting for Spatiotemporal Local Discontinuities