🤖 AI Summary
This study addresses the challenge of disentangling spatiotemporal confounding from nonlinear exposure–response relationships in air pollution health effect estimation. We propose a Bayesian dynamic generalized linear model featuring a novel spatiotemporal confounding decomposition mechanism that separates confounding effects into fine- and coarse-scale components, coupled with spline-based modeling to capture the nonlinear association between PM₂.₅ and mortality. Empirical analysis using nationwide Italian data reveals significant seasonal heterogeneity in PM₂.₅-related mortality risk, with the highest risk observed in summer. Simulation studies demonstrate that the proposed method effectively corrects for confounding bias and accurately recovers the true exposure–response curve. This work establishes a new paradigm for causal inference in environmental epidemiology under high-dimensional spatiotemporal confounding.
📝 Abstract
Spatial confounding, often regarded as a major concern in epidemiological studies, relates to the difficulty of recovering the effect of an exposure on an outcome when these variables are associated with unobserved factors. This issue is particularly challenging in spatio-temporal analyses, where it has been less explored so far. To study the effects of air pollution on mortality in Italy, we argue that a model that simultaneously accounts for spatio-temporal confounding and for the non-linear form of the effect of interest is needed. To this end, we propose a Bayesian dynamic generalized linear model, which allows for a non-linear association and for a decomposition of the exposure effect into two components. This decomposition accommodates associations with the outcome at fine and coarse temporal and spatial scales of variation. These features, when combined, allow reducing the spatio-temporal confounding bias and recovering the true shape of the association, as demonstrated through simulation studies. The results from the real-data application indicate that the exposure effect seems to have different magnitudes in different seasons, with peaks in the summer. We hypothesize that this could be due to possible interactions between the exposure variable with air temperature and unmeasured confounders.