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.

Learning Network-Structured Dependence from Non-Stationary Multivariate Point Process Data

Statistical Inference for Networks of High-Dimensional Point Processes

Efficient Estimation of Directed Connectivity in Nonlinear and Nonstationary Spiking Neuron Networks.

Learning Networks from Wide-Sense Stationary Stochastic Processes

An Empirical Study: Extensive Deep Temporal Point Process

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

Continuous-time Graph Representation with Sequential Survival Process

Tracking Dynamic Point Processes on Networks

Neural Point Process for Learning Spatiotemporal Event Dynamics

Learning Time Series Associated Event Sequences with Recurrent Point Process Networks

Deep graph kernel point processes

Learning continuous network emerging dynamics from scarce observations via data-adaptive stochastic processes

Decoupled Marked Temporal Point Process using Neural Ordinary Differential Equations

Online Estimation and Community Detection of Network Point Processes for Event Streams

Marked Neural Spatio-Temporal Point Process Involving a Dynamic Graph Neural Network

Temporal Point Process Graphical Models

Assessment of synchrony in multiple neural spike trains using loglinear point process models

Modeling The Intensity Function Of Point Process Via Recurrent Neural Networks

Learning general temporal point processes based on dynamic weight generation

Learning Registered Point Processes from Idiosyncratic Observations

A Multivariate Point Process Model for Simultaneously Recorded Neural Spike Trains