##REPLICATION LOEHLE Data=array(NA,dim=c(1995,18)); proxy_sm=array(NA,dim=c(1995,18)) for (i in 1:18) { test=proxy[[i]] test=test[!is.na(test$t),] f=approxfun(test[,1],test[,2]) Data[,i]=f(1:1995) proxy_sm[,i]=filter(Data[,i],rep(1/30,30) ) } proxy_sm=scale(proxy_sm,scale=FALSE); dimnames(proxy_sm)[[2]]=id annual=ts(apply(proxy_sm,1,mean,na.rm=TRUE),start=1) count=ts(apply(!is.na(Data),1,sum,na.rm=TRUE),start=1) ##LOEHLE VERSION url<-"http://www.ncasi.org/programs/areas/climate/LoehleE&E2007.csv" test=read.csv(url,skip=2) recon=ts(test[,2],start=test[1,1]) #17 to 1980 mean(recon[(1902:1980)-16]) # 0.01586732 ##PLOT par(mar=c(3,3,2,1)) ts.plot(annual) lines(recon,col="red") ##12 CRU NH #potential target 1961-1990 url<-"http://www.cru.uea.ac.uk/cru/data/temperature/hadcrut3nh.txt" #1850 2007 h<-scan(url) start0<-h[1] N<-length(h) h<-array(h,dim=c(27,N/27)) h<-t(h) CRU<-ts(h[,14],start=start0,end=2006) # 1850 2006 delta=mean(CRU[(1902:1980)-1849]);delta # -0.1592278