Multivariate pattern analysis for MEG: A comparison of dissimilarity measures

&NA; Multivariate pattern analysis (MVPA) methods such as decoding and representational similarity analysis (RSA) are growing rapidly in popularity for the analysis of magnetoencephalography (MEG) data. However, little is known about the relative performance and characteristics of the specific dissimilarity measures used to describe differences between evoked activation patterns. Here we used a multisession MEG data set to qualitatively characterize a range of dissimilarity measures and to quantitatively compare them with respect to decoding accuracy (for decoding) and between‐session reliability of representational dissimilarity matrices (for RSA). We tested dissimilarity measures from a range of classifiers (Linear Discriminant Analysis – LDA, Support Vector Machine – SVM, Weighted Robust Distance – WeiRD, Gaussian Naïve Bayes – GNB) and distances (Euclidean distance, Pearson correlation). In addition, we evaluated three key processing choices: 1) preprocessing (noise normalisation, removal of the pattern mean), 2) weighting decoding accuracies by decision values, and 3) computing distances in three different partitioning schemes (non‐cross‐validated, cross‐validated, within‐class‐corrected). Four main conclusions emerged from our results. First, appropriate multivariate noise normalization substantially improved decoding accuracies and the reliability of dissimilarity measures. Second, LDA, SVM and WeiRD yielded high peak decoding accuracies and nearly identical time courses. Third, while using decoding accuracies for RSA was markedly less reliable than continuous distances, this disadvantage was ameliorated by decision‐value‐weighting of decoding accuracies. Fourth, the cross‐validated Euclidean distance provided unbiased distance estimates and highly replicable representational dissimilarity matrices. Overall, we strongly advise the use of multivariate noise normalisation as a general preprocessing step, recommend LDA, SVM and WeiRD as classifiers for decoding and highlight the cross‐validated Euclidean distance as a reliable and unbiased default choice for RSA. HighlightsWe provide Python and MATLAB tutorials for key analysis steps.We compared dissimilarity measures and preprocessing choices for MEG MVPA.Multivariate noise normalisation is a key preprocessing step.LDA, SVM and WeiRD are recommended classifiers for decoding.The cross‐validated Euclidean distance is a reliable and unbiased choice for RSA.

[1]  Nikolaus Kriegeskorte,et al.  Frontiers in Systems Neuroscience Systems Neuroscience , 2022 .

[2]  Chih-Jen Lin,et al.  LIBSVM: A library for support vector machines , 2011, TIST.

[3]  Dimitrios Pantazis,et al.  Multivariate pattern analysis of MEG and EEG: a comparison of representational structure in time and space , 2016 .

[4]  Radoslaw Martin Cichy,et al.  Probing principles of large‐scale object representation: Category preference and location encoding , 2013, Human brain mapping.

[5]  Guillaume A. Rousselet,et al.  Does Filtering Preclude Us from Studying ERP Time-Courses? , 2012, Front. Psychology.

[6]  Michael Eickenberg,et al.  Machine learning for neuroimaging with scikit-learn , 2014, Front. Neuroinform..

[7]  Frank Tong,et al.  Representational Dynamics of Facial Viewpoint Encoding , 2017, Journal of Cognitive Neuroscience.

[8]  K. Seeliger,et al.  CNN-based Encoding and Decoding of Visual Object Recognition in Space and Time , 2017, bioRxiv.

[9]  Philipp Sterzer,et al.  WeiRD - a fast and performant multivoxel pattern classifier , 2016, 2016 International Workshop on Pattern Recognition in Neuroimaging (PRNI).

[10]  Dimitrios Pantazis,et al.  Multivariate pattern analysis of MEG and EEG: A comparison of representational structure in time and space , 2016, NeuroImage.

[11]  Antonio Torralba,et al.  Comparison of deep neural networks to spatio-temporal cortical dynamics of human visual object recognition reveals hierarchical correspondence , 2016, Scientific Reports.

[12]  Susan G. Wardle,et al.  Decoding Dynamic Brain Patterns from Evoked Responses: A Tutorial on Multivariate Pattern Analysis Applied to Time Series Neuroimaging Data , 2016, Journal of Cognitive Neuroscience.

[13]  Radoslaw Martin Cichy,et al.  The representational dynamics of task and object processing in humans , 2018, eLife.

[14]  Nikolaus Kriegeskorte,et al.  Neural dynamics of real-world object vision that guide behaviour , 2017, bioRxiv.

[15]  Ken Nakayama,et al.  The consequences of subtracting the mean pattern in fMRI multivariate correlation analyses , 2013, Front. Neurosci..

[16]  G. Rees,et al.  Predicting the orientation of invisible stimuli from activity in human primary visual cortex , 2005, Nature Neuroscience.

[17]  Nikolaus Kriegeskorte,et al.  Relating Population-Code Representations between Man, Monkey, and Computational Models , 2009, Front. Neurosci..

[18]  Hans P. Op de Beeck,et al.  Against hyperacuity in brain reading: Spatial smoothing does not hurt multivariate fMRI analyses? , 2010, NeuroImage.

[19]  Nikolaus Kriegeskorte,et al.  Deep Supervised, but Not Unsupervised, Models May Explain IT Cortical Representation , 2014, PLoS Comput. Biol..

[20]  R. Goebel,et al.  Human Object-Similarity Judgments Reflect and Transcend the Primate-IT Object Representation , 2013, Front. Psychol..

[21]  Francesca Pizzorni Ferrarese,et al.  Low-frequency oscillations employ a general coding of the spatio-temporal similarity of dynamic faces , 2017, NeuroImage.

[22]  N. Kriegeskorte,et al.  Author ' s personal copy Representational geometry : integrating cognition , computation , and the brain , 2013 .

[23]  Emiliano Ricciardi,et al.  Beyond sensory images: Object-based representation in the human ventral pathway. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Olivier Ledoit,et al.  A well-conditioned estimator for large-dimensional covariance matrices , 2004 .

[25]  N. Kanwisher,et al.  Feedback of pVisual Object Information to Foveal Retinotopic Cortex , 2008, Nature Neuroscience.

[26]  Radoslaw Martin Cichy,et al.  Resolving human object recognition in space and time , 2014, Nature Neuroscience.

[27]  Jason Weston,et al.  Gene Selection for Cancer Classification using Support Vector Machines , 2002, Machine Learning.

[28]  Jörn Diedrichsen,et al.  Reliability of dissimilarity measures for multi-voxel pattern analysis , 2016, NeuroImage.

[29]  Y.S. Hung,et al.  Gene selection for Brain Cancer Classification , 2006, 2006 International Conference of the IEEE Engineering in Medicine and Biology Society.

[30]  David D. Cox,et al.  Functional magnetic resonance imaging (fMRI) “brain reading”: detecting and classifying distributed patterns of fMRI activity in human visual cortex , 2003, NeuroImage.

[31]  L. Tyler,et al.  Decoding the Cortical Dynamics of Sound-Meaning Mapping , 2017, The Journal of Neuroscience.

[32]  Li Su,et al.  Spatiotemporal Searchlight Representational Similarity Analysis in EMEG Source Space , 2012, 2012 Second International Workshop on Pattern Recognition in NeuroImaging.

[33]  A. Ishai,et al.  Distributed and Overlapping Representations of Faces and Objects in Ventral Temporal Cortex , 2001, Science.

[34]  Dimitrios Pantazis,et al.  Dynamics of scene representations in the human brain revealed by magnetoencephalography and deep neural networks , 2015 .

[35]  Weihong Qian Low-Frequency Oscillations , 2017 .

[36]  K. Grill-Spector,et al.  fMRI-adaptation and category selectivity in human ventral temporal cortex: regional differences across time scales. , 2010, Journal of neurophysiology.

[37]  Alessandro Presacco,et al.  Robust decoding of selective auditory attention from MEG in a competing-speaker environment via state-space modeling , 2016, NeuroImage.

[38]  Keiji Tanaka,et al.  Object category structure in response patterns of neuronal population in monkey inferior temporal cortex. , 2007, Journal of neurophysiology.

[39]  S. Taulu,et al.  Spatiotemporal signal space separation method for rejecting nearby interference in MEG measurements , 2006, Physics in medicine and biology.

[40]  K. Grill-Spector,et al.  Differential development of high-level visual cortex correlates with category-specific recognition memory , 2007, Nature Neuroscience.

[41]  F. Tong,et al.  Decoding the visual and subjective contents of the human brain , 2005, Nature Neuroscience.

[42]  N. Kanwisher,et al.  Only some spatial patterns of fMRI response are read out in task performance , 2007, Nature Neuroscience.

[43]  S. Taulu,et al.  Suppression of Interference and Artifacts by the Signal Space Separation Method , 2003, Brain Topography.

[44]  Karl J. Friston,et al.  Comparing the similarity and spatial structure of neural representations: A pattern-component model , 2011, NeuroImage.

[45]  Jörn Diedrichsen,et al.  Representational models: A common framework for understanding encoding, pattern-component, and representational-similarity analysis , 2017, bioRxiv.

[46]  Dimitrios Pantazis,et al.  Similarity-based fusion of MEG and fMRI reveals spatio-temporal dynamics in human cortex during visual object recognition , 2015 .

[47]  M. Forshaw What Are the Consequences , 2006 .

[48]  Keiji Tanaka,et al.  Matching Categorical Object Representations in Inferior Temporal Cortex of Man and Monkey , 2008, Neuron.

[49]  Venkatapavani Pallavi Punugu Machine Learning in Neuroimaging , 2017 .

[50]  Quanzheng Li,et al.  Decoding the orientation of contrast edges from MEG evoked and induced responses , 2017 .

[51]  Nikolaus Kriegeskorte,et al.  Perceptual similarity of visual patterns predicts dynamic neural activation patterns measured with MEG , 2016, NeuroImage.