title: Bayesian estimation of phase response curves.
abstract: Phase response curve (PRC) of an oscillatory neuron describes the response of the neuron to external perturbation. The PRC is useful to predict synchronized dynamics of neurons; hence, its measurement from experimental data attracts increasing interest in neural science. This paper introduces a Bayesian method for estimating PRCs from data, which allows for the correlation of errors in explanatory and response variables of the PRC. The method is implemented with a replica exchange Monte Carlo technique; this avoids local minima and enables efficient calculation of posterior averages. A test with artificial data generated by the noisy Morris-Lecar equation shows that the proposed method outperforms conventional regression that ignores errors in the explanatory variable. Experimental data from the pyramidal cells in the rat motor cortex is also analyzed with the method; a case is found where the result with the proposed method is considerably different from that obtained by conventional regression.
authors.affiliation: Department of Statistical Science, School of Multidisciplinary Sciences, The Graduate University for Advanced Studies, 10-3 Midori-Machi, Tachikawa, Tokyo, Japan. nakae@ism.ac.jp
authors.forename: Ken
authors.initials: K
authors.lastname: Nakae
authors.forename: Yukito
authors.initials: Y
authors.lastname: Iba
authors.forename: Yasuhiro
authors.initials: Y
authors.lastname: Tsubo
authors.forename: Tomoki
authors.initials: T
authors.lastname: Fukai
authors.forename: Toshio
authors.initials: T
authors.lastname: Aoyagi
mesh_terms: D000200:Action Potentials / Q000502:physiology; D000818:Animals; D001499:Bayes Theorem; D001683:Biological Clocks / Q000502:physiology; D003198:Computer Simulation; D009010:Monte Carlo Method; D009044:Motor Cortex / Q000502:physiology; D009415:Nerve Net / Q000502:physiology*; D016571:Neural Networks, Computer; D009474:Neurons / Q000502:physiology*; D051381:Rats; D011930:Reaction Time / Q000502:physiology*; D012815:Signal Processing, Computer-Assisted / Q000295:instrumentation