2,880 research outputs found
Disentangling causal webs in the brain using functional Magnetic Resonance Imaging: A review of current approaches
In the past two decades, functional Magnetic Resonance Imaging has been used
to relate neuronal network activity to cognitive processing and behaviour.
Recently this approach has been augmented by algorithms that allow us to infer
causal links between component populations of neuronal networks. Multiple
inference procedures have been proposed to approach this research question but
so far, each method has limitations when it comes to establishing whole-brain
connectivity patterns. In this work, we discuss eight ways to infer causality
in fMRI research: Bayesian Nets, Dynamical Causal Modelling, Granger Causality,
Likelihood Ratios, LiNGAM, Patel's Tau, Structural Equation Modelling, and
Transfer Entropy. We finish with formulating some recommendations for the
future directions in this area
Serial Correlations in Single-Subject fMRI with Sub-Second TR
When performing statistical analysis of single-subject fMRI data, serial
correlations need to be taken into account to allow for valid inference.
Otherwise, the variability in the parameter estimates might be under-estimated
resulting in increased false-positive rates. Serial correlations in fMRI data
are commonly characterized in terms of a first-order autoregressive (AR)
process and then removed via pre-whitening. The required noise model for the
pre-whitening depends on a number of parameters, particularly the repetition
time (TR). Here we investigate how the sub-second temporal resolution provided
by simultaneous multislice (SMS) imaging changes the noise structure in fMRI
time series. We fit a higher-order AR model and then estimate the optimal AR
model order for a sequence with a TR of less than 600 ms providing whole brain
coverage. We show that physiological noise modelling successfully reduces the
required AR model order, but remaining serial correlations necessitate an
advanced noise model. We conclude that commonly used noise models, such as the
AR(1) model, are inadequate for modelling serial correlations in fMRI using
sub-second TRs. Rather, physiological noise modelling in combination with
advanced pre-whitening schemes enable valid inference in single-subject
analysis using fast fMRI sequences
Structure Learning in Coupled Dynamical Systems and Dynamic Causal Modelling
Identifying a coupled dynamical system out of many plausible candidates, each
of which could serve as the underlying generator of some observed measurements,
is a profoundly ill posed problem that commonly arises when modelling real
world phenomena. In this review, we detail a set of statistical procedures for
inferring the structure of nonlinear coupled dynamical systems (structure
learning), which has proved useful in neuroscience research. A key focus here
is the comparison of competing models of (ie, hypotheses about) network
architectures and implicit coupling functions in terms of their Bayesian model
evidence. These methods are collectively referred to as dynamical casual
modelling (DCM). We focus on a relatively new approach that is proving
remarkably useful; namely, Bayesian model reduction (BMR), which enables rapid
evaluation and comparison of models that differ in their network architecture.
We illustrate the usefulness of these techniques through modelling
neurovascular coupling (cellular pathways linking neuronal and vascular
systems), whose function is an active focus of research in neurobiology and the
imaging of coupled neuronal systems
Adaptive Gaussian Markov Random Fields with Applications in Human Brain Mapping
Functional magnetic resonance imaging (fMRI) has become the standard technology in human brain mapping. Analyses of the massive spatio-temporal fMRI data sets often focus on parametric or nonparametric modeling of the temporal component, while spatial smoothing is based on Gaussian kernels or random fields. A weakness of Gaussian spatial smoothing is underestimation of activation peaks or blurring of high-curvature transitions between activated and non-activated brain regions. In this paper, we introduce a class of inhomogeneous Markov random fields (MRF) with spatially adaptive interaction weights in a space-varying coefficient model for fMRI data. For given weights, the random field is conditionally Gaussian, but marginally it is non-Gaussian. Fully Bayesian inference, including estimation of weights and variance parameters, is carried out through efficient MCMC simulation. An application to fMRI data from a visual stimulation experiment demonstrates the performance of our approach in comparison to Gaussian and robustified non-Gaussian Markov random field models
Physiological Gaussian Process Priors for the Hemodynamics in fMRI Analysis
Background: Inference from fMRI data faces the challenge that the hemodynamic
system that relates neural activity to the observed BOLD fMRI signal is
unknown.
New Method: We propose a new Bayesian model for task fMRI data with the
following features: (i) joint estimation of brain activity and the underlying
hemodynamics, (ii) the hemodynamics is modeled nonparametrically with a
Gaussian process (GP) prior guided by physiological information and (iii) the
predicted BOLD is not necessarily generated by a linear time-invariant (LTI)
system. We place a GP prior directly on the predicted BOLD response, rather
than on the hemodynamic response function as in previous literature. This
allows us to incorporate physiological information via the GP prior mean in a
flexible way, and simultaneously gives us the nonparametric flexibility of the
GP.
Results: Results on simulated data show that the proposed model is able to
discriminate between active and non-active voxels also when the GP prior
deviates from the true hemodynamics. Our model finds time varying dynamics when
applied to real fMRI data.
Comparison with Existing Method(s): The proposed model is better at detecting
activity in simulated data than standard models, without inflating the false
positive rate. When applied to real fMRI data, our GP model in several cases
finds brain activity where previously proposed LTI models does not.
Conclusions: We have proposed a new non-linear model for the hemodynamics in
task fMRI, that is able to detect active voxels, and gives the opportunity to
ask new kinds of questions related to hemodynamics.Comment: 18 pages, 14 figure
Multimodal imaging of human brain activity: rational, biophysical aspects and modes of integration
Until relatively recently the vast majority of imaging and electrophysiological studies of human brain activity have relied on single-modality measurements usually correlated with readily observable or experimentally modified behavioural or brain state patterns. Multi-modal imaging is the concept of bringing together observations or measurements from different instruments. We discuss the aims of multi-modal imaging and the ways in which it can be accomplished using representative applications. Given the importance of haemodynamic and electrophysiological signals in current multi-modal imaging applications, we also review some of the basic physiology relevant to understanding their relationship
Time Series Analysis of fMRI Data: Spatial Modelling and Bayesian Computation
Time series analysis of fMRI data is an important area of medical statistics
for neuroimaging data. The neuroimaging community has embraced mean-field
variational Bayes (VB) approximations, which are implemented in Statistical
Parametric Mapping (SPM) software. While computationally efficient, the quality
of VB approximations remains unclear even though they are commonly used in the
analysis of neuroimaging data. For reliable statistical inference, it is
important that these approximations be accurate and that users understand the
scenarios under which they may not be accurate.
We consider this issue for a particular model that includes spatially-varying
coefficients. To examine the accuracy of the VB approximation we derive
Hamiltonian Monte Carlo (HMC) for this model and conduct simulation studies to
compare its performance with VB. As expected we find that the computation time
required for VB is considerably less than that for HMC. In settings involving a
high or moderate signal-to-noise ratio (SNR) we find that the two approaches
produce very similar results suggesting that the VB approximation is useful in
this setting. On the other hand, when one considers a low SNR, substantial
differences are found, suggesting that the approximation may not be accurate in
such cases and we demonstrate that VB produces Bayes estimators with larger
mean squared error (MSE). A real application related to face perception is also
carried out. Overall, our work clarifies the usefulness of VB for the
spatiotemporal analysis of fMRI data, while also pointing out the limitation of
VB when the SNR is low and the utility of HMC in this case
- …