\begin{Schunk}
\begin{Sinput}
 bushideoglm4 <- glm(bushvote ~ ideoln * partyid, family=binomial, data=nes2002) 
 summary(bushideoglm4)
\end{Sinput}
\begin{Soutput}
Call:
glm(formula = bushvote ~ ideoln * partyid, family = binomial, 
    data = nes2002)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.8414  -0.4345   0.2624   0.3635   2.5001  

Coefficients:
                Estimate Std. Error z value Pr(>|z|)    
(Intercept)      -3.9410     0.6768  -5.823 5.79e-09 ***
ideoln            1.0171     0.1619   6.283 3.31e-10 ***
partyidD          0.6046     0.9363   0.646 0.518440    
partyidR          3.2869     1.1228   2.927 0.003419 ** 
ideoln:partyidD  -0.7612     0.2233  -3.409 0.000651 ***
ideoln:partyidR  -0.3495     0.2457  -1.423 0.154873    
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 1074.69  on 780  degrees of freedom
Residual deviance:  507.24  on 775  degrees of freedom
  (730 observations deleted due to missingness)
AIC: 519.24

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