tests/testthat/helper-model-mixture3p.R

In bmm: Easy and Accessible Bayesian Measurement Models Using 'brms'

# the cases the shared tests run for mixture3p; see helper-model-cases.R
model_cases_mixture3p <- function() {
  model <- mixture3p("dev_rad", nt_features = paste0("col_nt", 1:7), set_size = "set_size")
  re_formula <- bmf(kappa ~ 1 + (1 | ID), thetat ~ 1 + (1 | ID), thetant ~ 1 + (1 | ID))
  list(
    # brms has no reason to keep set_size in the first case's frame; the second
    # names it, so brms keeps the real column instead
    stored_frame = list(
      mixture3p = list(model = model, formula = bmf(kappa ~ 1, thetat ~ 1, thetant ~ 1),
                       data = small_data()),
      mixture3p_set_size = list(model = model,
                                formula = bmf(kappa ~ 1, thetat ~ 1, thetant ~ 0 + set_size),
                                data = small_data())
    ),
    sd_default = list(
      mixture3p = list(
        model = model, formula = re_formula, data = small_data(),
        sd = c(kappa = "exponential(1)", thetat = "exponential(1)", thetant = "exponential(1)")
      )
    ),
    init = list(mixture3p = list(model = model, formula = re_formula, data = small_data())),
    links = list(mixture3p = list(model = model, settable = character(0)))
  )
}