Causal Inference in Biomedical Imaging via Functional Linear Structural Equation Models
Ting Li, Ethan Fan, Tengfei Li, Hongtu Zhu
TL;DR
This work addresses causal inference for the effect of endogenous, infinite-dimensional functional imaging exposures on scalar clinical outcomes under unobserved confounding. The authors formulate a Functional Linear Structural Equation Model ($FLSEM$) with scalar covariates and functional instruments, establishing identifiability through an injective integral operator $\mathcal{K}$ and Mercer's decomposition of the kernel $K(s,t)$. Estimation combines a three-step, $L_0$-penalized, RKHS-based approach with a Functional Group Support Detection and Root Finding ($\text{FGSDAR}$) algorithm, implemented via a region-based divide-and-conquer scheme, and a nullity test for the functional coefficient $B(t)$ with a derived null distribution. Theoretical guarantees include selection consistency, nonasymptotic error bounds for both scalar and functional components, and a practical $S_n$ statistic approximated by a scaled $\chi^2$ distribution. Empirical evaluation in simulations and a UK Biobank application demonstrates robust performance in detecting causal relationships between genetics, brain imaging, and cognition, offering a scalable framework for causal imaging genetics research.
Abstract
Understanding the causal effects of organ-specific features from medical imaging on clinical outcomes is essential for biomedical research and patient care. We propose a novel Functional Linear Structural Equation Model (FLSEM) to capture the relationships among clinical outcomes, functional imaging exposures, and scalar covariates like genetics, sex, and age. Traditional methods struggle with the infinite-dimensional nature of exposures and complex covariates. Our FLSEM overcomes these challenges by establishing identifiable conditions using scalar instrumental variables. We develop the Functional Group Support Detection and Root Finding (FGS-DAR) algorithm for efficient variable selection, supported by rigorous theoretical guarantees, including selection consistency and accurate parameter estimation. We further propose a test statistic to test the nullity of the functional coefficient, establishing its null limit distribution. Our approach is validated through extensive simulations and applied to UK Biobank data, demonstrating robust performance in detecting causal relationships from medical imaging.
