Model of heart rate variability including the effect of sleep stages

This model of heart rate variability is based on the two main assumptions: the dominating component of heart rate variability is random and that sleep architecture has a significant effect on HRV. We retained the stochastic properties of RR interval time series of the model by Kantelhardt et al.—the starting point in creating our model [1]. Human sleep architecture (the hypnogram) was simulated using our own method based on 30 real EEG recordings (database of the Faculty of Psychiatry, Warsaw University of Medicine, Warsaw, Poland) considering two properties of sleep architecture: the probabilities of transitions between all possible pairs of sleep stages and the probability distributions of sleep stage durations. The matrices of transitions probabilities with the percentage of the occurrence of the sleep stages in the complete time series are shown in Table 1.

We considered several combinations of additional components introducing into the model, e.g., the strength of the effect of breathing on the heart rate simulated using the inverse Fourier transform of a series obtained from a pre-specified power spectrum (based on McSharry et. al model [2]), the variability of the model parameters, patterns of the shortening of the RR interval lengths during exercise periods during sleep, and smoothing out the transitions between sleep stages. That whole set of elements, often very easy to implement, provides to generate synthetic signals which are indistinguishable from the real data using either standard or nonlinear (including Multiscale Multifractal Analysis method [3]) heart rate variability measurements. The model was tested by using 34 RR time series from Institute of Cardiology, Warsaw, Poland each extracted from the nighttime part of Holter recordings [4]. Examples of real and generated by our model RR time series are shown in Figure 1.

Figure 1: RR time series obtained for the real data (a) and from our model (b). The characteristic fuzzy structure of the upper envelope, visible in the signals from our model and in real data, reflects human sleep architecture.

One of the strengths of this model is the easy to interpret parameters that define the sleep architecture (i.e., the probability of transitions between the sleep phases). Variation of their values allows to simulate cases of sleep disorders (e.g., sleep apnea and insomnia) and even impossible or rare structures of sleep architecture, e.g., sleep without REM or deep stage, sleep with many transitions to wake stage, etc., as well.

The Matlab code, NightTimeHRVmodel.m, in this directory contains function generates RR time interval series. Function name: NightTimeHRVmodel .

Usage

Usage example:

NightTimeHRVmodel(n,nos,ifplot)

where n is number of samples of signal (default 25000), nos is number of signals to generate (default 1) and ifplot parameter enable/disable plotting of generated signals (signals are plotted when ifplot =1). The output as .txt files included RR time series and (optionaly) .png plots are storage at a folder “. /RR” placed in the same directory as Matlab script.

Example data

The .zip archive HRVmodelExample.zip contains 20 RR time series generated by the model. Each signal is 25,000 samples length. The values of RR time intervals (in ms unit) are stored in .txt files . Additionally, .png picture with tachogram and hypnogram is joined to each recording . The files are created according to the default output option of the Matlab script NightTimeHRVmodel.m

Acknowledgements

The authors would like to thank their colleagues from the Cardiovascular Physics Group at the Faculty of Physics of Warsaw University of Technology for support. Special thanks go to Gari Clifford from Emory University and Georgia Institute of Technology for sharing uploaded version of HRV model by McSharry and Clifford. Thanks also to the following physicians for providing the medical data: Rafał Baranowski, Aleksander Prejbisz, Marek Kabat, and Andrzej Januszewicz from The Cardinal Stefan Wyszyński Institute of Cardiology, Warsaw, Poland, and Anna Piotrowska and Waldemar Szelenberg from the Faculty of Psychiatry, Warsaw University of Medicine, Poland.

References:

 Name                    Last modified      Size  Description
 Parent Directory                             -   
 local.css               04-Mar-2016 23:22  3.1K  C# source file
 NightTimeHRVmodel.m     04-Mar-2016 23:22   16K  
 HRVmodelExample.zip     04-Mar-2016 23:22  1.9M