Nothing
#
# Example file for creating a simple MNL design
#
rm(list = ls(all = TRUE))
# library(spdesign)
# Define the list of utility functions ----
utility <- list(
sq = "b_sq[0.1] * sq[1] ",
alt1 = " b_new_wind_mills_dummy[c(-0.2, -0.4)] * new_wind_mills[c(1500, 750, 500)] + b_distance[0.15] * distance_km[c(0.25, 0.50, 0.75, 1, 1.25)] + b_fee[-0.25] * fee[c(5, 10, 15, 20, 25)]",
alt2 = " b_new_wind_mills_dummy * new_wind_mills + b_distance * distance_km + b_fee * fee"
)
# Generate designs ----
design <- generate_design(utility, rows = 20,
model = "mnl", efficiency_criteria = "d-error",
algorithm = "federov", draws = "scrambled-sobol",
control = list(
max_iter = 15000,
efficiency_threshold = 0.01
))
# Add a blocking variable to the design with 2 blocks.
design <- block(design, 2)
summary(design)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.