Many systems arising in biological applications are subject to periodic forcing. In these systems the forcing parameter is not only time-varying but also known to have a periodic structure. We present an approach to estimating periodic, time-varying parameters that imposes periodic structure by treating the time-varying parameter as a piecewise function with unknown coefficients. This method allows the resulting parameter estimate more flexibility in shape than prescribing a specific functional form (e.g. sinusoidal) to model its behavior, while still maintaining periodicity. We employ nonlinear filtering, more specifically, a version of the augmented ensemble Kalman filter (EnKF), to estimate the unknown coefficients comprising the piecewise approximation of the periodic, time-varying parameter. This allows for straightforward comparison of the proposed method with an EnKF-based parameter tracking algorithm, where periodicity is not guaranteed. We demonstrate the effectiveness of the proposed approach on two bio- logical examples: a synthetic example with data generated from the nonlinear FitzHugh–Nagumo system, modeling the excitability of a nerve cell, to estimate the external voltage parameter, and a case study using reported measles incidence data from three locations during the pre-vaccine era to estimate the seasonal transmission parameter. The formulation of the proposed approach also allows for simultaneous estimation of initial conditions and other static system parameters, such as the reporting probability of measles cases, which is vital for predicting under-reported incidence data.
The classic Hodgkin-Huxley model is widely used for understanding the electrophysiological dynamics of a single neuron. While applying a low-amplitude constant current to the system results in a single voltage spike, it is possible to produce multiple voltage spikes by applying time-varying currents, which may not be experimentally measurable. The aim of this work is to estimate time-varying applied currents of different deterministic forms given noisy voltage data. In particular, we utilize an augmented ensemble Kalman filter with parameter tracking to estimate four different time-varying applied current parameters and associated Hodgkin-Huxley model states, along with uncertainty bounds in each case. We test the efficiency of the parameter tracking algorithm in this setting by analyzing the effects of changing the standard deviation of the parameter drift and the frequency of data available on the resulting time-varying applied current estimates and related uncertainty.
28 January 2022 in International Journal for Numerical Methods in Biomedical Engineering