我们如何在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