On applicability of PCA, voxel-wise variance normalization and dimensionality assumptions for sliding temporal window sICA in resting-state fMRI.

Subject-level resting-state fMRI (RS-fMRI) spatial independent component analysis (sICA) may provide new ways to analyze the data when performed in the sliding time window. However, whether principal component analysis (PCA) and voxel-wise variance normalization (VN) are applicable pre-processing procedures in the sliding-window context, as they are for regular sICA, has not been addressed so far. Also model order selection requires further studies concerning sliding-window sICA. In this paper we have addressed these concerns. First, we compared PCA-retained subspaces concerning overlapping parts of consecutive temporal windows to answer whether in-window PCA and VN can confound comparisons between sICA analyses in consecutive windows. Second, we compared the PCA subspaces between windowed and full data to assess expected comparability between windowed and full-data sICA results. Third, temporal evolution of dimensionality estimates in RS-fMRI data sets was monitored to identify potential challenges in model order selection in a sliding-window sICA context. Our results illustrate that in-window VN can be safely used, in-window PCA is applicable with most window widths and that comparisons between windowed and full data should not be performed from a subspace similarity point of view. In addition, our studies on dimensionality estimates demonstrated that there are sustained, periodic and very case-specific changes in signal-to-noise ratio within RS-fMRI data sets. Consequently, dimensionality estimation is needed for well-founded model order determination in the sliding-window case. The observed periodic changes correspond to a frequency band of ≤0.1 Hz, which is commonly associated with brain activity in RS-fMRI and become on average most pronounced at window widths of 80 and 60 time points (144 and 108 s, respectively). Wider windows provided only slightly better comparability between consecutive windows, and 60 time point or shorter windows also provided the best comparability with full-data results. Further studies are needed to determine the cause for dimensionality variations.

[1]  Rainer Goebel,et al.  Real-time independent component analysis of fMRI time-series , 2003, NeuroImage.

[2]  Rajesh Nandy,et al.  Estimation of the intrinsic dimensionality of fMRI data , 2006, NeuroImage.

[3]  Tülay Adali,et al.  Estimating the number of independent components for functional magnetic resonance imaging data , 2007, Human brain mapping.

[4]  Xiaoping Xie,et al.  Spatiotemporal nonlinearity in resting-state fMRI of the human brain , 2008, NeuroImage.

[5]  Vesa Kiviniemi,et al.  A Sliding Time-Window ICA Reveals Spatial Variability of the Default Mode Network in Time , 2011, Brain Connect..

[6]  O. Tervonen,et al.  Neuroglial Plasticity at Striatal Glutamatergic Synapses in Parkinson's Disease , 2011, Front. Syst. Neurosci..

[7]  Samuel Kaski,et al.  Dependencies between stimuli and spatially independent fMRI sources: Towards brain correlates of natural stimuli , 2009, NeuroImage.

[8]  Catie Chang,et al.  Time–frequency dynamics of resting-state brain connectivity measured with fMRI , 2010, NeuroImage.

[9]  I Daubechies,et al.  Independent component analysis for brain fMRI does not select for independence , 2009 .

[10]  Christopher M. Bishop,et al.  Mixtures of Probabilistic Principal Component Analyzers , 1999, Neural Computation.

[11]  Tom Minka,et al.  Automatic Choice of Dimensionality for PCA , 2000, NIPS.

[12]  G. Glover,et al.  Resting-State Functional Connectivity in Major Depression: Abnormally Increased Contributions from Subgenual Cingulate Cortex and Thalamus , 2007, Biological Psychiatry.

[13]  Xiangyu Long,et al.  Functional segmentation of the brain cortex using high model order group PICA , 2009, Human brain mapping.

[14]  Mark W. Woolrich,et al.  Advances in functional and structural MR image analysis and implementation as FSL , 2004, NeuroImage.

[15]  Stephen M. Smith,et al.  Probabilistic independent component analysis for functional magnetic resonance imaging , 2004, IEEE Transactions on Medical Imaging.

[16]  Jean-Baptiste Poline,et al.  A group model for stable multi-subject ICA on fMRI datasets , 2010, NeuroImage.

[17]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[18]  O. Tervonen,et al.  The effect of model order selection in group PICA , 2010, Human brain mapping.

[19]  M. Fox,et al.  Spontaneous fluctuations in brain activity observed with functional magnetic resonance imaging , 2007, Nature Reviews Neuroscience.

[20]  Xi-Nian Zuo,et al.  Reliable intrinsic connectivity networks: Test–retest evaluation using ICA and dual regression approach , 2010, NeuroImage.

[21]  Christian F. Beckmann,et al.  Modelling with independent components , 2012, NeuroImage.

[22]  Riitta Hari,et al.  Towards natural stimulation in fMRI—Issues of data analysis , 2007, NeuroImage.

[23]  Stephen M Smith,et al.  Correspondence of the brain's functional architecture during activation and rest , 2009, Proceedings of the National Academy of Sciences.

[24]  Fabian J. Theis,et al.  Spatial ICA of fMRI data in time windows , 2004 .

[25]  Stephen M. Smith,et al.  Investigations into resting-state connectivity using independent component analysis , 2005, Philosophical Transactions of the Royal Society B: Biological Sciences.

[26]  Thomas Kailath,et al.  Detection of signals by information theoretic criteria , 1985, IEEE Trans. Acoust. Speech Signal Process..

[27]  David Borsook,et al.  A Simple and Objective Method for Reproducible Resting State Network (RSN) Detection in fMRI , 2011, PloS one.

[28]  Li Yao,et al.  An Empirical Comparison of Information-Theoretic Criteria in Estimating the Number of Independent Components of fMRI Data , 2011, PloS one.

[29]  Xu Chen,et al.  Dimensionality estimation for optimal detection of functional networks in BOLD fMRI data , 2011, NeuroImage.

[30]  Christian Windischberger,et al.  Toward discovery science of human brain function , 2010, Proceedings of the National Academy of Sciences.

[31]  Olli Silvén,et al.  Mapping Transient Hyperventilation Induced Alterations with Estimates of the Multi-Scale Dynamics of BOLD Signal , 2009, Front. Neuroinform..

[32]  Xiaoping Xie,et al.  Estimating intrinsic dimensionality of fMRI dataset incorporating an AR(1) noise model with cubic spline interpolation , 2009, Neurocomputing.

[33]  J. Rissanen,et al.  Modeling By Shortest Data Description* , 1978, Autom..