Skip to content


  • Poster presentation
  • Open Access

A maximum likelihood estimator of neural network synaptic weights

  • 1 and
  • 1
BMC Neuroscience201314 (Suppl 1) :P59

  • Published:


  • Maximum Likelihood Estimator
  • Reverse Engineering
  • Synaptic Weight
  • Gibbs Distribution
  • Collective Dynamic

The statistics of spikes in a neuronal network is constrained on one hand by the stimulus and shared noise, and on the other hand by neuron interactions and collective dynamics. The join spike statistics and its spatio-temporal correlations can be explicitly computed in conductance-based Integrate-and-Fire models [1, 2]. The probability distribution of spike is a Gibbs distribution (in its most general definition allowing to consider non-stationarity) which encompasses existing statistical models such as Maximum Entropy models or Generalized-Linear Models.

Moreover, the dependence of spike statistics in network parameters such as synaptic weights and stimulus is explicit.

Here, we address the following reverse engineering problem. Given a conductance-based Integrate-and-Fire model as above where the spike statistics dependence on synaptic weights is known, can one reconstruct this network of synaptic weights from the observation of a raster plot generated by the network ? We have solved this inverse problem using an explicit expression of a maximum likelihood estimator based on the Newton-Raphson method. This estimator employs analytically computed gradients and Hessian of the likelihood function given by the product of conditional probabilities. The explicit form of these conditional probabilities can be found in [1]. Our results show that this method allows to estimate the set of connections weights knowing the input, the noise distribution and the leak function. Moreover, we found that, in a log scale scheme, the estimation mean percentage error Err decreases linearly with observation time T (Figure 1).
Figure 1
Figure 1

A: The mean percentage error is calculated over the estimated weights of 10 fully connected neurons (100 weights) using for each point a different bloc raster of same or different sizes. B. The error bars correspond to the variation of the estimated weights function of the real weights for an observation time T = 500.

This estimator is based on a plausible probabilistic model of spiking activity, and not a Poisson likelihood processing. So, it offers a flexible framework that should allow better statistical analysis of real data.



This work was supported by INRIA, ERC-NERVI number 227747, KEOPS ANR-CONICYT and European Union Project # FP7-269921 (BrainScales), Renvision grant agreement N 600847 and Mathemacs FP7-ICT_2011.9.7.

Authors’ Affiliations

NeuroMathComp team (INRIA, UNSA LJAD), Sophia Antipolis, France


  1. Cessac B: Statistics of spike trains in conductance-based neural networks: Rigorous results. The journal of Mathematical Neuroscience. 2011, 1: 8-10.1186/2190-8567-1-8.View ArticlePubMedGoogle Scholar
  2. Cofré R, Cessac B: Dynamics and spike trains statistics in conductance-based Integrate-and-Fire neural networks with chemical and electric synapses, to appear in. Chaos, Solitons and Fractals. 2013Google Scholar


© Taouali and Cessac; licensee BioMed Central Ltd. 2013

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.