如何在rstanarm中为每次迭代找到不同的值？

Question

我确实用rstanarm编写了一个基本的GLM代码。当我执行以下代码时，它给出了所执行的所有迭代的预测变量的平均值：

  summary(TestGLM) 
  print(TestGLM)  

我希望在我的解决方案的每一次迭代中都要预测变量的值。有人能帮我解决这个问题吗？

Answer 1

可以提取每个系数的绘图，如下所示：

require(rstanarm)

TestGLM <- stan_glm(mpg ~ wt, data = mtcars)
draws <- as.data.frame(TestGLM)
print(colnames(draws))
str(draws)

# See http://mc-stan.org/rstanarm/reference/as.matrix.stanreg.html