tests/testthat/test-log.R

In TOSTER: Two One-Sided Tests (TOST) Equivalence Testing

# need hush function to run print through examples

hush = function(code) {
  sink("NUL") # use /dev/null in UNIX
  tmp = code
  sink()
  return(tmp)
}


test_that("Run examples for two sample", {

  set.seed(2525894)
  dat1 = mtcars
  expect_error(log_TOST(
    x = dat1$mpg,
    var.equal = TRUE,
    eqb = 1.25
  ))
  test1 = boot_log_TOST(data = dat1,
                        mpg ~ am,
                        var.equal = TRUE,
                        eqb = 1.25,
                        R = 199)

  test1_t = log_TOST(data = dat1,
                        mpg ~ am,
                        var.equal = TRUE,
                        eqb = 1.25)

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)

  test1 = boot_log_TOST(data = dat1,
                        mpg ~ am,
                        var.equal = TRUE,
                        eqb = 1.75,
                        R = 199)

  test1_t = log_TOST(data = dat1,
                     mpg ~ am,
                     var.equal = TRUE,
                     eqb = 1.75)

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)

  test1 = boot_log_TOST(data = dat1,
                        mpg ~ am,
                        var.equal = TRUE,
                        hypothesis = "MET",
                        eqb = 1.25,
                        R = 199)

  test1_t = log_TOST(data = dat1,
                     mpg ~ am,
                     var.equal = TRUE,
                     hypothesis = "MET",
                     eqb = 1.25)

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)

  test2 = boot_log_TOST(data = dat1,
                        mpg ~ am,
                        var.equal = FALSE,
                        eqb = 1.25,
                        R = 199)

  test2_t = log_TOST(data = dat1,
                     mpg ~ am,
                     var.equal = FALSE,
                     eqb = 1.25)

  expect_equal(test2$effsize$estimate,
               test2_t$effsize$estimate)

  test2 = boot_log_TOST(data = dat1,
                        mpg ~ am,
                        var.equal = FALSE,
                        eqb = 1.25,
                        R = 199,
                        hypothesis = "MET")

  test2_t = log_TOST(data = dat1,
                     mpg ~ am,
                     var.equal = FALSE,
                     eqb = 1.25,
                     hypothesis = "MET")

  expect_equal(test2$effsize$estimate,
               test2_t$effsize$estimate)

})


test_that("Run examples for paired samples", {

  set.seed(922287)
  sleep2 = sleep
  sleep2$sleep = sleep2$extra + 4
  expect_error(boot_log_TOST(
    data = sleep2,
    extra ~ group,
    paired = TRUE,
    eqb = 1.25,
    R = 199
  ))
  expect_error(log_TOST(
    data = sleep2,
    extra ~ group,
    paired = TRUE,
    eqb = 1.25,
    R = 199
  ))

  expect_error(boot_log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.25,
    null = -1,
    R = 199
  ))
  expect_error(log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.25,
    null = -1,
    R = 199
  ))

  expect_error(boot_log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.25,
    alpha = -1,
    R = 199
  ))
  expect_error(log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.25,
    alpha = -1,
    R = 199
  ))

  test1 = boot_log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.25,
    R = 199
  )

  test1_t = log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.25
  )

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)

  test1 = boot_log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    hypothesis = "MET",
    eqb = 1.25,
    R = 199
  )

  test1_t = log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    hypothesis = "MET",
    eqb = 1.25
  )

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)

  test1 = boot_log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    hypothesis = "MET",
    eqb = 1.05,
    R = 199
  )

  test1_t = log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    hypothesis = "MET",
    eqb = 1.05
  )

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)


  test1 = boot_log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.9,
    R = 199
  )

  test1_t = log_TOST(
    data = sleep2,
    sleep ~ group,
    paired = TRUE,
    eqb = 1.9
  )

  expect_equal(test1$effsize$estimate,
               test1_t$effsize$estimate)

})