diff options
author | Ken Kellner <ken@kenkellner.com> | 2023-12-10 20:26:24 -0500 |
---|---|---|
committer | Ken Kellner <ken@kenkellner.com> | 2023-12-10 20:26:24 -0500 |
commit | 9d26eb95b11731ccdf9828d94d9681694df2e003 (patch) | |
tree | a3aa03acddb37814243d06add04920684cc749a2 | |
parent | 68cca855fdbf51906a2c255e24b5f0df67434f97 (diff) |
Fix tests
-rw-r--r-- | R/mixedModelTools.R | 4 | ||||
-rw-r--r-- | tests/testthat/test_mixed_models.R | 107 |
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", { |