tests/testthat/test-da-with-lmm.R

In dominanceanalysis: Dominance Analysis

context("Dominance analysis with lmer")

library(lme4)

lmer.npk.1<-lmer(yield~N+P+K+(1|block),npk)
lmer.npk.0<-lmer(yield~1+(1|block),npk)
da.lmer<-dominanceAnalysis(lmer.npk.1,null.model=lmer.npk.0)

test_that("contribution by level are correct", {
  cbyl<-contributionByLevel(da.lmer,fit.functions = "rb.r2.1")
  expect_equal(sum(cbyl$rb.r2.1$N>cbyl$rb.r2.1$P),3L)
  expect_equal(sum(cbyl$rb.r2.1$K>cbyl$rb.r2.1$P),3L)
})

test_that("contribution.average are correct",{
  ca<-averageContribution(da.lmer,fit.functions = "rb.r2.1")
  expect_equal(ca$rb.r2.1,c(N=0.34,P=-0.03, K=.15),tolerance=0.005)
})


performance.available<-requireNamespace("performance", quietly=TRUE)

if(performance.available) {
  test_that("if performance package is available, related functions are incorporated", {
    ca<-averageContribution(da.lmer)
    expect_true(all(c("n.marg","n.cond") %in% names(ca)))
  })

}