tests/testthat/test_contrasts.R

In dmenne/dmisc2: Dieter Menne's Misc Functions

context("Reading Excel contrast tables")
excelfile = system.file("extdata", "contrasts.xlsx", package = "Dmisc2")
cname = "peripostinterval"

test_that("Excel test file exists", {
  expect_true(file.exists(excelfile))
})

test_that("Reading valid file with readContrast returns contrast table", {
  ct = readContrasts(cname,excelfile)
  expect_equal(nrow(ct),6)
  expect_equal(ncol(ct),9)
  expect_equal(sum(ct[,-1]) ,22)
})

test_that("Reading invalid file name throws", {
  expect_error(readContrasts(cname,str_c(excelfile,"xxxx")),"not found")
  expect_error(getContrasts(cname,str_c(excelfile,"xxxx")),"not found")
})

test_that("Reading invalid range name throws", {
  expect_error(readContrasts("blub",excelfile),"Sheet blub not found")
  expect_error(getContrasts("blub",excelfile),"Sheet blub not found")
})

test_that("Reading valid file with getContrast returns contrast table with attributes", {
  ct = getContrasts(cname,excelfile)
  expect_equal(attr(ct,"varnames"),c("peri","post","interval"))
  vars = attr(ct,"vars")
  expect_equal(levels(vars$peri),c("Tea","Wine"))
  expect_equal(levels(vars$post),c("Kirsch","Water"))
  expect_equal(levels(vars$interval),c("Post","Pre"))
})


test_that("Using gmodels::estimable gives contrasts and confidence intervals", {
  options(digits = 3)
  set.seed(4711)
  d = expand.grid(subject = LETTERS[1:8],
                  peri = c("Wine","Tea"),
                  post = c("Water","Kirsch"),
                  interval = c("Pre","Post"))
  d$vol = round(rnorm(nrow(d),10,2),1)
  d.lme = nlme::lme(vol ~ interval + peri + post + peri:interval + interval:post,
              data = d, random = ~1|subject)
  ct = getContrasts(cname,excelfile)
  est = estimable(d.lme,ct,conf.int = 0.95)
  expect_equal(nrow(est), 8)
  expect_equal(ncol(est), 7)
  expect_equal(round(sum(est)), 765)
})