Как нам найти разные значения для каждой итерации в rstanarm?

Я написал основной код GLM в rstanarm, Когда я выполняю следующий код, он дает среднее значение переменных-предикторов всех выполненных итераций:

  summary(TestGLM) 
  print(TestGLM)  

Я хотел получить значение переменных-предикторов на каждой итерации для моего решения. Может кто-нибудь, пожалуйста, помогите мне решить эту проблему?

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
Другие вопросы по тегам