tests/testthat/test_noether.R

In happma/WMWssp: Wilcoxon-Mann-Whitney Sample Size Planning

context("Noether Formula")

# data
x <- c(315,375,356,374,412,418,445,403,431,410,391,475,379)
y <- x - 20
x_ties <- c(315, 315, 320, 234, 234, 439)

# calculate sample size, true result
result_noetherN <- 114.74467
result_noetherN2 <- 114.0657
result_noetherN3 <- 105.1826


test_that("function Noether", {
  expect_equivalent(WMWssp::WMWssp_noether(alpha = 0.05, power = 0.8, t =1/2, p = 0.349)$result[4, ], result_noetherN, tolerance=1e-4)
  expect_equivalent(WMWssp::WMWssp_noether(alpha = 0.05, power = 0.8, t = 1/2, p = 0.349, x = x, ties = TRUE)$result[4, ], result_noetherN2, tolerance=1e-4)
  expect_equivalent(WMWssp::WMWssp_noether(alpha = 0.05, power = 0.8, t = 1/2, p = 0.349, x = x_ties, ties = TRUE)$result[4, ], result_noetherN3, tolerance=1e-4)
  expect_warning(WMWssp::WMWssp_noether(alpha = 0.05, power = 0.8, t = 1/2, p = 0.349, x = x, ties = FALSE))
})