\begin{Schunk}
\begin{Sinput}
 options(digits=3)
 logitsummary <- summary(bushideolrm5)
 probitsummary <- summary(bushideoProbit5)
 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.832 1.35e-06   -4.696 2.65e-06  1.029
V023010       8.814 1.21e-18    9.138 6.39e-20  0.965
V023022       3.344 8.25e-04    3.518 4.35e-04  0.951
repub         3.558 3.74e-04    3.633 2.81e-04  0.980
democ        -8.470 2.45e-17   -8.929 4.31e-19  0.949
V023027      -0.349 7.27e-01   -0.262 7.93e-01  1.328
V023131      -1.273 2.03e-01   -1.812 7.00e-02  0.702
\end{Soutput}
\end{Schunk}
