Site icon R-bloggers

Code: extended model support for mtable

[This article was first published on Left Censored » R, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.

I finally got around to organizing and packaging my complete set of extended model support for mtable in Martin Elff’s memisc library. Here is a list of the models supported:

Here’s the file: mtable-ext

It’s extremely easy to use, though I am sure there are some bugs and not all model features get reported. Here is an example of how to compare three Cox proportional hazard models as estimated with coxph (yes, I know, the Breslow and Exact models are the same; Exact was refusing to converge so I just stuck the Breslow model in to get the example done):

library(survival)
library(memisc)
source("mtable-ext.R")

> Efron <- coxph(Surv(start, end, switch) ~ Poll + Government + OldParty +
+                Experience + Right + Members, robust = TRUE, data = Sejm,
+                method = "efron")

> Breslow <- coxph(Surv(start, end, switch) ~ Poll + Government + OldParty +
+                  Experience + Right + Members, robust = TRUE, data = Sejm,
+                  method = "breslow")

> Exact <- coxph(Surv(start, end, switch) ~ Poll + Government + OldParty +
+                Experience + Right + Members, robust = TRUE, data = Sejm,
+                method = "exact")

mtable(Efron, Breslow, Exact)

[... some model call info ...]

=============================================
                  Efron    Breslow    Exact
---------------------------------------------
Poll            -2.120*** -1.931*** -1.931***
                (0.477)   (0.430)   (0.430)
Government      -1.603*   -1.461*   -1.461*
                (0.706)   (0.684)   (0.684)
OldParty        -0.186    -0.279    -0.279
                (0.243)   (0.209)   (0.209)
Experience      -0.195    -0.198    -0.198
                (0.167)   (0.156)   (0.156)
Right            1.427***  1.290***  1.290***
                (0.198)   (0.178)   (0.178)
Members         -0.002    -0.003    -0.003
                (0.007)   (0.007)   (0.007)
---------------------------------------------
Log-likelihood   -818.487  -832.625  -832.625
AIC              1648.973  1677.250  1677.250
BIC              1701.297  1729.574  1729.574
N               45280     45280     45280
=============================================

To convert the mtable object into a table, you just need to apply the toLatex function:

toLatex(mtable(Efron, Breslow, Exact))

Of course, many formatting options are available. See the documentation for details.

To leave a comment for the author, please follow the link and comment on their blog: Left Censored » R.

R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.