Gaussian-Process Factor Analysis for Low-Dimensional Single-Trial Analysis of Neural Population Activity

• Published in: Journal of neurophysiology 2009-07, Vol.102 (1), p.614-635
• Main Authors: Yu, Byron M, Cunningham, John P, Santhanam, Gopal, Ryu, Stephen I, Shenoy, Krishna V, Sahani, Maneesh
• Format: Article
• Language: English
• Subjects: Action Potentials - physiology; Animals; Innovative Methodology; Models, Neurological; Nerve Net - physiology; Neural Networks, Computer; Neurons - physiology; Nonlinear Dynamics; Normal Distribution; Principal Component Analysis; Reaction Time - physiology; Signal Processing, Computer-Assisted; Time Factors
• ISSN: 0022-3077; 1522-1598
• DOI: 10.1152/jn.90941.2008

1 Department of Electrical Engineering, 2 Neurosciences Program, and 3 Department of Neurosurgery, Stanford University, Stanford, California; and 4 Gatsby Computational Neuroscience Unit, University College London, London, United Kingdom Submitted 19 August 2008; accepted in final form 24 March 2009 We consider the problem of extracting smooth, low-dimensional neural trajectories that summarize the activity recorded simultaneously from many neurons on individual experimental trials. Beyond the benefit of visualizing the high-dimensional, noisy spiking activity in a compact form, such trajectories can offer insight into the dynamics of the neural circuitry underlying the recorded activity. Current methods for extracting neural trajectories involve a two-stage process: the spike trains are first smoothed over time, then a static dimensionality-reduction technique is applied. We first describe extensions of the two-stage methods that allow the degree of smoothing to be chosen in a principled way and that account for spiking variability, which may vary both across neurons and across time. We then present a novel method for extracting neural trajectories—Gaussian-process factor analysis (GPFA)—which unifies the smoothing and dimensionality-reduction operations in a common probabilistic framework. We applied these methods to the activity of 61 neurons recorded simultaneously in macaque premotor and motor cortices during reach planning and execution. By adopting a goodness-of-fit metric that measures how well the activity of each neuron can be predicted by all other recorded neurons, we found that the proposed extensions improved the predictive ability of the two-stage methods. The predictive ability was further improved by going to GPFA. From the extracted trajectories, we directly observed a convergence in neural state during motor planning, an effect that was shown indirectly by previous studies. We then show how such methods can be a powerful tool for relating the spiking activity across a neural population to the subject's behavior on a single-trial basis. Finally, to assess how well the proposed methods characterize neural population activity when the underlying time course is known, we performed simulations that revealed that GPFA performed tens of percent better than the best two-stage method. Address for reprint requests and other correspondence: M. Sahani, Gatsby Computational Neuroscience Unit, UCL, 17 Queen Square, London, WC1N 3AR, UK (E-mail: manees