Data and MATLAB Code for Reproducing RPs and RQA of Westerhold et al, Science, 369, 2020

The data file `CENOGRID_Loess_20.txt` contains the astronomically tuned deep-sea benthic foraminifer carbon (δ¹³C) and oxygen (δ¹⁸O) isotope reference records uniformly covering the entire Cenozoic. The first column is the tuned age in Ma, the second column the δ¹³C, and the third column the δ¹⁸O record.

The original calculations were performed using the CRP Toolbox for MATLAB. In order to avoid installing the toolbox and for better performance, the functions for calculating RP and RQA were here reimplemented, providing identical result.

To reproduce the RPs in Fig. 2, use the script `perform_rp.m`, for reproducing the determinism values and upper confidence bounds, use the script `perform_rqa.m`.

Published in
T. Westerhold, N. Marwan, A. J. Drury, D. Liebrand, C. Agnini, E. Anagnostou, J. S. K. Barnet, S. M. Bohaty, D. De Vleeschouwer, F. Florindo, T. Frederichs, D. A. Hodell, A. E. Holbourn, D. Kroon, V. Lauretano, K. Littler, L. J. Lourens, M. Lyle, H. Pälike, U. Röhl, J. Tian, R. H. Wilkens, P. A. Wilson, J. C. Zachos: An astronomically dated record of Earth's climate and its predictability over the last 66 million years, Science, 369(6509), 1383–1387 (2020). DOI:10.1126/science.aba6853