tests/testthat/test-f.R

In olsrr: Tools for Building OLS Regression Models

context("f_test")

model <- lm(mpg ~ disp + hp + wt + drat + qsec, data = mtcars)

test_that("when fitted.values == TRUE, fitted values from the regression\n\tare used for the test", {
  k <- ols_test_f(model)

  expect_equal(round(k$f, 3), 1.239)
  expect_equal(round(k$p, 3), 0.275)
  expect_equal(k$numdf, 1)
  expect_equal(k$dendf, 30)

  expect_true(k$fv)
  expect_false(k$rhs)

  expect_null(k$vars)

  expect_match(k$resp, "mpg")
  expect_equivalent(k$preds, c("disp", "hp", "wt", "drat", "qsec"))
})


test_that("when fitted_values == TRUE and rhs == TRUE, predictors from the\n\tmodel are used for the test", {
  k <- ols_test_f(model, fitted_values = TRUE, rhs = TRUE)

  expect_equal(round(k$f, 3), 0.444)
  expect_equal(round(k$p, 3), 0.814)
  expect_equal(k$numdf, 5)
  expect_equal(k$dendf, 26)

  expect_false(k$fv)
  expect_true(k$rhs)

  expect_null(k$vars)

  expect_match(k$resp, "mpg")
  expect_equivalent(k$preds, c("disp", "hp", "wt", "drat", "qsec"))
})


test_that("when vars != NULL, variables specified from the are\n\tused for the test", {
  k <- ols_test_f(model, vars = c("disp", "hp"))

  expect_equal(round(k$f, 3), 0.453)
  expect_equal(round(k$p, 3), 0.64)
  expect_equal(k$numdf, 2)
  expect_equal(k$dendf, 29)

  expect_false(k$fv)
  expect_false(k$rhs)

  expect_equivalent(k$vars, c("disp", "hp"))
  expect_match(k$resp, "mpg")
  expect_equivalent(k$preds, c("disp", "hp", "wt", "drat", "qsec"))
})


test_that("when vars != NULL and rhs == TRUE, predictors in the model are\n\tused for the test", {
  k <- ols_test_f(model, rhs = TRUE, vars = c("disp", "hp"))

  expect_equal(round(k$f, 3), 0.444)
  expect_equal(round(k$p, 3), 0.814)
  expect_equal(k$numdf, 5)
  expect_equal(k$dendf, 26)

  expect_false(k$fv)
  expect_true(k$rhs)

  expect_equivalent(k$vars, c("disp", "hp"))
  expect_match(k$resp, "mpg")
  expect_equivalent(k$preds, c("disp", "hp", "wt", "drat", "qsec"))
})