#DArrigo 2006 url="ftp://ftp.ncdc.noaa.gov/pub/data/paleo/treering/reconstructions/n_hem_temp/nhtemp-darrigo2006.txt" #"YearAD NHLandObs STD recon RCS recon" #[1229] " 1856 -0.55 -0.55 -0.69" #[1230] " 1857 -0.38 -0.39 -0.48" #[1231] " 1858 -0.3 -0.19 -0.3" darrigo=read.fwf(url,skip=85,fill=TRUE,widths=c(6,14,15,13)) names(darrigo)=c("year","instr","STD","RCS") temp=!is.na(darrigo[,"RCS"]) darrigo=ts(darrigo[temp,"RCS"],start=min(darrigo[temp,1]) ) #darrigo=ts(darrigo[!is.na(darrigo)],start=tsp(darrigo)[1])