rstanarm: визуализировать / извлечь данные о накоплении доказательств (сдвиг от предыдущего к заднему) - PullRequest
3 голосов
/ 25 мая 2019

Я хотел бы визуализировать сдвиг от предыдущего к последнему, когда данные собираются (для каждой новой точки данных).

Я могу визуализировать как предшествующий, так и задний:

library(ggplot2)
library(rstanarm)
library(insight)
library(bayestestR)

model <- stan_glm(mpg ~ wt, data = mtcars, chains = 2, iter = 500)

priors <- insight::get_priors(model)
priors <- priors[priors$parameter == "wt", ]

prior <- bayestestR::distribution_normal(n = 500,
                                        mean = priors$location,
                                        sd = priors$adjusted_scale)
prior <- data.frame(x = prior, type = "prior")

posterior <- data.frame(x = insight::get_parameters(model)$wt, type = "posterior")
df <- rbind(prior, posterior)


ggplot(df, aes(x=x, fill=type)) +
  geom_density(alpha=0.5)

Создано в 2019-05-25 пакетом Представить (v0.3.0)

Однако я быхотелось бы иметь все «промежуточные» шаги (прогрессивный сдвиг от предыдущего к заднему для каждой новой точки данных).Есть ли способ извлечь это из модели?Спасибо!

...