tests/testthat/test-eztoptab.R

In jdreyf/ezlimma: Streamlines and extends limma package

context('eztoptab')

test_that("adjust.method none", {
  eztt.none <- ezlimma:::eztoptab(fit, adjust.method = 'none', prefix = 'abc')
  expect_false(any(grepl('adj\\.P\\.Val', colnames(eztt.none))))
  expect_true(all(grepl('abc', colnames(eztt.none))))
  expect_equal(eztt.none$abc.p, eztt.none$abc.none)
})

test_that("adjust.method holm", {
  eztt.holm <- ezlimma:::eztoptab(fit, adjust.method = 'holm', prefix = 'a-b')
  expect_false(any(grepl('adj\\.P\\.Val', colnames(eztt.holm))))
  expect_true(all(grepl('a-b', colnames(eztt.holm), fixed = TRUE)))
  expect_true(all(eztt.holm$`a-b.p` <= eztt.holm$`a-b.holm`))
})

test_that("z", {
  eztt <- ezlimma:::eztoptab(fit, prefix = 'abc', cols=c("z", "t"))
  expect_true( all(abs(eztt$abc.z) <= abs(eztt$abc.t)) )
  
  eztt2 <- ezlimma:::eztoptab(fit, prefix = 'abc', cols=c("z"), adjust.method="BH", coef=NULL)
  expect_false("abc.t" %in% colnames(eztt2))
  
  eztt3 <- ezlimma:::eztoptab(fit, prefix = 'abc', cols=c("z", "logFC"), adjust.method="BH", coef=NULL)
  expect_false("abc.t" %in% colnames(eztt3))
  expect_gte(ncol(eztt3), 2)
})