test-generic.R
In expirest: Expiry Estimation Procedures

context("Generic summary, print and plot functions")

test_that("print_and_thus_summary.expirest_osle_succeeds", {
  re1 <-
    expirest_osle(data = exp1[exp1$Batch %in% c("b2", "b5", "b7"), ],
                  response_vbl = "Potency", time_vbl = "Month",
                  batch_vbl = "Batch", sl = 95, sl_sf = 2,
                  srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
                  xform = c("no", "no"), shift = c(0, 0), sf_option = "tight",
                  ivl = "confidence", ivl_type = "one.sided",
                  ivl_side = "lower")
  re2 <-
    expirest_osle(data = exp1[exp1$Batch %in% c("b4", "b5", "b8"), ],
                  response_vbl = "Potency", time_vbl = "Month",
                  batch_vbl = "Batch", sl = 95, sl_sf = 2,
                  srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
                  xform = c("no", "no"), shift = c(0, 0), sf_option = "tight",
                  ivl = "confidence", ivl_type = "one.sided",
                  ivl_side = "lower")
  re3 <- suppressWarnings(
    expirest_osle(
      data = exp3[exp3$Batch == "b1", ], response_vbl = "Moisture",
      time_vbl = "Month", batch_vbl = "Batch", sl = 1.5, sl_sf = 2,
      srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
      xform = c("no", "no"), shift = c(0, 0), sf_option = "tight",
      ivl = "prediction", ivl_type = "one.sided",
      ivl_side = "lower"))

  # <-><-><-><->

  expect_s3_class(expect_output(print(re1)), "expirest_osle")
  expect_output(print(re1), "acronym: cics")
  expect_output(print(re1), "batch: NA")

  expect_output(print(re1, digits = 5), "intercept: 100.57")
  expect_output(print(re1, digits = 5), "model: 25.996")
  expect_output(print(re1, digits = 5),
                "cics         100.57 25.996 lower    NA")
  expect_output(print(re1, digits = 5),
                "dics         100.36 24.567 lower    b2")
  expect_output(print(re1, digits = 5),
                "dids.pmse    100.78  23.47 lower    b5")
  expect_output(print(re1, digits = 5),
                "dids         100.78 23.148 lower    b5")

  expect_s3_class(expect_output(print(re2)), "expirest_osle")
  expect_output(print(re2), "acronym: dids")
  expect_output(print(re2), "batch: b8")

  expect_output(print(re2, digits = 5), "intercept: 101.26")
  expect_output(print(re2, digits = 5), "model: 15.845")
  expect_output(print(re2, digits = 5),
                "cics         101.55 27.925 lower    NA")
  expect_output(print(re2, digits = 5),
                "dics         100.49 22.267 lower    b8")
  expect_output(print(re2, digits = 5),
                "dids.pmse    101.26 15.606 lower    b8")
  expect_output(print(re2, digits = 5),
                "dids         101.26 15.845 lower    b8")

  expect_s3_class(expect_output(print(re3)), "expirest_osle")
  expect_output(print(re3), "acronym: n.a.")
  expect_output(print(re3), "NA intercepts and NA slopes")

  expect_output(print(re3, digits = 5), "intercept: NA")
  expect_output(print(re3, digits = 5), "model: NA")
  expect_output(print(re3, digits = 5),
                "cics             NA     NA lower    NA")
  expect_output(print(re3, digits = 5),
                "dics             NA     NA lower    NA")
  expect_output(print(re3, digits = 5),
                "dids.pmse        NA     NA lower    NA")
  expect_output(print(re3, digits = 5),
                "dids             NA     NA lower    NA")
})

test_that("print_and_plot.plot_expirest_osle_succeeds", {
  re1 <-
    expirest_osle(
      data = exp2, response_vbl = "Related", time_vbl = "Month",
      batch_vbl = "Batch", sl = 0.3, sl_sf = 1, srch_range = c(0, 500),
      alpha = 0.05, alpha_pool = 0.25, xform = c("no", "no"),
      shift = c(0, 0), sf_option = "loose", ivl = "confidence",
      ivl_type = "one.sided", ivl_side = "upper")

  # <-><-><-><->

  ggre1_1 <-
    plot_expirest_osle(
      model = re1, show_grouping = "yes", response_vbl_unit = "%",
      x_range = NULL, y_range = c(-0.01, 0.5), plot_option = "full",
      ci_app = "line")
  ggre1_2 <- expect_invisible(suppressWarnings(print(x = ggre1_1)))

  # <-><-><-><->

  expect_s3_class(ggre1_2, "plot_expirest_osle")
  expect_length(ggre1_2, 7)
  expect_s3_class(ggre1_2$Graph, c("gg", "ggplot"))
  expect_equal(
    ggre1_2$Graph$layers[[7]]$aes_params$label, c("USL: 0.34%", "19.0"))
})

test_that("print_and_thus_summary.expirest_wisle_succeeds", {
  t_dat <- exp1[exp1$Batch %in% c("b4", "b5", "b8"), ]

  # <-><-><-><->

  re1 <-
    expirest_wisle(
      data = exp3, response_vbl = "Moisture", time_vbl = "Month",
      batch_vbl = "Batch", rl = 2.0, rl_sf = 3, sl = 0.5, sl_sf = 1,
      srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
      xform = c("no", "no"), shift = c(0, 0), sf_option = "loose",
      ivl = "confidence", ivl_type = "one.sided", ivl_side = "lower")
  re2 <-
    expirest_wisle(
      data = t_dat, response_vbl = "Potency", time_vbl = "Month",
      batch_vbl = "Batch", rl = 98, rl_sf = 2, sl = 95, sl_sf = 2,
      srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
      xform = c("no", "no"), shift = c(0, 0), sf_option = "loose",
      ivl = "confidence", ivl_type = "one.sided", ivl_side = "lower")
  re3 <-
    expirest_wisle(
      data = t_dat, response_vbl = "Potency", time_vbl = "Month",
      batch_vbl = "Batch", rl = seq(98, 99, 0.5), rl_sf = rep(3, 3), sl = 95,
      sl_sf = 2, srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
      xform = c("no", "no"), shift = c(0, 0), sf_option = "loose",
      ivl = "confidence", ivl_type = "one.sided", ivl_side = "lower")

  # <-><-><-><->

  expect_s3_class(expect_output(print(re1)), "expirest_wisle")
  expect_output(print(re1), "cics")
  expect_output(print(re1, digits = 5), "2    NA    2.4568")
  expect_output(print(re1, digits = 5), "0.5  2 92.999 118.19")

  expect_s3_class(expect_output(print(re2)), "expirest_wisle")
  expect_output(print(re2), "dids")
  expect_output(print(re2, digits = 5), "98    b8    101.26")
  expect_output(print(re2, digits = 5),
                "95 98 7.6197       7.4832 17.039      16.777")

  expect_output(print(re3, digits = 5),
                "95 98.0  8.8631        8.725 17.039      16.777")
  expect_output(print(re3, digits = 5),
                "95 98.5 10.1748       11.221 17.039      16.777")
  expect_output(print(re3, digits = 5),
                "95 99.0 11.4405       11.277 17.039      16.777")
})

test_that("print_and_plot.plot_expirest_wisle_succeeds", {
  re1 <-
    expirest_wisle(
      data = exp2, response_vbl = "Related", time_vbl = "Month",
      batch_vbl = "Batch", rl = 0.15, rl_sf = 2, sl = 0.3, sl_sf = 1,
      srch_range = c(0, 500), alpha = 0.05, alpha_pool = 0.25,
      xform = c("no", "no"), shift = c(0, 0), sf_option = "loose",
      ivl = "confidence", ivl_type = "one.sided", ivl_side = "upper")

  # <-><-><-><->

  ggre1_1 <-
    plot_expirest_wisle(
      model = re1, rl_index = 1, show_grouping = "yes",
      response_vbl_unit = "%", x_range = NULL, y_range = c(-0.01, 0.50),
      scenario = "standard", plot_option = "full", ci_app = "line")
  ggre1_2 <- expect_invisible(suppressWarnings(print(x = ggre1_1)))

  # <-><-><-><->

  expect_s3_class(ggre1_2, "plot_expirest_wisle")
  expect_length(ggre1_2, 9)
  expect_s3_class(ggre1_2$Graph, c("gg", "ggplot"))
  expect_equal(
    ggre1_2$Graph$layers[[6]]$aes_params$label,
    c("USL: 0.34%", "0.298% ", "0.112% ", "15.7\n(worst case\nscenario)",
      "19.0\n(standard\nscenario)", "URL: 0.154%"))
})