Spatially varying distributed lag non-linear models using Laplacian P-splines

2026-04-10

AI Summary

This study addresses the limitation of traditional distributed lag nonlinear models (DLNMs), which assume spatially constant exposure–response relationships and thus struggle to capture spatial heterogeneity in sparse count data at fine geographical scales. To overcome this, we propose a computationally efficient Bayesian single-stage spatially varying coefficient DLNM that innovatively integrates Laplace P-splines with Laplace approximation as an alternative to conventional MCMC methods. This approach substantially reduces computational cost while preserving modeling flexibility. The framework accommodates four distinct spatial dependence structures and spline specifications, accompanied by corresponding information criteria for model selection. Simulation studies demonstrate the method’s strong performance, and its practical utility is illustrated through an application to temperature–mortality association analysis at the municipal level in Sicily, Italy, confirming its effectiveness with real-world sparse data.

Technology Category

• Machine Learning: Bayesian Learning
• Data Mining & Knowledge Management: Mining of Spatial, Temporal or Spatio-Temporal Data
• Constraint Satisfaction and Optimization: Distributed CSP/Optimization

Application Category

• Graph Algorithms and Modeling for the Web: Algorithms and analysis for incomplete, noisy, or partially observed Web-related graphs
• User Modeling, Personalization and Recommendation: Explainable and interpretable methods for personalization
• Web Mining and Content Analysis: Models for Web evolution

Abstract

Although distributed lag non-linear models (DLNMs) are commonly used to quantify delayed and non-linear exposure-response relationships, most existing applications assume that these relationships are constant across space. However, in many geographical and environmental studies, local characteristics vary substantially across areas, making a spatially varying effect more realistic. Extending DLNMs to allow for spatial heterogeneity remains challenging, and only a limited number of modelling strategies have been proposed in literature. The most popular extension is a two-stage meta-analysis approach, which requires sufficiently large sample sizes at each location. Therefore, its usefulness is limited when working with sparse count data in small area data analyses. Although a number of alternative one-stage approaches have been introduced, their computational burden restricts their applicability in real-life data applications. In this paper, we introduce a computationally efficient Bayesian one-stage spatially-varying DLNM for count data. We define four model variants, differing in the assumed spatial dependence structure and the flexibility of the DLNM spline specification. To address the computational burden typically associated with these flexible models, we use Laplace approximations, offering an efficient alternative to classically used Markov Chain Monte Carlo (MCMC) approaches. Model comparison criteria are provided to facilitate the selection of a suitable model in a real-life data application. The proposed methods are evaluated through simulation studies, and their practical usefulness is illustrated through a real-life data application, investigating the temperature-mortality relationship in every municipality of Sicily, Italy.

Problem

• distributed lag non-linear models
• spatial heterogeneity
• count data
• small area estimation
• exposure-response relationship

Innovation

• spatially varying DLNM
• Laplacian P-splines
• Bayesian one-stage model
• Laplace approximation
• small area count data