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

Question

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

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

这是我的R代码：

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

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

Answer 1

探索例如关于结果的vs（在您的情况下为mpg），您可以分别在posterior_predict和vs == 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中提供了更多详细信息。