tests/testthat/test-predict.R

In sommer: Solving Mixed Model Equations in R

# test_that("Interactions with * work properly", {
#     data(DT_yatesoats)
#     DT <- DT_yatesoats
# 
#     m1 <- mmer(fixed=Y ~ V+N+V:N, random = ~ B + B:MP, rcov=~units, data = DT, verbose = F)
#     m2 <- mmer(fixed=Y ~ V*N, random = ~ B + B:MP, rcov=~units, data = DT, verbose = F)
# 
#     Dt <- m2$Dtable; Dt
#     # first fixed effect include and average
#     Dt[1,"include"] = TRUE
#     Dt[1,"average"] = TRUE
#     # second fixed effect include and average
#     Dt[2,"include"] = TRUE
#     Dt[2,"average"] = TRUE
#     # third fixed effect include and average
#     Dt[3,"include"] = TRUE
#     Dt[3,"average"] = TRUE
#     Dt[4,"include"] = TRUE
#     Dt[4,"average"] = TRUE
#     Dt
# 
#     p3 <- predict.mmer(object=m1, D = "V:N", verbose = F, Dtable = Dt)
#     p4 <- predict.mmer(object=m2, D = "V:N", verbose = F, Dtable = Dt)
#     # expect_equal(p1$pvals$predicted.value, p2$pvals$predicted.value)
#     # expect_equal(p1$hypertable$average, p2$hypertable$average)
#     expect_equal(p3$pvals$predicted.value, p4$pvals$predicted.value)
#     expect_equal(p3$hypertable$average, p4$hypertable$average)
# })