微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

呈现 rstanarm 输出结果的最佳方式?

如何解决呈现 rstanarm 输出结果的最佳方式?

我想知道是否有人对显示 rstanarm 输出的最佳方式有任何建议。我通常喜欢 texreg::screenreg 和 stargazer,但都不接受 rstanarm 对象。如果它也需要一个常规的 lm 对象来直接比较具有相似外观的表的两个输出,那就额外奖励了!谢谢!

解决方法

modelsummary package for R 支持由 rstanarm 包开箱即用生成的模型(免责声明:我是作者)。

例如,如果您想估计两个模型并在并排回归表中显示结果,您可以这样做:

library(rstanarm)
library(modelsummary)

mod1 <- stan_glm(hp ~ mpg,data = mtcars)
mod2 <- stan_glm(hp ~ mpg + drat,data = mtcars)

models <- list(mod1,mod2)

modelsummary(models,statistic = "conf.int")

请注意,我必须手动将 statistic 定义为置信区间,因为此模型类型不会产生标准误差。

生成的 HTML 表格如下所示(其他可用格式:LaTeX、JPG、Word):

enter image description here

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。