tests/testthat/test-varrankfile.R

In HegemanLab/w4mclassfilter: W4M Class Filter

# test w4mclassfilter::w4m_filter_imputation

library(base)
library(testthat)
library(w4mclassfilter)

test_that("variance by rank or file test", {
  m <- matrix(
    c(
      1, 2, 3,
      5, 7, 11,
      13, 17, 19
    )
  , nrow = 3
  , ncol = 3
  , byrow = TRUE
  )
  rowvars <- w4m__var_by_rank_or_file(m = m, dim = 1)
  expecteds <- c(var(c(1, 2, 3)), var(c(5, 7, 11)), var(c(13, 17, 19)))
  all.equal(rowvars, expecteds)
  expect_equivalent(rowvars, expecteds, info = "variance by row validation")

  colvars <- w4m__var_by_rank_or_file(m = m, dim = 2)
  expecteds <- c(var(c(1, 5, 13)), var(c(2, 7, 17)), var(c(3, 11, 19)))
  all.equal(colvars, expecteds)
  expect_equivalent(colvars, expecteds, info = "variance by column validation")
})

test_that("variance one-column", {
  m <- matrix(
    c(
      1,
      5,
      13
    )
  , nrow = 3
  , ncol = 1
  , byrow = TRUE
  )
  rowvars <- w4m__var_by_rank_or_file(m = m, dim = 1)
  expecteds <- c(NaN, NaN, NaN)
  all.equal(rowvars, expecteds)
  expect_equivalent(rowvars, expecteds, info = "variance by row validation")

  colvars <- w4m__var_by_rank_or_file(m = m, dim = 2)
  expecteds <- var(c(1, 5, 13))
  all.equal(colvars, expecteds)
  expect_equivalent(colvars, expecteds, info = "variance by column validation")
})

test_that("variance one-row", {
  m <- matrix(
    c( 1, 5, 13 )
  , nrow = 1
  , ncol = 3
  , byrow = FALSE
  )
  colvars <- w4m__var_by_rank_or_file(m = m, dim = 2)
  expecteds <- c(NaN, NaN, NaN)
  all.equal(colvars, expecteds)
  expect_equivalent(colvars, expecteds, info = "variance by column validation")

  rowvars <- w4m__var_by_rank_or_file(m = m, dim = 1)
  expecteds <- var(c(1, 5, 13))
  all.equal(rowvars, expecteds)
  expect_equivalent(rowvars, expecteds, info = "variance by row validation")
})