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.

Mesoscopic modeling of hidden spiking neurons

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

Mesoscopic population equations for spiking neural networks with synaptic short-term plasticity

Deep Spiking Neural Network Using Spatio-temporal Backpropagation with Variable Resistance.

Efficient Point-Process Modeling of Spiking Neurons for Neuroprosthesis

Reconstruction of recurrent synaptic connectivity of thousands of neurons from simulated spiking activity

Toward Efficient Processing and Learning with Spikes: New Approaches for Multispike Learning

A Simple Probabilistic Spiking Neuron Model with Hebbian Learning Rules.

Generative Modeling of Neural Dynamics via Latent Stochastic Differential Equations

Decoding Spiking Mechanism with Dynamic Learning on Neuron Population

Silicon Modeling of Spiking Neurons with Diverse Dynamic Behaviors

Modeling Neuromorphic Persistent Firing Networks

Common-input models for multiple neural spike-train data

Recurrent Switching Dynamical Systems Models for Multiple Interacting Neural Populations

Memristive leaky integrate-and-fire neuron and learnable straight-through estimator in spiking neural networks

Next generation neural population models

Characterizing the nonlinear structure of shared variability in cortical neuron populations using latent variable models

Efficient Inference on a Network of Spiking Neurons using Deep Learning

Latent Diffusion for Neural Spiking Data

Toward stochastic neural computing