Measuring the signal-to-noise ratio of a neuron

Significance Neurons represent both signal and noise in binary electrical discharges termed action potentials. Hence, the standard signal-to-noise ratio (SNR) definition of signal amplitude squared and divided by the noise variance does not apply. We show that the SNR estimates a ratio of expected prediction errors. Using point process generalized linear models, we extend the standard definition to one appropriate for single neurons. In analyses of four neural systems, we show that single neuron SNRs range from −29 dB to −3 dB and that spiking history is often a more informative predictor of spiking propensity than the signal or stimulus activating the neuron. By generalizing the standard SNR metric, we make explicit the well-known fact that individual neurons are highly noisy information transmitters. The signal-to-noise ratio (SNR), a commonly used measure of fidelity in physical systems, is defined as the ratio of the squared amplitude or variance of a signal relative to the variance of the noise. This definition is not appropriate for neural systems in which spiking activity is more accurately represented as point processes. We show that the SNR estimates a ratio of expected prediction errors and extend the standard definition to one appropriate for single neurons by representing neural spiking activity using point process generalized linear models (PP-GLM). We estimate the prediction errors using the residual deviances from the PP-GLM fits. Because the deviance is an approximate χ2 random variable, we compute a bias-corrected SNR estimate appropriate for single-neuron analysis and use the bootstrap to assess its uncertainty. In the analyses of four systems neuroscience experiments, we show that the SNRs are −10 dB to −3 dB for guinea pig auditory cortex neurons, −18 dB to −7 dB for rat thalamic neurons, −28 dB to −14 dB for monkey hippocampal neurons, and −29 dB to −20 dB for human subthalamic neurons. The new SNR definition makes explicit in the measure commonly used for physical systems the often-quoted observation that single neurons have low SNRs. The neuron’s spiking history is frequently a more informative covariate for predicting spiking propensity than the applied stimulus. Our new SNR definition extends to any GLM system in which the factors modulating the response can be expressed as separate components of a likelihood function.

[1]  N. Nagelkerke,et al.  A note on a general definition of the coefficient of determination , 1991 .

[2]  Trevor Hastie,et al.  A Closer Look at the Deviance , 1987 .

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

[4]  Sean P. MacEvoy,et al.  Macaque V1 activity during natural vision: effects of natural scenes and saccades. , 2008, Journal of neurophysiology.

[5]  Steven Kay,et al.  Fundamentals Of Statistical Signal Processing , 2001 .

[6]  Vasilis Z. Marmarelis,et al.  Nonlinear Dynamic Modeling of Physiological Systems: Marmarelis/Nonlinear , 2004 .

[7]  Thippur V. Sreenivas,et al.  Signal-to-noise ratio estimation using higher-order moments , 2006, Signal Process..

[8]  Jacob Cohen,et al.  A power primer. , 1992, Psychological bulletin.

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

[10]  Emery N. Brown,et al.  Using Point Process Models to Compare Neural Spiking Activity in the Subthalamic Nucleus of Parkinson's Patients and a Healthy Primate , 2010, IEEE Transactions on Biomedical Engineering.

[11]  Y. Pawitan In all likelihood , 2001 .

[12]  L. Frank,et al.  Single Neurons in the Monkey Hippocampus and Learning of New Associations , 2003, Science.

[13]  Uri T Eden,et al.  A point process framework for relating neural spiking activity to spiking history, neural ensemble, and extrinsic covariate effects. , 2005, Journal of neurophysiology.

[14]  William Clouser Boyd,et al.  Fundamentals of Immunology , 1958 .

[15]  Ehsan S. Soo Principal Information Theoretic Approaches , 2000 .

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

[17]  Erik G. Larsson,et al.  Measuring the signal-to-noise ratio in magnetic resonance imaging: a caveat , 2004, Signal Process..

[18]  L. Paninski Maximum likelihood estimation of cascade point-process neural encoding models , 2004, Network.

[19]  Emery N. Brown,et al.  A Signal-to-Noise Ratio Estimator for Generalized Linear Model Systems , 2008 .

[20]  L. Optican,et al.  Temporal encoding of two-dimensional patterns by single units in primate inferior temporal cortex. III. Information theoretic analysis. , 1987, Journal of neurophysiology.

[21]  F. Windmeijer,et al.  An R-squared measure of goodness of fit for some common nonlinear regression models , 1997 .

[22]  P. McCullagh,et al.  Generalized Linear Models , 1984 .

[23]  Jakob H. Macke,et al.  Frontiers in Computational Neuroscience , 2022 .

[24]  Yves Rosseel,et al.  On the Definition of Signal-To-Noise Ratio and Contrast-To-Noise Ratio for fMRI Data , 2013, PloS one.

[25]  William Bialek,et al.  Spikes: Exploring the Neural Code , 1996 .

[26]  W. Newsome,et al.  The Variable Discharge of Cortical Neurons: Implications for Connectivity, Computation, and Information Coding , 1998, The Journal of Neuroscience.

[27]  Ken Kelley,et al.  Confidence Intervals for Standardized Effect Sizes: Theory, Application, and Implementation , 2007 .

[28]  J. Kent Information gain and a general measure of correlation , 1983 .

[29]  William R. Softky,et al.  The highly irregular firing of cortical cells is inconsistent with temporal integration of random EPSPs , 1993, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[30]  Geert Molenberghs,et al.  Prentice's Approach and the Meta‐Analytic Paradigm: A Reflection on the Role of Statistics in the Evaluation of Surrogate Endpoints , 2004, Biometrics.

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

[32]  David J. Anderson,et al.  Auditory cortical responses to electrical stimulation of the inferior colliculus: implications for an auditory midbrain implant. , 2006, Journal of neurophysiology.

[33]  Martina Mittlböck,et al.  Adjustments for R 2 -measures for Poisson regression models , 2000 .

[34]  E. N. Brown,et al.  nSTAT: Open-source neural spike train analysis toolbox for Matlab , 2012, Journal of Neuroscience Methods.

[35]  F. Collins,et al.  Policy: NIH plans to enhance reproducibility , 2014, Nature.

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

[37]  E. Brown Course 14 - Theory of Point Processes for Neural Systems , 2005 .

[38]  Simona Temereanca,et al.  Rapid Changes in Thalamic Firing Synchrony during Repetitive Whisker Stimulation , 2008, The Journal of Neuroscience.

[39]  Jennifer J. Richler,et al.  Effect size estimates: current use, calculations, and interpretation. , 2012, Journal of experimental psychology. General.

[40]  Vasilis Z. Marmarelis,et al.  Nonlinear Dynamic Modeling of Physiological Systems , 2004 .

[41]  Gary Simon,et al.  Additivity of Information in Exponential Family Probability Laws , 1973 .

[42]  Norman C. Beaulieu,et al.  Maximum likelihood estimation of SNR using digitally modulated signals , 2007, IEEE Transactions on Wireless Communications.

[43]  Uri T Eden,et al.  Analysis of between-trial and within-trial neural spiking dynamics. , 2008, Journal of neurophysiology.

[44]  Emery N. Brown,et al.  A Point Process Model for Auditory Neurons Considering Both Their Intrinsic Dynamics and the Spectrotemporal Properties of an Extrinsic Signal , 2011, IEEE Transactions on Biomedical Engineering.