Estimating the number of independent components for functional magnetic resonance imaging data

Abstract

Multivariate analysis methods such as independent component analysis (ICA) have been applied to the analysis of functional magnetic resonance imaging (fMRI) data to study brain function. Because of the high dimensionality and high noise level of the fMRI data, order selection, i.e., estimation of the number of informative components, is critical to reduce over/underfitting in such methods. Dependence among fMRI data samples in the spatial and temporal domain limits the usefulness of the practical formulations of information-theoretic criteria (ITC) for order selection, since they are based on likelihood of independent and identically distributed (i.i.d.) data samples. To address this issue, we propose a subsampling scheme to obtain a set of effectively i.i.d. samples from the dependent data samples and apply the ITC formulas to the effectively i.i.d. sample set for order selection. We apply the proposed method on the simulated data and show that it significantly improves the accuracy of order selection from dependent data. We also perform order selection on fMRI data from a visuomotor task and show that the proposed method alleviates the over-estimation on the number of brain sources due to the intrinsic smoothness and the smooth preprocessing of fMRI data. We use the software package ICASSO (Himberg et al. [ 2004]: Neuroimage 22:1214-1222) to analyze the independent component (IC) estimates at different orders and show that, when ICA is performed at overestimated orders, the stability of the IC estimates decreases and the estimation of task related brain activations show degradation.

Treatment with Olanzapine is Associated with Modulation of the Default Mode Network in Patients with Schizophrenia

“…Given that in spatial ICA, the individual subject TC is assumed to be constant across the entire brain, the voxel-wise estimated signal deviations from the average network TC reflects local changes in the strength of functional connectivity within a given network. The minimum description length criteria were used to estimate the order selection that is the number of ICs from the smoothed data sets after taking into account the spatial and temporal correlation of the fMRI data (Li et al, 2007). These components were then used to back-reconstruct to individual ICs using the aggregate mixing matrix created during the reduction steps of dimensionality data.…”

EEG microstates as a tool for studying the temporal dynamics of whole-brain neuronal networks: A review

“…(Li et al, 2007)). Originally, Pascual-Marqui and colleagues (Pascual-Marqui et al, 1995) proposed a cross-validation criterion for selecting the optimal number of cluster maps, which optimizes the ratio between the global explained variance and the degrees of freedom for a given set of cluster maps.…”

Functional Genetic Variation of the Cannabinoid Receptor 1 and Cannabis Use Interact on Prefrontal Connectivity and Related Working Memory Behavior

“…fMRI datasets were split into a final set of 28 spatially independent components (ICs), using a modified minimum description length algorithm (Li et al, 2007). This is a stochastic process and the end results are not always identical.…”