\begin{Schunk}
\begin{Sinput}
 df1 <- as.data.frame(mi.inference(betas, se.glm))
 df2 <- cbind(df1[,1:2], ldbeta = coef(ldglm), ldse = sqrt(diag(vcov(ldglm))))
 df2
\end{Sinput}
\begin{Soutput}
              est std.err ldbeta   ldse
(Intercept)  4.03   0.532  4.371 0.4709
pclass2nd   -1.44   0.258 -1.462 0.2847
pclass3rd   -2.71   0.308 -2.744 0.3403
sexmale     -2.53   0.175 -3.077 0.2414
age         -0.05   0.011 -0.046 0.0088
\end{Soutput}
\end{Schunk}
