Non-parametric estimation of the spiking rate in systems of interacting neurons

2016/04/25 by Hodara, Pierre, Krell, Nathalie, Löcherbach, Eva
#60J75 #62G05 #62M05 #FOS: Mathematics #Statistics Theory (math.ST)

paper · doi:10.48550/arxiv.1604.07300

Abstract

We consider a model of interacting neurons where the membrane potentials of the neurons are described by a multidimensional piecewise deterministic Markov process (PDMP) with values in \mathbb R^N, where N is the number of neurons in the network. A deterministic drift attracts each neuron's membrane potential to an equilibrium potential m. When a neuron jumps, its membrane potential is reset to 0, while the other neurons receive an additional amount of potential (1)/(N). We are interested in the estimation of the jump (or spiking) rate of a single neuron based on an observation of the membrane potentials of the N neurons up to time t. We study a Nadaraya-Watson type kernel estimator for the jump rate and establish its rate of convergence in L² . This rate of convergence is shown to be optimal for a given Hölder class of jump rate functions. We also obtain a central limit theorem for the error of estimation. The main probabilistic tools are the uniform ergodicity of the process and a fine study of the invariant measure of a single neuron.

Related

• Nonparametric estimation of the jump rate in mean field interacting systems of neurons
• On a toy model of interacting neurons
• Strong error bounds for the convergence to its mean field limit for systems of interacting neurons in a diffusive scaling
• Fluctuations for mean field limits of interacting systems of spiking neurons
• Stochastic dynamics of a finite-size spiking neural network