| 1. |
- Abramian, David, 1992-, et al.
(författare)
-
Diffusion-Informed Spatial Smoothing of fMRI Data in White Matter Using Spectral Graph Filters
- 2021
-
Ingår i: NeuroImage. - : Elsevier. - 1053-8119 .- 1095-9572. ; 237
-
Tidskriftsartikel (refereegranskat)abstract
- Brain activation mapping using functional magnetic resonance imaging (fMRI) has been extensively studied in brain gray matter (GM), whereas in large disregarded for probing white matter (WM). This unbalanced treatment has been in part due to controversies in relation to the nature of the blood oxygenation level-dependent (BOLD) contrast in WM and its detachability. However, an accumulating body of studies has provided solid evidence of the functional significance of the BOLD signal in WM and has revealed that it exhibits anisotropic spatio-temporal correlations and structure-specific fluctuations concomitant with those of the cortical BOLD signal. In this work, we present an anisotropic spatial filtering scheme for smoothing fMRI data in WM that accounts for known spatial constraints on the BOLD signal in WM. In particular, the spatial correlation structure of the BOLD signal in WM is highly anisotropic and closely linked to local axonal structure in terms of shape and orientation, suggesting that isotropic Gaussian filters conventionally used for smoothing fMRI data are inadequate for denoising the BOLD signal in WM. The fundamental element in the proposed method is a graph-based description of WM that encodes the underlying anisotropy observed across WM, derived from diffusion-weighted MRI data. Based on this representation, and leveraging graph signal processing principles, we design subject-specific spatial filters that adapt to a subject’s unique WM structure at each position in the WM that they are applied at. We use the proposed filters to spatially smooth fMRI data in WM, as an alternative to the conventional practice of using isotropic Gaussian filters. We test the proposed filtering approach on two sets of simulated phantoms, showcasing its greater sensitivity and specificity for the detection of slender anisotropic activations, compared to that achieved with isotropic Gaussian filters. We also present WM activation mapping results on the Human Connectome Project’s 100-unrelated subject dataset, across seven functional tasks, showing that the proposed method enables the detection of streamline-like activations within axonal bundles.
|
|
| 2. |
- Eklund, Anders, et al.
(författare)
-
A Bayesian Heteroscedastic GLM with Application to fMRI Data with Motion Spikes
- 2017
-
Ingår i: NeuroImage. - : Elsevier. - 1053-8119 .- 1095-9572. ; 155, s. 354-369
-
Tidskriftsartikel (refereegranskat)abstract
- We propose a voxel-wise general linear model with autoregressive noise and heteroscedastic noise innovations (GLMH) for analyzing functional magnetic resonance imaging (fMRI) data. The model is analyzed from a Bayesian perspective and has the benefit of automatically down-weighting time points close to motion spikes in a data-driven manner. We develop a highly efficient Markov Chain Monte Carlo (MCMC) algorithm that allows for Bayesian variable selection among the regressors to model both the mean (i.e., the design matrix) and variance. This makes it possible to include a broad range of explanatory variables in both the mean and variance (e.g., time trends, activation stimuli, head motion parameters and their temporal derivatives), and to compute the posterior probability of inclusion from the MCMC output. Variable selection is also applied to the lags in the autoregressive noise process, making it possible to infer the lag order from the data simultaneously with all other model parameters. We use both simulated data and real fMRI data from OpenfMRI to illustrate the importance of proper modeling of heteroscedasticity in fMRI data analysis. Our results show that the GLMH tends to detect more brain activity, compared to its homoscedastic counterpart, by allowing the variance to change over time depending on the degree of head motion.
|
|
| 3. |
- Eklund, Anders, et al.
(författare)
-
Does Parametric fMRI Analysis with SPM Yield Valid Results? - An Empirical Study of 1484 Rest Datasets
- 2012
-
Ingår i: NeuroImage. - : Elsevier. - 1053-8119 .- 1095-9572. ; 61:3, s. 565-578
-
Tidskriftsartikel (refereegranskat)abstract
- The validity of parametric functional magnetic resonance imaging (fMRI) analysis has only been reported for simulated data.Recent advances in computer science and data sharing make it possible to analyze large amounts of real fMRI data. In this study,1484 rest datasets have been analyzed in SPM8, to estimate true familywise error rates. For a familywise significance threshold of5%, significant activity was found in 1% - 70% of the 1484 rest datasets, depending on repetition time, paradigm and parametersettings. This means that parametric significance thresholds in SPM both can be conservative or very liberal. The main reason forthe high familywise error rates seems to be that the global AR(1) auto correlation correction in SPM fails to model the spectra ofthe residuals, especially for short repetition times. The findings that are reported in this study cannot be generalized to parametricfMRI analysis in general, other software packages may give different results. By using the computational power of the graphicsprocessing unit (GPU), the 1484 rest datasets were also analyzed with a random permutation test. Significant activity was thenfound in 1% - 19% of the datasets. These findings speak to the need for a better model of temporal correlations in fMRI timeseries.
|
|
| 4. |
- Sidén, Per, et al.
(författare)
-
Fast Bayesian whole-brain fMRI analysis with spatial 3D priors
- 2017
-
Ingår i: NeuroImage. - : Elsevier. - 1053-8119 .- 1095-9572. ; 146, s. 211-225
-
Tidskriftsartikel (refereegranskat)abstract
- Spatial whole-brain Bayesian modeling of task-related functional magnetic resonance imaging (fMRI) is a great computational challenge. Most of the currently proposed methods therefore do inference in subregions of the brain separately or do approximate inference without comparison to the true posterior distribution. A popular such method, which is now the standard method for Bayesian single subject analysis in the SPM software, is introduced in Penny et al. (2005b). The method processes the data slice-by-slice and uses an approximate variational Bayes (VB) estimation algorithm that enforces posterior independence between activity coefficients in different voxels. We introduce a fast and practical Markov chain Monte Carlo (MCMC) scheme for exact inference in the same model, both slice-wise and for the whole brain using a 3D prior on activity coefficients. The algorithm exploits sparsity and uses modern techniques for efficient sampling from high-dimensional Gaussian distributions, leading to speed-ups without which MCMC would not be a practical option. Using MCMC, we are for the first time able to evaluate the approximate VB posterior against the exact MCMC posterior, and show that VB can lead to spurious activation. In addition, we develop an improved VB method that drops the assumption of independent voxels a posteriori. This algorithm is shown to be much faster than both MCMC and the original VB for large datasets, with negligible error compared to the MCMC posterior.
|
|
| 5. |
- Sjölund, Jens, 1987-, et al.
(författare)
-
Bayesian uncertainty quantification in linear models for diffusion MRI
- 2018
-
Ingår i: NeuroImage. - : Elsevier BV. - 1053-8119 .- 1095-9572. ; 175, s. 272-285
-
Tidskriftsartikel (refereegranskat)abstract
- Diffusion MRI (dMRI) is a valuable tool in the assessment of tissue microstructure. By fitting a model to the dMRI signal it is possible to derive various quantitative features. Several of the most popular dMRI signal models are expansions in an appropriately chosen basis, where the coefficients are determined using some variation of least-squares. However, such approaches lack any notion of uncertainty, which could be valuable in e.g. group analyses. In this work, we use a probabilistic interpretation of linear least-squares methods to recast popular dMRI models as Bayesian ones. This makes it possible to quantify the uncertainty of any derived quantity. In particular, for quantities that are affine functions of the coefficients, the posterior distribution can be expressed in closed-form. We simulated measurements from single- and double-tensor models where the correct values of several quantities are known, to validate that the theoretically derived quantiles agree with those observed empirically. We included results from residual bootstrap for comparison and found good agreement. The validation employed several different models: Diffusion Tensor Imaging (DTI), Mean Apparent Propagator MRI (MAP-MRI) and Constrained Spherical Deconvolution (CSD). We also used in vivo data to visualize maps of quantitative features and corresponding uncertainties, and to show how our approach can be used in a group analysis to downweight subjects with high uncertainty. In summary, we convert successful linear models for dMRI signal estimation to probabilistic models, capable of accurate uncertainty quantification. © 2018 Elsevier Inc.
|
|
