我们如何在 rstanarm 中找到每次迭代的不同值?

How do we find the different values for each iteration in rstanarm?

我确实在 rstanarm 中编写了基本的 GLM 代码。当我执行以下代码时,它给出了所有已完成迭代的预测变量的平均值:

  summary(TestGLM) 
  print(TestGLM)  

我想要我的解决方案在每次迭代时的预测变量值。谁能帮我解决这个问题?

您可以按如下方式提取每个系数的绘图:

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