tests/testthat/test-delta_ind_t.R

In MOTE: Effect Size and Confidence Interval Calculator

test_that("delta_ind_t returns correctly structured list", {
  res <- delta_ind_t(
    m1  = 17.75,
    m2  = 23,
    sd1 = 3.30,
    sd2 = 2.16,
    n1  = 4,
    n2  = 4,
    a   = 0.05
  )

  expect_type(res, "list")

  # Legacy names
  expect_true(all(c(
    "d", "dlow", "dhigh",
    "M1", "sd1", "se1", "M1low", "M1high",
    "M2", "sd2", "se2", "M2low", "M2high",
    "spooled", "sepooled",
    "n1", "n2", "df",
    "t", "p",
    "estimate", "statistic"
  ) %in% names(res)))

  # Snake_case aliases
  expect_true(all(c(
    "d_value", "d_lower_limit", "d_upper_limit",
    "mean1_value", "mean1_lower", "mean1_upper",
    "mean2_value", "mean2_lower", "mean2_upper",
    "sample_sd1", "sample_sd2",
    "sample_se1", "sample_se2",
    "pooled_sd", "pooled_se",
    "sample_size1", "sample_size2",
    "t_value", "p_value"
  ) %in% names(res)))
})


test_that("delta_ind_t computes correct d and CI for known example", {
  m1  <- 17.75
  m2  <- 23
  sd1 <- 3.30
  sd2 <- 2.16
  n1  <- 4
  n2  <- 4
  a   <- 0.05

  res <- delta_ind_t(m1, m2, sd1, sd2, n1, n2, a)

  # Expected d (direction matches function: (m1 - m2) / sd1)
  d_expected <- (m1 - m2) / sd1
  expect_equal(res$d, d_expected, tolerance = 1e-6)

  # Expected dlow and dhigh from known output
  expect_equal(res$dlow,  -3.56801,  tolerance = 1e-5)
  expect_equal(res$dhigh, -0.1021864, tolerance = 1e-6)

  # CI should bracket d
  expect_true(res$dlow <= res$d)
  expect_true(res$dhigh >= res$d)

  # P-value should be between 0 and 1
  expect_true(res$p >= 0 && res$p <= 1)
})


test_that("APA-style strings are generated correctly", {
  res <- delta_ind_t(17.75, 23, 3.30, 2.16, 4, 4, .05)

  expect_type(res$estimate, "character")
  expect_true(grepl("d", res$estimate))
  expect_true(grepl("CI", res$estimate))

  expect_type(res$statistic, "character")
  expect_true(grepl("t", res$statistic))
})