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.

Modelling spine locations on dendrite trees using inhomogeneous Cox point processes

Reaction–Diffusion Model-Based Research on Formation Mechanism of Neuron Dendritic Spine Patterns

Hotspots of dendritic spine turnover facilitate clustered spine addition and learning and memory

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

The cylindrical K-function and Poisson line cluster point processes

Anomalous transport and nonlinear reactions in spiny dendrites

A Neurocomputational Method for Fully Automated 3D Dendritic Spine Detection and Segmentation of Medium-Sized Spiny Neurons

Competitive processes shape multi-synapse plasticity along dendritic segments

Models and simulation of 3D neuronal dendritic trees using Bayesian networks

Comprehensive Analysis of Human Dendritic Spine Morphology and Density: Exploring Diversity of Human Dendritic Spines

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

A New Algorithm For 3d Dendritic Spine Detection

Bayesian inference for Matérn repulsive processes

Optimizing a Neuron for Reliable Dendritic Subunit Pooling

Tracking the Morphological Evolution of Neuronal Dendrites by First-Passage Analysis

A Multivariate Point Process Model for Simultaneously Recorded Neural Spike Trains

Modelling Point Referenced Spatial Count Data: A Poisson Process Approach

Electro-diffusive modeling and the role of spine geometry on action potential propagation in neurons

The spine decomposition of branching processes and their applications

Multi-contact synapses for stable networks: a spike-timing dependent model of dendritic spine plasticity and turnover