An integrated neighborhood correlation and hierarchical clustering approach of functional MRI

A decomposition model and voxel selection framework for fMRI analysis to predict neural response of visual stimuli

This paper presents an fMRI signal analysis methodology using geometric mean curve decomposition (GMCD) and mutual information-based voxel selection framework. Previously, the fMRI signal analysis has been conducted using empirical mean curve decomposition (EMCD) model and voxel selection on raw fMRI signal. The erstwhile methodology loses frequency component, while the latter methodology suffers from signal redundancy. Both challenges are addressed by our methodology in which the frequency component is considered by decomposing the raw fMRI signal using geometric mean rather than arithmetic mean and the voxels are selected from EMCD signal using GMCD components, rather than raw fMRI signal. The proposed methodologies are adopted for predicting the neural response. Experimentations are conducted in the openly available fMRI data of six subjects, and comparisons are made with existing decomposition models and voxel selection frameworks. Subsequently, the effect of degree of selected voxels and the selection constraints are analyzed. The comparative results and the analysis demonstrate the superiority and the reliability of the proposed methodology.

Combining self-organizing mapping and supervised affinity propagation clustering approach to investigate functional brain networks involved in motor imagery and execution with fMRI measurements

Clustering analysis methods have been widely applied to identifying the functional brain networks of a multitask paradigm. However, the previously used clustering analysis techniques are computationally expensive and thus impractical for clinical applications. In this study a novel method, called SOM-SAPC that combines self-organizing mapping (SOM) and supervised affinity propagation clustering (SAPC), is proposed and implemented to identify the motor execution (ME) and motor imagery (MI) networks. In SOM-SAPC, SOM was first performed to process fMRI data and SAPC is further utilized for clustering the patterns of functional networks. As a result, SOM-SAPC is able to significantly reduce the computational cost for brain network analysis. Simulation and clinical tests involving ME and MI were conducted based on SOM-SAPC, and the analysis results indicated that functional brain networks were clearly identified with different response patterns and reduced computational cost. In particular, three activation clusters were clearly revealed, which include parts of the visual, ME and MI functional networks. These findings validated that SOM-SAPC is an effective and robust method to analyze the fMRI data with multitasks.

Quantification of the power changes in BOLD signals using Welch spectrum method during different single-hand motor imageries

Motor imagery is an experimental paradigm implemented in cognitive neuroscience and cognitive psychology. To investigate the asymmetry of the strength of cortical functional activity due to different single-hand motor imageries, functional magnetic resonance imaging (fMRI) data from right handed normal subjects were recorded and analyzed during both left-hand and right-hand motor imagery processes. Then the average power of blood oxygenation level-dependent (BOLD) signals in temporal domain was calculated using the developed tool that combines Welch power spectrum and the integral of power spectrum approach of BOLD signal changes during motor imagery. Power change analysis results indicated that cortical activity exhibited a stronger power in the precentral gyrus and medial frontal gyrus with left-hand motor imagery tasks compared with that from right-hand motor imagery tasks. These observations suggest that right handed normal subjects mobilize more cortical nerve cells for left-hand motor imagery. Our findings also suggest that the approach based on power differences of BOLD signals is a suitable quantitative analysis tool for quantification of asymmetry of brain activity intensity during motor imagery tasks.

Unsupervised spatiotemporal analysis of fMRI data using graph-based visualizations of self-organizing maps

We present novel graph-based visualizations of self-organizing maps for unsupervised functional magnetic resonance imaging (fMRI) analysis. A self-organizing map is an artificial neural network model that transforms high-dimensional data into a low-dimensional (often a 2-D) map using unsupervised learning. However, a postprocessing scheme is necessary to correctly interpret similarity between neighboring node prototypes (feature vectors) on the output map and delineate clusters and features of interest in the data. In this paper, we used graph-based visualizations to capture fMRI data features based upon 1) the distribution of data across the receptive fields of the prototypes (density-based connectivity); and 2) temporal similarities (correlations) between the prototypes (correlation-based connectivity). We applied this approach to identify task-related brain areas in an fMRI reaction time experiment involving a visuo-manual response task, and we correlated the time-to-peak of the fMRI responses in these areas with reaction time. Visualization of self-organizing maps outperformed independent component analysis and voxelwise univariate linear regression analysis in identifying and classifying relevant brain regions. We conclude that the graph-based visualizations of self-organizing maps help in advanced visualization of cluster boundaries in fMRI data enabling the separation of regions with small differences in the timings of their brain responses.

Analysis of FMRI data using an integrated principal component analysis and supervised affinity propagation clustering approach

Clustering analysis is a promising data-driven method for analyzing functional magnetic resonance imaging (fMRI) time series data. The huge computational load, however, creates practical difficulties for this technique. We present a novel approach, integrating principal component analysis (PCA) and supervised affinity propagation clustering (SAPC). In this method, fMRI data are initially processed by PCA to obtain a preliminary image of brain activation. SAPC is then used to detect different brain functional activation patterns. We used a supervised Silhouette index to optimize clustering quality and automatically search for the optimal parameter p in SAPC, so that the basic affinity propagation clustering is improved by applying SAPC. Four simulation studies and tests with three in vivo fMRI datasets containing data from both block-design and event-related experiments revealed that functional brain activation was effectively detected and different response patterns were distinguished using our integrated method. In addition, the improved SAPC method was superior to the k -centers clustering and hierarchical clustering methods in both block-design and event-related fMRI data, as measured by the average squared error. These results suggest that our proposed novel integrated approach will be useful for detecting brain functional activation in both block-design and event-related experimental fMRI data.

Analysis of activity in fMRI data using affinity propagation clustering

Clustering analysis is a promising data-driven method for the analysis of functional magnetic resonance imaging (fMRI) data. The huge computation load, however, makes it difficult for the practical use. We use affinity propagation clustering (APC), a new clustering algorithm especially for large data sets to detect brain functional activation from fMRI. It considers all data points as possible exemplars through the minimisation of an energy function and message-passing architecture, and obtains the optimal set of exemplars and their corresponding clusters. Four simulation studies and three in vivo fMRI studies reveal that brain functional activation can be effectively detected and that different response patterns can be distinguished using this method. Our results demonstrate that APC is superior to the k-centres clustering, as revealed by their performance measures in the weighted Jaccard coefficient and average squared error. These results suggest that the proposed APC will be useful in detecting brain functional activation from fMRI data.

Analysis of functional networks involved in motor execution and motor imagery using combined hierarchical clustering analysis and independent component analysis

Cognitive experiments involving motor execution (ME) and motor imagery (MI) have been intensively studied using functional magnetic resonance imaging (fMRI). However, the functional networks of a multitask paradigm which include ME and MI were not widely explored. In this article, we aimed to investigate the functional networks involved in MI and ME using a method combining the hierarchical clustering analysis (HCA) and the independent component analysis (ICA). Ten right-handed subjects were recruited to participate a multitask experiment with conditions such as visual cue, MI, ME and rest. The results showed that four activation clusters were found including parts of the visual network, ME network, the MI network and parts of the resting state network. Furthermore, the integration among these functional networks was also revealed. The findings further demonstrated that the combined HCA with ICA approach was an effective method to analyze the fMRI data of multitasks.

Convolution power spectrum analysis for FMRI data based on prior image signal

Functional MRI (fMRI) data-processing methods based on changes in the time domain involve, among other things, correlation analysis and use of the general linear model with statistical parametric mapping (SPM). Unlike conventional fMRI data analysis methods, which aim to model the blood-oxygen-level-dependent (BOLD) response of voxels as a function of time, the theory of power spectrum (PS) analysis focuses completely on understanding the dynamic energy change of interacting systems. We propose a new convolution PS (CPS) analysis of fMRI data, based on the theory of matched filtering, to detect brain functional activation for fMRI data. First, convolution signals are computed between the measured fMRI signals and the image signal of prior experimental pattern to suppress noise in the fMRI data. Then, the PS density analysis of the convolution signal is specified as the quantitative analysis energy index of BOLD signal change. The data from simulation studies and in vivo fMRI studies, including block-design experiments, reveal that the CPS method enables a more effective detection of some aspects of brain functional activation, as compared with the canonical PS SPM and the support vector machine methods. Our results demonstrate that the CPS method is useful as a complementary analysis in revealing brain functional information regarding the complex nature of fMRI time series.

A connectivity-based method for defining regions-of-interest in fMRI data

In this paper, we describe a new methodology for defining brain regions-of-interset (ROIs) in functional magnetic resonance imaging (fMRI) data. The ROIs are defined based on their functional connectivity to other ROIs, i.e., ROIs are defined as sets of voxels with similar connectivity patterns to other ROIs. The method relies on 1) a spatially regularized canonical correlation analysis for identifying maximally correlated signals, which are not due to correlated noise; 2) a test for merging ROIs which have similar connectivity patterns to the other ROIs; and 3) a graph-cuts optimization for assigning voxels to ROIs. Since our method is fully connectivity-based, the extracted ROIs and their corresponding time signals are ideally suited for a subsequent brain connectivity analysis.

Analysis of fMRI data using improved self-organizing mapping and spatio-temporal metric hierarchical clustering

The self-organizing mapping (SOM) and hierarchical clustering (HC) methods are integrated to detect brain functional activation; functional magnetic resonance imaging (fMRI) data are first processed by SOM to obtain a primary merged neural nodes image, and then by HC to obtain further brain activation patterns. The conventional Euclidean distance metric was replaced by the correlation distance metric in SOM to improve clustering and merging of neural nodes. To improve the use of spatial and temporal information in fMRI data, a new spatial distance (node coordinates in the 2-D lattice) and temporal correlation (correlation degree of each time course in the exemplar matrix) are introduced in HC to merge the primary SOM results. Two simulation studies and two in vivo fMRI data that both contained block-design and event-related experiments revealed that brain functional activation can be effectively detected and that different response patterns can be distinguished using these methods. Our results demonstrate that the improved SOM and HC methods are clearly superior to the statistical parametric mapping (SPM), independent component analysis (ICA), and conventional SOM methods in the block-design, especially in the event-related experiment, as revealed by their performance measured by receiver operating characteristic (ROC) analysis. Our results also suggest that the proposed new integrated approach could be useful in detecting block-design and event-related fMRI data.

Obstructive sleep apnea detection using clustering classification of nonlinear features from nocturnal oximetry

This study is focused on the classification of patients suspected of suffering from obstructive sleep apnea (OSA) by means of cluster analysis. We assessed the diagnostic ability of three clustering algorithms: k-means, hierarchical and fuzzy c-means (FCM). Nonlinear features of blood oxygen saturation (SaO2) from nocturnal oximetry were used as inputs to the clustering methods. Three nonlinear methods were used: approximate entropy (ApEn), central tendency measure (CTM) and Lempel-Ziv (LZ) complexity. A population of 74 subjects (44 OSA positive and 30 OSA negative) was studied. 90.5%, 87.8% and 86.5% accuracies were reached with k-means, hierarchical and FCM algorithms, respectively. The diagnostic accuracy values improved those obtained with each nonlinear method individually. Our results suggest that nonlinear analysis and clustering classification could provide useful information to help in the diagnosis of OSA syndrome.

fMRI data analysis with nonstationary noise models: a Bayesian approach

The assumption of noise stationarity in the functional magnetic resonance imaging (fMRI) data analysis may lead to the loss of crucial dynamic features of the data and thus result in inaccurate activation detection. In this paper, a Bayesian approach is proposed to analyze the fMRI data with two nonstationary noise models (the time-varying variance noise model and the fractional noise model). The covariance matrices of the time-varying variance noise and the fractional noise after wavelet transform are diagonal matrices. This property is investigated under the Bayesian framework. The Bayesian estimator not only gives an accurate estimate of the weights in general linear model, but also provides posterior probability of activation in a voxel and, hence, avoids the limitations (i.e., using only hypothesis testing) in the classical methods. The performance of the proposed Bayesian methods (under the assumption of different noise models) are compared with the ordinary least squares (OLS) and the weighted least squares (WLS) methods. Results from the simulation studies validate the superiority of the proposed approach to the OLS and WLS methods considering the complex noise structures in the fMRI data.