Multisubject activation detection in fMRI by testing correlation of data with a signal subspace

In this article, a generalized likelihood ratio test is proposed to assess the correlation between multisubject functional MRI (fMRI) time series and bases of a signal subspace for detecting the existence of group activation in each voxel of the brain. The signal subspace is generated by a design matrix using the time series of the desired effects. The proposed method leads to testing the product of eigenvalues of a specific matrix. The eigenvector corresponding to the largest eigenvalue is the weighting vector for the linear combination of time series of various subjects that has the maximum correlation with the signal subspace. In another method, namely, canonical correlation analysis, the largest eigenvalue of the above matrix is tested for activation detection. Surrogate data on resting state (no activation) are generated by randomization and used to estimate the statistical distribution of these parameters under the null hypothesis condition. A postprocessing step is applied to prevent false detection of voxels that are not sufficiently active (among subjects) by defining a minimum ratio for the active population. The proposed methods are applied on simulated and experimental fMRI data, and the results are compared with those of the general linear model (GLM; using the SPM and FMRISTAT toolboxes). The proposed methods showed higher detection sensitivity as compared with the GLM for activation detection in simulated data. Similarly, they detected more activated regions than did the GLM from multisubject experimental fMRI data on a visual (sensorimotor) event-related task.