tests/testthat/test-rank_group.R

In LieberInstitute/TREG: Tools for finding Total RNA Expression Genes in single nucleus RNA-seq data

## Input checks
test_that("Bad Inputs Throw Error", {
    bad_input <- data.frame(row.names = c("g1", "g2", "g3"), c1 = c(0, 1, 2))
    ## wrong input format
    expect_error(rank_group(sce = bad_input))

    ## bad column name
    expect_error(rank_group(sce = sce_zero_test, group_col = "NOT_THERE"))
})


test_that("Empty Levels are Dropped", {
    sce_zero_test$cellType <- factor(sce_zero_test$cellType, levels = c("A", "B", "C"))
    n_non_empty <- sum(table(sce_zero_test$cellType) != 0)

    expect_warning(rank_group(sce_zero_test))
    # expect_equal(length(rank_group(sce_zero_test)), n_non_empty)
})