Spike trains as (in)homogeneous Poisson processes or Hawkes processes: non-parametric adaptive estimation and goodness-of-fit tests

This article aims to propose new non-parametric adaptive estimation methods and to adapt other recent similar results to the setting of spike trains analysis. After briefly recalling main features of the homogeneous Poisson model, we focus on two main generalizations of this process : the inhomogeneous Poisson model, which is non-stationary, and the Hawkes model, which can take into account interactions. Goodness-of-fit tests are also proposed and are proved to be of prescribed asymptotical level. They enable us to test these non-parametric models. Various simulations show good performance of the estimation and test procedures. A complete analysis is also performed with these tools on single unit activity recorded on a monkey during a sensory-motor task. We can show that the homogeneous Poisson process hypothesis is always rejected and that the inhomogeneous Poisson process hypothesis is rarely accepted. The Hawkes model seems to fit most of the data.

[1]  B. Laurent,et al.  Adaptive tests of homogeneity for a Poisson process , 2009, 0905.0989.

[2]  P. Reynaud-Bouret,et al.  Near optimal thresholding estimation of a Poisson intensity on the real line , 2008, 0810.5204.

[3]  Shigeru Shinomoto,et al.  Kernel bandwidth optimization in spike rate estimation , 2009, Journal of Computational Neuroscience.

[4]  Ole Winther,et al.  Multivariate Hawkes process models of the occurrence of regulatory elements , 2010, BMC Bioinformatics.

[5]  Stefan Rotter,et al.  Recurrent interactions in spiking networks with arbitrary topology. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  Christine Tuleau-Malot,et al.  A model of Poissonian interactions and detection of dependence , 2015, Stat. Comput..

[7]  D. Vere-Jones,et al.  Some examples of statistical estimation applied to earthquake data , 1982 .

[8]  P. Brémaud.  Point Processes and Queues , 1981 .

[9]  Alexandre B. Tsybakov,et al.  Introduction to Nonparametric Estimation , 2008, Springer series in statistics.

[10]  J. Klein,et al.  Statistical Models Based On Counting Process , 1994 .

[11]  Franck Grammont,et al.  Multiple Tests Based on a Gaussian Approximation of the Unitary Events Method with Delayed Coincidence Count , 2012, Neural Computation.

[12]  Emery N. Brown,et al.  The Time-Rescaling Theorem and Its Application to Neural Spike Train Data Analysis , 2002, Neural Computation.

[13]  M. Rudemo.  Empirical Choice of Histograms and Kernel Density Estimators , 1982 .

[14]  W. Wildman,et al.  Theoretical Neuroscience , 2014 .

[15]  P. Brémaud,et al.  STABILITY OF NONLINEAR HAWKES PROCESSES , 1996 .

[16]  F. Papangelou.  Integrability of expected increments of point processes and a related random change of scale , 1972 .

[17]  Sonja Grün,et al.  Significance of joint-spike events based on trial-shuffling by efficient combinatorial methods , 2003, Complex..

[18]  A. P. Georgopoulos,et al.  Neuronal population coding of movement direction. , 1986, Science.

[19]  M. Duflo,et al.  Probability and statistics , 2012, Actuarial Mathematics for Life Contingent Risks.

[20]  Alexa Riehle,et al.  Spike synchronization and firing rate in a population of motor cortical neurons in relation to movement direction and reaction time , 2003, Biological Cybernetics.

[21]  Stefan Rotter,et al.  How Structure Determines Correlations in Neuronal Networks , 2011, PLoS Comput. Biol..

[22]  Sonja Grün,et al.  Non-parametric significance estimation of joint-spike events by shuffling and resampling , 2003, Neurocomputing.

[23]  R. Tibshirani.  Regression Shrinkage and Selection via the Lasso , 1996 .

[24]  Sonja Grün,et al.  Unitary joint events in multiple neuron spiking activity: detection, significance, and interpretation , 1996 .

[25]  Daryl J. Daley,et al.  An Introduction to the Theory of Point Processes , 2013 .

[26]  Sonja Grün,et al.  Dynamical changes and temporal precision of synchronized spiking activity in monkey motor cortex during movement preparation , 2000, Journal of Physiology-Paris.

[27]  M K Habib,et al.  Dynamics of neuronal firing correlation: modulation of "effective connectivity". , 1989, Journal of neurophysiology.

[28]  Yosihiko Ogata,et al.  Statistical Models for Earthquake Occurrences and Residual Analysis for Point Processes , 1988 .

[29]  G. Rizzolatti,et al.  The mirror-neuron system. , 2004, Annual review of neuroscience.

[30]  N. Bingham.  EMPIRICAL PROCESSES WITH APPLICATIONS TO STATISTICS (Wiley Series in Probability and Mathematical Statistics) , 1987 .

[31]  Christophe Pouzat,et al.  On Goodness of Fit Tests For Models of Neuronal Spike Trains Considered as Counting Processes , 2009, 0909.2785.

[32]  Shy Shoham,et al.  Frontiers in Computational Neuroscience Correlation-based Analysis and Generation of Multiple Spike Trains Using Hawkes Models with an Exogenous Input , 2022 .

[33]  P. Reynaud-Bouret,et al.  Adaptive estimation for Hawkes processes; application to genome analysis , 2009, 0903.2919.

[34]  Marie Doumic,et al.  Nonparametric Estimation of the Division Rate of a Size-Structured Population , 2011, SIAM J. Numer. Anal..

[35]  Khashayar Pakdaman,et al.  Dynamics of a structured neuron population , 2009 .

[36]  Lasso and probabilistic inequalities for multivariate point processes , 2015, 1208.0570.

[37]  T. Rydén.  An EM algorithm for estimation in Markov-modulated Poisson processes , 1996 .

[38]  Alexa Riehle,et al.  Cancellation of a planned movement in monkey motor cortex , 2006, Neuroreport.

[39]  Magalie Fromont,et al.  Adaptive goodness-of-fit tests in a density model , 2006, math/0607013.

[40]  Antoine Chaffiol,et al.  Automatic spike train analysis and report generation. An implementation with R, R2HTML and STAR , 2009, Journal of Neuroscience Methods.

[41]  A. Goldenshluger,et al.  Bandwidth selection in kernel density estimation: Oracle inequalities and adaptive minimax optimality , 2010, 1009.1016.

[42]  Patricia Reynaud-Bouret,et al.  Adaptive estimation of the intensity of inhomogeneous Poisson processes via concentration inequalities , 2003 .

[43]  Robert D. Nowak,et al.  Multiscale Poisson Intensity and Density Estimation , 2007, IEEE Transactions on Information Theory.

[44]  Sophie Schbath,et al.  FADO: A Statistical Method to Detect Favored or Avoided Distances between Occurrences of Motifs using the Hawkes' Model , 2005, Statistical applications in genetics and molecular biology.

[45]  C. D. Kemp,et al.  Density Estimation for Statistics and Data Analysis , 1987 .

[46]  Larry Wasserman,et al.  All of Statistics: A Concise Course in Statistical Inference , 2004 .

[47]  Sonja Grün,et al.  Analysis of Parallel Spike Trains , 2010 .

[48]  Harold B. Sackrowitz,et al.  Evaluating tests for increasing intensity of a Poisson process , 1993 .

[49]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[50]  P. Reynaud-Bouret,et al.  Adaptive density estimation: a curse of support? ✩ , 2009, 0907.1794.

[51]  S. Holm.  A Simple Sequentially Rejective Multiple Test Procedure , 1979 .

[52]  A. Hawkes.  Point Spectra of Some Mutually Exciting Point Processes , 1971 .