\begin{Schunk}
\begin{Sinput}
 glmpred1 <- predict(glmthermo1, newdata=data.frame(numV023010 =seq(0,100)), type="response")
 plot(nes2002$V023010, nes2002$bushvote, ylim=c(-0.2,1.2), xlim=c(0,100), xlab="Bush Thermometer",ylab="Predic6.1.2.5ted Probability of Bush Vote")
 lines(seq(0,100),glmpred1)
 abline(olsthermo1,lty=2)
\end{Sinput}
\end{Schunk}
