# NOT RUN {
# Rstanarm ----------------------
if (require("rstanarm")) {
  # Simple regressions
  model <- stan_glm(Sepal.Length ~ Petal.Length * Species,
    data = iris, chains = 2, refresh = 0
  )
  bayesian_as_frequentist(model)
  model <- stan_glm(vs ~ mpg,
    family = "binomial",
    data = mtcars, chains = 2, refresh = 0
  )
  bayesian_as_frequentist(model)
}
# }
# NOT RUN {
# }
# NOT RUN {
if (require("rstanarm")) {
  # 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