tests/testthat/test_testLibrary.R

In simr: Power Analysis for Generalised Linear Mixed Models by Simulation

context("Test Library")

test_that("tests run", {

    t1 <- doTest(fm1)
    expect_equal(c(t1), 0.0007336556, tolerance=1e-5)
    expect_output(print(t1), "Kenward Roger \\(package pbkrtest\\)")

    t2 <- doTest(fm1, fixed("x", "lr"))
    expect_equal(c(t2), 0.0005147106, tolerance=1e-5)
    expect_output(print(t2), "Likelihood ratio")

    t3 <- doTest(fm1, random())
    expect_equal(c(t3), 0, tolerance=1e-5)
    expect_output(print(t3), "Exact restricted LRT \\(package RLRsim\\)")

    t4 <- doTest(fm1, compare(~ (1|g)))
    expect_equal(c(t4), 0.0005147106, tolerance=1e-5)
    expect_output(print(t4), "Likelihood ratio")

    t5 <- doTest(fm1, fcompare(~ 1))
    expect_equal(c(t5), 0.0005147106, tolerance=1e-5)
    expect_output(print(t5), "Likelihood ratio")

    t6 <- doTest(fm3, rcompare(~ (1|g)))
    expect_equal(c(t6), 0.6029132, tolerance=1e-5)
    expect_output(print(t6), "Likelihood ratio")

    set.seed(333) # pb test is random

    t7 <- suppressWarnings(doTest(fm3, rcompare(~ (1|g), "pb")))
    expect_equal(c(t7), 1/3, tolerance=1e-5)
    expect_output(print(t7), "Parametric bootstrap \\(package pbkrtest\\)")


})

test_that("lmer tests work", {

    tst <- doTest(fm1, fixed("x"))
    expect_equal(c(tst), 0.0007336556, tolerance=1e-8)
    expect_output(print(tst), "Kenward Roger")

    #tst <- doTest(fm1, fixed("x", "z"))
    #expect_equal(c(tst), 0.0005147106, tolerance=1e-5)
    #expect_output(print(tst), "Likelihood ratio")

    tst <- doTest(fm1, fixed("x", "kr"))
    expect_equal(c(tst), 0.0007336556, tolerance=1e-8)
    expect_output(print(tst), "Kenward Roger")

    tst <- doTest(fm1, fixed("x", "sa"))
    expect_equal(c(tst), 0.0007336547, tolerance=1e-8)
    expect_output(print(tst), "Satterthwait")

})