\begin{Schunk}
\begin{Sinput}
 ndf <- data.frame(V043250 = range(mydta1$V043250))
 ndf$fit <- predict(mod1age, newdata = ndf )
 head(ndf)
\end{Sinput}
\begin{Soutput}
  V043250       fit
1      18 -3.524611
2      90  9.742260
\end{Soutput}
\begin{Sinput}
 plot(th.bush.kerry ~ V043250, data = mydta1, xlab = "V403250: Age", ylab = "th.bush.kerry = Bush-Kerry thermometers",  col = gray(.7)) ## Scatterplot
 lines(fit ~ V043250, data = ndf, lwd = 1.5)
\end{Sinput}
\end{Schunk}
