test_that("RegEnRF output is correct", {
  x <- matrix(rnorm(100 * 20), 100, 20)
  y <- rnorm(100)
  mod <- RegEnRF(x, y, 0.1)

  expect_s3_class(mod, "RegEnRF")
})