A critique of Tensor Probabilistic Independent Component Analysis: Implications and recommendations for multi-subject fMRI data analysis

• Published in: Journal of neuroscience methods Vol. 213; no. 2; pp. 263 - 273
• Main Authors: Helwig, Nathaniel E., Hong, Sungjin
• Format: Journal Article
• Language: English
• Published: Netherlands Elsevier B.V 15.03.2013
• Subjects: Humans; Image Interpretation, Computer-Assisted - methods; Magnetic Resonance Imaging; Neuroimage data analysis; Parafac; Parallel Factor Analysis; PICA; Tensor PICA; Tensor Probabilistic Independent Component Analysis; Neuroimage data analysis; Tensor PICA; Tensor Probabilistic Independent Component Analysis; PICA; Parallel Factor Analysis; Parafac
• ISSN: 0165-0270; 1872-678X; 1872-678X
• DOI: 10.1016/j.jneumeth.2012.12.009

Tensor Probabilistic Independent Component Analysis (TPICA) is a popular tool for analyzing multi-subject fMRI data (voxels×time×subjects) because of TPICA's supposed robustness. In this paper, we show that TPICA is not as robust as its authors claim. Specifically, we discuss why TPICA's overall objective is questionable, and we present some flaws related to the iterative nature of the TPICA algorithm. To demonstrate the relevance of these issues, we present a simulation study that compares TPICA versus Parallel Factor Analysis (Parafac) for analyzing simulated multi-subject fMRI data. Our simulation results demonstrate that TPICA produces a systematic bias that increases with the spatial correlation between the true components, and that the quality of the TPICA solution depends on the chosen ICA algorithm and iteration scheme. Thus, TPICA is not robust to small-to-moderate deviations from the model's spatial independence assumption. In contrast, Parafac produces unbiased estimates regardless of the spatial correlation between the true components, and Parafac with orthogonality-constrained voxel maps produces smaller biases than TPICA when the true voxel maps are moderately correlated. As a result, Parafac should be preferred for the analysis multi-subject fMRI data where the underlying components may have spatially overlapping voxel activation patterns.