#
#   Generate am example rstanarm logit fitted model object
#

if (RUN_STAN && requireNamespace("rstanarm", quietly = TRUE) &&
    requireNamespace("carData", quietly = TRUE)) {
  
  library(rstanarm)
  
  df <- carData::Cowles
  df$female <- (as.numeric(df$sex) - 2) * (-1)
  df$volunteer <- as.numeric(df$volunteer) - 1
  df$extraversion <- (df$extraversion - mean(df$extraversion)) / (2 * sd(df$extraversion))
  df$neuroticism <- (df$neuroticism - mean(df$neuroticism)) / (2 * sd(df$neuroticism))
  
  rstanarm_logit <- fit_or_warn("rstanarm", stan_glm(volunteer ~ female + neuroticism + extraversion, 
                                                     data = df, family = binomial(link = "logit"),
                                                     prior = normal(0, 3),
                                                     prior_intercept = normal(0, 3),
                                                     chains = 4, 
                                                     iter = 2000,
                                                     seed = 123))
  
  if (!is.null(rstanarm_logit)) {
    #saveRDS(rstanarm_logit, "tests/testdata/rstanarm-logit.rds")
    saveRDS(rstanarm_logit, file.path(TESTDATA_DIR, "rstanarm_logit.rds"))
  }
  
}