# Set seed for this example
set.seed(11001101)
# Generate data for 4 binary methods
my_sim <- generate_multimethod_data(
"binary",
n_obs = 75,
n_method = 4,
se = c(0.87, 0.92, 0.79, 0.95),
sp = c(0.85, 0.93, 0.94, 0.80),
method_names = c("alpha", "beta", "gamma", "delta"))
# Bootstrap ML results
boot_ex <- boot_ML(
"binary",
data = my_sim$generated_data,
n_boot = 20)
# Aggregate Bootstrap ML results by statistic
aggregate_boot_ML(boot_ex)
# Plot Bootstrap ML estimate distributions
plot(boot_ex)
Run the code above in your browser using DataLab