Model-Based Human Circadian Phase Estimation Using a Particle Filter

We present a method for tracking an individual's circadian phase that integrates dynamic models of circadian physiology with physiological measurements in a Bayesian statistical framework. A model of the circadian pacemaker's response to light exposure is transformed into a nonlinear state-space model with a circadian phase state. The probability distribution of the circadian phase is estimated by a particle filter that predicts changes over time based on the model, and performs updates with information gained from physiological measurements. Simulations demonstrate how probability distributions allow flexible initialization of model states and enable statistical quantification of entrainment and divergence properties of the circadian pacemaker. The combined use of sleep-wake scheduling data and physiological measurements is demonstrated in a case study highlighting advantages for addressing the challenge of noninvasive ambulatory monitoring of circadian physiology.

[1]  D. Skene,et al.  An action spectrum for melatonin suppression: evidence for a novel non‐rod, non‐cone photoreceptor system in humans , 2001, The Journal of physiology.

[2]  A. Caufriez,et al.  Immediate effects of an 8h advance shift of the rest-activity cycle on 24h profiles of cortisol , 2002 .

[3]  Branko Ristic,et al.  Beyond the Kalman Filter: Particle Filters for Tracking Applications , 2004 .

[4]  Nando de Freitas,et al.  Sequential Monte Carlo Methods in Practice , 2001, Statistics for Engineering and Information Science.

[5]  D. Boivin,et al.  Phase-Dependent Effect of Room Light Exposure in a 5-h Advance of the Sleep-Wake Cycle: Implications for Jet Lag , 2002, Journal of biological rhythms.

[6]  T. Reilly,et al.  A COMPARISON OF SOME DIFFERENT METHODS FOR PURIFYING CORE TEMPERATURE DATA FROM HUMANS , 2000, Chronobiology international.

[7]  Jürgen Kurths,et al.  The Unscented Kalman Filter, a Powerful Tool for Data Analysis , 2004, Int. J. Bifurc. Chaos.

[8]  Francis Lévi,et al.  From circadian rhythms to cancer chronotherapeutics , 2002, Chronobiology international.

[9]  Markus Meister,et al.  Individual neurons dissociated from rat suprachiasmatic nucleus express independently phased circadian firing rhythms , 1995, Neuron.

[10]  Athanasios Papoulis,et al.  Probability, Random Variables and Stochastic Processes , 1965 .

[11]  C A Czeisler,et al.  Linear Demasking Techniques Are Unreliable for Estimating the Circadian Phase of Ambulatory Temperature Data , 1999, Journal of biological rhythms.

[12]  S. Folkard,et al.  Shift work, safety and productivity. , 2003, Occupational medicine.

[13]  A. Wirz-Justice,et al.  Circadian Rhythms and Depression: Clinical and Experimental Models , 1997 .

[14]  Charles A. Czeisler,et al.  The Timing of the Human Circadian Clock Is Accurately Represented by the Core Body Temperature Rhythm following Phase Shifts to a Three-Cycle Light Stimulus Near the Critical Zone , 2000, Journal of biological rhythms.

[15]  Richard E. Kronauer,et al.  Quantifying Human Circadian Pacemaker Response to Brief, Extended, and Repeated Light Stimuli over the Phototopic Range , 1999, Journal of biological rhythms.

[16]  Daniel B. Forger,et al.  Revised Limit Cycle Oscillator Model of Human Circadian Pacemaker , 1999, Journal of biological rhythms.

[17]  Hans P. A. Van Dongen,et al.  Chapter 35 – Circadian Rhythms in Sleepiness, Alertness, and Performance , 2005 .

[18]  D. Skene,et al.  Physiology and pharmacology of melatonin in relation to biological rhythms , 2009, Pharmacological reports : PR.

[19]  B. Anderson,et al.  Optimal Filtering , 1979, IEEE Transactions on Systems, Man, and Cybernetics.

[20]  Richard E Kronauer,et al.  Comparisons of the Variability of Three Markers of the Human Circadian Pacemaker , 2002, Journal of biological rhythms.

[21]  Richard E. Kronauer,et al.  Erratum: Quantifying human circadian pacemaker response to brief, extended, and repeated light stimuli over the photopic range (Journal of Biological Rhythms) , 2000 .

[22]  Neil J. Gordon,et al.  Editors: Sequential Monte Carlo Methods in Practice , 2001 .

[23]  D. Boivin,et al.  Circadian clock genes oscillate in human peripheral blood mononuclear cells. , 2003, Blood.

[24]  E N Brown,et al.  The Statistical Analysis of Circadian Phase and Amplitude in Constant-Routine Core-Temperature Data , 1992, Journal of biological rhythms.

[25]  H. Saunders,et al.  Probability, Random Variables and Stochastic Processes (2nd Edition) , 1989 .

[26]  Neil J. Gordon,et al.  A tutorial on particle filters for online nonlinear/non-Gaussian Bayesian tracking , 2002, IEEE Trans. Signal Process..

[27]  M. Huzmezan,et al.  Modifying the human circadian pacemaker using model based predictive control , 2003, Proceedings of the 2003 American Control Conference, 2003..

[28]  R. Kronauer,et al.  Stability, precision, and near-24-hour period of the human circadian pacemaker. , 1999, Science.

[29]  Jiří Zelinka,et al.  Kernel Density Estimation Toolbox for Matlab , 2011 .

[30]  R. Kronauer,et al.  Bright light induction of strong (type 0) resetting of the human circadian pacemaker. , 1989, Science.

[31]  D. Boivin,et al.  Working on atypical schedules. , 2007, Sleep medicine.

[32]  Simon J. Godsill,et al.  On sequential Monte Carlo sampling methods for Bayesian filtering , 2000, Stat. Comput..

[33]  G. Brainard,et al.  Action Spectrum for Melatonin Regulation in Humans: Evidence for a Novel Circadian Photoreceptor , 2001, The Journal of Neuroscience.