Inferential framework for nonstationary dynamics. II. Application to a model of physiological signaling.

The problem of how to reconstruct the parameters of a stochastic nonlinear dynamical system when they are time-varying is considered in the context of online decoding of physiological information from neuron signaling activity. To model the spiking of neurons, a set of FitzHugh-Nagumo (FHN) oscillators is used. It is assumed that only a fast dynamical variable can be detected for each neuron, and that the monitored signals are mixed by an unknown measurement matrix. The Bayesian framework introduced in paper I (immediately preceding this paper) is applied both for reconstruction of the model parameters and elements of the measurement matrix, and for inference of the time-varying parameters in the nonstationary system. It is shown that the proposed approach is able to reconstruct unmeasured (hidden) slow variables of the FHN oscillators, to learn to model each individual neuron, and to track continuous, random, and stepwise variations of the control parameter for each neuron in real time.

[1]  A. McCulloch,et al.  A collocation-Galerkin finite element model of cardiac action potential propagation , 1994, IEEE Transactions on Biomedical Engineering.

[2]  Riccardo Mannella,et al.  Integration Of Stochastic Differential Equations On A Computer , 2002 .

[3]  A. Duggento,et al.  Inferential framework for nonstationary dynamics. I. Theory. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[4]  S Akselrod,et al.  Bifurcation in a simple model of the cardiovascular system. , 2000, Methods of information in medicine.

[5]  R. FitzHugh.  Impulses and Physiological States in Theoretical Models of Nerve Membrane. , 1961, Biophysical journal.

[6]  E. Hill.  Journal of Theoretical Biology , 1961, Nature.

[7]  A. Winfree.  The geometry of biological time , 1991 .

[8]  Ki H. Chon,et al.  Nonlinear autoregressive and nonlinear autoregressive moving average model parameter estimation by minimizing hypersurface distance , 2003, IEEE Trans. Signal Process..

[9]  S. Abboud,et al.  Simulation of cardiac activity and the ECG using a heart model with a reaction-diffusion action potential. , 1996, Medical engineering & physics.

[10]  S. Grossberg,et al.  Cortical networks for control of voluntary arm movements under variable force conditions. , 1998, Cerebral cortex.

[11]  James P. Keener,et al.  Mathematical physiology , 1998 .

[12]  Anthony Papagiannis.  Intern , 2010, BMJ : British Medical Journal.

[13]  S. Yoshizawa,et al.  An Active Pulse Transmission Line Simulating Nerve Axon , 1962, Proceedings of the IRE.

[14]  A V Panfilov,et al.  Modeling of heart excitation patterns caused by a local inhomogeneity. , 1996, Journal of theoretical biology.

[15]  M Congedo,et al.  A review of classification algorithms for EEG-based brain–computer interfaces , 2007, Journal of neural engineering.

[16]  A Stefanovska,et al.  Nonlinear statistical modeling and model discovery for cardiorespiratory data. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[17]  D. Wilkin,et al.  Neuron , 2001, Brain Research.

[18]  Andrew B. Schwartz,et al.  Brain-Controlled Interfaces: Movement Restoration with Neural Prosthetics , 2006, Neuron.

[19]  R J Cohen,et al.  A forward model-based validation of cardiovascular system identification. , 2001, American journal of physiology. Heart and circulatory physiology.

[20]  E. N. Best,et al.  Null space in the Hodgkin-Huxley Equations. A critical test. , 1979, Biophysical journal.