显示rstanarm输出结果的最佳方式？

Question

我想知道是否有人对显示rstanarm输出的最佳方式有什么建议。我通常喜欢texreg::screenreg和stargazer，但它们都不接受rstanarm对象。如果它是一个常规的lm对象，用于直接比较两个外观相似的表的输出，那就更好了！谢谢!

Answer 1

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定义为置信区间，因为此模型类型不会产生标准错误。

生成的超文本标记语言表格如下所示(其他可用格式: LaTeX、JPG、Word)：

​

​