\begin{Schunk}
\begin{Sinput}
 summary(bush1)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = pres04 ~ partyid + sex + owngun, family = binomial(link = logit), 
    data = dat)

Deviance Residuals: 
   Min      1Q  Median      3Q     Max  
-2.941  -0.488   0.163   0.390   2.683  

Coefficients:
                        Estimate Std. Error z value
(Intercept)              -3.5712     0.3934   -9.08
partyidDem.               1.9103     0.3972    4.81
partyidInd. Near Dem.     1.4559     0.4348    3.35
partyidIndependent        3.4642     0.4105    8.44
partyidInd. Near Repub.   5.4677     0.5073   10.78
partyidRepub.             6.0307     0.4502   13.39
partyidStrong Repub.      7.1908     0.6213   11.57
sexFemale                 0.0488     0.1928    0.25
owngunYES                 0.6424     0.1937    3.32
                        Pr(>|z|)    
(Intercept)              < 2e-16 ***
partyidDem.              1.5e-06 ***
partyidInd. Near Dem.    0.00081 ***
partyidIndependent       < 2e-16 ***
partyidInd. Near Repub.  < 2e-16 ***
partyidRepub.            < 2e-16 ***
partyidStrong Repub.     < 2e-16 ***
sexFemale                0.80006    
owngunYES                0.00091 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 1721.9  on 1242  degrees of freedom
Residual deviance:  764.0  on 1234  degrees of freedom
  (3267 observations deleted due to missingness)
AIC: 782

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