pdf(pdffile, width=8, height=5)
par(oma=rep(0,4), mar=c(4,4,.7,.1))


tmpin <- as.numeric(names(yhat_orig[[29]])) 
tmpout <- setdiff(x_in, tmpin)
vv <- mean(c( max(tmpout), min(tmpin) ))

plot( x_in,yhat_correct[[1]], type="n", xlab="Year", ylab="Temperature Anomaly (C)", ylim=c(-1.0,0.75), col="grey", main=title1, cex.main=0.8)

for(i in 2:length(yhat_correct)){
  lines( x_in, yhat_correct[[i]], col="grey" )
}
 # OLS PC1:
 i <- 2
 lines( x_in, yhat_correct[[i]], col=2,lwd=1)

 # OLS PC4:
# i <- 3
# lines( x_in, yhat_correct[[i]], col=5,lwd=1)

 # OLS PC10:
 i <- 4
 lines( x_in, yhat_correct[[i]], col=3,lwd=1)

 # OLS G5 PC5:
 i <- 17
 lines( x_in, yhat_correct[[i]], col=4,lwd=1)

 # EIV from M08
# lines(as.numeric(names(yhat_orig[[28]])), yhat_orig[[28]], col=6, lwd=2)

 # Add CRU:
 lines(as.numeric(names(yhat_orig[[29]])), yhat_orig[[29]], col=1, lwd=1)

 abline(v=vv)
 legend(1480,0.75,c("CRU","OLS PC1","OLS PC10","OLS G5 PC5"),col=c("black","red","green","blue"),lwd=c(2,2,2,2))

dev.off()