\begin{Schunk}
\begin{Sinput}
 #start adding variables 
 bushideoglm2 <- glm(bushvote~V023022+repub+democ,family=binomial, data=nes2002) 
 summary(bushideoglm2)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = bushvote ~ V023022 + repub + democ, family = binomial, 
    data = nes2002)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.8301  -0.4134   0.2648   0.4978   2.9631  

Coefficients:
            Estimate Std. Error z value Pr(>|z|)    
(Intercept) -2.37298    0.37868  -6.266 3.69e-10 ***
V023022      0.65359    0.08681   7.529 5.12e-14 ***
repub        1.78430    0.28728   6.211 5.26e-10 ***
democ       -2.65819    0.27611  -9.627  < 2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 1106.83  on 805  degrees of freedom
Residual deviance:  550.29  on 802  degrees of freedom
  (705 observations deleted due to missingness)
AIC: 558.29

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