r - 基于来自`rstanarm`包中`stan_glm（）`的分组变量的后验预测？

Question

stan_glm()我想知道如何根据包中的分组变量获得后验预测rstanarm？

例如，如果我的数据中有一个二进制(0, 1)编码的分组变量"vs"（基本 R 数据：）mtcars，我如何获得何时vs == 0和何时的预测vs == 1？

这是我的 R 代码：

library(rstanarm)
fit <- stan_glm(mpg ~., data = mtcars)

posterior_predict(fit, newdata = WHAT SHOULD BE HERE?)

Answer 1

要探索 egvs对结果的影响（在您的情况下mpg），您可以分别posterior_predict在 wherevs == 0和子集上使用vs == 1：

posterior_predict(fit, newdata = subset(mtcars[1:10, ], vs == 0));

和

posterior_predict(fit, newdata = subset(mtcars[1:10, ], vs == 1));

更多细节在?rstanarm::posterior_predict.