Applications of the expectation-maximization algorithm to quantal analysis of postsynaptic potentials

The expectation-maximization (EM) algorithm is a robust method for maximum likelihood estimation of the parameters of an incompletely sampled distribution. It has been used to resolve the trial-to-trial amplitude fluctuations of postsynaptic potentials, when these are recorded in the presence of noise. Its use has however been limited by the need for different recursion equations for each set of conditions defined by the signal and noise processes. These equations are derived for the following conditions which arise in studies of synaptic transmission: non-gaussian noise process; quantal fluctuation; quantal variability. In addition, a constraint can be incorporated to accommodate simple and compound binomial models of transmitter release. Some advantages of these methods are illustrated by Monte Carlo simulations.

[1]  D. A. Sprott.  Estimating the Parameters of a Convolution by Maximum Likelihood , 1983 .

[2]  P. Harrison,et al.  Monosynaptic EPSPs in cat lumbosacral motoneurones from group Ia afferents and fibres descending in the spinal cord. , 1989, The Journal of physiology.

[3]  A. W. Kemp,et al.  Analysis with Standard Contagious Distributions , 1982 .

[4]  B. Walmsley,et al.  Nonuniform release probabilities underlie quantal synaptic transmission at a mammalian excitatory central synapse. , 1988, Journal of neurophysiology.

[5]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[6]  P. Andersen,et al.  Amplitude fluctuations in small EPSPs recorded from CA1 pyramidal cells in the guinea pig hippocampal slice , 1989, The Journal of neuroscience : the official journal of the Society for Neuroscience.

[7]  K. Pearson.  Contributions to the Mathematical Theory of Evolution , 1894 .

[8]  Brian Everitt,et al.  Maximum Likelihood Estimation of the Parameters in a Mixture of Two Univariate Normal Distributions; a Comparison of Different Algorithms , 1984 .

[9]  D. Kullmann,et al.  Reduction by general anaesthetics of group Ia excitatory postsynaptic potentials and currents in the cat spinal cord. , 1989, The Journal of physiology.

[10]  J. Jack,et al.  The components of synaptic potentials evoked in cat spinal motoneurones by impulses in single group Ia afferents. , 1981, The Journal of physiology.

[11]  Samuel Kotz,et al.  A dictionary and bibliography of discrete distributions , 1968 .

[12]  Geoffrey J. McLachlan,et al.  Mixture models : inference and applications to clustering , 1989 .

[13]  E. McLachlan.  The statistics of transmitter release at chemical synapses. , 1978, International review of physiology.

[14]  V. Hasselblad.  Estimation of parameters for a mixture of normal distributions , 1966 .

[15]  B. Katz,et al.  Spontaneous subthreshold activity at motor nerve endings , 1952, The Journal of physiology.

[16]  S. Redman.  Quantal analysis of synaptic potentials in neurons of the central nervous system. , 1990, Physiological reviews.

[17]  M. Agha,et al.  Maximum Likelihood Estimation of Mixtures of Distributions , 1984 .

[18]  L. Ling,et al.  Recovering the parameters of finite mixtures of normal distributions from a noisy record: an empirical comparison of different estimating procedures , 1983, Journal of Neuroscience Methods.