tests/testthat/test.data.interim.R

In MikeKSmith/MSToolkit: The MSToolkit Library For Clinical Trial Design

test_that("test.data.interim.1", {
    expect_error(createInterims(10, method = "a"),
                 regexp = "method must be `Sample` or `Proportion`")
    expect_error(createInterims(-10),
                 regexp = "subjects must be positive")
    expect_error(createInterims(5, proportion = ".1,1.2"),
                 regexp = "Proportion greater than 1")
    expect_error(createInterims(5, proportion = ".4,.3"),
                 regexp = "Proportion not increasing")
})

test_that("test.data.interim.2", {
    a <- createInterims(10, proportion = ".1,.3,.6")
    expect_equal(10, nrow(a))
    expect_equal(c("SUBJ", "INTERIM"), names(a))
    expect_equal(1:10, a[, 1])
    expect_error(createInterims(10, proportion = ".1,.3,.6", idCol = ".1"),
                 regexp = "invalid R name")
    expect_error(createInterims(10, proportion = ".1,.3,.6", interimCol = ".1"),
                 regexp = "invalid R name")
    expect_error(createInterims(10, proportion = ".1,.3,.6",
                                interimCol = "XX",
                                idCol = "XX"),
                 regexp = "`idCol` and `interimCol` should be different")
    expect_equal("IDCOL", names(createInterims(10,
                                               proportion = ".1,.3,.6",
                                               idCol = "IDCOL"))[1])
})

test_that("test.data.interim.3", {
    # test that the proportion of subjects allocated matches expectation
    b <- createInterims(1000, proportion = ".1", method = "Sample")
    expect_equal(1:2, sort(unique(b[, 2])))
    expect_true(binom.test(table(b[, 2]),
                           n = 1,
                           p = 0.1)$p.value > 0.05)
})

test_that("test.data.interim.5", {
    a <- createInterims(10)
    expect_true(all(a[, 2] == 1))
})