Voxel selection in fMRI data analysis based on sparse representation

Multivariate pattern analysis approaches toward detection of brain regions from fMRI data have been gaining attention recently. In this study, we introduce an iterative sparse-representation-based algorithm for detection of voxels in functional MRI (fMRI) data with task relevant information. In each iteration of the algorithm, a linear programming problem is solved and a sparse weight vector is subsequently obtained. The final weight vector is the mean of those obtained in all iterations. The characteristics of our algorithm are as follows: 1) the weight vector (output) is sparse; 2) the magnitude of each entry of the weight vector represents the significance of its corresponding variable or feature in a classification or regression problem; and 3) due to the convergence of this algorithm, a stable weight vector is obtained. To demonstrate the validity of our algorithm and illustrate its application, we apply the algorithm to the Pittsburgh Brain Activity Interpretation Competition 2007 functional fMRI dataset for selecting the voxels, which are the most relevant to the tasks of the subjects. Based on this dataset, the aforementioned characteristics of our algorithm are analyzed, and a comparison between our method with the univariate general-linear-model-based statistical parametric mapping is performed. Using our method, a combination of voxels are selected based on the principle of effective/sparse representation of a task. Data analysis results in this paper show that this combination of voxels is suitable for decoding tasks and demonstrate the effectiveness of our method

Auto detection in autism

Autism is a neurobiological disorder in which, certain regions of the brain are affected. The main features of autism are impairment in communication, social interaction, language and deficit in imitation and theory of mind. Using Functional Magnetic Resonance Imaging (fMRI), haemodynamic responses during a bilateral finger tapping task are analyzed for both autistic subjects and normal control subjects. fMRI is a noninvasive technique to image the activity of the brain related to a specific task. Generally, the active voxels in the IMRI images are detected using parametric or non-parametric statistical methods in which the fMRI response is assumed to have a model. Such methods are not applicable to detect the active voxels when the fMRI response is unknown. The data driven methods are also used for analyzing the fMRI data. The data driven methods are computationally expensive. In this study, a method for detecting activated voxels without using prior knowledge of the input stimulus is presented. The assumption in this method is that the activation typically involves larger region comprising of several voxels and that these neighboring activated voxels are also temporally correlated. To validate the accuracy of this method, Principal component Analysis and Independent Component Analysis are also performed. A significant overlap in the sensorimotor cortex is found between the various methods suggesting that the automatic detecting method presented does provide accurate detection and localization

Ultrahigh Field Functional Magnetic Resonance Electrical Impedance Tomography (fMREIT) in Neural Activity Imaging

abstract: A direct Magnetic Resonance (MR)-based neural activity mapping technique with high spatial and temporal resolution may accelerate studies of brain functional organization. The most widely used technique for brain functional imaging is functional Magnetic Resonance Image (fMRI). The spatial resolution of fMRI is high. However, fMRI signals are highly influenced by the vasculature in each voxel and can be affected by capillary orientation and vessel size. Functional MRI analysis may, therefore, produce misleading results when voxels are nearby large vessels. Another problem in fMRI is that hemodynamic responses are slower than the neuronal activity. Therefore, temporal resolution is limited in fMRI. Furthermore, the correlation between neural activity and the hemodynamic response is not fully understood. fMRI can only be considered an indirect method of functional brain imaging. Another MR-based method of functional brain mapping is neuronal current magnetic resonance imaging (ncMRI), which has been studied over several years. However, the amplitude of these neuronal current signals is an order of magnitude smaller than the physiological noise. Works on ncMRI include simulation, phantom experiments, and studies in tissue including isolated ganglia, optic nerves, and human brains. However, ncMRI development has been hampered due to the extremely small signal amplitude, as well as the presence of confounding signals from hemodynamic changes and other physiological noise. Magnetic Resonance Electrical Impedance Tomography (MREIT) methods could have the potential for the detection of neuronal activity. In this technique, small external currents are applied to a body during MR scans. This current flow produces a magnetic field as well as an electric field. The altered magnetic flux density along the main magnetic field direction caused by this current flow can be obtained from phase images. When there is neural activity, the conductivity of the neural cell membrane changes and the current paths around the neurons change consequently. Neural spiking activity during external current injection, therefore, causes differential phase accumulation in MR data. Statistical analysis methods can be used to identify neuronal-current-induced magnetic field changes.Dissertation/ThesisDoctoral Dissertation Biomedical Engineering 201

Change detection in children with autism: an auditory event-related fMRI study

Autism involves impairments in communication and social interaction, as well as high levels of repetitive, stereotypic and ritualistic behaviours, and extreme resistance to change. This latter dimension, whilst required for a diagnosis, has received less research attention. We hypothesise that this extreme resistance to change in autism is rooted in atypical processing of unexpected stimuli. We tested this using auditory event-related fMRI to determine regional brain activity associated with passive detection of infrequently occurring frequency-deviant and complex novel sounds in a no-task condition. Participants were twelve 10 to 15-year-old children with autism, and a group of 12 age- and sex-matched healthy controls

fMRI activation detection with EEG priors

The purpose of brain mapping techniques is to advance the understanding of the relationship between structure and function in the human brain in so-called activation studies. In this work, an advanced statistical model for combining functional magnetic resonance imaging (fMRI) and electroencephalography (EEG) recordings is developed to fuse complementary information about the location of neuronal activity. More precisely, a new Bayesian method is proposed for enhancing fMRI activation detection by the use of EEG-based spatial prior information in stimulus based experimental paradigms. I.e., we model and analyse stimulus influence by a spatial Bayesian variable selection scheme, and extend existing high-dimensional regression methods by incorporating prior information on binary selection indicators via a latent probit regression with either a spatially-varying or constant EEG effect. Spatially-varying effects are regularized by intrinsic Markov random field priors. Inference is based on a full Bayesian Markov Chain Monte Carlo (MCMC) approach. Whether the proposed algorithm is able to increase the sensitivity of mere fMRI models is examined in both a real-world application and a simulation study. We observed, that carefully selected EEG--prior information additionally increases sensitivity in activation regions that have been distorted by a low signal-to-noise ratio

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