tests/testthat/test-wilcoxExt.R

In Wancen/airpart: Differential cell-type-specific allelic imbalance

library(S4Vectors)

test_that("wrong inputs", {

  # missed input arrays
  set.seed(2021)
  sce <- makeSimulatedData()
  expect_error(wilcoxExt(sce))

  # require genecluster identified
  sce <- preprocess(sce)
  expect_error(wilcoxExt(sce), "No gene cluster number")

  # no cluster column
  expect_error(wilcoxExt(sce, genecluster = 1))

  # wrong cell state column names
  names(colData(sce)) <- "ct"
  expect_error(wilcoxExt(sce, genecluster = 1))
})

test_that("basic wilcox analyses", {
  set.seed(2021)
  sce <- makeSimulatedData()
  sce <- preprocess(sce)
  sce <- geneCluster(sce, G = seq_len(4), plot = FALSE)
  sce_sub <- wilcoxExt(sce, genecluster = 1)
  expect_true(is.numeric(metadata(sce_sub)$threshold))
})