# NOT RUN {
library(bayestestR)
if (require("rstanarm")) {
model <- stan_glm(
Sepal.Length ~ Petal.Width * Species,
data = iris,
chains = 2, iter = 200, refresh = 0
)
x <- rope(model)
plot(x)
x <- hdi(model)
plot(x) + theme_modern()
data <- rnorm(1000, 1)
x <- p_direction(data)
plot(x)
x <- p_direction(model)
plot(x)
model <- stan_glm(
mpg ~ wt + gear + cyl + disp,
chains = 2,
iter = 200,
refresh = 0,
data = mtcars
)
x <- equivalence_test(model)
plot(x)
}
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab