# Make graph:

tmpin <- as.numeric(names(yy_in))
tmpout <- setdiff(as.numeric(rownames(ppc_out)), tmpin)
vv <- mean(c( max(tmpout), min(tmpin) ))

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

plot( as.numeric(rownames(ppc_out)), yhat[[1]], type="n", xlab="Year", ylab="Temperature Anomaly (C)", ylim=c(-1,0.75), col="grey", main="Reconstruction with lowpassed instrumental and proxies", cex.main=0.8)
abline(v=vv)

 # Add OLS PC1:
  i <- 2
  lines( as.numeric(rownames(ppc_out)), yhat[[i]], col=2,lwd=2)
  lines( as.numeric(rownames(ppc_out2)), yhat2[[i]], col=2,lwd=2,lty=2)

  # Add OLS PC4:
  i <- 3
  lines( as.numeric(rownames(ppc_out)), yhat[[i]], col=5,lwd=2)
  lines( as.numeric(rownames(ppc_out2)), yhat2[[i]], col=5,lwd=2,lty=2)

  # Add OLS PC10:
  i <- 4
  lines( as.numeric(rownames(ppc_out)), yhat[[i]], col=3,lwd=2)
  lines( as.numeric(rownames(ppc_out2)), yhat2[[i]], col=3,lwd=2,lty=2)

# add EIV (20 yr smooth component)
  nheiv <- read.table("eiv_crunhglno7_froz_20yrsmooth", nrows=1996) 
  lines(nheiv[1:851,1],nheiv[1:851,2], col=6,lwd=2)

# Add CRU:
  lines( lowpass.nh[,1], lowpass.nh[,2], col=1, lwd=2)

legend(1330,0.75,c("CRU","OLS PC1","OLS PC4","OLS PC10","M08 EIV"),col=c("black","red","5","green","6"),lwd=c(2,2,2,2,2),title="Calibration 1850-1980 (dashed)/1995 (solid)")

dev.off()