\begin{Schunk}
\begin{Sinput}
 nd <- data.frame(V023010 =seq(0,100))
 nd$glmpred1 <- predict(glmthermo1, newdata=nd, type="response")
 nd$olspred1 <- predict(olsthermo1, newdata=nd)
 plot(nes2002$V023010, nes2002$bushvote, ylim=c(-0.2,1.2), xlim=c(0,100), xlab="Bush Thermometer", ylab="Predicted Probability of Bush Vote")
 lines(glmpred1 ~ V023010, data = nd, lty = 1)
 lines(olspred1 ~ V023010, data = nd, lty = 2)
 legend("left", legend = c("OLS", "Logistic"), lty = c(2,1))
\end{Sinput}
\end{Schunk}
