期刊
JOURNAL OF COMPUTATIONAL AND GRAPHICAL STATISTICS
卷 32, 期 2, 页码 413-433出版社
TAYLOR & FRANCIS INC
DOI: 10.1080/10618600.2022.2104289
关键词
Bayesian methods; Computationally intensive methods; Expectation-maximization; Neuroimaging; Spatial analysis
Independent Component Analysis (ICA) is commonly used in extracting functional brain networks from fMRI data. However, single-subject ICA often produces noisy results. This study proposes a spatial template ICA (stICA) method that incorporates spatial priors into the template ICA framework for greater estimation efficiency. The results show that stICA produces more accurate and reliable estimates compared to benchmark approaches, and identifies larger and more reliable areas of engagement.
Independent component analysis is commonly applied to functional magnetic resonance imaging (fMRI) data to extract independent components (ICs) representing functional brain networks. While ICA produces reliable group-level estimates, single-subject ICA often produces noisy results. Template ICA is a hierarchical ICA model using empirical population priors to produce more reliable subject-level estimates. However, this and other hierarchical ICA models assume unrealistically that subject effects are spatially independent. Here, we propose spatial template ICA (stICA), which incorporates spatial priors into the template ICA framework for greater estimation efficiency. Additionally, the joint posterior distribution can be used to identify brain regions engaged in each network using an excursions set approach. By leveraging spatial dependencies and avoiding massive multiple comparisons, stICA has high power to detect true effects. We derive an efficient expectation-maximization algorithm to obtain maximum likelihood estimates of the model parameters and posterior moments of the latent fields. Based on analysis of simulated data and fMRI data from the Human Connectome Project, we find that stICA produces estimates that are more accurate and reliable than benchmark approaches, and identifies larger and more reliable areas of engagement. The algorithm is computationally tractable, achieving convergence within 12 hr for whole-cortex fMRI analysis. Supplementary materials for this article are available online.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据