tests/testthat/test-groupFisher.R

In YunhuiQi/TestPMD: Post Selection Inference on Results from Penalized Matrix Decomposition (SCCA only)

X <- matrix(rnorm(200), nrow = 10, ncol = 20)
loadings <- seq(0.1, 2, 0.1)
qvalues <- seq(0.05, 1, 0.05)
feature_groups <- c(rep("A", 10), rep("B", 10))
loadings1 <- 0.1
qvalues1 <- 0.1
feature_groups1 <- c("A", "B")
test_that("Can identify correct input", {
  expect_error(groupFisher(data = X, loadings = loadings1, qvalues = qvalues, feature_groups = feature_groups, alpha = 0.05))
  expect_error(groupFisher(data = X, loadings = loadings, qvalues = qvalues1, feature_groups = feature_groups, alpha = 0.05))
  expect_error(groupFisher(data = X, loadings = loadings, qvalues = qvalues, feature_groups = feature_groups1, alpha = 0.05))
  expect_error(groupFisher(data = X, loadings = loadings, qvalues = qvalues, feature_groups = feature_groups, alpha = 2))
  expect_error(groupFisher(data = X, loadings = loadings, qvalues = qvalues, feature_groups = feature_groups, alpha = "2"))
})


test_that("Can output correct values", {
  expect_length(groupFisher(data = X, loadings = loadings, qvalues = qvalues, feature_groups = feature_groups, alpha = 0.05), 9)
  expect_length(groupFisher(data = X, loadings = loadings, qvalues = qvalues, feature_groups = feature_groups, alpha = 0.05, draw = F), 2)


  })