Methods for computing the maximum performance of computational models of fMRI responses

• Published in: PLoS computational biology Vol. 15; no. 3; p. e1006397
• Main Authors: Lage-Castellanos, Agustin, Valente, Giancarlo, Formisano, Elia, De Martino, Federico
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 01.03.2019; Public Library of Science (PLoS)
• Subjects: Accuracy; Biology and Life Sciences; Brain; Brain - physiology; Brain mapping; Cognitive ability; Computation; Computational neuroscience; Computer and Information Sciences; Computer Simulation; Data acquisition; Diagnostic imaging; Exact solutions; Functional magnetic resonance imaging; Humans; Magnetic resonance imaging; Magnetic Resonance Imaging - methods; Mathematical analysis; Medical imaging; Medicine and Health Sciences; Methods; Model accuracy; Monte Carlo Method; Monte Carlo methods; Monte Carlo simulation; Neuroimaging; Neurology; Neurosciences; NMR; Noise; Noise measurement; Nuclear magnetic resonance; Physical Sciences; Physiology; Regularization; Reproducibility of Results; Research and Analysis Methods; Supervision; Time series; United States > US; Netherlands
• ISSN: 1553-7358; 1553-734X; 1553-7358
• DOI: 10.1371/journal.pcbi.1006397

Computational neuroimaging methods aim to predict brain responses (measured e.g. with functional magnetic resonance imaging [fMRI]) on the basis of stimulus features obtained through computational models. The accuracy of such prediction is used as an indicator of how well the model describes the computations underlying the brain function that is being considered. However, the prediction accuracy is bounded by the proportion of the variance of the brain response which is related to the measurement noise and not to the stimuli (or cognitive functions). This bound to the performance of a computational model has been referred to as the noise ceiling. In previous fMRI applications two methods have been proposed to estimate the noise ceiling based on either a split-half procedure or Monte Carlo simulations. These methods make different assumptions over the nature of the effects underlying the data, and, importantly, their relation has not been clarified yet. Here, we derive an analytical form for the noise ceiling that does not require computationally expensive simulations or a splitting procedure that reduce the amount of data. The validity of this analytical definition is proved in simulations, we show that the analytical solution results in the same estimate of the noise ceiling as the Monte Carlo method. Considering different simulated noise structure, we evaluate different estimators of the variance of the responses and their impact on the estimation of the noise ceiling. We furthermore evaluate the interplay between regularization (often used to estimate model fits to the data when the number of computational features in the model is large) and model complexity on the performance with respect to the noise ceiling. Our results indicate that when considering the variance of the responses across runs, computing the noise ceiling analytically results in similar estimates as the split half estimator and approaches the true noise ceiling under a variety of simulated noise scenarios. Finally, the methods are tested on real fMRI data acquired at 7 Tesla.