A man is looking for a bar in Antarctica when he comes across a sign outside a pub :-

“A pint

A pie

and a Friendly Word”

This seems fine so in he goes.

Customer: “A pint please”

Customer: “I’ll also have a pie”

Silence

Customer: “Landlord what about the Friendly Word”

Landlord: “Don’t eat the pie”

]]>(quote from: “Real Programmers Don’t Do Pascal”)

]]>]]>id=names(Model[[1]])

test=sapply(id,function(x) strsplit(x,”/”)[[1]])

Info.Model=data.frame(model=test[1,],run=test[2,])

row.names(Info.Model)=1:49h=function(A) {

temp=(time(A)>=1979)& (time(A)<2000)

fm=lm(A[temp,"TROP"]~A[temp,"RTIME"])

h=fm$coef[2];

h}Info.Model$trend_T2=sapply(Model[[1]],h)

Info.Model$trend_T2LT=sapply(Model[["T2LT"]],h)Info=data.frame(model=unique(Info.Model$model))

Info$trend_T2=tapply(Info.Model$trend_T2,Info.Model$model,mean)

Info$trend_T2LT=tapply(Info.Model$trend_T2LT,Info.Model$model,mean)SANTER TABLE 1 #1. Santer 2008 Table 1

loc=”http://data.climateaudit.org/data/models/santer_2008_table1.dat”

santer=read.table(loc,skip=1)

names(santer)=c(“item”,”layer”,”trend”,”se”,”sd”,”r1″,”neff”)

row.names(santer)=paste(santer[,1],santer[,2],sep=”_”)

santer=santer[,3:ncol(santer)]rbind(10*apply(Info[,2:3],2,sd),

santer[c(16,12),1])

# trend_T2 trend_T2LT

#[1,] 0.09778996 0.09208094

#[2,] 0.09800000 0.09200000rbind(10*apply(Info[,2:3],2,mean),

santer[c(15,11),1])

# trend_T2 trend_T2LT

#[1,] 0.1992421 0.2147754

#[2,] 0.1990000 0.2150000

The groups are given (using the notation and group order in the paper in Figure 3A) by the script:

modgrp = c(12,1,1,1,1,1,13,14,15,15,15,2,2,2,3,3,3,17,17,4,4,4,4,4,5,

5,5,5,5,10,11,16,16,16,18,19,7,6,6,6,8,8,8,8,8,9,9,9,9)labs =c(“CCSM3.0″,”GFDL-CM2.0″,”GFDL-CM2.1″,”GISS-EH”,”GISS-ER”,

“MIROC3.2(medres)”,”MIROC3.2(hires)”,”MRI-CGCM2.3.2″,”PCM”,

“UKMO-HadCM3″,”UKMO-HadGEM1″,”CCCma-CGCM3.1(T47)”,”CNRM-CM3″,

“CSIRO-Mk3.0″,”ECHAM5/MPI-OM”,”FGOALS-g1.0″,”GISS-AOM”,”INM-CM3.0″,”IPSL-CM4″)modgrp = factor(modgrp,labels=labs)

My scripts for extracting 1979 to 1999 and getting the regression coefficients were a bit clumsy so I won’t post them.

]]>