Surrogate Data with Correlations, Trends, and Nonstationarities

These data were contributed by Plamen Ch. Ivanov, Zhi Chen and Kun Hu, who used them in:

Hu K, Ivanov PCh, Chen Z, Carpena P, Stanley HE. Effects of trends on detrended fluctuation analysis. Phys Rev E 2001; 64:011114.
Chen Z, Ivanov PCh, Hu K, Stanley HE. Effects of nonstationarities on detrended fluctuation analysis. Phys Rev E 2002; 65:041107.

Please cite these publications when referencing this material, and also 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. Circulation 101(23):e215-e220 [Circulation Electronic Pages;]; 2000 (June 13).

The data in this collection include: (1) 6 surrogate stationary signals with different correlations; (2) 7 surrogate correlated signals with linear, sinusoidal and power-law trends; and (3) 15 surrogate correlated signals with different types of nonstationarities. Each data file contains one column of data in ASCII format. Results on correlated signals with trends are discussed in Physical Review E 64, 011114 (2001). Results on correlated signals with different types of nonstationarities are discussed in Physical Review E 65, 041107 (2002). The parameter "alpha" (see below) is an exponent measuring the degree of correlations in a signal, and Nmax is the signal length. A detailed description of these signals can be found in the original articles.

Correlations in these signals can be quantified using Detrended Fluctuation Analysis (DFA). Limitations of the DFA method are discussed in the articles cited above. In particular, the second paper notes that

... for anti-correlated signals, the scaling exponent obtained from the DFA method overestimates the true correlations at small scales. To avoid this problem, one needs first to integrate the original anti-correlated signal and then apply the DFA method. The correct scaling exponent can thus be obtained from the relation between n [the DFA box length] and F(n)/n instead of F(n) ... In order to provide a more accurate estimate of F(n), the largest box size n we use is Nmax/10, where Nmax is the total number of points in the signal.

Since these files are quite large, they are provided as gzip-compressed text.

1. Correlated stationary signals

2. Surrogate signals with trends

    2a) Signals with linear trends

    2b) Signals with sinusoidal trends
    2c) Signals with power-law trends

3. Surrogate nonstationary signals

    3a) Signals with cutout segments (discontinuities)

    3b) Signals with spikes     3c) Signals with different local standard deviation     3d) Signals with different local correlations

Address for correspondence:
Plamen Ch. Ivanov, Ph.D.
Room 247, Dept. of Physics
Boston Univeristy
590 Commonwealth Avenue
Boston, MA 02215, USA

Questions and Comments

If you would like help understanding, using, or downloading content, please see our Frequently Asked Questions.

If you have any comments, feedback, or particular questions regarding this page, please send them to the webmaster.

Comments and issues can also be raised on PhysioNet's GitHub page.

Updated Tuesday, 4 October 2016 at 16:25 BRT

PhysioNet is supported by the National Institute of General Medical Sciences (NIGMS) and the National Institute of Biomedical Imaging and Bioengineering (NIBIB) under NIH grant number 2R01GM104987-09.