Inferring neural firing rates from spike trains using Gaussian processes

Neural spike trains present challenges to analytical efforts due to their noisy, spiking nature. Many studies of neuroscientic and neural prosthetic importance rely on a smoothed, denoised estimate of the spike train’s underlying ring rate. Current techniques to nd time-varying ring rates require ad hoc choices of parameters, offer no condence intervals on their estimates, and can obscure potentially important single trial variability. We present a new method, based on a Gaussian Process prior, for inferring probabilistically optimal estimates of ring rate functions underlying single or multiple neural spike trains. We test the performance of the method on simulated data and experimentally gathered neural spike trains, and we demonstrate improvements over conventional estimators.