Activated region fitting: a robust high-power method for fMRI analysis using parameterized regions of activation

3D spatially-adaptive canonical correlation analysis: Local and global methods

Local spatially-adaptive canonical correlation analysis (local CCA) with spatial constraints has been introduced to fMRI multivariate analysis for improved modeling of activation patterns. However, current algorithms require complicated spatial constraints that have only been applied to 2D local neighborhoods because the computational time would be exponentially increased if the same method is applied to 3D spatial neighborhoods. In this study, an efficient and accurate line search sequential quadratic programming (SQP) algorithm has been developed to efficiently solve the 3D local CCA problem with spatial constraints. In addition, a spatially-adaptive kernel CCA (KCCA) method is proposed to increase accuracy of fMRI activation maps. With oriented 3D spatial filters anisotropic shapes can be estimated during the KCCA analysis of fMRI time courses. These filters are orientation-adaptive leading to rotational invariance to better match arbitrary oriented fMRI activation patterns, resulting in improved sensitivity of activation detection while significantly reducing spatial blurring artifacts. The kernel method in its basic form does not require any spatial constraints and analyzes the whole-brain fMRI time series to construct an activation map. Finally, we have developed a penalized kernel CCA model that involves spatial low-pass filter constraints to increase the specificity of the method. The kernel CCA methods are compared with the standard univariate method and with two different local CCA methods that were solved by the SQP algorithm. Results show that SQP is the most efficient algorithm to solve the local constrained CCA problem, and the proposed kernel CCA methods outperformed univariate and local CCA methods in detecting activations for both simulated and real fMRI episodic memory data.

Characterizing cross-subject spatial interaction patterns in functional magnetic resonance imaging studies: A two-stage point-process model

We develop a two-stage spatial point process model that introduces new characterizations of activation patterns in multisubject functional Magnetic Resonance Imaging (fMRI) studies. Conventionally multisubject fMRI methods rely on combining information across subjects one voxel at a time in order to identify locations of peak activation in the brain. The two-stage model that we develop here addresses shortcomings of standard methods by explicitly modeling the spatial structure of functional signals and recognizing that corresponding cross-subject functional signals can be spatially misaligned. In our first stage analysis, we introduce a marked spatial point process model that captures the spatial features of the functional response and identifies a configuration of activation units for each subject. The locations of these activation units are used as input for the second stage model. The point process model of the second stage analysis is developed to characterize multisubject activation patterns by estimating the strength of cross-subject interactions at different spatial ranges. The model uses spatial neighborhoods to account for the cross-subject spatial misalignment in corresponding functional units. We applied our methods to an fMRI study of 21 individuals who performed an attention test. We identified four brain regions that are involved in the test and found that our model results agree well with our understanding of how these regions engage with the tasks performed during the attention test. Our results highlighted that cross-subject interactions are stronger in brain areas that have a more specific function in performing the experimental tasks than in other areas.

Multiple testing corrections, nonparametric methods, and random field theory

I provide a selective review of the literature on the multiple testing problem in fMRI. By drawing connections with the older modalities, PET in particular, and how software implementations have tracked (or lagged behind) theoretical developments, my narrative aims to give the methodological researcher a historical perspective on this important aspect of fMRI data analysis.

Optimizing the performance of local canonical correlation analysis in fMRI using spatial constraints

The benefits of locally adaptive statistical methods for fMRI research have been shown in recent years, as these methods are more proficient in detecting brain activations in a noisy environment. One such method is local canonical correlation analysis (CCA), which investigates a group of neighboring voxels instead of looking at the single voxel time course. The value of a suitable test statistic is used as a measure of activation. It is customary to assign the value to the center voxel for convenience. The method without constraints is prone to artifacts, especially in a region of localized strong activation. To compensate for these deficiencies, the impact of different spatial constraints in CCA on sensitivity and specificity are investigated. The ability of constrained CCA (cCCA) to detect activation patterns in an episodic memory task has been studied. This research shows how any arbitrary contrast of interest can be analyzed by cCCA and how accurate P-values optimized for the contrast of interest can be computed using nonparametric methods. Results indicate an increase of up to 20% in detecting activation patterns for some of the advanced cCCA methods, as measured by ROC curves derived from simulated and real fMRI data.

Image analysis and statistical inference in neuroimaging with R

R is a language and environment for statistical computing and graphics. It can be considered an alternative implementation of the S language developed in the 1970s and 1980s for data analysis and graphics (Becker and Chambers, 1984; Becker et al., 1988). The R language is part of the GNU project and offers versions that compile and run on almost every major operating system currently available. We highlight several R packages built specifically for the analysis of neuroimaging data in the context of functional MRI, diffusion tensor imaging, and dynamic contrast-enhanced MRI. We review their methodology and give an overview of their capabilities for neuroimaging. In addition we summarize some of the current activities in the area of neuroimaging software development in R.

Functional connectivity analysis of fMRI data using parameterized regions-of-interest

Connectivity analysis of fMRI data requires correct specification of regions-of-interest (ROIs). Selection of ROIs based on outcomes of a GLM analysis may be hindered by conservativeness of the multiple comparison correction, while selection based on brain anatomy may be biased due to inconsistent structure-to-function mapping. To alleviate these problems we propose a method to define functional ROIs without the need for a stringent multiple comparison correction. We extend a flexible framework for fMRI analysis (Activated Region Fitting, Weeda et al. 2009) to connectivity analysis of fMRI data. This method describes an entire fMRI data volume by regions of activation defined by a limited number of parameters. Therefore a less stringent multiple comparison procedure is required. The regions of activation from this analysis can be directly used to estimate functional connectivity. Simulations show that Activated Region Fitting can recover the connectivity of brain regions. An application to real data of a Go/No-Go experiment highlights the advantages of the method.

Robust and unbiased variance of GLM coefficients for misspecified autocorrelation and hemodynamic response models in fMRI

As a consequence of misspecification of the hemodynamic response and noise variance models, tests on general linear model coe cients are not valid. Robust estimation of the variance of the general linear model (GLM) coecients in fMRI time series is therefore essential. In this paper an alternative method to estimate the variance of the GLM coe cients accurately is suggested and compared to other methods. The alternative, referred to as the sandwich, is based primarily on the fact that the time series are obtained from multiple exchangeable stimulus presentations. The analytic results show that the sandwich is unbiased. Using this result, it is possible to obtain an exact statistic which keeps the 5% false positive rate. Extensive Monte Carlo simulations show that the sandwich is robust against misspeci cation of the autocorrelations and of the hemodynamic response model. The sandwich is seen to be in many circumstances robust, computationally efficient, and flexible with respect to correlation structures across the brain. In contrast, the smoothing approach can be robust to a certain extent but only with specific knowledge of the circumstances for the smoothing parameter.