# add in M08 EIV and CRU to the yhat/yhat_sm arrays for convenience

 # EIV from M08 (1000 AD network, no Tiljander)
 nheiv <- read.table("eiv_crunhglno7_froz", nrows=996)
 nheiv_tmp<- nheiv[,2]
 years=1000:1995
 yhat[[28]] = nheiv_tmp
 names(yhat[[28]]) <- years 
 tmplo <- loess(nheiv_tmp ~years, span=.05)
 tmpy <- predict(tmplo, years, se=TRUE) 
 yhat_sm[[28]]=tmpy$fit
 names(yhat_sm[[28]]) <- years 

 # Add CRU:  (different loess smooth since timeseries is shorter) 
 # (go up to 2006)
 yhat[[29]] = instrument[,2]  
 names(yhat[[29]]) = as.character(instrument[,1])
 tmplo <- loess(instrument[,2] ~instrument[,1], span=.33)
 tmpy <- predict(tmplo, instrument[,1], se=TRUE) 
 yhat_sm[[29]]=tmpy$fit
 names(yhat_sm[[29]]) = as.character(instrument[,1]) 

 # EIV from M08 (1000 AD network, w/Tiljander)
 nheiv2 <- read.table("eiv_crunhgl_froz", nrows=996)
 nheiv2_tmp<- nheiv2[,1]
 years=1000:1995
 yhat[[30]] = nheiv2_tmp
 names(yhat[[30]]) <- years 
 tmplo <- loess(nheiv2_tmp ~years, span=.05)
 tmpy <- predict(tmplo, years, se=TRUE) 
 yhat_sm[[30]]=tmpy$fit
 names(yhat_sm[[30]]) <- years