How machine learning is shaping cognitive neuroimaging
GigaScience volume 3, Article number: 28 (2014)
Functional brain images are rich and noisy data that can capture indirect signatures of neural activity underlying cognition in a given experimental setting. Can data mining leverage them to build models of cognition? Only if it is applied to well-posed questions, crafted to reveal cognitive mechanisms. Here we review how predictive models have been used on neuroimaging data to ask new questions, i.e., to uncover new aspects of cognitive organization. We also give a statistical learning perspective on these progresses and on the remaining gaping holes.
Functional neuroimaging has opened the door to quantitative yet non invasive experiments on brain function. These experiments contribute to bridging the gap between cognitive sciences and neuroscience: the former analyse thought and mind while the latter probes the nervous system at various spatial and temporal scales. To study high-level aspects of human cognition, the two modalities of choice are functional Magnetic Resonance Imaging (fMRI) and electro-and magneto-encephalography (EEG/MEG), both can be used to observe brain activity with good spatial resolution for fMRI and temporal temporal resolution for EEG/MEG. The concurrent progress of scanners and experimental paradigms has made it possible to accumulate very rich imaging data that quantify specific correlates of brain function in an uncountable variety of cognitive tasks and processes. In parallel, the advent of machine learning has brought huge progress to data processing of large datasets. But these techniques are geared towards well-posed predictive tasks. The key question is then; how can they be leveraged to push forward understanding of the brain, beyond merely predicting a numerical signal?
This paper presents a subjective view on the work that has been done combining machine learning with functional neuroimaging to advance the understanding of brain function. It dwells mostly on modeling considerations: how and what do the predictive models teach us about the brain? But it also touches upon machine learning and statistical issues. This review focuses on fMRI in humans, that represents most of the accumulated functional neuroimaging data; however, most of the concepts carry to other imaging modalities. FMRI provides images of the brain at the mm scale, however it is only sensitive to the metabolic counterpart of neural activity and suffers from a poor temporal resolution. The first two sections of this paper discuss supervised learning, used first to model brain activity from the stimuli, then to predict the task performed from the evoked activity. The last section reviews the use of unsupervised learning to extract relevant structures in functional images: the interaction structure that underlies brain function, or their natural spatial organization.
Encoding: richer models of evoked activity
The keystone to the use of fMRI in cognitive neuroscience is the standard mass-univariate analysis framework. It consists of modeling the brain response evoked via an experimental paradigm as the linear combination of different experimental conditions [1, 2]. A statistical test is performed at each voxel to delineate regions recruited differently by the various conditions. The art of fMRI experiment design and analysis then consists in crafting the succession of conditions so that, when properly contrasted, they reveal the neural support corresponding to the cognitive function of interest. With regards to brain function, this statistical analysis answers naturally a “where” question, but to a lesser extent a “how” question. Indeed the tests for differences between experimental conditions are statistically well-posed, but not very expressive to refine cognitive models.
In contrast, the study of neural coding, lead historically via intra-cellular recordings of neural activity, has opened the door to breaking down many cognitive functions into atomic steps implemented by ensembles of neurons. The seminal work of Hubel and Wiesel  showed that neurons in the primary visual cortex have receptive fields tuned to a variety of image features, from simple cells sensitive to local orientation in an image, to more complex cells capturing in addition, motion and length of local image features. Progress on uncovering the link between stimuli and neural response revealed neurons tuned to richer and higher-level descriptions of the stimulus, such as receptive fields specific to complex shapes , but also a richer description of neural responses, in particular coding distributed across a population of neurons .
Beyond individual neurons, at the spatial scales probed in fMRIa, and high-level cognition arises from functional integration of multiple specialized brain regions .
The stepping stones of this line of work are to find the right features of the stimuli and neuronal population that can be matched closely. How well the former explains the latter gives a natural figure of merit of these models, in a setting known as encoding. Given models that explain neural responses at the spatial scales captured by fMRI [9, 10] rather than at the neural level, encoding research can be lead with fMRI data, which benefits from full-brain coverage. Technically, designing an encoding model is not different from specifying the design matrix in a standard fMRI analysis and can be seen as model-based fMRI. However relinquishing the methodology of contrasts for more diverse, albeit indirect, statistical tests opens the door to richer modeling. In particular, it is possible to tackle more complex stimuli, such as natural stimuli , very high-level and diverse descriptions of the stimuli , or a cognitive model of the observed behavior .
This increase in model complexity is the driving force behind the use of machine learning in encoding. First it entails fitting many parameters on limited data, and thus conventional in-sample statistical testing is thorny. For this reason, goodness of fit of the encoding model is best assessed via its cross-validated ability to predict brain signals . Similarly, the predictive engine that links stimuli features to brain signal is best chosen amongst machine learning tools, that balance modeling flexibility and regularization, such as a naive Bayes predictor , sparse  or ridge  regression. Finally, the computational models that derive encoding features from the stimuli often draw from the domain-specific feature extraction techniques developed in applied machine learning research. These provide simple quantitative proxies for the cognitive features of interest. For instance, to map semantic concepts  and  used natural language processing techniques: word co-occurrence or an ontology on words. The ties between brain science and machine learning are strikingly close in the study of vision: computer vision, i.e., the use of computers and machine learning to analyze and interpret images, has built upon, but also fostered our understanding of the brain visual system. David Marr’s seminal work  formalized the idea of hierarchical levels of representation that tie together the receptive fields observed in visual cortex, but is also reflected in modern state-of-the-art computer vision architecture based on convolutional networks . Very recently, Yamins et al. have shown a striking correspondence between 296 neural recordings in the infero-temporal cortex of the monkey and intermediate layers of computer-vision convolutional networks. This work is a quintessential example of machine learning in encoding models: a predictive engine performs the same task as the brain system under study; machine learning is used to fit its parameters on a set of stimuli and the final architecture matches neural data.
Transferring such results to fMRI would open doors to studying the full complete brain of healthy human subjects rather than 296 neurons in implanted monkeys. However, it poses significant challenges. Indeed, fMRI is an indirect and noisy measurement of brain activity, that captures the average effect of many spikes and does not resolve cortical columns, let alone individual neurons. The concept of population receptive field is sometimes used to refer to the aggregate properties of neurons in one voxel. Thus, encoding models need to be adapted to the resulting structured noise and signal convolutions. Model evaluation and selection is in itself often a major roadblock.
Decoding: towards principled reverse inference
In the study of neural recordings, decoding models reconstruct stimuli or behavior from the neural code [5, 19]. More generally, the decoding approach can be seen as solving the inverse problem to the encoding model, even when applied on fMRI signals that do not capture individual neural firing [20, 21].
Since a decoding model often predicts quantities that are directly observable, it can provide a very useful tool to validate an encoding model. Indeed, decoding performance is an omnibus test of goodness of fit: it tests the overall significance of the model, but does not test which variables have a significant contribution to the model. As an omnibus test, decoding is used with explicit sophisticated encodings [8, 11–14], but also with simple fMRI analysis to perform an omnibus test at the region level  or on a wide family of regions as in searchlight analysis . Interestingly, an early fMRI study  on neural representation hypothesized that “objects are represented by a relatively widely distributed activity of functional modules”, but considered this statement to be insufficiently quantitative to allow a statistical test. Nowadays this study would probably be formulated in an encoding/decoding framework , using a multivariate predictor to provide evidence for the author’s hypothesis, as in . It is often considered that multi-voxel analysis, as used in decoding, provides an increase in sensitivity compared to standard mass-univariate analysis ; however, we stress that it does not correspond to an increase in statistical power, but rather to a different test performed: decoding performs a global (omnibus) test of the model, while voxel-level tests are useful to delineate regions, but are subject to corrections for multiple comparisons.
As noted in , decoding analysis provides a good framework to interpret overlapping activation patterns. Brain maps in encoding and decoding settings carry actually a different meaning. An inference with an encoding model, or in the fMRI standard analysis framework, is a statement on whether or not the signal in a brain region is well explained by the model that we have of the task: we can conclude that the task implies this brain activation, and we say that the region is recruited by the task. A decoding analysis tells us that if we observe a certain brain activity, we can deduce properties of the task or the stimulus. Such a conclusion is the converse implication of the encoding settings, sometimes dubbed reverse inference. Reverse inference, i.e., drawing conclusions on behavior and mental processes from the brain activations, answers natural questions in cognitive neuroimaging, e.g.,: what is the function of neural sub-system? But reverse inferences drawn from maps, estimated using encoding models, are a logical fallacy . On the other hand, decoding models provide a path to principled reverse inferences . However, it is important to keep in mind that, in general, a decoding experiment does not tell us anything about tasks and cognitive processes that it did not probe. For example, an experiment studying brain regions discriminating images of faces from images of houses  does not inform us on how these regions are related to recognizing letters.
The appealing idea of inferring brain processes from brain activation only carries meaning if the decoding model has captured a large variety of brain processes. Beyond interpretation of brain images, the basic neuroscience questions at stakes here are that of functional specificity. For instance, while many brain regions are more activated under physical pain, a decoding analysis including many different aspects of pain showed that a network comprising parts of the thalamus, the insulae, and the somatosensory cortex was specific of physical pain . At the spatial scale probed by fMRI, the multiplicity of regions needed to come to precise conclusions on the cognitive function recruited is consistent with the modern view that high-level cognitive processes arise from distributed networks. This calls for multivariate decoding engines.
Going beyond a specific cognitive domain, such as vision or pain, and studying functional specialization in a broad sense require probing more functions than can be addressed in one experimental paradigm. For this reason, investigators have turned to accumulating data across experiments. Using 8 studies, covering 22 different cognitive concepts, Poldrack et al. were able to predict the concepts involved from activation images in unseen subjects. The use of a variety of studies, with different experimental protocols, can overcome the idiosyncrasies of each protocol that are not relevant to cognitive concepts of interest; for instance, to study high-level decision mechanisms independently of the modality used to present stimuli –visual or auditory. However, in , the train set contained images from the same protocols as the test set; thus, the hypothesis that the decoder was actually detecting protocols rather than cognitive concepts cannot be ruled out. To generalize to unseen protocols, the challenge is to describe them in terms that are common enough to be shared across many protocols, but also sufficiently rich to capture their cognitive content. Schwartz et al. used an ontology of experimental paradigms and multi-label classification: labeling 83 different conditions, from 19 studies, with a set of different terms from the ontology. The resulting predicting engine can not only describe the content of an unseen experiment from the corresponding brain activation, but also give brain maps associated with each term in a reverse inference. Covering more cognitive concepts requires accumulating many brain images. Sharing data across institutions is a practical means to this end, for instance relying on the OpenfMRI project  that hosts to this day 24 different fMRI studies. Another interesting alley is to collect from the literature the coordinates, in standard brain space, of observed activation foci, as in the Neurosynth project .
Although decoding gives a principled methodological framework for reverse inference, there are some tough statistical challenges. Indeed, the discriminant brain maps extracted may be the most relevant information captured by the model from a neuroscience perspective. However, decoders solve a high-dimensional multivariate statistical estimation problem that is very ill-posed  given the typical small sample size. Many different brain maps will give rise to similar predictive performance. Worst yet, minimizing a prediction risk does not lead to any control on the brain maps. For instance, if two neighboring voxels carry the same information but one is less noisy than the other, a decoder might favor selecting only that one. For related reasons, sparse models can only capture a subset of relevant voxels . Injecting priors –or regularization– in the estimation makes it well-posed and shapes the brain maps extracted. Capturing large-scale brain systems calls for spatial regularization such as sparsity and spatial smoothness  or total-variation (TV) for piecewise smooth maps . In particular TV- ℓ 1 regularization, combining sparsity and total-variation, selects well the predictive regions . Unlike widespread belief, multivariate tools used commonly, such as support vector machines or searchlight, seem to do a worse job at selecting predictive regions than univariate tools .
Encoding and decoding models explore the two directions linking brain activation to stimuli and cognitive processes  (see Figure 1). Both of these methodologies do not form credible models of how the brain creates cognition. They are rather experimental devices to test hypotheses and retrieve brain maps, where the critical modeling work goes in the formal description of the cognitive concepts associated with the brain signals under study. This description is most often a non-trivial transformation of the stimuli, non-linear  or calling for concept ontologies [14, 29, 30]. Following the concepts of neural coding and Marr’s vision that good representations give rise to powerful computational processing , encoding and decoding models are often understood as revealing a representational space, distributed representations in the cortex that reflect fundamental concepts [9, 38]. However, the combination of the lack of temporal dynamics in fMRI and the linear models that we rely upon naturally create such an understanding of the data in terms of representations, while for some functions studied, the actual neural implementation may be closer to processes dynamically sustained information, as in theories of conscious processing . In this light, the use of linear models for decoding may be criticized as too simple to capture non-linear interactions. However, from the neuroscience point-of-view they lead to probing well-posed questions  and from the statistical learning point of view, they can be relatively well-behaved even in very high dimensional settings with the typical small sample sizes faced by fMRI .
Finding hidden structure: parcellations and connectomes
In machine learning applications, it is often easier to accumulate unlabeled data than labeled data. This is also the case in neuroimaging, as controlling the cognitive state of the subject is very challenging and calls for careful experimental design and analysis. Data collection speed is limited by the timescale of psychological experiments. On the opposite, accumulating so-called resting-state acquisitions, in which the mental state of the subject is not controlled, is much more tractable , and is applicable to diminished populations . The interest of resting-state data for cognitive neuroimaging is not immediate, as it does not carry explicit cognitive information; however, they reflect on-going activity, which is an important part of brain function. Indeed, on-going activity shares a common substrate with activity explicitly evoked by controlled cognition, priming cognition but also shaped by task . Unsupervised learning on resting-state scans holds the promise of extracting intrinsic brain structures .
Capturing brain interactions
The brain is a heavily interacting system. Mapping its connections in the form of a connectome can help to understand the flow of information in the brain. Fluctuations in brain activity, for example, observed at rest, reveal functional interactions and thus can be used to estimate a functional connectome. From a statistical learning perspective, the estimation of a connectome can be formalized as extracting the conditional independence structure from observed correlations . The challenge here is that of the paucity of data, and can be tackled with graph estimators that have good small-sample structure recovery properties, such as sparse covariance models [47, 48].
Learning functional units
Interest in resting-state data arose originally from the observation that voxel-based signals observed at rest could be used to segment spatial structures known from task studies . Subsequently, researchers realized that these could exhibit some additional functional structures . What spatial distributed brain networks are modulated during rest? This question can be formulated as that of blind source separation, and independent component analysis (ICA) provides a good algorithm to recover these networks [51, 52]. Datasets of increasing size and quality lead to extracting more networks, that break up in a set of smaller regions, paving the brain in a parcellation . Breaking down the brain into homogeneous units is a long quest in neuroscience that can be traced back to Brodmann areas. Such parcellations have been historically driven by anatomical features. Resting-state fMRI provides valuable data to learn a functional parcellation, as it gives a spatially-resolved window into intrinsic brain function. Indeed, functionally-homogeneous regions can be extracted by clustering voxels with similar fMRI time-series [54, 55]. The unmixing model underlying ICA can be adapted to extracting regions by formulating it in the more general framework of dictionary learning  and adding sparsity-inducing penalty that also favor clustered spatial components, thus yielding region segmentations . While identifying intrinsic functional brain modules is crucial from a basic neuroscience point of view, brain parcellation can also provide useful data reduction even if they don’t capture true functional units [21, 34]. These different purposes give rise to different methodological trade-offs . Beyond resting-state data, applying similar methods to databases of evoked activity exploring a large variety of cognitive concepts can have the additional benefit of appending cognitive labels to the spatial units extracted [59, 60].
However, care must be exercised when applying the brain-parcellation techniques. By construction, such methods will return a parcellation, even if there is little to no structure in the data. They do not build upon well-posed statistical hypothesis testing. The methods can often be unstable, with a small modification of the input data leading to large changes in the results. This unstability can be explained by, on one hand the lack of explicit noise model, and on the other hand the fact that unsupervised learning is an intrinsically hard problem from the statistical standpoint. Validation of the functional units is very challenging beyond a simple confirmation bias that boils down to checking for known structures, the variability of which is unknown and uncontrolled. Some researchers have explored quantifying variability of the patterns [55, 57, 58, 61] or controlling how well they explain the data [57, 58] but these are weak proxys of the neuroscientific questions on brain organization.
Practical considerations: methods and implementations matter
The focus of this review is not on methodological details, but on general concepts and approaches that further our understanding of brain function. However, it is important to stress that many of the roadblocks to the use of machine-learning-related techniques in cognitive neuroimaging lie in the methods. From a theoretical point of view, the statistical control is seldom warranted by the models used [34, 37]. On the empirical side of things, best practices are not established. The high-dimensionality of the statistical models and the plurality of methods considered mean that, at the level of the literature, machine-learning techniques probably give rise to more variability, although they do come with more expressiveness and power.
A final critical aspect, all too often overlooked, is that of software. The standard GUI-based fMRI data processing environments, such as SPM, FSL  or AFNI , do not implement most of the modeling approaches described in this review. FSL and AFNI do provide some methods tailored to fMRI uses (respectively ICA  and basic decoding ). There is progress on dedicated tools such as PyMVPA , but these require the practitioners to learn new skills, in particular some understanding of machine learning and basic programming. The challenges of a good environment for machine-learning on fMRI data is that it should be simple enough to be within reach of the practitioner, yet leverage a powerful machine-learning toolkit, such as the popular scikit-learn package in Python , and offer flexibility to assemble new models, encoding, decoding, or unsupervised .
The goals of cognitive neurosciences are to link cognition with its neural basis. FMRI gives a noisy and incomplete window on neural mechanisms. Nevertheless, to map effects at a large scale, it is priceless, as it can be applied massively on healthy human subjects, and thus enables the systematic study of high-level cognition. Machine learning tools are instrumental in making the most of this data, as they do not require a precise mechanistic understanding of the signal, but rather to frame a prediction problem that captures some relevant aspects of brain function, as in encoding or decoding. However, for progress in neuroscience, black-box prediction engines do not suffice as the key to understanding brain function lies in the properties of the signal used for prediction. For these reasons, the statistics aspects in statistical learning cannot be neglected: different methods give rise to different results and the figure of merit does not simply boil down to predictive power.
a It is unlikely that standard fMRI acquisitions, even after analysis with powerful multivariate methods, capture information at the level of the cortical column .
Functional magnetic resonnance imaging
Independent component analysis
Graphical User Interface.
Penny WD, Friston KJ, Ashburner JT, Kiebel SJ, Nichols TE: Statistical Parametric Mapping: The Analysis of Functional Brain mages. 2007, London: Academic Press
Poldrack RA, Mumford JA, Nichols TE: Handbook of Functional MRI Data Analysis. 2011, Cambridge: University Press
Hubel DH, Wiesel TN: Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex. J Physiol. 1962, 160: 106.
Logothetis NK, Pauls J, Poggio T: Shape representation in the inferior temporal cortex of monkeys. Curr Biol. 1995, 5: 552-10.1016/S0960-9822(95)00108-4.
Georgopoulos AP, Schwartz AB, Kettner RE: Neuronal population coding of movement direction. Science. 1986, 233: 1416-10.1126/science.3749885.
Freeman J, Brouwer GJ, Heeger DJ, Merriam EP: Orientation decoding depends on maps, not columns. J Neurosci. 2011, 31: 4792-10.1523/JNEUROSCI.5160-10.2011.
Tononi G, Edelman GM, Sporns O: Complexity and coherency: integrating information in the brain. Trends Cogn Sci. 1998, 2: 474-10.1016/S1364-6613(98)01259-5.
Naselaris T, Kay KN, Nishimoto S, Gallant JL: Encoding and decoding in fMRI. Neuroimage. 2011, 56: 400-10.1016/j.neuroimage.2010.07.073.
Edelman S, Grill-Spector K, Kushnir T, Malach R: Toward direct visualization of the internal shape representation space by fMRI. Psychobiology. 1998, 26: 309.
O’Doherty JP, Hampton A, Kim H: Model-based fMRI and its application to reward learning and decision making. Ann N Y Acad Sci. 2007, 1104: 35-10.1196/annals.1390.022.
Kay KN, Naselaris T, Prenger RJ, Gallant JL: Identifying natural images from human brain activity. Nature. 2008, 452: 352-10.1038/nature06713.
Mitchell TM, Shinkareva SV, Carlson A, Chang K-M, Malave VL, Mason RA, Just MA: Predicting human brain activity associated with the meanings of nouns. Science. 2008, 320: 1191-10.1126/science.1152876.
Nishimoto S, Vu AT, Naselaris T, Benjamini Y, Yu B, Gallant JL: Reconstructing visual experiences from brain activity evoked by natural movies. Curr Biol. 2011, 21: 1641-10.1016/j.cub.2011.08.031.
Huth AG, Nishimoto S, Vu AT, Gallant JL: A continuous semantic space describes the representation of thousands of object and action categories across the human brain. Neuron. 2012, 76: 1210-10.1016/j.neuron.2012.10.014.
Marr D: Vision: A Computational Investigation Into the Human Representation and Processing of Visual Information. 1982, Cambridge: The MIT Press
LeCun Y, Kavukcuoglu K, Farabet C: Convolutional networks and applications in vision. Circuits and Systems (ISCAS). 2010, IEEE, 253-253.
Yamins DL, Hong H, Cadieu CF, Solomon EA, Seibert D, DiCarlo JJ: Performance-optimized hierarchical models predict neural responses in higher visual cortex. Proc Natl Acad Sci. 2014, 111: 201403112.
Dumoulin SO, Wandell BA: Population receptive field estimates in human visual cortex. Neuroimage. 2008, 39 (2): 647-660. 10.1016/j.neuroimage.2007.09.034.
Abbott L: Decoding neuronal firing and modelling neural networks. Q Rev Biophys. 1994, 27: 291-10.1017/S0033583500003024.
Dehaene S, Le Clec’H G, Cohen L, Poline J-B, van de Moortele P-F, Le Bihan D: Inferring behavior from functional brain images. Nat Neurosci. 1998, 1: 549-10.1038/2785.
Thirion B, Flandin G, Pinel P, Roche A, Ciuciu P, Poline JB: Dealing with the shortcomings of spatial normalization: Multi-subject parcellation of fMRI datasets. Hum Brain Map. 2006, 27: 678-10.1002/hbm.20210.
Haxby JV, Gobbini IM, Furey ML, Ishai A, Schouten JL, Pietrini P: Distributed and overlapping representations of faces and objects in ventral temporal cortex. Science. 2001, 293 (5539): 2425-10.1126/science.1063736.
Kriegeskorte N, Goebel R, Bandettini P: Information-based functional brain mapping. Proc Natl Acad Sci U S A. 2006, 103: 3863-10.1073/pnas.0600244103.
Norman KA, Polyn SM, Detre GJ, Haxby JV: Beyond mind-reading: multi-voxel pattern analysis of fMRI data. Trends Cogn Sci. 2006, 10: 424-10.1016/j.tics.2006.07.005.
Peelen MV, Downing PE: Using multi-voxel pattern analysis of fMRI data to interpret overlapping functional activations. Trends Cogn Sci. 2007, 11: 4-10.1016/j.tics.2006.10.009.
Poldrack R: Can cognitive processes be inferred from neuroimaging data?. Trends Cogn Sci. 2006, 10: 59-10.1016/j.tics.2005.12.004.
Poldrack RA: Inferring mental states from neuroimaging data: from reverse inference to large-scale decoding. Neuron. 2011, 72: 692-10.1016/j.neuron.2011.11.001.
Wager TD, Atlas LY, Lindquist MA, Roy M, Woo C-W, Kross E: An fMRI-based neurologic signature of physical pain. N Engl J Med. 2013, 368: 1388-10.1056/NEJMoa1204471.
Poldrack RA, Halchenko YO, Hanson SJ: Decoding the large-scale structure of brain function by classifying mental states across individuals. Psychol Sci. 2009, 20: 1364-10.1111/j.1467-9280.2009.02460.x.
Schwartz Y, Thirion B, Varoquaux G: Mapping cognitive ontologies to and from the brain. NIPS. 2013, Curran Associates, Inc.
Poldrack RA, Barch DM, Mitchell JP, Wager TD, Wagner AD, Devlin JT, Cumba C, Koyejo O, Milham MP: Toward open sharing of task-based fMRI data: the OpenfMRI project. Frontiers Neuroinformatics. 2013,,7,
Yarkoni T, Poldrack RA, Nichols TE, Van Essen DC, Wager TD: Large-scale automated synthesis of human functional neuroimaging data. Nat Methods. 2011, 8: 665-10.1038/nmeth.1635.
Haufe S, Meinecke F, Görgen K, Däne S, Haynes J-D, Blankertz B, Bießmann F: On the interpretation of weight vectors of linear models in multivariate neuroimaging. Neuroimage. 2014, 87: 96-110.
Varoquaux G, Gramfort A, Thirion B: Small-sample brain mapping: sparse recovery on spatially correlated designs with randomization and clustering. ICML. 2012, Omnipress, 1375-1375.
Grosenick L, Klingenberg B, Katovich K, Knutson B, Taylor JE: Interpretable whole-brain prediction analysis with graphnet. NeuroImage. 2013, 72: 304.
Michel V, Gramfort A, Varoquaux G, Eger E, Thirion B: Total variation regularization for fMRI-based prediction of behavior. Med Imaging IEEE Trans. 2011, 30: 1328.
Gramfort A, Thirion B, Varoquaux G: Identifying predictive regions from fMRI with TV-L1 prior. Pattern Recognition in Neuroimaging (PRNI). 2013, IEEE, 17-17.
Haxby JV, Connolly AC, Guntupalli JS: Decoding neural representational spaces using multivariate pattern analysis. Annu Rev Neurosci. 2014, 37: 435-456. 10.1146/annurev-neuro-062012-170325.
Davis T, Poldrack RA: Measuring neural representations with fMRI: practices and pitfalls. Ann N Y Acad Sci. 2013, 1296: 108-10.1111/nyas.12156.
Dehaene S, Naccache L: Towards a cognitive neuroscience of consciousness: basic evidence and a workspace framework. Cognition. 2001, 79: 1-10.1016/S0010-0277(00)00123-2.
Biswal BB, Mennes M, Zuo XN, Gohel S, Kelly C, Smith SM, Beckmann CF, Adelstein JS, Buckner RL, Colcombe S, Dogonowski AM, Ernst M, Fair D, Hampson M, Hoptman MJ, Hyde JS, Kiviniemi VJ, Kötter R, Li SJ, Lin CP, Lowe MJ, Mackay C, Madden DJ, Madsen KH, Margulies DS, Mayberg HS, McMahon K, Monk CS, Mostofsky SH, Nagel BJ: Toward discovery science of human brain function. Proc Ntl Acad Sci. 2010, 107: 4734-10.1073/pnas.0911855107.
Greicius M: Resting-state functional connectivity in neuropsychiatric disorders. Curr Opin Neurol. 2008, 21: 424.
Sadaghiani S, Hesselmann G, Friston KJ, Kleinschmidt A: The relation of ongoing brain activity, evoked neural responses, and cognition. Frontiers Syst Neuroscience. 2010,,4,
Sporns O, Tononi G, Kotter R: The human connectome: a structural description of the human brain. PLoS Comput Biol. 2005, 1: 42-10.1371/journal.pcbi.0010042.
Varoquaux G, Craddock RC: Learning and comparing functional connectomes across subjects. NeuroImage. 2013, 80: 405.
Varoquaux G, Gramfort A, Poline JB, Thirion B: Markov models for fMRI correlation structure: is brain functional connectivity small world, or decomposable into networks?. J Physiol - Paris. 2012, 106: 212-10.1016/j.jphysparis.2012.01.001.
Varoquaux G, Gramfort A, Poline JB, Thirion B: Brain covariance selection: better individual functional connectivity models using population prior. NIPS. 2010, Curran Associates, Inc.
Smith SM, Miller KL, Salimi-Khorshidi G, Webster M, Beckmann CF, Nichols TE, Ramsey JD, Woolrich MW: Network modelling methods for fMRI. Neuroimage. 2011, 54: 875-10.1016/j.neuroimage.2010.08.063.
Biswal B, Zerrin Yetkin F, Haughton VM, Hyde JS: Functional connectivity in the motor cortex of resting human brain using echo-planar MRI. Magn Reson Med. 1995, 34: 53719.
Raichle ME, MacLeod AM, Snyder AZ, Powers WJ, Gusnard DA, Shulman GL: A default mode of brain function. Proc Natl Acad Sci U S A. 2001, 98: 676-10.1073/pnas.98.2.676.
Kiviniemi V, Kantola JH, Jauhiainen J, Hyvärinen A, Tervonen O: Independent component analysis of nondeterministic fmri signal sources. Neuroimage. 2003, 19: 253-10.1016/S1053-8119(03)00097-1.
Beckmann CF, DeLuca M, Devlin JT, Smith SM: Investigations into resting-state connectivity using independent component analysis. Philos Trans R Soc Lond B. 2005, 360: 1001-10.1098/rstb.2005.1634.
Kiviniemi V, Starck T, Remes J, Long X, Nikkinen J, Haapea M, Veijola J, Moilanen I, Isohanni M, Zang YF, Tervonen O: Functional segmentation of the brain cortex using high model order group PICA. Hum Brain Map. 2009, 30: 3865-10.1002/hbm.20813.
Yeo BTT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, Roffman JL, Smoller JW, Zöllei L, Polimeni JR, Fischl B, Liu H, Buckner RL: The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysio. 2011, 106: 1125-10.1152/jn.00338.2011.
Craddock RC, James GA, Holtzheimer PE, Hu XP, Mayberg HS: A whole brain fMRI atlas generated via spatially constrained spectral clustering. Human Brain Mapping. 2012, 33: 1914-10.1002/hbm.21333.
Varoquaux G, Gramfort A, Pedregosa F, Michel V, Thirion B: Multi-subject dictionary learning to segment an atlas of brain spontaneous activity. Inf Proc Med Imag. 2011, IEEE, 562-562.
Abraham A, Dohmatob E, Thirion B, Samaras D, Varoquaux G: Extracting brain regions from rest fMRI with total-variation constrained dictionary learning. MICCAI. 2013, Frontiers Media S.A., 607-607.
Thirion B, Varoquaux G, Dohmatob E, Poline J: Which fMRI clustering gives good brain parcellations?. Name: Frontiers Neurosci. 2014, 8: 167.
Smith SM, Fox PT, Miller KL, Glahn DC, Fox PM, Mackay CE, Filippini N, Watkins KE, Toro R, Laird AR, Beckmann CF: Correspondence of the brain’s functional architecture during activation and rest. Proc Natl Acad Sci. 2009, 106: 13040-10.1073/pnas.0905267106.
Varoquaux G, Schwartz Y, Pinel P, Thirion B: Cohort-level brain mapping: learning cognitive atoms to single out specialized regions. Information Processing in Medical Imaging. 2013, Springer, 438-438.
Varoquaux G, Sadaghiani S, Pinel P, Kleinschmidt A, Poline JB, Thirion B: A group model for stable multi-subject ICA on fMRI datasets. NeuroImage. 2010, 51: 288-10.1016/j.neuroimage.2010.02.010.
Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM: Advances in functional and structural mr image analysis and implementation as FSL. Neuroimage. 2004, 23: 208-219.
Cox RW: AFNI: software for analysis and visualization of functional magnetic resonance neuroimages. Comput Biomed Res. 1996, 29: 162-10.1006/cbmr.1996.0014.
LaConte SM: Decoding fMRI brain states in real-time. Neuroimage. 2011, 56: 440-10.1016/j.neuroimage.2010.06.052.
Hanke M, Halchenko YO, Sederberg PB, Hanson SJ, Haxby JV, Pollmann S: PyMVPA: A python toolbox for multivariate pattern analysis of fmri data. Neuroinformatics. 2009, 7: 37-10.1007/s12021-008-9041-y.
Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, Blondel M, Prettenhofer P, Weiss R, Dubourg V, Vanderplas J, Passos A, Cournapeau D, Brucher M, Perrot M, Duchesnay E: Scikit-learn: machine learning in Python. J Mach Learn Res. 2011, 12: 2825.
Abraham A, Pedregosa F, Eickenberg M, Gervais P, Mueller A, Kossaifi J, Gramfort A, Thirion B, Varoquaux G: Machine learning for neuroimaging with scikit-learn. Frontiers Neuroinformatics. 2014,,8,
The authors acknowledge fruitful discussions with Catherine Wacongne and thank the two reviewers, Krzysztof Gorgolewski and Tal Yarkoni, for their review which improved the manuscript.
The authors declare that they have no competing interests.
GV and BT carried out the research and drafted the manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.