rm(list = ls()) devtools::load_all() designpath <- system.file("extdata", "SE_AGRI", package = "simulateDCE") destype <- "ngene" resps <- 360 # number of respondents nosim <- 2 # number of simulations to run (about 500 is minimum) # betacoefficients should not include "-" bcoeff <- list( basc = 4.2, ## very high asc bprof = 0.3, bexp = 0.3, bdomestic = 0.3, bforeign = 0.3, bdamage = 0.6, bprice = 0.2 ) manipulations <- list( alt1.professional = rlang::expr(alt1.initiator == 1), alt2.professional = rlang::expr(alt2.initiator == 1), alt1.expert = rlang::expr(alt1.initiator == 2), alt2.expert = rlang::expr(alt2.initiator == 2), alt1.domestic = rlang::expr(alt1.funding == 1), alt2.domestic = rlang::expr(alt2.funding == 1), alt1.foreign = rlang::expr(alt1.funding == 2), alt2.foreign = rlang::expr(alt2.funding == 2) ) # place your utility functions here ul <- list( u1 = list( v1 = V.1 ~ bprof * alt1.professional + bexp * alt1.expert + bdomestic * alt1.domestic + bforeign * alt1.foreign + bdamage * alt1.damage + bprice * alt1.compensation, v2 = V.2 ~ bprof * alt2.professional + bexp * alt2.expert + bdomestic * alt2.domestic + bforeign * alt2.foreign + bdamage * alt2.damage + bprice * alt2.compensation, v3 = V.3 ~ basc ) ) seagri <- sim_all( nosim = nosim, resps = resps, designtype = destype, designpath = designpath, u = ul, bcoeff = bcoeff, manipulations = manipulations, utility_transform_type = "exact" )