Comparison of Picard Versions for Analyzing functional Magnetic Resonance Imaging Data

Abstract

Independent Component Analysis (ICA) is a popular method that uses statistical principles to separate the mixture into statistically independent non-Gaussian sources. It has been well used in functional Magnetic Resonance Imaging (fMRI) data. However, real fMRI data can rarely be accurately modeled as mixtures of independent components, the convergence of ICA may be impaired. This paper is based on the idea of preconditioned ICA for real data (Picard), which involves a preprocessing L-BFGS strategy based on orthogonal matrix sets. In this study, we designed an experiment to validate the idea that Picard can improve ICA algorithms such as Infomax, Extended-Infomax, and FastICA, respectively named Picard 1, Picard 2, and Picard 3, for fMRI data analysis. Three Picard versions were performed on the simulated and noisy fMRI mixtures to verify the ability to separate independent sources. Experimental results showed that Picard 3 outperformed Picard 1 and Picard 2 on both noiseless and noisy simulated fMRI data, which implied the priority of Picard 3 in fMRI data analysis.