Code: extended model support for mtable
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:
coxph
,survreg
– Cox proportional hazards models and parametric survival models from thesurvival
library.aftreg
,phreg
,weibreg
– parametric AFT, proportional hazards, and Weibull models available from theeha
library.mer
– multilevel mixed effects models as provided by thelme4
package. The model objects produced in the betalme4a
package are not 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.
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.