r语言 - 基于“rstanarm”包中“stan_glm()”分组变量的后验预测

Question

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

例如，如果我的数据中有一个名为 "vs" 的二进制(0, 1)编码分组变量(基本 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

探索例如 vs结果(在您的mpg情况下(，您可以分别对vs == 0和vs == 1的子集使用posterior_predict：

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

和

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

更多详细信息在?rstanarm::posterior_predict中给出。