design <- setup_fujikawa(k = 3, p0 = 0.2)
scenarios <- get_scenarios(design, p1 = 0.5)
# Without simulated data
opt_design(design, n = 20, alpha = 0.05, design_params =
list(epsilon = c(1, 2), tau = c(0, 0.5)), scenarios = scenarios,
prec_digits = 3)
# With simulated data
scenario_list <- as.list(data.frame(scenarios))
data_list <- lapply(scenario_list,
function(x) get_data(k = 3, n = 20, p = x, iter = 1000))
opt_design(design, n = 20, alpha = 0.05, design_params =
list(epsilon = c(1, 2), tau = c(0, 0.5)), scenarios = scenarios,
prec_digits = 3, data = data_list)
Run the code above in your browser using DataLab