An information-geometric framework for statistical inferences in the neural spike train space

Statistical inferences are essentially important in analyzing neural spike trains in computational neuroscience. Current approaches have followed a general inference paradigm where a parametric probability model is often used to characterize the temporal evolution of the underlying stochastic processes. To directly capture the overall variability and distribution in the space of the spike trains, we focus on a data-driven approach where statistics are defined and computed in the function space in which spike trains are viewed as individual points. To this end, we at first develop a parametrized family of metrics that takes into account different warpings in the time domain and generalizes several currently used spike train distances. These new metrics are essentially penalized Lp norms, involving appropriate functions of spike trains, with penalties associated with time-warping. The notions of means and variances of spike trains are then defined based on the new metrics when p = 2 (corresponding to the “Euclidean distance”). Using some restrictive conditions, we present an efficient recursive algorithm, termed Matching-Minimization algorithm, to compute the sample mean of a set of spike trains with arbitrary numbers of spikes. The proposed metrics as well as the mean spike trains are demonstrated using simulations as well as an experimental recording from the motor cortex. It is found that all these methods achieve desirable performance and the results support the success of this novel framework.

[1]  Benjamin Schrauwen,et al.  Linking non-binned spike train kernels to several existing spike train metrics , 2006, ESANN.

[2]  Anuj Srivastava,et al.  Towards statistical summaries of spike train data , 2011, Journal of Neuroscience Methods.

[3]  G. Seber Multivariate observations / G.A.F. Seber , 1983 .

[4]  Conor J. Houghton,et al.  A New Multineuron Spike Train Metric , 2008, Neural Computation.

[5]  G. P. Moore,et al.  Neuronal spike trains and stochastic point processes. I. The single spike train. , 1967, Biophysical journal.

[6]  G. P. Moore,et al.  Neuronal spike trains and stochastic point processes. II. Simultaneous spike trains. , 1967, Biophysical journal.

[7]  D. Mumford,et al.  A Metric on Shape Space with Explicit Geodesics , 2007, 0706.4299.

[8]  D. Aronov Fast algorithm for the metric-space analysis of simultaneous responses of multiple single neurons , 2003, Journal of Neuroscience Methods.

[9]  Anuj Srivastava,et al.  Analysis of planar shapes using geodesic paths on shape spaces , 2004, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[10]  J. D. Hunter,et al.  Amplitude and frequency dependence of spike timing: implications for dynamic regulation. , 2003, Journal of neurophysiology.

[11]  Mark C. W. van Rossum,et al.  A Novel Spike Distance , 2001, Neural Computation.

[12]  D. Mumford,et al.  An overview of the Riemannian metrics on spaces of curves using the Hamiltonian approach , 2006, math/0605009.

[13]  Jonathan D. Victor,et al.  Dynamic programming algorithms for comparing multineuronal spike trains via cost-based metrics and alignments , 2007, Journal of Neuroscience Methods.

[14]  Douglas Curran-Everett,et al.  Guidelines for reporting statistics in journals published by the American Physiological Society. , 2004, American journal of physiology. Endocrinology and metabolism.

[15]  José Carlos Príncipe,et al.  A comparison of binless spike train measures , 2010, Neural Computing and Applications.

[16]  G. Laurent,et al.  Who reads temporal information contained across synchronized and oscillatory spike trains? , 1998, Nature.

[17]  Conor J. Houghton,et al.  Studying spike trains using a van Rossum metric with a synapse-like filter , 2009, Journal of Computational Neuroscience.

[18]  Uri T Eden,et al.  A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. , 2005, Journal of neurophysiology.

[19]  H. Karcher Riemannian center of mass and mollifier smoothing , 1977 .

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

[21]  R E Kass,et al.  Recursive bayesian decoding of motor cortical signals by particle filtering. , 2004, Journal of neurophysiology.

[22]  Jonathan D Victor,et al.  Non-Euclidean properties of spike train metric spaces. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  Paul H. E. Tiesinga,et al.  A New Correlation-Based Measure of Spike Timing Reliability , 2002, Neurocomputing.

[24]  Peter Dayan,et al.  Theoretical Neuroscience: Computational and Mathematical Modeling of Neural Systems , 2001 .

[25]  Robert E Kass,et al.  Statistical issues in the analysis of neuronal data. , 2005, Journal of neurophysiology.

[26]  Antonio Politi,et al.  Measuring spike train synchrony , 2007, Journal of Neuroscience Methods.

[27]  Jonathan D. Victor,et al.  Metric-space analysis of spike trains: theory, algorithms and application , 1998, q-bio/0309031.

[28]  F. Mechler,et al.  Neural coding of spatial phase in V1 of the macaque monkey. , 2003, Journal of neurophysiology.

[29]  J. Victor,et al.  Nature and precision of temporal coding in visual cortex: a metric-space analysis. , 1996, Journal of neurophysiology.

[30]  José Carlos Príncipe,et al.  A Reproducing Kernel Hilbert Space Framework for Spike Train Signal Processing , 2009, Neural Computation.

[31]  Anuj Srivastava,et al.  Riemannian Analysis of Probability Density Functions with Applications in Vision , 2007, 2007 IEEE Conference on Computer Vision and Pattern Recognition.

[32]  Dukhwan Lim,et al.  Measurement of temporal regularity of spike train responses in auditory nerve fibers of the green treefrog , 1994, Journal of Neuroscience Methods.

[33]  Robert E. Kass,et al.  A Spike-Train Probability Model , 2001, Neural Computation.

[34]  Marcelo O. Magnasco,et al.  A Fast p Spike Alignment Metric , 2010, Neural Computation.

[35]  William Bialek,et al.  Spikes: Exploring the Neural Code , 1996 .

[36]  R. Kass,et al.  Geometrical Foundations of Asymptotic Inference: Kass/Geometrical , 1997 .

[37]  R Quian Quiroga,et al.  Event synchronization: a simple and fast method to measure synchronicity and time delay patterns. , 2002, Physical review. E, Statistical, nonlinear, and soft matter physics.

[38]  N. Čencov Statistical Decision Rules and Optimal Inference , 2000 .

[39]  J. S. Hunter,et al.  Statistics for Experimenters: An Introduction to Design, Data Analysis, and Model Building. , 1979 .

[40]  Nicholas G Hatsopoulos,et al.  Template-based spike pattern identification with linear convolution and dynamic time warping. , 2007, Journal of neurophysiology.

[41]  Dimitri P. Bertsekas,et al.  Dynamic Programming and Optimal Control, Two Volume Set , 1995 .

[42]  John W. Tukey,et al.  Exploratory Data Analysis , 1980, ACM SIGSPATIAL International Workshop on Advances in Geographic Information Systems.