title: Gaussian process based independent analysis for temporal source separation in fMRI.
abstract: Functional Magnetic Resonance Imaging (fMRI) gives us a unique insight into the processes of the brain, and opens up for analyzing the functional activation patterns of the underlying sources. Task-inferred supervised learning with restrictive assumptions in the regression set-up, restricts the exploratory nature of the analysis. Fully unsupervised independent component analysis (ICA) algorithms, on the other hand, can struggle to detect clear classifiable components on single-subject data. We attribute this shortcoming to inadequate modeling of the fMRI source signals by failing to incorporate its temporal nature. fMRI source signals, biological stimuli and non-stimuli-related artifacts are all smooth over a time-scale compatible with the sampling time (TR). We therefore propose Gaussian process ICA (GPICA), which facilitates temporal dependency by the use of Gaussian process source priors. On two fMRI data sets with different sampling frequency, we show that the GPICA-inferred temporal components and associated spatial maps allow for a more definite interpretation than standard temporal ICA methods. The temporal structures of the sources are controlled by the covariance of the Gaussian process, specified by a kernel function with an interpretable and controllable temporal length scale parameter. We propose a hierarchical model specification, considering both instantaneous and convolutive mixing, and we infer source spatial maps, temporal patterns and temporal length scale parameters by Markov Chain Monte Carlo. A companion implementation made as a plug-in for SPM can be downloaded from https://github.com/dittehald/GPICA.
authors.affiliation: DTU Compute B321, Technical University of Denmark, DK-2800 Lyngby, Denmark. Electronic address: ditha@dtu.dk.
authors.forename: Ditte Høvenhoff
authors.initials: DH
authors.lastname: Hald
authors.affiliation: DTU Compute B321, Technical University of Denmark, DK-2800 Lyngby, Denmark; Electrical and Computer Engineering, Duke University, Durham, NC 27708, USA.
authors.forename: Ricardo
authors.initials: R
authors.lastname: Henao
authors.affiliation: DTU Compute B321, Technical University of Denmark, DK-2800 Lyngby, Denmark.
authors.forename: Ole
authors.initials: O
authors.lastname: Winther
mesh_terms: D000465:Algorithms; D016477:Artifacts; D001921:Brain / Q000502:physiology*; D001931:Brain Mapping / Q000379:methods*; D006801:Humans; D007089:Image Enhancement; D008279:Magnetic Resonance Imaging; D009010:Monte Carlo Method; D016011:Normal Distribution; D012815:Signal Processing, Computer-Assisted
keywords: Bayesian inference; Convolutive mixing; FMRI; Gaussian processes; Independent component analysis; Source separation