Hockeystick for Matlab
Here’s the version 1.1: hockeystick1.txt
UPD Jan 2010: change
urlwrite(’http://www.climateaudit.org/data/mbh99/proxy.txt’,'proxy.txt’)
to
urlwrite(’http://www.climateaudit.info/data/mbh99/proxy.txt’,'proxy.txt’)
Some notes:
- Download to empty folder and rename to hockeystick.m
- Program downloads necessary data from the web (once), uses urlwrite.m (newish Matlab needed)
- It’s a script
- Shows what PC1_fixed does
- Only one file is downloaded from CA (AD1000 proxies), sorry RC, but I don’t know where to find morc014 elsewhere..
- Pl. tell me if it works or not, uc_edit at yahoo.com !
Updated to Ver 1.1, added cooling trends:
Spence_UK:
Worked okay for me, MATLAB 6.5 at home. I could try it on later versions at work but I suspect the firewall may disagree with me.
This is what it produced:
http://i27.tinypic.com/309pt02.gif
It also claimed a verification R2 of 0.998 for the 1400AD step …
… just kidding
3 July 2008, 1:41 pmuced:
So, now you know that Hockey Stick is 600 lines (+ code for U and V vectors + Proxy PCs) of data formatting, plus incorrectly implemented Classical Calibration (lines 595-606) with incorrectly computed uncertainties. And astronomical cooling is due to fixed PC1.
4 July 2008, 3:27 pm