Author Archives: uc00

Code: The Trick

load MBHsmooths1.txt % http://www.climateaudit.info/wp-content/uploads/2009/11/mbhsmooths1.txt % MBHsmooths1 = [ Year AnnualRecon Instrumental MB98SmoothTrick MBH98SmoothNoTrick MBH99SmoothTrick MBH99SmoothNoTrick ] [B98,A98]=butter(10,2/50); % ’50 year lowpass’ in98=MBHsmooths1(401:981,2); % Annual Recon MBH98, %source: % http://www.nature.com/nature/journal/v430/n6995/extref/FigureData/nhmean.txt , % col 2 Trick=MBHsmooths1(982:996,3); % pad with instrumental (1981..1995) .. % source: % http://www.nature.com/nature/journal/v430/n6995/extref/FigureData/nhmean.txt % col 3 out98t=flipud(filter(B98,A98,flipud([filter(B98,A98,[in98;Trick;zeros(100,1)])]))); %..then smooth % Note that filter.m initializes […]

The Trick Timeline

Date: 16 Nov 1999, Phil I’ve just completed Mike’s Nature trick of adding in the real temps to each series for the last 20 years (ie from 1981 onwards) amd from 1961 for Keith’s to hide the decline. Date: 22 Dec 2004, mike   No researchers in this field have ever, to our knowledge, “grafted […]

Mann 2008 – Replication II

Let’s continue Mann 2008 – Replication with EIV. To run regreclow.m and regrechigh.m you’ll need files climate, eofnumb, and proxyN, where N runs from 1 to 19. I’ve run prepinputforrecon.m with required folder structure (C:\holocene\s1\zuz10\work1\temann\ etc.) in my computer. After that I did run regrechigh.m in folder C:\holocene\s1\zuz10\work1\temann\zzrecon1209\nhnhscrihad\highf and regreclow.m in folder C:\holocene\s1\zuz10\work1\temann\zzrecon1209\nhnhscrihad\lowf The regem […]

UC on Mann Smoothing

UC writes in that there’s another Mannian problem: And GRL08 ( http://www.meteo.psu.edu/~mann/shared/articles/MannGRL08.pdf ) is the prologue. Tried this new smoother ( http://www.meteo.psu.edu/~mann/smoothing08/lowpassadaptive.m ) with HadCRUT monthly and f=0.0104, [smoothedbest,w0,w1,w2,msebest] = lowpassadaptive(HadM, 0.0104); got this figure, But the weights are [w0 w1 w2] ans = 0.4100 0 0.5500 Surprisingly, this does not agree with the text, […]

Code: The Hockey Stick

%%%% Matlab version of Hockey Stick % %Mann, M. E., R. S. Bradley, and M. K. Hughes (1998), Global-scale %temperature patterns and climate forcing over the past six centuries, %Nature, 392, 779– 787. % %Mann, M. E., R. S. Bradley, and M. K. Hughes (1999), Northern %Hemisphere temperatures during the past millennium: Inferences, uncertainties, %and […]

UC on CCE

I then compared verification statistics for the different reconstructions as shown below. OLS yielded much the “best” fit in the calibration period, but the worst fit in the verification period. If OLS is equivalent to ICE, it actually finds the best fit (minimizes calibration residuals), and in proxy-temperature case makes the most obvious overfit. Let […]

Multivariate Calibration

In calibration problem we have accurately known data values (X) and a responses to those values (Y). Responses are scaled and contaminated by noise (E), but easier to obtain. Given the calibration data (X,Y), we want to estimate new data values (X’) when we observe response Y’. Using Brown’s (Brown 1982) notation, we have a […]

Follow

Get every new post delivered to your Inbox.

Join 3,262 other followers