Comparison of a Bayesian and a regression model for stimulus classification
© Köpcke et al; licensee BioMed Central Ltd. 2011
Published: 18 July 2011
The central nervous system relies on spikes of retinal ganglion cells (RGC) as the only source of information about the visual environment. Therefore, the RGC response must contain information about e.g. velocity and luminance of an object. A common approach to find out which spike train features encode different stimulus properties is to estimate stimuli by assigning the observed spike trains to stimulus classes according to their features . Here, we compare two statistical methods – a Bayesian and a categorical regression model – for stimulus estimation.
Bayesian stimulus estimation has been used in neuroscientific literature for many years . Its basic idea is to determine the most probable stimulus given the observed response by using prior knowledge about the probability of stimulus occurrence combined with the statistics of response properties, conditioned on the stimulus. The second approach, a regression model based on categorical ordered data, is new to the field of neuronal stimulus estimation. The so called bivariate cumulative probit model  assumes that the bivariate ordinal response variable (e.g. the combination of ordered light intensity and velocity) is a discretization of a continuous variable that cannot be observed directly but can be described by a standard linear model with spike train features as covariates.
Both approaches led to similar results both for the classification based on the population response combining all 114 cells and for almost all individual cell responses. On average, the Bayesian model yielded slightly higher classification performances than the regression model for the population estimation and most of the single cells. Additionally, less CPU time was needed for the Bayesian analysis. The regression model, however, yielded on average a slightly lower standard deviation of correctly classified spike trains. So both methods proved to be well suited to analyze neural coding.
Comparing both models, the Bayesian approach has the advantage that it does not require a specific order of the response and can include prior knowledge of stimulus occurrence. On the other hand, the cumulative probit model allows for a generalization to higher dimensional responses without much computational effort. It also has the advantage to be intuitively comprehensible and flexible with respect to the number and choice of covariates.
- Rieke F, Warland D, de Ruyter van Steveninck R, Bialek W: Spikes – Exploring the Neural Code. 1999, The MIT PressGoogle Scholar
- Zhang K, Ginzburg I, McNaughton BL, Sejnowski TJ: Interpreting neuronal population activity by reconstruction: Unified framework with application to hippocampal place cells. Journal of Neurophysiology. 1998, 79: 1017-1044.PubMedGoogle Scholar
- Kim K: A Bivariate Cumulative Probit Regression Model for Ordered Categorical Data. Statistics in Medicine. 1995, 14: 1341-1352. 10.1002/sim.4780141207.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.