READ ME

SiZer_changepoints v1.0



Code archived for publication:

Abram, N.J., McGregor, H.V., Tierney, J.E., Evans, M.N., McKay, N.P., Kaufmann, D. and the PAGES 2k Consortium (2016). Early onset of industrial-era warming across the oceans and continents. Nature, doi:10.1038/nature19082.

Code v1.0 prepared by Nerilie Abram, The Australian National University, email: nerilie.abram@anu.edu.au, June 2016.  Tested with Matlab R2014a on Mac OSX 10.10 (Yosemite); also tested on matlab R2012b (64-bit, glnxa64) and Debian Linux 8.4 (Jessie), April 2016.

Please see publication for a general description of the analysis methods.  This code package was used for detecting the onset of the most recent sustained, significant trend in annually resolved paleoclimate data.



This zipfile three matlab routines:


1. calculate_changepoints_usingsizer.m: this is the main matlab file that runs the analysis routine.


2. sizer_run.m: is called in the main matlab file to run the sizer analysis.


3. sizer_changepoints.m: is a function called in the main matlab file to calculate climate change points from the matlab output



The zipfile also includes four folders:

1. input/ contains an example dataset to demonstrate how the program works. Example data set uses the primary PAGES2k continental and tropical ocean temperature reconstructions (as shown in Figure 2a and Extended Data Figure 2a in Abram et al., and the PAGES 2k Consortium, 2016).


2. output/ is where output results are written.


3. sizer/ contains the subroutines for the SiZer analysis and is based on original code by S.J. Marron with some modification by Nerilie Abram to suit the assessment of industrial-era climate trends. Edited subroutines flagged by *_edit, with changes within files marked by comments initialled with NJA. 

Please cite: 
Hannig, J., & Marron, J. S. (2006). Advanced distribution theory for SiZer. Journal of the American Statistical Association, 101(474), 484-499.

Chaudhuri, P. & Marron, J. S. SiZer for exploration of structure in curves. Journal of the American Statistical Association 94, 807-823 (1999), 

and refer to http://www.unc.edu/~marron/marron_software.html for original, unmodified sizer package.


4. working/ is where interim data is saved during the analysis, prior to final collation and saving to the output folder.



To use this code to calculate the onset of sustained, significant climate trends some modifications may be needed to suit the input data.



To run the code and see additional documentation, please open:
calculate_changepoints_usingsizer.m







