##BRIFFA JGR 2001 SPAGHETTI DIAGRAM #this is scaled differently than original sources for non-Briffa 2001 versions url<-"ftp://ftp.ncdc.noaa.gov/pub/data/paleo/contributions_by_author/briffa1998/briffa2001jgr3.txt" #this has 20N version from Rev Geophys 1999 Apr-Sept #1961-990 mean Briffa<-read.table(url,skip=24,fill=TRUE) temp<-(Briffa< (-999)) Briffa[temp]<-NA dimnames(Briffa)[[2]]<-c("Year","Jones98","MBH99","Briffa01","Briffa00","Overpeck97","Crowley00","CRU99") Briffa<-ts(Briffa,start=Briffa[1,1]) #this gives series 1, 12,3, NA,NA,4,NA. #starts 1000; 1000; 1402; 1000; 1600;1000;1871 #ends 1991; 1980; 1960; 1993; 1990; 1987; 1997 #array from 1000 to 1999 #Briffa00 gores to 1993; Briffa 01 truncated in 1960 (I've emulated extension) briffa01<-ts(Briffa[!is.na(Briffa[,4]),4],start=1402)