Construction of Point Process Adaptive Filter Algorithms for Neural Systems Using Sequential Monte Carlo Methods

The stochastic state point process filter (SSPPF) and steepest descent point process filter (SDPPF) are adaptive filter algorithms for state estimation from point process observations that have been used to track neural receptive field plasticity and to decode the representations of biological signals in ensemble neural spiking activity. The SSPPF and SDPPF are constructed using, respectively, Gaussian and steepest descent approximations to the standard Bayes and Chapman-Kolmogorov (BCK) system of filter equations. To extend these approaches for constructing point process adaptive filters, we develop sequential Monte Carlo (SMC) approximations to the BCK equations in which the SSPPF and SDPPF serve as the proposal densities. We term the two new SMC point process filters SMC-PPFS and SMC-PPFD , respectively. We illustrate the new filter algorithms by decoding the wind stimulus magnitude from simulated neural spiking activity in the cricket cercal system. The SMC-PPFS and SMC-PPFD provide more accurate state estimates at low number of particles than a conventional bootstrap SMC filter algorithm in which the state transition probability density is the proposal density. We also use the SMC-PPFS algorithm to track the temporal evolution of a spatial receptive field of a rat hippocampal neuron recorded while the animal foraged in an open environment. Our results suggest an approach for constructing point process adaptive filters using SMC methods

[1]  M. Pettet,et al.  Dynamic changes in receptive-field size in cat primary visual cortex. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

[2]  R E Kass,et al.  Recursive bayesian decoding of motor cortical signals by particle filtering. , 2004, Journal of neurophysiology.

[3]  J. Kaas,et al.  The reorganization of somatosensory cortex following peripheral nerve damage in adult and developing mammals. , 1983, Annual review of neuroscience.

[4]  Frederic Paik Schoenberg,et al.  Multidimensional Residual Analysis of Point Process Models for Earthquake Occurrences , 2003 .

[5]  Matthew A. Wilson,et al.  Dynamic Analyses of Information Encoding in Neural Ensembles , 2004, Neural Computation.

[6]  Yong Rui,et al.  Better proposal distributions: object tracking using unscented particle filter , 2001, Proceedings of the 2001 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. CVPR 2001.

[7]  B L McNaughton,et al.  Dynamics of the hippocampal ensemble code for space. , 1993, Science.

[8]  Michael Isard,et al.  Contour Tracking by Stochastic Propagation of Conditional Density , 1996, ECCV.

[9]  M. Wilson,et al.  Analyzing Functional Connectivity Using a Network Likelihood Model of Ensemble Neural Spiking Activity , 2005, Neural Computation.

[10]  L. Frank,et al.  Behavioral/Systems/Cognitive Hippocampal Plasticity across Multiple Days of Exposure to Novel Environments , 2022 .

[11]  G. Jacobs,et al.  Direction sensitivity of the filiform hair population of the cricket cereal system , 1995, Journal of Comparative Physiology A.

[12]  Hong Wang,et al.  Voice source localization for automatic camera pointing system in videoconferencing , 1997, Proceedings of 1997 Workshop on Applications of Signal Processing to Audio and Acoustics.

[13]  Nando de Freitas,et al.  Sequential Monte Carlo Methods in Practice , 2001, Statistics for Engineering and Information Science.

[14]  J. P. Miller,et al.  Stimulus-response properties of cricket cereal filiform receptors , 1995, Journal of Comparative Physiology A.

[15]  Emery N. Brown,et al.  Dynamic Analysis of Neural Encoding by Point Process Adaptive Filtering , 2004, Neural Computation.

[16]  E. S. Chornoboy,et al.  Maximum likelihood identification of neural point process systems , 1988, Biological Cybernetics.

[17]  Emery N. Brown,et al.  Estimating a State-space Model from Point Process Observations Emery N. Brown , 2022 .

[18]  Petar M. Djurić.  Sequential Estimation of Signals under Model Uncertainty , 2001, Sequential Monte Carlo Methods in Practice.

[19]  J. Mendel.  Lessons in Estimation Theory for Signal Processing, Communications, and Control , 1995 .

[20]  Emery N. Brown,et al.  The Time-Rescaling Theorem and Its Application to Neural Spike Train Data Analysis , 2002, Neural Computation.

[21]  Dawn M. Taylor,et al.  Direct Cortical Control of 3D Neuroprosthetic Devices , 2002, Science.

[22]  Nicholas G. Hatsopoulos,et al.  Brain-machine interface: Instant neural control of a movement signal , 2002, Nature.

[23]  N. Weinberger.  Learning-induced changes of auditory receptive fields , 1993, Current Opinion in Neurobiology.

[24]  Hong Wang,et al.  Voice source localization for automatic camera pointing system in videoconferencing , 1997, 1997 IEEE International Conference on Acoustics, Speech, and Signal Processing.

[25]  Miguel A. L. Nicolelis,et al.  Real-time control of a robot arm using simultaneously recorded neurons in the motor cortex , 1999, Nature Neuroscience.

[26]  Jun S. Liu,et al.  Sequential Monte Carlo methods for dynamic systems , 1997 .

[27]  David C. Hoaglin,et al.  Applications, basics, and computing of exploratory data analysis , 1983 .

[28]  Wei Wu,et al.  Bayesian Population Decoding of Motor Cortical Activity Using a Kalman Filter , 2006, Neural Computation.

[29]  Emery N. Brown,et al.  Computational Neuroscience: A Comprehensive Approach , 2022 .

[30]  V. Solo,et al.  Contrasting Patterns of Receptive Field Plasticity in the Hippocampus and the Entorhinal Cortex: An Adaptive Filtering Approach , 2002, The Journal of Neuroscience.

[31]  E N Brown,et al.  A Statistical Paradigm for Neural Spike Train Decoding Applied to Position Prediction from Ensemble Firing Patterns of Rat Hippocampal Place Cells , 1998, The Journal of Neuroscience.

[32]  E. Brown,et al.  A point-process model of human heartbeat intervals: new definitions of heart rate and heart rate variability. , 2005, American journal of physiology. Heart and circulatory physiology.

[33]  William Bialek,et al.  Spikes: Exploring the Neural Code , 1996 .

[34]  M. Quirk,et al.  Experience-Dependent Asymmetric Shape of Hippocampal Receptive Fields , 2000, Neuron.

[35]  A. Doucet,et al.  Monte Carlo Smoothing for Nonlinear Time Series , 2004, Journal of the American Statistical Association.

[36]  E N Brown,et al.  An analysis of neural receptive field plasticity by point process adaptive filtering , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[37]  M. Cynader,et al.  Somatosensory cortical map changes following digit amputation in adult monkeys , 1984, The Journal of comparative neurology.

[38]  Simon Haykin,et al.  Adaptive filter theory (2nd ed.) , 1991 .

[39]  E. Bizzi,et al.  Neuronal Correlates of Motor Performance and Motor Learning in the Primary Motor Cortex of Monkeys Adapting to an External Force Field , 2001, Neuron.

[40]  G. Kitagawa.  Smoothness priors analysis of time series , 1996 .

[41]  R. Andersen,et al.  Cognitive Control Signals for Neural Prosthetics , 2004, Science.

[42]  C. I. Connolly,et al.  Building neural representations of habits. , 1999, Science.

[43]  Jerald D. Kralik,et al.  Real-time prediction of hand trajectory by ensembles of cortical neurons in primates , 2000, Nature.

[44]  J. Donoghue.  Plasticity of adult sensorimotor representations , 1995, Current Opinion in Neurobiology.

[45]  Teresa H. Y. Meng,et al.  Model-based neural decoding of reaching movements: a maximum likelihood approach , 2004, IEEE Transactions on Biomedical Engineering.

[46]  J. Fritz,et al.  Rapid task-related plasticity of spectrotemporal receptive fields in primary auditory cortex , 2003, Nature Neuroscience.