# NOT RUN {
# Rstanarm ----------------------
if (require("rstanarm")) {
# Simple regressions
model <- stan_glm(Sepal.Length ~ Species,
data = iris, chains = 2, refresh = 0
)
bayesian_as_frequentist(model)
}
# }
# NOT RUN {
# }
# NOT RUN {
if (require("rstanarm")) {
model <- stan_glm(vs ~ mpg,
family = "binomial",
data = mtcars, chains = 2, refresh = 0
)
bayesian_as_frequentist(model)
# Mixed models
model <- stan_glmer(Sepal.Length ~ Petal.Length + (1 | Species),
data = iris, chains = 2, refresh = 0
)
bayesian_as_frequentist(model)
model <- stan_glmer(vs ~ mpg + (1 | cyl),
family = "binomial",
data = mtcars, chains = 2, refresh = 0
)
bayesian_as_frequentist(model)
}
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab