tests/testthat/test-matlm.R

In variani/matlm: Linear models in matrix form

context("matlm")

test_that("matlm: basics", {
  N <- 100
  M <- 5
  
  simdat <- matlm_sim_randpred(seed = 10, N = N, M = M) 
  
  assoc <- matlm(simdat$form, simdat$dat, pred = simdat$pred)
  assoc_mod <- matlm_mod(simdat$form, simdat$dat, pred = simdat$pred)
  
  expect_equal(assoc_mod$tab$pval, assoc$tab$pval, tol = 1e-3)
})

test_that("matlm: interaction", {
  N <- 500
  M <- 5
  
  simdat <- matlm_sim_randpred(seed = 1, N = N, M = M) 

  assoc <- matlm(simdat$form_int, simdat$dat, pred = simdat$pred, int = "d")
  assoc_mod <- matlm_mod(simdat$form_int, simdat$dat, pred = simdat$pred, int = "d")
  
  expect_equal(assoc_mod$tab$pval, assoc$tab$pval, tol = 1e-3)
})