\begin{Schunk}
\begin{Sinput}
 library(mix)
 se.glm <- MIextract(allimplogreg, fun = function(x){sqrt(diag(vcov(x)))}) 
 as.data.frame(mi.inference(betas, se.glm))
\end{Sinput}
\begin{Soutput}
              est std.err  df  signif  lower  upper    r fminf
(Intercept)  4.03   0.532  23 1.2e-07  2.926  5.128 1.72  0.66
pclass2nd   -1.44   0.258 129 1.3e-07 -1.951 -0.931 0.36  0.27
pclass3rd   -2.71   0.308  31 6.3e-10 -3.339 -2.081 1.17  0.57
sexmale     -2.53   0.175 300 0.0e+00 -2.872 -2.184 0.21  0.18
age         -0.05   0.011  18 2.5e-04 -0.073 -0.027 2.39  0.73
\end{Soutput}
\end{Schunk}
