Predictive sparse modeling of fMRI data for improved classification, regression, and visualization using the k-support norm

We explore various sparse regularization techniques for analyzing fMRI data, such as the ℓ1 norm (often called LASSO in the context of a squared loss function), elastic net, and the recently introduced k-support norm. Employing sparsity regularization allows us to handle the curse of dimensionality, a problem commonly found in fMRI analysis. In this work we consider sparse regularization in both the regression and classification settings. We perform experiments on fMRI scans from cocaine-addicted as well as healthy control subjects. We show that in many cases, use of the k-support norm leads to better predictive performance, solution stability, and interpretability as compared to other standard approaches. We additionally analyze the advantages of using the absolute loss function versus the standard squared loss which leads to significantly better predictive performance for the regularization methods tested in almost all cases. Our results support the use of the k-support norm for fMRI analysis and on the clinical side, the generalizability of the I-RISA model of cocaine addiction.

[1]  Anna Rose Childress,et al.  Dopamine transporter genotype modulation of neural responses to smoking cues: confirmation in a new cohort , 2011, Addiction biology.

[2]  Andreas Bartels,et al.  Augmenting Feature-driven fMRI Analyses: Semi-supervised learning and resting state activity , 2009, NIPS.

[3]  Dimitris Samaras,et al.  Oral methylphenidate normalizes cingulate activity in cocaine addiction during a salient cognitive task , 2010, Proceedings of the National Academy of Sciences.

[4]  Rita Z. Goldstein,et al.  Drug addiction and its underlying neurobiological basis: neuroimaging evidence for the involvement of the frontal cortex. , 2002, The American journal of psychiatry.

[5]  Matthew B. Blaschko,et al.  FMRI analysis of cocaine addiction using k-support sparsity , 2013, 2013 IEEE 10th International Symposium on Biomedical Imaging.

[6]  Bertrand Thirion,et al.  Multiscale Mining of fMRI Data with Hierarchical Structured Sparsity , 2012, SIAM J. Imaging Sci..

[7]  Rita Z. Goldstein,et al.  Anterior cingulate cortex hypoactivations to an emotionally salient task in cocaine addiction , 2009, Proceedings of the National Academy of Sciences.

[8]  Jin Fan,et al.  Common and distinct networks underlying reward valence and processing stages: A meta-analysis of functional neuroimaging studies , 2011, Neuroscience & Biobehavioral Reviews.

[9]  R. Tibshirani Regression Shrinkage and Selection via the Lasso , 1996 .

[10]  Andreas Bartels,et al.  The chronoarchitecture of the human brain—natural viewing conditions reveal a time-based anatomy of the brain , 2004, NeuroImage.

[11]  Marc Teboulle,et al.  A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems , 2009, SIAM J. Imaging Sci..

[12]  D. Tomasi,et al.  Dopaminergic involvement during mental fatigue in health and cocaine addiction , 2012, Translational Psychiatry.

[13]  A. Laird,et al.  The Neural Basis of Drug Stimulus Processing and Craving: An Activation Likelihood Estimation Meta-Analysis , 2011, Biological Psychiatry.

[14]  Rita Z. Goldstein,et al.  Role of the anterior cingulate and medial orbitofrontal cortex in processing drug cues in cocaine addiction , 2007, Neuroscience.

[15]  Junzhou Huang,et al.  Efficient MR Image Reconstruction for Compressed MR Imaging , 2010, MICCAI.

[16]  Matthew B. Blaschko A Note on k-support Norm Regularized Risk Minimization , 2013, ArXiv.

[17]  H. Zou,et al.  Regularization and variable selection via the elastic net , 2005 .

[18]  Andreas Bartels,et al.  Semi-supervised kernel canonical correlation analysis with application to human fMRI , 2011, Pattern Recognit. Lett..

[19]  Nathan Srebro,et al.  Sparse Prediction with the $k$-Support Norm , 2012, NIPS.

[20]  Trevor Hastie,et al.  Regularization Paths for Generalized Linear Models via Coordinate Descent. , 2010, Journal of statistical software.

[21]  Junzhou Huang,et al.  Composite splitting algorithms for convex optimization , 2011, Comput. Vis. Image Underst..

[22]  Jennifer Bramen,et al.  Effect of bupropion treatment on brain activation induced by cigarette-related cues in smokers. , 2011, Archives of general psychiatry.

[23]  Matthew B. Blaschko,et al.  Predicting cross-task behavioral variables from fMRI data using the k-support norm , 2014 .

[24]  Andreas Bartels,et al.  Brain dynamics during natural viewing conditions—A new guide for mapping connectivity in vivo , 2005, NeuroImage.

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

[26]  Dimitris Samaras,et al.  Can a Single Brain Region Predict a Disorder? , 2012, IEEE Transactions on Medical Imaging.

[27]  Janaina Mourão Miranda,et al.  Unsupervised analysis of fMRI data using kernel canonical correlation , 2007, NeuroImage.

[28]  Li Yao,et al.  Comparative Study of SVM Methods Combined with Voxel Selection for Object Category Classification on fMRI Data , 2011, PloS one.

[29]  N. Logothetis,et al.  Natural vision reveals regional specialization to local motion and to contrast-invariant, global flow in the human brain. , 2008, Cerebral cortex.