A least angle regression method for fMRI activation detection in phase-encoded experimental designs

This paper presents a new regression method for functional magnetic resonance imaging (fMRI) activation detection. Unlike general linear models (GLM), this method is based on selecting models for activation detection adaptively which overcomes the limitation of requiring a predefined design matrix in GLM. This limitation is because GLM designs assume that the response of the neuron populations will be the same for the same stimuli, which is often not the case. In this work, the fMRI hemodynamic response model is selected from a series of models constructed online by the least angle regression (LARS) method. The slow drift terms in the design matrix for the activation detection are determined adaptively according to the fMRI response in order to achieve the best fit for each fMRI response. The LARS method is then applied along with the Moore-Penrose pseudoinverse (PINV) and fast orthogonal search (FOS) algorithm for implementation of the selected model to include the drift effects in the design matrix. Comparisons with GLM were made using 11 normal subjects to test method superiority. This paper found that GLM with fixed design matrix was inferior compared to the described LARS method for fMRI activation detection in a phased-encoded experimental design. In addition, the proposed method has the advantage of increasing the degrees of freedom in the regression analysis. We conclude that the method described provides a new and novel approach to the detection of fMRI activation which is better than GLM based analyses.

[1]  Jody Tanabe,et al.  See Blockindiscussions, Blockinstats, Blockinand Blockinauthor Blockinprofiles Blockinfor Blockinthis Blockinpublication Comparison Blockinof Blockindetrending Blockinmethods Blockinfor Optimal Blockinfmri Blockinpreprocessing , 2022 .

[2]  Colin L. Mallows,et al.  Some Comments on Cp , 2000, Technometrics.

[3]  G. Glover,et al.  Retinotopic organization in human visual cortex and the spatial precision of functional MRI. , 1997, Cerebral cortex.

[4]  Y. Benjamini,et al.  THE CONTROL OF THE FALSE DISCOVERY RATE IN MULTIPLE TESTING UNDER DEPENDENCY , 2001 .

[5]  M. Korenberg Identifying nonlinear difference equation and functional expansion representations: The fast orthogonal algorithm , 2006, Annals of Biomedical Engineering.

[6]  E. Bullmore,et al.  Statistical methods of estimation and inference for functional MR image analysis , 1996, Magnetic resonance in medicine.

[7]  Michael Elad,et al.  Optimally sparse representation in general (nonorthogonal) dictionaries via ℓ1 minimization , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[8]  X. Huo,et al.  When do stepwise algorithms meet subset selection criteria , 2007, 0708.2149.

[9]  Tim Hesterberg,et al.  Least Angle Regression and LASSO for Large Datasets , 2009, Stat. Anal. Data Min..

[10]  S. Weisberg Applied Linear Regression , 1981 .

[11]  Behzad Mansouri,et al.  Cortical deficits in human amblyopia: their regional distribution and their relationship to the contrast detection deficit. , 2007, Investigative ophthalmology & visual science.

[12]  C. L. Mallows Some comments on C_p , 1973 .

[13]  Thomas E. Nichols,et al.  Controlling the familywise error rate in functional neuroimaging: a comparative review , 2003, Statistical methods in medical research.

[14]  Alan C. Evans,et al.  A general statistical analysis for fMRI data , 2000, NeuroImage.

[15]  Jafar A. Khan,et al.  Robust Linear Model Selection Based on Least Angle Regression , 2007 .

[16]  C. Mallows Some Comments on Cp , 2000, Technometrics.

[17]  J W Belliveau,et al.  Borders of multiple visual areas in humans revealed by functional magnetic resonance imaging. , 1995, Science.

[18]  R. Tibshirani,et al.  Least angle regression , 2004, math/0406456.

[19]  A. Dale,et al.  Tonotopic organization in human auditory cortex revealed by progressions of frequency sensitivity. , 2004, Journal of neurophysiology.

[20]  George A. F. Seber,et al.  Linear regression analysis , 1977 .

[21]  Behzad Mansouri,et al.  The fidelity of the cortical retinotopic map in human amblyopia , 2007, The European journal of neuroscience.

[22]  J. Duyn,et al.  Investigation of Low Frequency Drift in fMRI Signal , 1999, NeuroImage.

[23]  C F Starmer,et al.  Computer computations of cardiac output using the gamma function. , 1970, Journal of applied physiology.

[24]  P. McCullagh,et al.  Generalized Linear Models , 1984 .

[25]  H. Akaike A new look at the statistical model identification , 1974 .

[26]  Clifford M. Hurvich,et al.  Regression and time series model selection in small samples , 1989 .

[27]  Chantal Delon-Martin,et al.  fMRI Retinotopic Mapping—Step by Step , 2002, NeuroImage.

[28]  R. Tibshirani,et al.  Forward stagewise regression and the monotone lasso , 2007, 0705.0269.

[29]  T. Hesterberg,et al.  Least angle and ℓ1 penalized regression: A review , 2008, 0802.0964.

[30]  G. Glover Deconvolution of Impulse Response in Event-Related BOLD fMRI1 , 1999, NeuroImage.

[31]  A. T. Smith,et al.  Estimating receptive field size from fMRI data in human striate and extrastriate visual cortex. , 2001, Cerebral cortex.

[32]  William H. Press,et al.  Numerical recipes , 1990 .

[33]  Ki H. Chon,et al.  Application of fast orthogonal search to linear and nonlinear stochastic systems , 2007, Annals of Biomedical Engineering.

[34]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[35]  D. Heeger,et al.  Linear Systems Analysis of Functional Magnetic Resonance Imaging in Human V1 , 1996, The Journal of Neuroscience.

[36]  C. Starmer,et al.  Indicator Transit Time Considered as a Gamma Variate , 1964, Circulation research.

[37]  Karl J. Friston,et al.  Statistical parametric maps in functional imaging: A general linear approach , 1994 .

[38]  S. Billings,et al.  Orthogonal parameter estimation algorithm for non-linear stochastic systems , 1988 .

[39]  E C Wong,et al.  Processing strategies for time‐course data sets in functional mri of the human brain , 1993, Magnetic resonance in medicine.