3,923 research outputs found
PEAR: PEriodic And fixed Rank separation for fast fMRI
In functional MRI (fMRI), faster acquisition via undersampling of data can
improve the spatial-temporal resolution trade-off and increase statistical
robustness through increased degrees-of-freedom. High quality reconstruction of
fMRI data from undersampled measurements requires proper modeling of the data.
We present an fMRI reconstruction approach based on modeling the fMRI signal as
a sum of periodic and fixed rank components, for improved reconstruction from
undersampled measurements. We decompose the fMRI signal into a component which
a has fixed rank and a component consisting of a sum of periodic signals which
is sparse in the temporal Fourier domain. Data reconstruction is performed by
solving a constrained problem that enforces a fixed, moderate rank on one of
the components, and a limited number of temporal frequencies on the other. Our
approach is coined PEAR - PEriodic And fixed Rank separation for fast fMRI.
Experimental results include purely synthetic simulation, a simulation with
real timecourses and retrospective undersampling of a real fMRI dataset.
Evaluation was performed both quantitatively and visually versus ground truth,
comparing PEAR to two additional recent methods for fMRI reconstruction from
undersampled measurements. Results demonstrate PEAR's improvement in estimating
the timecourses and activation maps versus the methods compared against at
acceleration ratios of R=8,16 (for simulated data) and R=6.66,10 (for real
data). PEAR results in reconstruction with higher fidelity than when using a
fixed-rank based model or a conventional Low-rank+Sparse algorithm. We have
shown that splitting the functional information between the components leads to
better modeling of fMRI, over state-of-the-art methods
Disambiguating the role of blood flow and global signal with partial information decomposition
Global signal (GS) is an ubiquitous construct in resting state functional magnetic resonance imaging (rs-fMRI), associated to nuisance, but containing by definition most of the neuronal signal. Global signal regression (GSR) effectively removes the impact of physiological noise and other artifacts, but at the same time it alters correlational patterns in unpredicted ways. Performing GSR taking into account the underlying physiology (mainly the blood arrival time) has been proven to be beneficial. From these observations we aimed to: 1) characterize the effect of GSR on network-level functional connectivity in a large dataset; 2) assess the complementary role of global signal and vessels; and 3) use the framework of partial information decomposition to further look into the joint dynamics of the global signal and vessels, and their respective influence on the dynamics of cortical areas. We observe that GSR affects intrinsic connectivity networks in the connectome in a non-uniform way. Furthermore, by estimating the predictive information of blood flow and the global signal using partial information decomposition, we observe that both signals are present in different amounts across intrinsic connectivity networks. Simulations showed that differences in blood arrival time can largely explain this phenomenon, while using hemodynamic and calcium mouse recordings we were able to confirm the presence of vascular effects, as calcium recordings lack hemodynamic information. With these results we confirm network-specific effects of GSR and the importance of taking blood flow into account for improving de-noising methods. Additionally, and beyond the mere issue of data denoising, we quantify the diverse and complementary effect of global and vessel BOLD signals on the dynamics of cortical areas
A simple and objective method for reproducible resting state network (RSN) detection in fMRI
Spatial Independent Component Analysis (ICA) decomposes the time by space
functional MRI (fMRI) matrix into a set of 1-D basis time courses and their
associated 3-D spatial maps that are optimized for mutual independence. When
applied to resting state fMRI (rsfMRI), ICA produces several spatial
independent components (ICs) that seem to have biological relevance - the
so-called resting state networks (RSNs). The ICA problem is well posed when the
true data generating process follows a linear mixture of ICs model in terms of
the identifiability of the mixing matrix. However, the contrast function used
for promoting mutual independence in ICA is dependent on the finite amount of
observed data and is potentially non-convex with multiple local minima. Hence,
each run of ICA could produce potentially different IC estimates even for the
same data. One technique to deal with this run-to-run variability of ICA was
proposed by Yang et al. (2008) in their algorithm RAICAR which allows for the
selection of only those ICs that have a high run-to-run reproducibility. We
propose an enhancement to the original RAICAR algorithm that enables us to
assign reproducibility p-values to each IC and allows for an objective
assessment of both within subject and across subjects reproducibility. We call
the resulting algorithm RAICAR-N (N stands for null hypothesis test), and we
have applied it to publicly available human rsfMRI data (http://www.nitrc.org).
Our reproducibility analyses indicated that many of the published RSNs in
rsfMRI literature are highly reproducible. However, we found several other RSNs
that are highly reproducible but not frequently listed in the literature.Comment: 54 pages, 13 figure
A group model for stable multi-subject ICA on fMRI datasets
Spatial Independent Component Analysis (ICA) is an increasingly used
data-driven method to analyze functional Magnetic Resonance Imaging (fMRI)
data. To date, it has been used to extract sets of mutually correlated brain
regions without prior information on the time course of these regions. Some of
these sets of regions, interpreted as functional networks, have recently been
used to provide markers of brain diseases and open the road to paradigm-free
population comparisons. Such group studies raise the question of modeling
subject variability within ICA: how can the patterns representative of a group
be modeled and estimated via ICA for reliable inter-group comparisons? In this
paper, we propose a hierarchical model for patterns in multi-subject fMRI
datasets, akin to mixed-effect group models used in linear-model-based
analysis. We introduce an estimation procedure, CanICA (Canonical ICA), based
on i) probabilistic dimension reduction of the individual data, ii) canonical
correlation analysis to identify a data subspace common to the group iii)
ICA-based pattern extraction. In addition, we introduce a procedure based on
cross-validation to quantify the stability of ICA patterns at the level of the
group. We compare our method with state-of-the-art multi-subject fMRI ICA
methods and show that the features extracted using our procedure are more
reproducible at the group level on two datasets of 12 healthy controls: a
resting-state and a functional localizer study
Modeling Covariate Effects in Group Independent Component Analysis with Applications to Functional Magnetic Resonance Imaging
Independent component analysis (ICA) is a powerful computational tool for
separating independent source signals from their linear mixtures. ICA has been
widely applied in neuroimaging studies to identify and characterize underlying
brain functional networks. An important goal in such studies is to assess the
effects of subjects' clinical and demographic covariates on the spatial
distributions of the functional networks. Currently, covariate effects are not
incorporated in existing group ICA decomposition methods. Hence, they can only
be evaluated through ad-hoc approaches which may not be accurate in many cases.
In this paper, we propose a hierarchical covariate ICA model that provides a
formal statistical framework for estimating and testing covariate effects in
ICA decomposition. A maximum likelihood method is proposed for estimating the
covariate ICA model. We develop two expectation-maximization (EM) algorithms to
obtain maximum likelihood estimates. The first is an exact EM algorithm, which
has analytically tractable E-step and M-step. Additionally, we propose a
subspace-based approximate EM, which can significantly reduce computational
time while still retain high model-fitting accuracy. Furthermore, to test
covariate effects on the functional networks, we develop a voxel-wise
approximate inference procedure which eliminates the needs of computationally
expensive covariance estimation. The performance of the proposed methods is
evaluated via simulation studies. The application is illustrated through an
fMRI study of Zen meditation.Comment: 36 pages, 5 figure
Nonparametric Modeling of Dynamic Functional Connectivity in fMRI Data
Dynamic functional connectivity (FC) has in recent years become a topic of
interest in the neuroimaging community. Several models and methods exist for
both functional magnetic resonance imaging (fMRI) and electroencephalography
(EEG), and the results point towards the conclusion that FC exhibits dynamic
changes. The existing approaches modeling dynamic connectivity have primarily
been based on time-windowing the data and k-means clustering. We propose a
non-parametric generative model for dynamic FC in fMRI that does not rely on
specifying window lengths and number of dynamic states. Rooted in Bayesian
statistical modeling we use the predictive likelihood to investigate if the
model can discriminate between a motor task and rest both within and across
subjects. We further investigate what drives dynamic states using the model on
the entire data collated across subjects and task/rest. We find that the number
of states extracted are driven by subject variability and preprocessing
differences while the individual states are almost purely defined by either
task or rest. This questions how we in general interpret dynamic FC and points
to the need for more research on what drives dynamic FC.Comment: 8 pages, 1 figure. Presented at the Machine Learning and
Interpretation in Neuroimaging Workshop (MLINI-2015), 2015 (arXiv:1605.04435
- …