1
|
Kostal L, Kovacova K. Estimation of firing rate from instantaneous interspike intervals. Neurosci Res 2024:S0168-0102(24)00085-3. [PMID: 38925356 DOI: 10.1016/j.neures.2024.06.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2024] [Revised: 05/27/2024] [Accepted: 06/19/2024] [Indexed: 06/28/2024]
Abstract
The rate coding hypothesis is the oldest and still one of the most accepted hypotheses of neural coding. Consequently, many approaches have been devised for the firing rate estimation, ranging from simple binning of the time axis to advanced statistical methods. Nonetheless the concept of firing rate, while informally understood, can be mathematically defined in several distinct ways. These definitions may yield mutually incompatible results unless implemented properly. Recently it has been shown that the notions of the instantaneous and the classical firing rates can be made compatible, at least in terms of their averages, by carefully discerning the time instant at which the neuronal activity is observed. In this paper we revisit the properties of instantaneous interspike intervals in order to derive several novel firing rate estimators, which are free of additional assumptions or parameters and their temporal resolution is 'locally self-adaptive'. The estimators are simple to implement and are numerically efficient even for very large sets of data.
Collapse
Affiliation(s)
- Lubomir Kostal
- Institute of Physiology of the Czech Academy of Sciences, Videnska 1083, Prague 4 14200, Czech Republic.
| | - Kristyna Kovacova
- Institute of Physiology of the Czech Academy of Sciences, Videnska 1083, Prague 4 14200, Czech Republic
| |
Collapse
|
2
|
Braun W, Matsuzaka Y, Mushiake H, Northoff G, Longtin A. Non-additive activity modulation during a decision making task involving tactic selection. Cogn Neurodyn 2022; 16:117-133. [PMID: 35116084 PMCID: PMC8807796 DOI: 10.1007/s11571-021-09702-0] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Revised: 06/28/2021] [Accepted: 07/14/2021] [Indexed: 12/01/2022] Open
Abstract
Human brain imaging has revealed that stimulus-induced activity does generally not simply add to the pre-stimulus activity, but rather builds in a non-additive way on this activity. Here we investigate this subject at the single neuron level and address the question whether and to what extent a strong form of non-additivity where activity drops post-cue is present in different areas of monkey cortex, including prefrontal and agranular frontal areas, during a perceptual decision making task involving action and tactic selection. Specifically we analyze spike train data recorded in vivo from the posterior dorsomedial prefrontal cortex (pmPFC), the supplementary motor area (SMA) and the presupplementary motor area (pre-SMA). For each neuron, we compute the ratio of the trial-averaged pre-stimulus spike count to the trial-averaged post-stimulus count. We also perform the ratio and averaging procedures in reverse order. We find that the statistics of these quantities behave differently across areas. pmPFC involved in tactic selection shows stronger non-additivity compared to the two other areas which more generically just increase their firing rate pos-stimulus. pmPFC behaved more similarly to pre-SMA, a likely consequence of the reciprocal connections between these areas. The trial-averaged ratio statistic was reproduced by a surrogate inhomogeneous Poisson process in which the measured trial-averaged firing rate for a given neuron is used as its time-dependent rate. Principal component analysis (PCA) of the trial-averaged firing rates of neuronal ensembles further reveals area-specific time courses of response to the stimulus, including latency to peak neural response, for the typical population activity. Our work demonstrates subtle forms of area-specific non-additivity based on the fine variability structure of pre- and post-stimulus spiking activity on the single neuron level. It also reveals significant differences between areas for PCA and surrogate analysis, complementing previous observations of regional differences based solely on post-stimulus responses. Moreover, we observe regional differences in non-additivity which are related to the monkey's successful tactic selection and decision making. Supplementary Information The online version contains supplementary material available at 10.1007/s11571-021-09702-0.
Collapse
Affiliation(s)
- Wilhelm Braun
- Institut für Genetik, Neural Network Dynamics and Computation, Universität Bonn, Kirschallee 1, 53115 Bonn, Germany.,Department of Physics and Centre for Neural Dynamics, University of Ottawa, 150 Louis-Pasteur Pvt, Ottawa, K1N 6N5 Canada
| | - Yoshiya Matsuzaka
- Division of Neuroscience, Faculty of Medicine, Tohoku Medical and Pharmaceutical University, 1-15-1 Fukumuro, Miyagino ward, Sendai, 983-8536 Japan
| | - Hajime Mushiake
- Department of Physiology, Graduate School of Medicine, Tohoku University, Aoba Ward, Sendai, 981-8558 Japan
| | - Georg Northoff
- Mind, Brain Imaging and Neuroethics Research Unit, University of Ottawa Institute of Mental Health Research, Royal Ottawa Mental Health Centre, 1145 Carling Avenue, Ottawa, K1Z 7K4 Canada
| | - André Longtin
- Department of Physics and Centre for Neural Dynamics, University of Ottawa, 150 Louis-Pasteur Pvt, Ottawa, K1N 6N5 Canada
| |
Collapse
|
3
|
Budzinski R, Lopes S, Masoller C. Symbolic analysis of bursting dynamical regimes of Rulkov neural networks. Neurocomputing 2021. [DOI: 10.1016/j.neucom.2020.05.122] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
4
|
Mankin R, Rekker A. Effects of transient subordinators on the firing statistics of a neuron model driven by dichotomous noise. Phys Rev E 2020; 102:012103. [PMID: 32794976 DOI: 10.1103/physreve.102.012103] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Accepted: 06/15/2020] [Indexed: 06/11/2023]
Abstract
The behavior of a stochastic perfect integrate-and-fire (PIF) model of neurons is considered. The effect of temporally correlated random activity of synaptic inputs is modeled as a combination of an asymmetric dichotomous noise and a random operation time in the form of an inverse strictly increasing Lévy-type subordinator. Using a first-passage-time formulation, we find exact expressions for the output interspike interval (ISI) statistics. Particularly, it is shown that at some parameter regimes the ISI density exhibits a multimodal structure. Moreover, it is demonstrated that the coefficient of variation, the serial correlation coefficient, and the Fano factor display a nonmonotonic dependence on the mean input current μ, i.e., the ISI's regularity is maximized at an intermediate value of μ. The features of spike statistics, analytically revealed in our study, are compared with previously obtained results for a perfect integrate-and-fire neuron model driven by dichotomous noise (without subordination).
Collapse
Affiliation(s)
- Romi Mankin
- School of Natural Sciences and Health, Tallinn University, 29 Narva Road, 10120 Tallinn, Estonia
| | - Astrid Rekker
- School of Natural Sciences and Health, Tallinn University, 29 Narva Road, 10120 Tallinn, Estonia
| |
Collapse
|
5
|
Bernardi D, Lindner B. Receiver operating characteristic curves for a simple stochastic process that carries a static signal. Phys Rev E 2020; 101:062132. [PMID: 32688497 DOI: 10.1103/physreve.101.062132] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/17/2020] [Accepted: 05/14/2020] [Indexed: 11/07/2022]
Abstract
The detection of a weak signal in the presence of noise is an important problem that is often studied in terms of the receiver operating characteristic (ROC) curve, in which the probability of correct detection is plotted against the probability for a false positive. This kind of analysis is typically applied to the situation in which signal and noise are stochastic variables; the detection problem emerges, however, also often in a context in which both signal and noise are stochastic processes and the (correct or false) detection is said to take place when the process crosses a threshold in a given time window. Here we consider the problem for a combination of a static signal which has to be detected against a dynamic noise process, the well-known Ornstein-Uhlenbeck process. We give exact (but difficult to evaluate) quadrature expressions for the detection rates for false positives and correct detections, investigate systematically a simple sampling approximation suggested earlier, compare to an approximation by Stratonovich for the limit of high threshold, and briefly explore the case of multiplicative signal; all theoretical results are compared to extensive numerical simulations of the corresponding Langevin equation. Our results demonstrate that the sampling approximation provides a reasonable description of the ROC curve for this system, and it clarifies limit cases for the ROC curve.
Collapse
Affiliation(s)
- Davide Bernardi
- Bernstein Center for Computational Neuroscience Berlin, 10115 Berlin, Germany and Physics Department of Humboldt University Berlin, 12489 Berlin, Germany
| | - Benjamin Lindner
- Bernstein Center for Computational Neuroscience Berlin, 10115 Berlin, Germany and Physics Department of Humboldt University Berlin, 12489 Berlin, Germany
| |
Collapse
|
6
|
Braun W, Longtin A. Interspike interval correlations in networks of inhibitory integrate-and-fire neurons. Phys Rev E 2019; 99:032402. [PMID: 30999498 DOI: 10.1103/physreve.99.032402] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2018] [Indexed: 11/07/2022]
Abstract
We study temporal correlations of interspike intervals, quantified by the network-averaged serial correlation coefficient (SCC), in networks of both current- and conductance-based purely inhibitory integrate-and-fire neurons. Numerical simulations reveal transitions to negative SCCs at intermediate values of bias current drive and network size. As bias drive and network size are increased past these values, the SCC returns to zero. The SCC is maximally negative at an intermediate value of the network oscillation strength. The dependence of the SCC on two canonical schemes for synaptic connectivity is studied, and it is shown that the results occur robustly in both schemes. For conductance-based synapses, the SCC becomes negative at the onset of both a fast and slow coherent network oscillation. We then show by means of offline simulations using prerecorded network activity that a neuron's SCC is highly sensitive to its number of presynaptic inputs. Finally, we devise a noise-reduced diffusion approximation for current-based networks that accounts for the observed temporal correlation transitions.
Collapse
Affiliation(s)
- Wilhelm Braun
- Neural Network Dynamics and Computation, Institut für Genetik, Universität Bonn, Kirschallee 1, 53115 Bonn, Germany.,Department of Physics and Centre for Neural Dynamics, University of Ottawa, 598 King Edward, Ottawa K1N 6N5, Canada
| | - André Longtin
- Department of Physics and Centre for Neural Dynamics, University of Ottawa, 598 King Edward, Ottawa K1N 6N5, Canada
| |
Collapse
|
7
|
Kostal L, Lansky P, Stiber M. Statistics of inverse interspike intervals: The instantaneous firing rate revisited. CHAOS (WOODBURY, N.Y.) 2018; 28:106305. [PMID: 30384662 DOI: 10.1063/1.5036831] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/19/2018] [Accepted: 07/20/2018] [Indexed: 06/08/2023]
Abstract
The rate coding hypothesis is the oldest and still one of the most accepted and investigated scenarios in neuronal activity analyses. However, the actual neuronal firing rate, while informally understood, can be mathematically defined in several different ways. These definitions yield distinct results; even their average values may differ dramatically for the simplest neuronal models. Such an inconsistency, together with the importance of "firing rate," motivates us to revisit the classical concept of the instantaneous firing rate. We confirm that different notions of firing rate can in fact be compatible, at least in terms of their averages, by carefully discerning the time instant at which the neuronal activity is observed. Two general cases are distinguished: either the inspection time is synchronised with a reference time or with the neuronal spiking. The statistical properties of the instantaneous firing rate, including parameter estimation, are analyzed, and compatibility with the intuitively understood concept is demonstrated.
Collapse
Affiliation(s)
- Lubomir Kostal
- Institute of Physiology of the Czech Academy of Sciences, Videnska 1083, 14220 Prague 4, Czech Republic
| | - Petr Lansky
- Institute of Physiology of the Czech Academy of Sciences, Videnska 1083, 14220 Prague 4, Czech Republic
| | - Michael Stiber
- Institute of Physiology of the Czech Academy of Sciences, Videnska 1083, 14220 Prague 4, Czech Republic
| |
Collapse
|
8
|
Sub-threshold signal encoding in coupled FitzHugh-Nagumo neurons. Sci Rep 2018; 8:8276. [PMID: 29844354 PMCID: PMC5974132 DOI: 10.1038/s41598-018-26618-8] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2017] [Accepted: 05/15/2018] [Indexed: 11/09/2022] Open
Abstract
Despite intensive research, the mechanisms underlying the neural code remain poorly understood. Recent work has focused on the response of a single neuron to a weak, sub-threshold periodic signal. By simulating the stochastic FitzHugh-Nagumo (FHN) model and then using a symbolic method to analyze the firing activity, preferred and infrequent spike patterns (defined by the relative timing of the spikes) were detected, whose probabilities encode information about the signal. As not individual neurons but neuronal populations are responsible for sensory coding and information transfer, a relevant question is how a second neuron, which does not perceive the signal, affects the detection and the encoding of the signal, done by the first neuron. Through simulations of two stochastic FHN neurons we show that the encoding of a sub-threshold signal in symbolic spike patterns is a plausible mechanism. The neuron that perceives the signal fires a spike train that, despite having an almost random temporal structure, has preferred and infrequent patterns which carry information about the signal. Our findings could be relevant for sensory systems composed by two noisy neurons, when only one detects a weak external input.
Collapse
|