##HEGERL ET AL NATURE 2006 #hegerl #this is archived in an Excel file ##source is "http://www.nature.com/nature/journal/v440/n7087/extref/nature04679-s2.xls" #I downloaded this and saved into txt file url<-"http://www.climateaudit.org/data/hegerl" fred<-scan(file.path(url,"nature04679-s2_orig.txt"),skip=1) #starts different at 1505 item 1017 dummy<-array(NA,dim=c(1960-1250,7)) dummy[1:(1504-1250),1:4]<-t(array(fred[1:1016],dim=c(4,1016/4))) dummy[(1505-1250):(1960-1250),]<-t(array(fred[1017:length(fred)],dim=c(7,(length(fred)-1016)/7))) hegerl<-ts(dummy[,2:7],start=1251) dimnames(hegerl)[[2]]<-c("long","plus","minus","blend","bplus","bminus") hegerl<-ts(hegerl[,1],start=1251) #1251 to 1960 ;different version at Juckes