Как нам найти разные значения для каждой итерации в 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