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")