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 )