test_that("check_scenarios() returns a plain data frame", {
  des <- ref_design("D1")
  sc <- check_scenarios(des, scen_c)
  expect_identical(sc, scen_c)
  # A tibble-like data frame is converted
  tb <- structure(scen_c, class = c("tbl_df", "tbl", "data.frame"))
  expect_identical(class(check_scenarios(des, tb)), "data.frame")
})

test_that("check_scenarios() validates the scenarios", {
  des <- ref_design("D1")
  expect_error(check_scenarios(des, scen_c[0, ]),
               "'scenarios' must be a data frame with at least one row.")
  expect_error(check_scenarios(des, scen_c, "sc", one_row = TRUE),
               "'sc' must have exactly one row.")
  expect_error(check_scenarios(des, data.frame(pi_c = 0.3)),
               "must have the column\\(s\\) pi_t, pi_c.")
  expect_error(check_scenarios(ref_design("D3"), data.frame(x = 0.3)),
               "must have the column\\(s\\) pi_t.")
  expect_error(check_scenarios(des, data.frame(pi_t = 0.3,
                                               pi_c = NA_real_)),
               "'scenarios\\$pi_c' must not contain missing values.")
  expect_error(check_scenarios(des, data.frame(pi_t = -0.1, pi_c = 0.2)),
               "'scenarios\\$pi_t' must be at least 0 and at most 1.")
  expect_error(check_scenarios(des, data.frame(pi_t = "a", pi_c = 0.2)),
               "'scenarios\\$pi_t' must be numeric.")
  expect_error(check_scenarios(des, data.frame(pi_t = 0.3, pi_c = 0.2,
                                               Miss = 0)),
               "must not have columns named Go, NoGo, Gray or Miss.")
})

test_that("check_scenarios() validates the scenarios of continuous designs", {
  des <- ref_design_cont("C1")
  expect_identical(check_scenarios(des, scen_cont_c), scen_cont_c)
  expect_error(check_scenarios(des, data.frame(mu_t = 1, sd_t = 1, mu_c = 0)),
               "must have the column\\(s\\) mu_t, sd_t, mu_c, sd_c.")
  expect_error(check_scenarios(des, data.frame(mu_t = 1, sd_t = 1, mu_c = 0,
                                               sd_c = -1)),
               "'scenarios\\$sd_c' must be greater than 0.")
  expect_error(check_scenarios(des, data.frame(mu_t = Inf, sd_t = 1,
                                               mu_c = 0, sd_c = 1)),
               "'scenarios\\$mu_t' must be finite.")
  u <- data.frame(mu_t = c(-5, 5), sd_t = c(1, 2))
  expect_identical(check_scenarios(ref_design_cont("C3"), u), u)
})

test_that("check_scenarios() validates scenarios for two binary endpoints", {
  des <- ref_design_bin2("B1")
  expect_identical(check_scenarios(des, scen_bin2), scen_bin2)
  bad <- scen_bin2[1, ]
  bad$pi_t1 <- 1
  expect_error(check_scenarios(des, bad),
               "'scenarios\\$pi_t1' must be greater than 0 and less than 1.")
  bad <- scen_bin2[1, ]
  bad$rho_c <- 1.5
  expect_error(check_scenarios(des, bad),
               "'scenarios\\$rho_c' must be at least -1 and at most 1.")
  # Infeasible correlations are accepted unless feasibility is required
  inf <- scen_bin2[1, ]
  inf$rho_t <- -0.9
  expect_identical(check_scenarios(des, inf), inf)
  expect_error(check_scenarios(des, inf, "sc", feasible = TRUE),
               "'sc\\$rho_t' is outside the feasible range")
  inf_c <- scen_bin2[1, ]
  inf_c$rho_c <- -0.9
  expect_error(check_scenarios(des, inf_c, "sc", feasible = TRUE),
               "'sc\\$rho_c' is outside the feasible range")
})

test_that("check_scenarios() validates scenarios of two continuous endpoints", {
  des <- ref_design_cont2("K1")
  sc <- scen_cont2("K1")
  expect_identical(check_scenarios(des, sc), sc)
  bad <- sc[1, ]
  bad$rho_t <- 1
  expect_error(check_scenarios(des, bad),
               "'scenarios\\$rho_t' must be greater than -1 and less than 1.")
  bad <- sc[1, ]
  bad$sd_c2 <- 0
  expect_error(check_scenarios(des, bad),
               "'scenarios\\$sd_c2' must be greater than 0.")
  expect_error(check_scenarios(des, sc[, 1:5]),
               "must have the column\\(s\\) mu_t1, mu_t2, sd_t1, sd_t2, rho_t")
})