Efficient modeling of neural activity using coupled renewal processes
-
1
Ecole Polytechnique Fédérale de Lausanne, Brain-Mind-Institute, Switzerland
Generalized Linear Models (GLMs) are stochastic models that have been successfully used to model neural activity of single cells and populations. Typically, this requires the spike train to be binned into a binary time series that is modeled as a Bernoulli-type GLM. Unless a loss in temporal precision is acceptable, the bin size has to be chosen on the order of milliseconds. This renders an application to long-term recordings infeasible for computational reasons.
We propose a new, binning-free application of the GLM framework to model neural activity by noting that spike trains can equally be characterized by the sequence of continuous-valued inter-spike intervals (ISIs). It was shown that single-neuron spike trains can often be sufficiently modeled as a non-stationary renewal process. Here, we focus on the Gamma distribution for the inter-spike intervals as it allows the accommodation of a relative refractory period and includes the inhomogeneous Poisson process as a special case. We generalize the idea of a modulated Gamma process to the population level by including cross-couplings between the neurons in case of recordings from multiple neurons in parallel.
In a Gamma-GLM, each observed inter-spike interval is modeled as a sample from a Gamma distribution whose mean can vary with time, e.g. through modulation by an external stimulus. Traditionally, ISIs are assumed to be independent. However, this renewal assumption can be relaxed by conditioning on the durations of previous intervals. An inter-neuron coupling is introduced by additionally conditioning on the times of spikes of the other neurons. The Gamma distribution and other canonical distributions for inter-spike intervals are part of the exponential family so that efficient maximum likelihood solutions for all model parameters can be obtained using the framework of Generalized Linear Models. Considering the ratio between bin sizes for time-discrete models and the typical length of cortical inter-spike intervals, the computational demand can be reduced by several orders of magnitude. Goodness-of-fit is assessed by adapting standard tests for point process models and evaluating spike prediction performance based on cross-validation.
We demonstrate the versatility of the coupled Gamma model for estimating functional connectivity from multi-electrode recordings. With a data set obtained from multi-electrode recordings of the visual system of the awake monkey, we show that the test power of the Gamma-GLM for detecting functional connectivity links is comparable to the performance of a classical discrete-time GLM albeit requiring much less computational time and memory resources. Using the two models on the same data set, we are able to estimate the unknown, underlying connection density of the network. Finally, we discuss general limitations of coupled Gamma-GLMs on the type of interactions that can be faithfully represented.
Acknowledgements
Felipe Gerhard thanks Gordon Pipa for useful discussions. The authors thank Sergio Neuenschwander and Bruss Lima for supplying the experimental data. Felipe Gerhard is supported by the Swiss National Science Foundation (SNSF) under grant number 200020-117975.
Keywords:
data analysis and machine learning,
functional connectivity,
Generalized Linear Models
Conference:
BC11 : Computational Neuroscience & Neurotechnology Bernstein Conference & Neurex Annual Meeting 2011, Freiburg, Germany, 4 Oct - 6 Oct, 2011.
Presentation Type:
Poster
Topic:
data analysis and machine learning (please use "data analysis and machine learning" as keyword)
Citation:
Gerhard
F and
Gerstner
W
(2011). Efficient modeling of neural activity using coupled renewal processes.
Front. Comput. Neurosci.
Conference Abstract:
BC11 : Computational Neuroscience & Neurotechnology Bernstein Conference & Neurex Annual Meeting 2011.
doi: 10.3389/conf.fncom.2011.53.00197
Copyright:
The abstracts in this collection have not been subject to any Frontiers peer review or checks, and are not endorsed by Frontiers.
They are made available through the Frontiers publishing platform as a service to conference organizers and presenters.
The copyright in the individual abstracts is owned by the author of each abstract or his/her employer unless otherwise stated.
Each abstract, as well as the collection of abstracts, are published under a Creative Commons CC-BY 4.0 (attribution) licence (https://creativecommons.org/licenses/by/4.0/) and may thus be reproduced, translated, adapted and be the subject of derivative works provided the authors and Frontiers are attributed.
For Frontiers’ terms and conditions please see https://www.frontiersin.org/legal/terms-and-conditions.
Received:
19 Aug 2011;
Published Online:
04 Oct 2011.
*
Correspondence:
Mr. Felipe Gerhard, Ecole Polytechnique Fédérale de Lausanne, Brain-Mind-Institute, Lausanne-EPFL, 1015, Switzerland, felipe.gerhard@epfl.ch