Unbiased Estimator of Shape Parameter for Spiking Irregularities under Changing Environments

We considered a gamma distribution of interspike intervals as a statistical model for neuronal spike generation. The model parameters consist of a time-dependent firing rate and a shape parameter that characterizes spiking irregularities of individual neurons. Because the environment changes with time, observed data are generated from the time-dependent firing rate, which is an unknown function. A statistical model with an unknown function is called a semiparametric model, which is one of the unsolved problem in statistics and is generally very difficult to solve. We used a novel method of estimating functions in information geometry to estimate the shape parameter without estimating the unknown function. We analytically obtained an optimal estimating function for the shape parameter independent of the functional form of the firing rate. This estimation is efficient without Fisher information loss and better than maximum likelihood estimation.

[1]  S. Amari,et al.  Estimation in the Presence of Infinitely many Nuisance Parameters--Geometry of Estimating Functions , 1988 .

[2]  Shigeru Shinomoto,et al.  Differences in Spiking Patterns Among Cortical Neurons , 2003, Neural Computation.

[3]  Shun-ichi Amari,et al.  Methods of information geometry , 2000 .

[4]  P. Bickel.  Efficient and Adaptive Estimation for Semiparametric Models , 1993 .

[5]  Henry C. Tuckwell,et al.  Introduction to theoretical neurobiology , 1988 .

[6]  S. Shinomoto,et al.  A measure of local variation of inter-spike intervals. , 2005, Bio Systems.

[7]  H. Tamura,et al.  Regional and laminar differences in in vivo firing patterns of primate cortical neurons. , 2005, Journal of neurophysiology.

[8]  S. Amari,et al.  Information geometry of estimating functions in semi-parametric statistical models , 1997 .

[9]  R N Lemon,et al.  Precise spatiotemporal repeating patterns in monkey primary and supplementary motor areas occur at chance levels. , 2000, Journal of neurophysiology.

[10]  James L. Powell,et al.  Estimation of semiparametric models , 1994 .

[11]  David R. Cox,et al.  The statistical analysis of series of events , 1966 .

[12]  William R. Softky,et al.  Comparison of discharge variability in vitro and in vivo in cat visual cortex neurons. , 1996, Journal of neurophysiology.

[13]  P. Bickel,et al.  Achieving Information Bounds in Non and Semiparametric Models , 1990 .

[14]  J. Neyman,et al.  Consistent Estimates Based on Partially Consistent Observations , 1948 .

[15]  Shun-ichi Amari,et al.  Information geometry on hierarchy of probability distributions , 2001, IEEE Trans. Inf. Theory.

[16]  Kazushi Ikeda.  Information Geometry of Interspike Intervals in Spiking Neurons , 2005, Neural Computation.

[17]  K. Do,et al.  Efficient and Adaptive Estimation for Semiparametric Models. , 1994 .

[18]  V. P. Godambe.  An Optimum Property of Regular Maximum Likelihood Estimation , 1960 .

[19]  Yutaka Sakai,et al.  Temporally correlated inputs to leaky integrate-and-fire models can reproduce spiking statistics of cortical neurons , 1999, Neural Networks.