aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorKen Kellner <ken@kenkellner.com>2023-12-10 20:26:24 -0500
committerKen Kellner <ken@kenkellner.com>2023-12-10 20:26:24 -0500
commit9d26eb95b11731ccdf9828d94d9681694df2e003 (patch)
treea3aa03acddb37814243d06add04920684cc749a2
parent68cca855fdbf51906a2c255e24b5f0df67434f97 (diff)
Fix tests
-rw-r--r--R/mixedModelTools.R4
-rw-r--r--tests/testthat/test_mixed_models.R107
2 files changed, 57 insertions, 54 deletions
diff --git a/R/mixedModelTools.R b/R/mixedModelTools.R
index d8dad59..985d06e 100644
--- a/R/mixedModelTools.R
+++ b/R/mixedModelTools.R
@@ -95,8 +95,8 @@ terms_in_bar <- function(bars, RHS=FALSE){
bars_sub <- bars[[2]][[2]]
if(RHS) bars_sub <- bars[[2]][[3]]
form <- formula(substitute(~X, list(X=bars_sub)))
- trms <- attr(terms(form), "term.labels")
- int <- attr(terms(form), "intercept")
+ trms <- attr(stats::terms(form), "term.labels")
+ int <- attr(stats::terms(form), "intercept")
if(int == 1 & !RHS) trms <- c("1", trms)
trms
}
diff --git a/tests/testthat/test_mixed_models.R b/tests/testthat/test_mixed_models.R
index f0f7d10..4cf3b88 100644
--- a/tests/testthat/test_mixed_models.R
+++ b/tests/testthat/test_mixed_models.R
@@ -2,68 +2,71 @@ context("mixed model tools")
test_that("get_reTrms matches lme4::mkReTrms", {
- skip_if(!requireNamespace("lme4", quietly=TRUE),
- "lme4 package unavailable")
-
+ #skip_if(!requireNamespace("lme4", quietly=TRUE),
+ # "lme4 package unavailable")
+ set.seed(123)
dat <- data.frame(x = rnorm(20), y = rnorm(20), z = factor(sample(letters[1:3], 20, replace=T)),
group = factor(sample(letters[4:6], 20, replace=T)),
id = factor(sample(letters[7:9], 20, replace=T)))
-
+
+ load('lme4_output.Rdata')
form1 <- ~x + (1|group)
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ #l1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
+ r1 <- get_reTrms(form1, dat)
+ expect_identical(r1$Z, Matrix::t(l1$Zt))
+ expect_identical(r1$cnms, l1$cnms)
+ attributes(l1$flist)$assign <- NULL
+ expect_identical(r1$flist, l1$flist)
- form1 <- ~x + (x||group)
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ form2 <- ~x + (x||group)
+ #l2 <- lme4::mkReTrms(lme4::findbars(form2), dat)
+ r2 <- get_reTrms(form2, dat)
+ expect_identical(r2$Z, Matrix::t(l2$Zt))
+ expect_identical(r2$cnms, l2$cnms)
+ attributes(l2$flist)$assign <- NULL
+ expect_identical(r2$flist, l2$flist)
- form1 <- ~x + (x||group) + (1|id)
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ form3 <- ~x + (x||group) + (1|id)
+ #l3 <- lme4::mkReTrms(lme4::findbars(form3), dat)
+ r3 <- get_reTrms(form3, dat)
+ expect_identical(r3$Z, Matrix::t(l3$Zt))
+ expect_identical(r3$cnms, l3$cnms)
+ attributes(l3$flist)$assign <- NULL
+ expect_identical(r3$flist, l3$flist)
- form1 <- ~x + (x||group) + (y||id)
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ form4 <- ~x + (x||group) + (y||id)
+ #l4 <- lme4::mkReTrms(lme4::findbars(form4), dat)
+ r4 <- get_reTrms(form4, dat)
+ expect_identical(r4$Z, Matrix::t(l4$Zt))
+ expect_identical(r4$cnms, l4$cnms)
+ attributes(l4$flist)$assign <- NULL
+ expect_identical(r4$flist, l4$flist)
- form1 <- ~x + (x*y||group) + (y||id)
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ form5 <- ~x + (x*y||group) + (y||id)
+ #l5 <- lme4::mkReTrms(lme4::findbars(form5), dat)
+ r5 <- get_reTrms(form5, dat)
+ expect_identical(r5$Z, Matrix::t(l5$Zt))
+ expect_identical(r5$cnms, l5$cnms)
+ attributes(l5$flist)$assign <- NULL
+ expect_identical(r5$flist, l5$flist)
- form1 <- ~(1|group)
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ form6 <- ~(1|group)
+ #l6 <- lme4::mkReTrms(lme4::findbars(form6), dat)
+ r6 <- get_reTrms(form6, dat)
+ expect_identical(r6$Z, Matrix::t(l6$Zt))
+ expect_identical(r6$cnms, l6$cnms)
+ attributes(l6$flist)$assign <- NULL
+ expect_identical(r6$flist, l6$flist)
- form1 <- ~(1|group) + x
- r1 <- lme4::mkReTrms(lme4::findbars(form1), dat)
- r2 <- get_reTrms(form1, dat)
- expect_identical(r2$Z, Matrix::t(r1$Zt))
- expect_identical(r1$cnms, r2$cnms)
- attributes(r1$flist)$assign <- NULL
- expect_identical(r1$flist, r2$flist)
+ form7 <- ~(1|group) + x
+ #l7 <- lme4::mkReTrms(lme4::findbars(form7), dat)
+ r7 <- get_reTrms(form7, dat)
+ expect_identical(r7$Z, Matrix::t(l7$Zt))
+ expect_identical(r7$cnms, l7$cnms)
+ attributes(l7$flist)$assign <- NULL
+ expect_identical(r7$flist, l7$flist)
+
+ #save(l1,l2,l3,l4,l5,l6,l7, file='lme4_output.Rdata')
})
test_that("get_reTrms errors correctly", {