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")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.