tests/testthat/test_runGradientDescent.R

In BEclear: Correction of batch effects in DNA methylation data

testthat::context("Testing the gradient descent for a multiple epochs")

testthat::test_that("quadratic D matrix, only integer values given", {
  res1 <- runGradientDescent(
    L = matrix(c(1, 2)), R = t(matrix(c(3, 4))),
    D = matrix(c(1, 1, 1, 1), nrow = 2),
    lambda = 0.2, is = c(1, 2, 1, 2), js = c(1, 1, 2, 2),
    gamma = 0.01, r = 1, epochs = 5, block = 1
  )

  testthat::expect_equal(res1$L, matrix(c(0.347, 0.393)), tolerance = .01)
  testthat::expect_equal(res1$R, t(matrix(c(2.61, 3.43))), tolerance = .01)
})

testthat::test_that("long D matrix, only integer values given", {
  res1 <- runGradientDescent(
    L = matrix(c(5, 6)), R = t(matrix(c(7, 4, 8))),
    D = matrix(c(7, 5, 3, 1, 3, 4), nrow = 2),
    lambda = 0.2, is = c(1, 2, 1, 2, 1, 2), js = c(1, 1, 2, 2, 3, 3),
    gamma = 0.01, r = 1, epochs = 5, block = 1
  )

  testthat::expect_equal(res1$L, matrix(c(-6.12, -7.76)), tolerance = .01)
  testthat::expect_equal(res1$R, t(matrix(c(-0.831, -0.266, -0.504))), tolerance = .01)
})

testthat::test_that("wide D matrix, only integer values given", {
  res1 <- runGradientDescent(
    L = matrix(c(5, 1, 2)), R = t(matrix(c(3, 3))),
    D = matrix(c(5, 6, 10, 11, 3, 1), nrow = 3),
    lambda = 0.2, is = c(1, 2, 3, 1, 2, 3), js = c(1, 1, 1, 2, 2, 2),
    gamma = 0.01, r = 1, epochs = 5, block = 1
  )

  testthat::expect_equal(res1$L, matrix(c(3.73, 1.68, 2.20)), tolerance = .01)
  testthat::expect_equal(res1$R, t(matrix(c(2.26, 2.26))), tolerance = .01)
})