A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects.

Multiple factors simultaneously affect the spiking activity of individual neurons. Determining the effects and relative importance of these factors is a challenging problem in neurophysiology. We propose a statistical framework based on the point process likelihood function to relate a neuron's spiking probability to three typical covariates: the neuron's own spiking history, concurrent ensemble activity, and extrinsic covariates such as stimuli or behavior. The framework uses parametric models of the conditional intensity function to define a neuron's spiking probability in terms of the covariates. The discrete time likelihood function for point processes is used to carry out model fitting and model analysis. We show that, by modeling the logarithm of the conditional intensity function as a linear combination of functions of the covariates, the discrete time point process likelihood function is readily analyzed in the generalized linear model (GLM) framework. We illustrate our approach for both GLM and non-GLM likelihood functions using simulated data and multivariate single-unit activity data simultaneously recorded from the motor cortex of a monkey performing a visuomotor pursuit-tracking task. The point process framework provides a flexible, computationally efficient approach for maximum likelihood estimation, goodness-of-fit assessment, residual analysis, model selection, and neural decoding. The framework thus allows for the formulation and analysis of point process models of neural spiking activity that readily capture the simultaneous effects of multiple covariates and enables the assessment of their relative importance.

[1]  Pia Veldt Larsen,et al.  In All Likelihood: Statistical Modelling and Inference Using Likelihood , 2003 .

[2]  A. P. Georgopoulos,et al.  Movement parameters and neural activity in motor cortex and area 5. , 1994, Cerebral cortex.

[3]  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.

[4]  M. Stone An Asymptotic Equivalence of Choice of Model by Cross‐Validation and Akaike's Criterion , 1977 .

[5]  B L McNaughton,et al.  Dynamics of the hippocampal ensemble code for space. , 1993, Science.

[6]  Yosihiko Ogata,et al.  On Lewis' simulation method for point processes , 1981, IEEE Trans. Inf. Theory.

[7]  B. Hille Ionic channels of excitable membranes , 2001 .

[8]  J. O’Keefe,et al.  Phase relationship between hippocampal place units and the EEG theta rhythm , 1993, Hippocampus.

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

[10]  D. V. van Essen,et al.  Neuronal responses to static texture patterns in area V1 of the alert macaque monkey. , 1992, Journal of neurophysiology.

[11]  H. Akaike,et al.  Information Theory and an Extension of the Maximum Likelihood Principle , 1973 .

[12]  G. Schwarz Estimating the Dimension of a Model , 1978 .

[13]  B. Efron The Estimation of Prediction Error , 2004 .

[14]  W. R. Buckland,et al.  Distributions in Statistics: Continuous Multivariate Distributions , 1973 .

[15]  D. G. Simpson,et al.  The Statistical Analysis of Discrete Data , 1989 .

[16]  Wilson Truccolo,et al.  Motor “Binding:” Do Functional Assemblies in Primary Motor Cortex Have a Role? , 2003, Neuron.

[17]  S. Scott The role of primary motor cortex in goal-directed movements: insights from neurophysiological studies on non-human primates , 2003, Current Opinion in Neurobiology.

[18]  M. Peruggia Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach (2nd ed.) , 2003 .

[19]  Thomas M. Cover,et al.  Elements of Information Theory , 2005 .

[20]  L. Paninski,et al.  Sequential movement representations based on correlated neuronal activity , 2003, Experimental Brain Research.

[21]  Mark Berman,et al.  Approximating Point Process Likelihoods with Glim , 1992 .

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

[23]  J. Mendel Lessons in Estimation Theory for Signal Processing, Communications, and Control , 1995 .

[24]  J. Csicsvari,et al.  Massively parallel recording of unit and local field potentials with silicon-based electrodes. , 2003, Journal of neurophysiology.

[25]  L. Paninski,et al.  Spatiotemporal tuning of motor cortical neurons for hand position and velocity. , 2004, Journal of neurophysiology.

[27]  J. Donoghue,et al.  Neural discharge and local field potential oscillations in primate motor cortex during voluntary movements. , 1998, Journal of neurophysiology.

[28]  D. Brillinger Maximum likelihood analysis of spike trains of interacting nerve cells , 2004, Biological Cybernetics.

[29]  A. Riehle,et al.  Precise spike synchronization in monkey motor cortex involved in preparation for movement , 1999, Experimental Brain Research.

[30]  R. Reid,et al.  Predicting Every Spike A Model for the Responses of Visual Neurons , 2001, Neuron.

[31]  L. Paninski,et al.  Information about movement direction obtained from synchronous activity of motor cortical neurons. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[32]  A. Jackson,et al.  Synchrony between Neurons with Similar Muscle Fields in Monkey Motor Cortex , 2003, Neuron.

[33]  Simon Grant,et al.  Reduced influence of the ipsilateral ear on spatial tuning of auditory neurons in the albino superior colliculus: a knock-on effect of anomalies of the acoustic chiasm? , 2003, Experimental Brain Research.

[34]  Artur Luczak,et al.  Multivariate receptive field mapping in marmoset auditory cortex , 2004, Journal of Neuroscience Methods.

[35]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[36]  J. Donoghue,et al.  Neuronal Interactions Improve Cortical Population Coding of Movement Direction , 1999, The Journal of Neuroscience.

[37]  B. McNaughton,et al.  Experience-dependent, asymmetric expansion of hippocampal place fields. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

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

[39]  T. Ebner,et al.  Temporal encoding of movement kinematics in the discharge of primate primary motor and premotor neurons. , 1995, Journal of neurophysiology.

[40]  V. Solo,et al.  Contrasting Patterns of Receptive Field Plasticity in the Hippocampus and the Entorhinal Cortex: An Adaptive Filtering Approach , 2002, The Journal of Neuroscience.

[41]  J. O'Keefe,et al.  The hippocampus as a spatial map. Preliminary evidence from unit activity in the freely-moving rat. , 1971, Brain research.

[42]  A B Schwartz,et al.  Motor cortical representation of speed and direction during reaching. , 1999, Journal of neurophysiology.

[43]  H. Wilson Simplified dynamics of human and mammalian neocortical neurons. , 1999, Journal of theoretical biology.

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

[45]  Emery N. Brown,et al.  Dynamic Analysis of Neural Encoding by Point Process Adaptive Filtering , 2004, Neural Computation.

[46]  Adrienne L. Fairhall,et al.  What Causes a Neuron to Spike? , 2003, Neural Computation.

[47]  John P. Donoghue,et al.  Connecting cortex to machines: recent advances in brain interfaces , 2002, Nature Neuroscience.

[48]  Adrienne L. Fairhall,et al.  Computation in a Single Neuron: Hodgkin and Huxley Revisited , 2002, Neural Computation.

[49]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[50]  H. Jones,et al.  Visual cortical mechanisms detecting focal orientation discontinuities , 1995, Nature.

[51]  Jianfeng Feng,et al.  Computational neuroscience , 1986, Behavioral and Brain Sciences.

[52]  M. Quirk,et al.  Experience-Dependent Asymmetric Shape of Hippocampal Receptive Fields , 2000, Neuron.

[53]  M. A. Tanner,et al.  Tools for Statistical Inference: Methods for the Exploration of Posterior Distributions and Likelihood Functions, 3rd Edition , 1998 .

[54]  E. S. Chornoboy,et al.  Maximum likelihood identification of neural point process systems , 1988, Biological Cybernetics.

[55]  M. Quirk,et al.  Construction and analysis of non-Poisson stimulus-response models of neural spiking activity , 2001, Journal of Neuroscience Methods.

[56]  Eric R. Ziegel,et al.  Generalized Linear Models , 2002, Technometrics.

[57]  Roger Sauter,et al.  In All Likelihood , 2002, Technometrics.

[58]  Jerald D. Kralik,et al.  Chronic, multisite, multielectrode recordings in macaque monkeys , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[59]  J L Gallant,et al.  Sparse coding and decorrelation in primary visual cortex during natural vision. , 2000, Science.

[60]  John A. Nelder,et al.  Generalized linear models. 2nd ed. , 1993 .