\begin{Schunk}
\begin{Sinput}
 nes2002$educFactor <- as.factor(nes2002$V023131)
 bushideolrm6 <- glm(bushvote ~ V023010+ V023022+ repub+democ+ V023027+ educFactor, data=nes2002, family=binomial(link=logit))
 summary(bushideolrm6)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = bushvote ~ V023010 + V023022 + repub + democ + 
    V023027 + educFactor, family = binomial(link = logit), data = nes2002)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.7569  -0.2847   0.1399   0.3475   3.3258  

Coefficients:
             Estimate Std. Error z value Pr(>|z|)    
(Intercept) -7.048866   1.394368  -5.055  4.3e-07 ***
V023010      0.068760   0.007817   8.796  < 2e-16 ***
V023022      0.368200   0.108234   3.402 0.000669 ***
repub        1.172697   0.333312   3.518 0.000434 ***
democ       -2.792379   0.330018  -8.461  < 2e-16 ***
V023027     -0.082868   0.125623  -0.660 0.509473    
educFactor2  1.393612   1.349233   1.033 0.301655    
educFactor3  2.176177   1.107965   1.964 0.049516 *  
educFactor4  2.459526   1.117383   2.201 0.027725 *  
educFactor5  1.516711   1.141784   1.328 0.184056    
educFactor6  1.308992   1.098129   1.192 0.233253    
educFactor7  1.872912   1.119012   1.674 0.094186 .  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 1098.89  on 799  degrees of freedom
Residual deviance:  409.45  on 788  degrees of freedom
  (711 observations deleted due to missingness)
AIC: 433.45

Number of Fisher Scoring iterations: 6
\end{Soutput}
\end{Schunk}
