Function Identification in Neuron Populations via Information Bottleneck

Sreenivas Kartik Buddha, Michael Christoph Gastpar
MDPI AG, 2013
Journal paper

Abstract

It is plausible to hypothesize that the spiking responses of certain neurons represent functions of the spiking signals of other neurons. A natural ensuing question concerns how to use experimental data to infer what kind of a function is being computed. Model-based approaches typically require assumptions on how information is represented. By contrast, information measures are sensitive only to relative behavior: information is unchanged by applying arbitrary invertible transformations to the involved random variables. This paper develops an approach based on the information bottleneck method that attempts to find such functional relationships in a neuron population. Specifically, the information bottleneck method is used to provide appropriate compact representations which can then be parsed to infer functional relationships. In the present paper, the parsing step is specialized to the case of remapped-linear functions. The approach is validated on artificial data and then applied to recordings from the motor cortex of a macaque monkey performing an arm-reaching task. Functional relationships are identified and shown to exhibit some degree of persistence across multiple trials of the same experiment.

Official source

https://infoscience.epfl.ch/record/186370?ln=en

About this result

This page is automatically generated and may contain information that is not correct, complete, up-to-date, or relevant to your search query. The same applies to every other page on this website. Please make sure to verify the information with EPFL's official sources.

Related concepts

Related publications

STDP in adaptive neurons gives close-to-optimal information transmission

Wulfram Gerstner, Guillaume Hennequin, Jean-Pascal Théodor Pfister

Spike-frequency adaptation is known to enhance the transmission of information in sensory spiking neurons by rescaling the dynamic range for input processing, matching it to the temporal statistics of the sensory stimulus. Achieving maximal information transmission has also been recently postulated as a role for spike-timing-dependent plasticity (STDP). However, the link between optimal plasticity and STDP in cortex remains loose, as does the relationship between STDP and adaptation processes. We investigate how STDP, as described by recent minimal models derived from experimental data, influences the quality of information transmission in an adapting neuron. We show that a phenomenological model based on triplets of spikes yields almost the same information rate as an optimal model specially designed to this end. In contrast, the standard pair-based model of STDP does not improve information transmission as much. This result holds not only for additive STDP with hard weight bounds, known to produce bimodal distributions of synaptic weights, but also for weight-dependent STDP in the context of unimodal but skewed weight distributions. We analyze the similarities between the triplet model and the optimal learning rule, and find that the triplet effect is an important feature of the optimal model when the neuron is adaptive. If STDP is optimized for information transmission, it must take into account the dynamical properties of the postsynaptic cell, which might explain the target-cell specificity of STDP. In particular, it accounts for the differences found in vitro between STDP at excitatory synapses onto principal cells and those onto fast-spiking interneurons.

Frontiers Research Foundation2010

Improved Similarity Measures for Small Sets of Spike Trains

Joao Emanuel Felipe Gerhard, Wulfram Gerstner, Skander Mensi, Richard Naud

Multiple measures have been developed to quantify the similarity between two spike trains. These measures have been used for the quantification of the mismatch between neuron models and experiments as well as for the classification of neuronal responses in neuroprosthetic devices and electrophysiological experiments. Frequently only a few spike trains are available in each class. We derive analytical expressions for the small-sample bias present when comparing estimators of the time-dependent firing intensity. We then exploit analogies between the comparison of firing intensities and previously used spike train metrics and show that improved spike train measures can be successfully used for fitting neuron models to experimental data, for comparisons of spike trains, and classification of spike train data. In classification tasks, the improved similarity measures can increase the recovered information. We demonstrate that when similarity measures are used for fitting mathematical models, all previous methods systematically underestimate the noise. Finally, we show a striking implication of this deterministic bias by reevaluating the results of the single-neuron prediction challenge.

Massachusetts Institute of Technology Press2011

The coordinated, collective spiking activity of neuronal populations encodes and processes information. One approach towards understanding such population based computation is to fit statistical models to simultaneously recorded spike trains and use these models to make inferences about correlations, functional connectivity, and so forth. Any statistical model must, prior to making inferences from it, be validated by an appropriate goodness of fit measure, but as yet there is no commonly agreed upon statistical sufficiency test for neuronal population models. For single neuron spike trains, the time rescaling theorem provides a goodness of fit test consistent with their point process nature. Interspike intervals (ISIs) are rescaled, as a function of the model's conditional intensity function (time and history dependent spike probability) to be independent and exponentially distributed if the model is accurate. A Kolmogorov Smirnov (KS) test between the rescaled ISIs and the exponential distribution is then used to check goodness of fit. The time rescaling theorem applies to models of single spike trains, however some researchers have used it to validate population models, by applying it separately to each neuron in the population. Here we demonstrate that this approach is insufficient and that tests based upon single neuron statistics can not validate a population model. If interactions between neurons are ignored univariate tests may be erroneously passed, even when the interactions are strong. To remedy this we present the multivariate version of the time rescaling theorem and develop a practical statistical test for implementing it. The test relies on the fact that if the conditional intensity functions for all neurons are modeled correctly then the rescaled ISIs should be independent not only within a single spike train, but also across spike trains. Testing for such independence across spike trains determines if the collective activity of the neuronal population is captured by the model. Using Generalized Linear Models (GLMs) fitted to both simulated spike trains and spike trains recorded experimentally in monkey V1 cortex we demonstrate that population models are much more reliably validated by the multivariate test. This suggests that neuronal populations must be considered collectively as important features of their dynamics can only be detected from a multivariate viewpoint.

2010