Software Open Access
Multiscale Poincare Plots
Published: April 12, 2016. Version: 1.0.0
Henriques TS, Mariani S, Burykin A, Rodrigues F, Silva TF, Goldberger AL, Multiscale Poincaré plots for visualizing the structure of heartbeat time series, BMC Medical Informatics and Decision Making. 2016 Feb 9;16(1):1.Please include the standard citation for PhysioNet:
Goldberger AL, Amaral LAN, Glass L, Hausdorff JM, Ivanov PCh, Mark RG, Mietus JE, Moody GB, Peng C-K, Stanley HE. PhysioBank, PhysioToolkit, and PhysioNet: Components of a New Research Resource for Complex Physiologic Signals (2003). Circulation. 101(23):e215-e220.
Poincaré delay maps are widely used in the analysis of cardiac interbeat interval (RR) dynamics. To facilitate visualization of the structure of these complex time series, we introduced a methodology called multiscale Poincaré (MSP) plots.
Briefly, the function
MsPplots.m (version 1.0) starts with the original time series and employs a coarse-graining procedure to create a set of time series, each of which represents the system’s dynamics on a different time scale. Next, Poincaré (delay) plots are constructed (using the function
dscatter.m by Robert Henson) for the original and the coarse-grained time series. Finally, as an optional adjunct, color can be added to each point to represent its normalized frequency.
The functions have been tested on Matlab R2010b, R2014b and R2015a on April 11 2016.
Inputs to the
MsPplots.m function are:
- The time series vector with one column (mandatory)
- The number of time scales (integer; default: scales = 12)
- Poincaré plots from scale 1 to scales by scstep. (integer; default: scstep = 1)
- The number of rows in the plot montage (integer; default: nrwo = 3)
- The number of columns in the plot montage (integer; default: ncol = 4)
- A option that allows to save the final figure (boolean; default: 1)
- The name of the figure to be saved (default: Figure1.png)
Also provided, are two example scripts that employ deidentified RR time series from two sources: one from the database for Chaos Is the Normal Heart Rate Chaotic? and the other from the Normal Sinus Rhythm RR Interval Database, both available at PhysioNet.
The first example script requires the text file
n2-filtered included in this directory, and the second one requires the WFDB Matlab toolbox to access data from Physiobank.
The helper function
dscatter.m is written by a separate author and is covered by the BSD License under Mathsworks. The rest of the software in this directory is licensed under the GNU Public License 3.
The Multiscale Poincare Plot (MSP) package was developed by Teresa S. Henriques, Sara Mariani and Anton Burykin. A detailed description can be found in the citation.
Anyone can access the files, as long as they conform to the terms of the specified license.
License (for files):
Open Data Commons Attribution License v1.0
Total uncompressed size: 0 B.