On the interpretation of weight vectors of linear models in multivariate neuroimaging

被引:878
作者
Haufe, Stefan [1 ,2 ]
Meinecke, Frank [1 ,3 ]
Goergen, Kai [4 ,5 ,6 ]
Daehne, Sven [1 ]
Haynes, John-Dylan [2 ,4 ,5 ]
Blankertz, Benjamin [2 ,6 ]
Biessgmann, Felix [1 ,7 ]
机构
[1] Tech Univ Berlin, Fachgebiet Maschinelles Lernen, Berlin, Germany
[2] Bernstein Focus Neurotechnol, Berlin, Germany
[3] Zalando GmbH, Berlin, Germany
[4] Charite, Bernstein Ctr Computat Neurosci, D-13353 Berlin, Germany
[5] Charite, Berlin Ctr Adv Neuroimaging, D-13353 Berlin, Germany
[6] Tech Univ Berlin, Fachgebiet Neurotechnol, Berlin, Germany
[7] Korea Univ, Seoul, South Korea
基金
新加坡国家研究基金会;
关键词
Neuroimaging; Multivariate; Univariate; fMRI; EEG; Forward/backward models; Generative/discriminative models; Encoding; Decoding; Activation patterns; Extraction filters; Interpretability; Regularization; Sparsity; HEMODYNAMIC SIGNALS; NEURAL ACTIVITY; EEG; CONNECTIVITY; LOCALIZATION; OSCILLATIONS; SELECTION;
D O I
10.1016/j.neuroimage.2013.10.067
中图分类号
Q189 [神经科学];
学科分类号
071006 ;
摘要
The increase in spatiotemporal resolution of neuroimaging devices is accompanied by a trend towards more powerful multivariate analysis methods. Often it is desired to interpret the outcome of these methods with respect to the cognitive processes under study. Here we discuss which methods allow for such interpretations, and provide guidelines for choosing an appropriate analysis for a given experimental goal: For a surgeon who needs to decide where to remove brain tissue it is most important to determine the origin of cognitive functions and associated neural processes. In contrast, when communicating with paralyzed or comatose patients via brain-computer interfaces, it is most important to accurately extract the neural processes specific to a certain mental state,These-equally-important but complementary objectives require differentanalysis methods. Determining the origin of neural processes in time or space from the parameters of a data-driven model requires what we call a forward model of the data; such a model explains how the measured data was generated from the neural sources. Examples are general linear models (GLMs). Methods for the extraction of neural information from data can be considered as backward models, as they attempt to reverse the data generating process. Examples are multivariate classifiers. Here we demonstrate that the parameters of forward models are neurophysiologically interpretable in the sense that significant nonzero weights are only observed at channels the activity of which is related to the brain process under study. In contrast, the interpretation of backward model parameters can lead to wrong conclusions regarding the spatial or temporal origin of the neural signals of interest, since significant nonzero weights may also be observed at channels the activity of which is statistically independent of the brain process under study. As a remedy for the linear case, we propose a procedure for transforming bacicward models into forward models. This procedure enables the neurophysiological interpretation of the parameters of linear backward models. We hope that this work raises awareness for an often encountered problem and provides a theoretical basis for conducting better interpretable multivariate neuroimaging analyses. (C) 2013 The Authors. Published by Elsevier Inc. All rights reserved.
引用
收藏
页码:96 / 110
页数:15
相关论文
共 65 条
[1]  
[Anonymous], 2012, Event-related potentials
[2]   Electromagnetic brain mapping [J].
Baillet, S ;
Mosher, JC ;
Leahy, RM .
IEEE SIGNAL PROCESSING MAGAZINE, 2001, 18 (06) :14-30
[3]  
Biessmann F., 2012, P 2 NIPS WO IN PRESS
[4]   Improved decoding of neural activity from fMRI signals using non-separable spatiotemporal deconvolutions [J].
Biessmann, Felix ;
Murayama, Yusuke ;
Logothetis, Nikos K. ;
Mueller, Klaus-Robert ;
Meinecke, Frank C. .
NEUROIMAGE, 2012, 61 (04) :1031-1042
[5]  
Biessmann Felix, 2011, IEEE Rev Biomed Eng, V4, P26, DOI 10.1109/RBME.2011.2170675
[6]   Temporal kernel CCA and its application in multimodal neuronal data analysis [J].
Biessmann, Felix ;
Meinecke, Frank C. ;
Gretton, Arthur ;
Rauch, Alexander ;
Rainer, Gregor ;
Logothetis, Nikos K. ;
Mueller, Klaus-Robert .
MACHINE LEARNING, 2010, 79 (1-2) :5-27
[7]  
Blankertz B, 2002, ADV NEUR IN, V14, P157
[8]   Optimizing spatial filters for robust EEG single-trial analysis [J].
Blankertz, Benjamin ;
Tomioka, Ryota ;
Lemm, Steven ;
Kawanabe, Motoaki ;
Mueller, Klaus-Robert .
IEEE SIGNAL PROCESSING MAGAZINE, 2008, 25 (01) :41-56
[9]   Single-trial analysis and classification of ERP components - A tutorial [J].
Blankertz, Benjamin ;
Lemm, Steven ;
Treder, Matthias ;
Haufe, Stefan ;
Mueller, Klaus-Robert .
NEUROIMAGE, 2011, 56 (02) :814-825
[10]   Jacobi angles for simultaneous diagonalization [J].
Cardoso, JF ;
Souloumiac, A .
SIAM JOURNAL ON MATRIX ANALYSIS AND APPLICATIONS, 1996, 17 (01) :161-164