Simultaneous Estimation of Population Receptive Field and Hemodynamic Parameters from Single Point BOLD Responses using Metropolis-Hastings Sampling

• Published in: bioRxiv
• Main Authors: Stanislaw Ryszard Adaszewski, Slater, David, Melie-Garcia, Lester, Draganski, Bogdan, Bogorodzki, Piotr
• Format: Paper
• Language: English
• Published: Cold Spring Harbor Cold Spring Harbor Laboratory Press 14.12.2017; Cold Spring Harbor Laboratory
• Edition: 1.1
• Subjects: Bayesian analysis; Heuristic; Mathematical models; Monte Carlo simulation; Neuroscience; Parameter estimation; Receptive field; Sampling; Bayesian estimation; Markov Chain Monte Carlo; Balloon-Windkessel hemodynamic model; Metropolis-Hastings; population receptive field; posterior sampling
• ISSN: 2692-8205; 2692-8205
• DOI: 10.1101/233619

We introduce a new approach to Bayesian Population Receptive Field model estimation using Markov Chain Monte Carlo sampling for simultaneous estimation of receptive field and hemodynamic parameters. To obtain high performance on commonly accessible hardware we present a novel heuristic consisting of interpolation between precomputed responses for predetermined stimuli and a large cross section of receptive field parameters. We investigate the validity of the proposed approach with respect to convergence, tuning and biases. We compare different combinations of receptive field (Compressive Spatial Summation, Dumoulin-Wandell) and hemodynamic (5-parameter and 3-parameter Balloon-Windkessel) models within our framework with and without the usage of the new heuristic. We evaluate estimation consistency and log probability across models. We perform as well a comparison of one model with and without lookup table within the RStan framework using its No-U-Turn Sampler. We present accelerated computation of whole-region-of-interest parameters for one subject. Finally, we discuss risks and limitations associated with the usage of the new heuristic as well as the means of resolving them. We found that the new algorithm is a valid sampling approach to joint receptive field/hemodynamic parameter estimation and that it exhibits very high performance.