Abstract:Abstract Background Distributed lag non-linear models (DLNMs) are the reference framework for modelling lagged non-linear associations. They are usually used in large-scale multi-location studies. Attempts to study these associations in small areas either did not include the lagged non-linear effects, did not allow for geographically-varying risks or downscaled risks from larger spatial units through socioeconomic and physical meta-predictors when the estimation of the risks was not feasible due to low statistical power. Methods Here we proposed spatial Bayesian DLNMs (SB-DLNMs) as a new framework for the estimation of reliable small-area lagged non-linear associations, and demonstrated the methodology for the case study of the temperature-mortality relationship in the 73 neighbourhoods of the city of Barcelona. We generalized location-independent DLNMs to the Bayesian framework (B-DLNMs), and extended them to SB-DLNMs by incorporating spatial models in a single-stage approach that accounts for the spatial dependence between risks. Results The results of the case study highlighted the benefits of incorporating the spatial component for small-area analysis. Estimates obtained from independent B-DLNMs were unstable and unreliable, particularly in neighbourhoods with very low numbers of deaths. SB-DLNMs addressed these instabilities by incorporating spatial dependencies, resulting in more plausible and coherent estimates and revealing hidden spatial patterns. In addition, the Bayesian framework enriches the range of estimates and tests that can be used in both large- and small-area studies. Conclusions SB-DLNMs account for spatial structures in the risk associations across small areas. By modelling spatial differences, SB-DLNMs facilitate the direct estimation of non-linear exposure-response lagged associations at the small-area level, even in areas with as few as 19 deaths. The manuscript includes an illustrative code to reproduce the results, and to facilitate the implementation of other case studies by other researchers.

What problem does this paper attempt to address?

The paper attempts to address the issues encountered when modeling the relationship between environmental exposure and health effects within small areas. Specifically, traditional Distributed Lag Non-linear Models (DLNMs) are typically used for large-scale multi-location studies, but their application within small areas faces the following problems: 1. **Ignoring lag non-linear effects**: Some studies do not include lag non-linear effects, which may lead to inaccurate models. 2. **Not allowing for geographic variation in risk**: Some studies do not consider the risk differences across different geographic locations, resulting in models that fail to reflect the true geographic distribution. 3. **Low statistical power**: Within small areas, due to the small sample size, the estimated risks may be unstable and unreliable. To overcome these issues, the authors propose the Spatial Bayesian Distributed Lag Non-linear Model (SB-DLNM), which improves estimation accuracy within small areas by introducing spatial dependency. This method combines the Bayesian framework and spatial models, allowing for the direct estimation of non-linear lag exposure-response relationships within small areas, even with small sample sizes, providing reliable results. ### Main Contributions of the Paper: 1. **Adaptation to Bayesian statistics**: Adapting the DLNM framework to Bayesian statistics provides a broader range of estimation and testing methods and intuitively explains uncertainty. 2. **Introduction of spatial models**: By introducing spatial models, the dependency between different geographic locations is considered, improving the stability and reliability of the estimates. 3. **Practical case study**: The study of the temperature-mortality relationship in 73 communities in Barcelona demonstrates the practical application of SB-DLNM. ### Highlights of the Results: - **Improved stability**: SB-DLNM provides more stable and reliable estimates in communities with small sample sizes. - **Revealing spatial patterns**: SB-DLNM reveals hidden spatial patterns, making the estimation results more reasonable and consistent. - **High statistical power**: The time series design of SB-DLNM shows higher statistical power in detecting increased risks. ### Conclusion: SB-DLNM, by modeling spatial differences and structures, can directly estimate non-linear lag exposure-response relationships within small areas, providing reliable results even with very small sample sizes. This method offers new tools and approaches for small-area environmental epidemiology research.

Spatial Bayesian distributed lag non-linear models (SB-DLNM) for small-area exposure-lag-response epidemiological modelling

Bayesian statistics in spatial epidemiology]

A Bayesian multivariate spatial approach for illness-death survival models

Bayesian modeling of spatial ordinal data from health surveys

Scalable Bayesian modeling for smoothing disease risks in large spatial data sets

A penalized distributed-lag non-linear model for modeling the joint delayed effect of two predictors: impact of minimum and maximum temperature on mortality

Treed distributed lag nonlinear models

A unifying modelling approach for hierarchical distributed lag models

Bayesian modelling for spatially misaligned health areal data: a multiple membership approach

Cluster detection and risk estimation for spatio-temporal health data

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

Bayesian Generalized Distributed Lag Regression with Variable Selection

High-dimensional order-free multivariate spatial disease mapping

Multiple exposure distributed lag models with variable selection

A Bayesian hierarchical small-area population model accounting for data source specific methodologies from American Community Survey, Population Estimates Program, and Decennial Census data

Spatial non-parametric Bayesian clustered coefficients

A scalable two-stage Bayesian approach accounting for exposure measurement error in environmental epidemiology

Bayesian nonparametric models for spatially indexed data of mixed type

Bayesian modelling of inseparable space-time variation in disease risk

A flexible hierarchical framework for improving inference in area-referenced environmental health studies

Using Integrated Nested Laplace Approximation for Modeling Spatial Healthcare Utilization