tests/testthat/test-lmm_f.R

In aosmith16/simpow: Simulation-based power analysis

test_that("error if means vector length doesn't match num trt", {
    expect_error(lmm_f(ntrt = 2, trtmeans = c(1, 2, 3)),
                 "You must provide a mean")
})

test_that("error if trt names vector length doesn't match num trt", {
    expect_error(lmm_f(ntrt = 2, trtmeans = c(1, 2),
                       trtnames = c("a")),
                 "You must provide a name")
})

test_that("error if nonconstant variance sd vector length doesn't match num trt", {
    expect_error(lmm_f(trtmeans = c(1, 2),
                       sd_resid = c(1, 2, 3), sd_eq = FALSE),
                 "You are allowing nonconstant variance")
})

test_that("mean to have constant variance", {
    expect_error(lmm_f(trtmeans = c(1, 2),
                       sd_block = 2,
                       sd_resid = c(1, 100)),
                 "You are allowing constant variance")
})

test_that("at least 3 blocking groups", {
    expect_error(lmm_f(ntrt = 2, trtmeans = c(1, 2),
                       nblock = 2, sd_resid = 2),
                 "The number of groups")
})

test_that("data length matches design parameters", {
    ntrt = 2
    nblock = 3
    nrep = 2
    res = lmm_f(test = "none",
                ntrt = ntrt,
                trtmeans = c(1, 25),
                nblock = nblock,
                nrep = nrep,
                sd_block = 2,
                sd_resid = 4,
                keep_data = TRUE)
    expect_equal(nrow(res$data[[1]]), ntrt*nblock*nrep)
})

test_that("data is ordered by blocks then treatments vector nrep times", {

    res1 = lmm_f(test = "none",
                ntrt = 2,
                trtmeans = c(1, 25),
                nblock = 5,
                nrep = 1,
                sd_block = 2,
                sd_resid = 4,
                keep_data = TRUE)$data[[1]]

    trt1 = rep(LETTERS[1:2], times = 5*1)
    block1 = rep(as.character(1:5), each = 2)

    expect_equal(res1[, 1], trt1)
    expect_equal(res1[, 2], block1)

    res2 = lmm_f(test = "none",
                 ntrt = 2,
                 trtmeans = c(1, 25),
                 nblock = 5,
                 nrep = 3,
                 sd_block = 2,
                 sd_resid = 4,
                 keep_data = TRUE)$data[[1]]

    trt2 = rep(LETTERS[1:2], times = 5*3)
    block2 = rep(as.character(1:5), each = 2*3)

    expect_equal(res2[, 1], trt2)
    expect_equal(res2[, 2], block2)
})