r - 呈现 rstanarm 输出结果的最佳方式？

Question

我想知道是否有人对显示 rstanarm 输出的最佳方式有任何建议。我通常喜欢 texreg::screenreg 和 stargazer，但没有一个接受 rstanarm 对象。如果它需要一个常规的 lm 对象以及直接将两个输出与外观相似的表进行比较，那么奖励！谢谢！

Answer 1

modelsummary开箱即用的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）：