Bayesian Active Learning of Neural Firing Rate Maps with Transformed Gaussian Process Priors

Abstract

A firing rate map, also known as a tuning curve, describes the nonlinear relationship between a neuron's spike rate and a low-dimensional stimulus (e.g., orientation, head direction, contrast, color). Here we investigate Bayesian active learning methods for estimating firing rate maps in closed-loop neurophysiology experiments. These methods can accelerate the characterization of such maps through the intelligent, adaptive selection of stimuli. Specifically, we explore the manner in which the prior and utility function used in Bayesian active learning affect stimulus selection and performance. Our approach relies on a flexible model that involves a nonlinearly transformed gaussian process (GP) prior over maps and conditionally Poisson spiking. We show that infomax learning, which selects stimuli to maximize the information gain about the firing rate map, exhibits strong dependence on the seemingly innocuous choice of nonlinear transformation function. We derive an alternate utility function that selects stimuli to minimize the average posterior variance of the firing rate map and analyze the surprising relationship between prior parameterization, stimulus selection, and active learning performance in GP-Poisson models. We apply these methods to color tuning measurements of neurons in macaque primary visual cortex.

Combined Social and Spatial Coding in a Descending Projection from the Prefrontal Cortex

“…To estimate spatial fields and their modulation by social investigation (Figs 4D–G, 5B–C and 6C–D), we used a Bayesian approach based on a Gaussian process (GP) regression model (Park et al 2014) for calcium fluorescence imaging data. This model provides a principled statistical framework for quantifying the relationship between spatial position and neural activity.…”

Dethroning the Fano Factor: a flexible, model-based approach to partitioning neural variability

“…tive optimal experimental design in statistics, seek to select stimuli based on the stimuli and responses observed so far during an experiment in order to characterize the neuron as quickly and efficiently as possible (Lindley, 1956;Bernardo, 1979;MacKay, 1992;Chaloner & Verdinelli, 1995;Cohn et al, 1996;Paninski, 2005). Adaptive stimulus selection is particularly useful in settings where data are limited or expensive to collect, and can substantially reduce in the number of trials needed for fitting an accurate model of neural responses (Paninski et al, 2007;Benda et al, 2007;Lewi et al, 2009;DiMattina & Zhang, 2011, 2013Bölinger & Gollisch, 2012;Park & Pillow, 2012;Park et al, 2014;Pillow & Park, 2016).…”

Bayesian Hypothesis Testing And Experimental Design For Two-Photon Imaging Data

“…We then used the model to predict neural activity for unseen parameter combinations and quantified how uncertain our predictions about the activity in response to these were [20][21][22]. Intuitively, the model should have the least uncertainty about stimulus parameters which had been observed.…”

“…For example, isolated mouse retinal tissue becomes unresponsive to light stimulation in a matter of hours, and single recording fields often bleach within half an hour of recording. To efficiently explore the space of possible stimulus features under severe time constraints is thus a critical problem, which GP models can be used to address [20,21].…”