Summary Independent component analysis (ICA) has become an important tool for analyzing data from functional magnetic resonance imaging (fMRI) studies. ICA has been successfully applied to single‐subject fMRI data. The extension of ICA to group inferences
in neuroimaging studies, however, is challenging due to the unavailability of a prespecified group design matrix and the uncertainty in between‐subjects variability in fMRI data. We present a general probabilistic ICA (PICA) model that can accommodate varying group structures of multisubject
spatiotemporal processes. An advantage of the proposed model is that it can flexibly model various types of group structures in different underlying neural source signals and under different experimental conditions in fMRI studies. A maximum likelihood (ML) method is used for estimating this
general group ICA model. We propose two expectation–maximization (EM) algorithms to obtain the ML estimates. The first method is an exact EM algorithm, which provides an exact E‐step and an explicit noniterative M‐step. The second method is a variational approximation EM
algorithm, which is computationally more efficient than the exact EM. In simulation studies, we first compare the performance of the proposed general group PICA model and the existing probabilistic group ICA approach. We then compare the two proposed EM algorithms and show the variational
approximation EM achieves comparable accuracy to the exact EM with significantly less computation time. An fMRI data example is used to illustrate application of the proposed methods.
No Supplementary Data
Document Type: Research Article
Department of Biostatistics and Bioinformatics, Rollins School of Public Health of Emory University, 1518 Clifton Road North East, Atlanta, Georgia 30322, U.S.A.
Publication date: 2011-12-01