Abstract:Objective. Temporal patterns in neuronal spiking encode stimulus uncertainty, and convey information about high-level functions such as memory and cognition. Estimating the associated information content and understanding how that evolves with time assume significance in the investigation of neuronal coding mechanisms and abnormal signaling. However, existing estimators of the entropy rate, a measure of information content, either ignore the inherent nonstationarity, or employ dictionary-based Lempel-Ziv (LZ) methods that converge too slowly for one to study temporal variations in sufficient detail. Against this backdrop, we seek estimates that handle nonstationarity, are fast converging, and hence allow meaningful temporal investigations.Approach. We proposed a homogeneous Markov model approximation of spike trains within windows of suitably chosen length and an entropy rate estimator based on empirical probabilities that converges quickly.Main results. We constructed mathematical families of nonstationary Markov processes with certain bi/multi-level properties (inspired by neuronal responses) with known entropy rates, and validated the proposed estimator against those. Further statistical validations were presented on data collected from hippocampal (and primary visual cortex) neuron populations in terms of single neuron behavior as well as population heterogeneity. Our estimator appears to be statistically more accurate and converges faster than existing LZ estimators, and hence well suited for temporal studies.Significance. The entropy rate analysis revealed not only informational and process memory heterogeneity among neurons, but distinct statistical patterns in neuronal populations (from two different brain regions) under basal and post-stimulus conditions. Taking inspiration, we envision future large-scale studies of different brain regions enabled by the proposed tool (estimator), potentially contributing to improved functional modeling of the brain and identification of statistical signatures of neurodegenerative diseases.

A multiple filter test for the detection of rate changes in renewal processes with varying variance

Online Neural Sequence Detection with Hierarchical Dirichlet Point Process

A K-Medoids Based Point-Process Modeling On Neural Spike Transformation Using Binless Kernel

A new family of statistical tests for neuronal spiking and autocorrelated timeseries data

Goodness-of-Fit Tests and Nonparametric Adaptive Estimation for Spike Train Analysis

Separation rates for the detection of synchronization of interacting point processes in a mean field frame. Application to neuroscience

Spike detection and sorting based on nonlinear energy operator and matched filter

Change Point Detection with Neural Online Density-Ratio Estimator

Particle-filtering approaches for nonlinear Bayesian decoding of neuronal spike trains

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

Change-point analysis of neuron spike train data

The "Twin Peaks" method of automated Spike-Wave detection: A two-step, two-criteria Matlab application

Evaluating stationarity via change-point alternatives with applications to fMRI data

Temporally Causal Discovery Tests for Discrete Time Series and Neural Spike Trains

Analysis and modelling of variability and covariability of population spike trains across multiple time scales

The Windowed Moments Change Test: A Novel Technique for Assessing Stationarity in Network Traffic

On Matched Filtering for Statistical Change Point Detection

A point process approach for the classification of noisy calcium imaging data

A doubly stochastic renewal framework for partitioning spiking variability

Robust entropy rate estimation for nonstationary neuronal calcium spike trains based on empirical probabilities