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’) , or use the new code page

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:

fig1.1.png

About these ads

3 Responses to “Hockeystick for Matlab”

  1. Spence_UK Says:

    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 ;)

  2. uced Says:

    :) 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.

  3. uc00 Says:

    The code is now moved to CA, http://climateaudit.org/2008/07/01/code-the-hockey-stick/ . Will close the comments here.

Comments are closed.


Follow

Get every new post delivered to your Inbox.

%d bloggers like this: