\begin{Schunk}
\begin{Sinput}
 options(digits=3)
 logitsummary <- summary(bushlogit5)
 probitsummary <- summary(bushprobit5)
 zcompare <- cbind(logitsummary$coefficients[, c(3,4)], probitsummary$coefficients[, c(3,4)])
 zcompare <- as.data.frame(zcompare)
 colnames(zcompare) <- c("logit-z","logit-p", "probit-z", "probit-p")
 zcompare$ratioz <- zcompare[ ,1] / zcompare[ ,3]
 zcompare
\end{Sinput}
\begin{Soutput}
            logit-z  logit-p probit-z probit-p ratioz
(Intercept)  -4.958 7.13e-07    -4.77 1.82e-06  1.039
V023010       8.464 2.58e-17     8.72 2.66e-18  0.970
ideoln        3.385 7.12e-04     3.52 4.25e-04  0.961
partyidD     -7.899 2.80e-15    -8.32 8.46e-17  0.949
partyidR      3.904 9.48e-05     4.00 6.44e-05  0.977
useconn      -0.444 6.57e-01    -1.49 1.37e-01  0.299
educn        -0.860 3.90e-01    -0.32 7.49e-01  2.686
\end{Soutput}
\end{Schunk}
