tests/testthat/test-model_parameters_std.R

In parameters: Processing of Model Parameters

skip_on_cran()

skip_if_not_installed("effectsize")
data(mtcars)
mtcars$am <- as.factor(mtcars$am)
d <- mtcars
model <- lm(mpg ~ wt * am, data = d)

test_that("model_parameters, standardize-refit", {
  params <- model_parameters(model, standardize = "refit")
  expect_identical(c(nrow(params), ncol(params)), c(4L, 9L))
  expect_equal(params$Coefficient, c(-0.14183, -0.61463, -0.35967, -0.86017), tolerance = 1e-3)
  expect_equal(params$SE, c(0.12207, 0.12755, 0.23542, 0.23454), tolerance = 1e-3)
  expect_equal(params$CI_high, c(0.10821, -0.35336, 0.12257, -0.37973), tolerance = 1e-3)
})

test_that("model_parameters, standardize-posthoc", {
  params <- model_parameters(model, standardize = "posthoc")
  expect_identical(c(nrow(params), ncol(params)), c(4L, 9L))
  expect_equal(params$Std_Coefficient, c(0, -0.61463, 2.46865, -0.87911), tolerance = 1e-3)
  expect_equal(params$SE, c(0, 0.12755, 0.7075, 0.23971), tolerance = 1e-3)
  expect_equal(params$CI_high, c(0, -0.35336, 3.91789, -0.38809), tolerance = 0.1)
})

test_that("model_parameters, standardize-basic", {
  params <- model_parameters(model, standardize = "basic")
  expect_identical(c(nrow(params), ncol(params)), c(4L, 9L))
  expect_equal(params$Std_Coefficient, c(0, -0.61463, 1.23183, -1.11016), tolerance = 1e-3)
  expect_equal(params$SE, c(0, 0.12755, 0.35303, 0.30271), tolerance = 1e-3)
  expect_equal(params$CI_high, c(0, -0.35336, 1.95499, -0.4901), tolerance = 0.1)
})

test_that("model_parameters, standardize-smart", {
  params <- model_parameters(model, standardize = "smart")
  expect_identical(c(nrow(params), ncol(params)), c(4L, 9L))
  expect_equal(params$Std_Coefficient, c(0, -0.61463, 2.41278, -0.85922), tolerance = 1e-3)
  expect_equal(params$SE, c(0, 0.12755, 0.69148, 0.23428), tolerance = 1e-3)
  expect_equal(params$CI_high, c(0, -0.35336, 3.82922, -0.37931), tolerance = 0.1)
})