tests/testthat/test-lda.R

In bimberlabinternal/CellMembrane: A package with high-level wrappers and pipelines for single-cell RNA-seq tools

context("scRNAseq")
library(Seurat)

test_that("LDA works as expected", {
    seuratObj <- suppressWarnings(Seurat::UpdateSeuratObject(readRDS('../testdata/seuratOutput.rds')))

    results <- CellMembrane::DoLdaParameterScan(seuratObj, ntopics = c(5, 10, 15))
    expect_equal(length(results), 3)
    expect_equal(length(results[[1]]$topic_sums), 5)
    expect_equal(length(results[[2]]$topic_sums), 10)
    expect_equal(length(results[[3]]$topic_sums), 15)
})