rm(list = ls())
devtools::load_all()

designpath <- system.file("extdata", "feedadditives", package = "simulateDCE")

# notes <- "This design consists of different heuristics. One group did not attend the methan attribute, another group only decided based on the payment"

notes <- "Three heuristics"

resps <- 396 # number of respondents
nosim <- 2 # number of simulations to run (about 500 is minimum)
destype <- "ngene"

# betacoefficients should not include "-"
bcoeff <- list(
  basc = 0.2,
  bcow = 0.3,
  badv = 0.3,
  bvet = 0.3,
  bfar = 0.3,
  bmet = 0.3,
  bbon = 0.3,
  bbon2 = 1.9,
  basc2 = 2
)

decisiongroups <- c(0, 0.3, 0.6, 0.8, 1)

# place your utility functions here
ul <- list(
  u1 = list(
    v1 = V.1 ~ bcow * alt1.cow + badv * alt1.adv + bvet * alt1.vet + bfar * alt1.far + bmet * alt1.met + bbon * alt1.bon,
    v2 = V.2 ~ bcow * alt2.cow + badv * alt2.adv + bvet * alt2.vet + bfar * alt2.far + bmet * alt2.met + bbon * alt2.bon,
    v3 = V.3 ~ basc
  ),
  u2 = list(
    v1 = V.1 ~ bcow * alt1.cow + badv * alt1.adv + bvet * alt1.vet + bfar * alt1.far + bbon * alt1.bon,
    v2 = V.2 ~ bcow * alt2.cow + badv * alt2.adv + bvet * alt2.vet + bfar * alt2.far + bbon * alt2.bon,
    v3 = V.3 ~ basc
  ),
  u3 = list(
    v1 = V.1 ~ bbon2 * alt1.bon,
    v2 = V.2 ~ bbon2 * alt2.bon,
    v3 = V.3 ~ basc
  ),
  u4 = list(
    v1 = V.1 ~ basc2 + bcow * alt1.cow + badv * alt1.adv + bvet * alt1.vet + bfar * alt1.far + bmet * alt1.met + bbon * alt1.bon,
    v2 = V.2 ~ bcow * alt2.cow + badv * alt2.adv + bvet * alt2.vet + bfar * alt2.far + bmet * alt2.met + bbon * alt2.bon,
    v3 = V.3 ~ basc
  )
)

feedadditives <- sim_all(
  nosim = nosim, resps = resps, designtype = destype,
  designpath = designpath, u = ul, bcoeff = bcoeff, decisiongroups = decisiongroups
)