Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2012:2012:578785.
doi: 10.1155/2012/578785. Epub 2012 Mar 6.

An efficient time-varying filter for detrending and bandwidth limiting the heart rate variability tachogram without resampling: MATLAB open-source code and Internet web-based implementation

A Eleuteri  1 , A C Fisher, D Groves, C J Dewhurst
Affiliations

An efficient time-varying filter for detrending and bandwidth limiting the heart rate variability tachogram without resampling: MATLAB open-source code and Internet web-based implementation

A Eleuteri et al. Comput Math Methods Med. 2012.

Abstract

The heart rate variability (HRV) signal derived from the ECG is a beat-to-beat record of RR intervals and is, as a time series, irregularly sampled. It is common engineering practice to resample this record, typically at 4 Hz, onto a regular time axis for analysis in advance of time domain filtering and spectral analysis based on the DFT. However, it is recognised that resampling introduces noise and frequency bias. The present work describes the implementation of a time-varying filter using a smoothing priors approach based on a Gaussian process model, which does not require data to be regular in time. Its output is directly compatible with the Lomb-Scargle algorithm for power density estimation. A web-based demonstration is available over the Internet for exemplar data. The MATLAB (MathWorks Inc.) code can be downloaded as open source.

PubMed Disclaimer

Figures

Figure 1
Figure 1
Weight functions (viz. D 2 operator).
Figure 2
Figure 2
Bode plot of theoretical transfer function of equivalent kernel filter.
Figure 3
Figure 3
Bode plot of discrete transfer function of equivalent kernel filter.
Figure 4
Figure 4
Synthetic and clinical HRV records band-pass filtered by sequential application of SGP: raw data vt 0 “smoothed” to give vt 1; vt 2 = vt 0 − vt 1 (not shown); vt 2 “smoothed” to give vt 3. Lomb Scargle Periodograms (LSPs) are for vt 0, vt 2, and vt 3.
Figure 5
Figure 5
Speed-up of SGP over SPA with increasing data set size.
See this image and copyright information in PMC

References

    1. Moody GB. Spectral analysis of heart rate without resampling. In: Proceedings of the IEEE Conference on Computers in Cardiology; September 1993; London, UK. pp. 715–718.
    1. Malik M, Camm AJ, Bigger JT, Jr., et al. Heart rate variability. Standards of measurement, physiological interpretation, and clinical use. European Heart Journal. 1996;17(3):354–381. - PubMed
    1. Clifford GD. ECG statistics, noise, artifacts, and missing data. In: Clifford GD, Azuaje F, McSharry PE, editors. Advanced Methods for ECG Analysis. Boston, Mass, USA: Artech-House; 2006. pp. 55–93.
    1. Litvack DA, Oberlander TF, Carney LH, Saul JP. Time and frequency domain methods for heart rate variability analysis: a methodological comparison. Psychophysiology. 1995;32(5):492–504. - PubMed
    1. Mitov IP. A method for assessment and processing of biomedical signals containing trend and periodic components. Medical Engineering and Physics. 1998;20(9):660–668. - PubMed

LinkOut - more resources