The Annals of Applied Statistics

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

Robert E. Kass, Ryan C. Kelly, and Wei-Liem Loh

Full-text: Open access


Neural spike trains, which are sequences of very brief jumps in voltage across the cell membrane, were one of the motivating applications for the development of point process methodology. Early work required the assumption of stationarity, but contemporary experiments often use time-varying stimuli and produce time-varying neural responses. More recently, many statistical methods have been developed for nonstationary neural point process data. There has also been much interest in identifying synchrony, meaning events across two or more neurons that are nearly simultaneous at the time scale of the recordings. A natural statistical approach is to discretize time, using short time bins, and to introduce loglinear models for dependency among neurons, but previous use of loglinear modeling technology has assumed stationarity. We introduce a succinct yet powerful class of time-varying loglinear models by (a) allowing individual-neuron effects (main effects) to involve time-varying intensities; (b) also allowing the individual-neuron effects to involve autocovariation effects (history effects) due to past spiking, (c) assuming excess synchrony effects (interaction effects) do not depend on history, and (d) assuming all effects vary smoothly across time. Using data from the primary visual cortex of an anesthetized monkey, we give two examples in which the rate of synchronous spiking cannot be explained by stimulus-related changes in individual-neuron effects. In one example, the excess synchrony disappears when slow-wave “up” states are taken into account as history effects, while in the second example it does not. Standard point process theory explicitly rules out synchronous events. To justify our use of continuous-time methodology, we introduce a framework that incorporates synchronous events and provides continuous-time loglinear point process approximations to discrete-time loglinear models.

Article information

Ann. Appl. Stat. Volume 5, Number 2B (2011), 1262-1292.

First available in Project Euclid: 13 July 2011

Permanent link to this document

Digital Object Identifier

Mathematical Reviews number (MathSciNet)

Zentralblatt MATH identifier

Discrete-time approximation loglinear model marked process nonstationary point process simultaneous events spike train synchrony detection


Kass, Robert E.; Kelly, Ryan C.; Loh, Wei-Liem. Assessment of synchrony in multiple neural spike trains using loglinear point process models. Ann. Appl. Stat. 5 (2011), no. 2B, 1262--1292. doi:10.1214/10-AOAS429.

Export citation


  • Brillinger, D. R. (1988). Maximum likelihood analysis of spike trains of interacting nerve cells. Biol. Cyber. 59 189–200.
  • Brillinger, D. R. (1992). Nerve cell spike train data analysis: A progression of technique. J. Amer. Statist. Assoc. 87 260–271.
  • Brown, E. N., Barbieri, R., Eden, U. T. and Frank L. M. (2004). Likelihood methods for neural spike train data analysis. In Computational Neuroscience: A Comprehensive Approach 253–286. CRC Press, London.
  • Brown, E. N., Barbieri, R., Ventura, V., Kass, R. E. and Frank L. M. (2001). A note on the time-rescaling theorem and its implications for neural data analysis. Neural Comput. 14 325–346.
  • Daley, D. J. and Vere-Jones, D. (2002). An Introduction to the Theory of Point Processes, Vol. 1: Elementary Theory and Methods, 2nd ed. Springer, New York.
  • Geman, S., Amarasingham, A., Harrison, M. and Hatsopoulos, N. (2010). The statistical analysis of temporal resolution in the nervous system. Unpublished manuscript.
  • Grün, S. (2009). Data-driven significance estimation for precise spike correlation. J. Neurophysiol. 101 1126–1140.
  • Grün, S., Diesmann, M. and Aertsen, A. (2002a). Unitary events in multiple single-neuron spiking activity: I. detection and significance. Neural Comput. 14 43–80.
  • Grün, S., Diesmann, M. and Aertsen, A. (2002b). Unitary events in multiple single-neuron spiking activity: II. nonstationary data. Neural Comput. 14 81–119.
  • Harris, K. D., Henze, D. A., Csicsvari, J., Hirase, H. and Buzsaki, G. (2000). Accuracy of tetrode spike separation as determined by simultaneous intracellular and extracellular measurements. J. Neurophysiol. 84 401–414.
  • Kass, R. E., Ventura, V. and Brown, E. N. (2005). Statistical issues in the analysis of neuronal data. J. Neurophysiol. 94 8–25.
  • Kelly, R. C., Smith, M. A., Kass, R. E. and Lee, T. S. (2010). Local field potentials indicate network state and account for neuronal response variability. J. Comput. Neurosci. 29 567–579.
  • Kelly, R. C., Smith, M. A., Samonds, J. M., Kohn, A., Bonds, A. B., Movshon, J. A. and Lee, T. S. (2007). Comparison of recordings from microelectrode arrays and single electrodes in the visual cortex. J. Neurosci. 27 261–264.
  • Martignon, L., Deco, G., Laskey, K., Diamond, M., Freiwald, W. and Vaadia, E. (2000). Neural coding: Higher-order temporal patterns in the neurostatistics of cell assemblies. Neural Comput. 12 2621–2653.
  • Paninski, L., Brown, E. N., Iyengar, S. and Kass, R. E. (2009). Statistical models of spike trains. In Stochastic Methods in Neuroscience (C. Liang and G. J. Lord, eds.) 278–303. Clarendon Press, Oxford.
  • Schneidman, E., Berry, M. J., Segev, R. and Bialek, W. (2006). Weak pairwise correlations imply strongly correlated network states in a neural population. Nature 440 1007–1012.
  • Solo, V. (2007). Likelihood functions for multivariate point processes with coincidences. In Proc. 46th IEEE Conf. on Decision and Control 4245–4250. IEEE, New Orleans.
  • Ventura, V. (2009). Traditional waveform based spike sorting yields biased rate code estimation. Proc. Natl. Acad. Sci. USA 106 6921–6926.
  • Ventura, V., Cai, C. and Kass, R. E. (2005). Statistical assessment of time-varying dependency between two neurons. J. Neurophysiol. 94 2940–2947.