\begin{Schunk}
\begin{Sinput}
 ndf1 <- data.frame(V043250 = range(mydta1$V043250))
 mod1pred <- predict(mod1age, newdata = ndf1, interval="confidence")
 ndf1 <- cbind(ndf1, mod1pred)
 lines(fit ~ V043250, data = ndf1, lwd = 3, lty = 1)
 lines(lwr ~ V043250, data = ndf1, lwd = 2, lty = 2)
 lines(upr ~ V043250, data = ndf1, lwd = 2, lty = 2)
 legend("topleft", legend = c("predicted","95% CI"), lty = c(1,2), lwd = c(3,2), bg = "white")
\end{Sinput}
\end{Schunk}
