\begin{Schunk}
\begin{Sinput}
 bushlogit5b <- glm(bushvote ~ V023010 + ideoln + educn, 
  data=model.frame(bushlogit5), family=binomial(link=logit), model=TRUE)
 anova(bushlogit5, bushlogit5b, test = "Chisq")
\end{Sinput}
\begin{Soutput}
Analysis of Deviance Table

Model 1: bushvote ~ V023010 + ideoln + partyid + educn + useconn
Model 2: bushvote ~ V023010 + ideoln + educn
  Resid. Df Resid. Dev Df Deviance  Pr(>Chi)    
1       768     404.64                          
2       771     556.54 -3   -151.9 < 2.2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
\end{Soutput}
\end{Schunk}
