Abstract:Neurons communicate nonlinearly through spike activities. Generalized linear models (GLMs) describe spike activities with a cascade of a linear combination across inputs, a static nonlinear function, and an inhomogeneous Bernoulli or Poisson process, or Cox process if a self-history term is considered. This structure considers the output nonlinearity in spike generation but excludes the nonlinear interaction among input neurons. Recent studies extend GLMs by modeling the interaction among input neurons with a quadratic function, which considers the interaction between every pair of input spikes. However, quadratic effects may not fully capture the nonlinear nature of input interaction. We therefore propose a staged point-process model to describe the nonlinear interaction among inputs using a few hidden units, which follows the idea of artificial neural networks. The output firing probability conditioned on inputs is formed as a cascade of two linear-nonlinear (a linear combination plus a static nonlinear function) stages and an inhomogeneous Bernoulli process. Parameters of this model are estimated by maximizing the log likelihood on output spike trains. Unlike the iterative reweighted least squares algorithm used in GLMs, where the performance is guaranteed by the concave condition, we propose a modified Levenberg-Marquardt(L-M) algorithm, which directly calculates the Hessian matrix of the log likelihood, for the nonlinear optimization in our model. The proposed model is tested on both synthetic data and real spike train data recorded from the dorsal premotor cortex and primary motor cortex of a monkey performing a center-out task. Performances are evaluated by discrete-time rescaled Kolmogorov-Smirnov tests, where our model statistically outperforms a GLM and its quadratic extension, with a higher goodness-of-fit in the prediction results. In addition, the staged point-process model describes nonlinear interaction among input neurons with fewer parameters than quadratic models, and the modified L-M algorithm also demonstrates fast convergence.

Neural Dynamics Discovery via Gaussian Process Recurrent Neural Networks

Modeling Latent Neural Dynamics with Gaussian Process Switching Linear Dynamical Systems

Linear Time GPs for Inferring Latent Trajectories from Neural Spike Trains

Identification of Recurrent Dynamics in Distributed Neural Populations

Extended Poisson Gaussian-Process Latent Variable Model for Unsupervised Neural Decoding

Nonlinear Modeling of Neural Interaction for Spike Prediction Using the Staged Point-Process Model

Latent Variable Double Gaussian Process Model for Decoding Complex Neural Data

LQResNet: A Deep Neural Network Architecture for Learning Dynamic Processes

Inferring stochastic low-rank recurrent neural networks from neural data

Tensor Programs I: Wide Feedforward or Recurrent Neural Networks of Any Architecture are Gaussian Processes

Efficient Learning Algorithms for Gaussian Processes

Continuous time Gaussian process dynamical models in gene regulatory network inference

Neural Latent Space Model for Dynamic Networks and Temporal Knowledge Graphs

Gaussian process learning of nonlinear dynamics

Probabilistic Decomposed Linear Dynamical Systems for Robust Discovery of Latent Neural Dynamics

Latent Variable Sequence Identification for Cognitive Models with Neural Network Estimators

Stochastic Deep Gaussian Processes over Graphs

Developing a Deep Learning Estimator to Learn Nonlinear Dynamic Systems

Learning Neural Point Processes with Latent Graphs

Generative Modeling of Neural Dynamics via Latent Stochastic Differential Equations