tests/testthat/test-do_CorrelationHeatmap.R

In enblacar/SCpubr: Generate Publication Ready Visualizations of Single Cell Transcriptomics Data

if (base::isFALSE(dep_check[["do_CorrelationHeatmap"]])){

  testthat::test_that("do_CorrelationHeatmap: CRAN essentials", {
    sample$orig.ident <- ifelse(sample$seurat_clusters %in% c("1", "2"), "A", "B")

    genes <- list("A" = Seurat::VariableFeatures(sample)[1:5],
                  "B" = Seurat::VariableFeatures(sample)[6:10],
                  "C" = Seurat::VariableFeatures(sample)[11:15])

    p <- SCpubr::do_CorrelationHeatmap(sample = sample, legend.position = "top")
    testthat::expect_true("ggplot" %in% class(p))
  })

  testthat::test_that("do_CorrelationHeatmap: PASS - normal", {

    testthat::skip_on_cran()
    sample$orig.ident <- ifelse(sample$seurat_clusters %in% c("1", "2"), "A", "B")

    p <- SCpubr::do_CorrelationHeatmap(sample = sample, legend.position = "top", group.by = "orig.ident")
    testthat::expect_true("ggplot" %in% class(p))

    p <- SCpubr::do_CorrelationHeatmap(sample = sample, legend.position = "right")
    testthat::expect_true("ggplot" %in% class(p))
    
    p <- SCpubr::do_CorrelationHeatmap(sample = sample, legend.position = "right", cluster = TRUE)
    testthat::expect_true("ggplot" %in% class(p))
    
    p <- SCpubr::do_CorrelationHeatmap(sample = sample, legend.position = "right", cluster = FALSE)
    testthat::expect_true("ggplot" %in% class(p))
  })
  
  testthat::test_that("do_CorrelationHeatmap: PASS - jaccard", {
    
    testthat::skip_on_cran()
    sample$orig.ident <- ifelse(sample$seurat_clusters %in% c("1", "2"), "A", "B")
    
    genes <- list("A" = rownames(sample)[1:5],
                  "B" = rownames(sample)[3:8],
                  "C" = rownames(sample)[5:13])
    
    p <- SCpubr::do_CorrelationHeatmap(input_gene_list = genes, mode = "jaccard", legend.position = "top", cluster = FALSE, use_viridis = TRUE)
    testthat::expect_true("ggplot" %in% class(p))
    
    p <- SCpubr::do_CorrelationHeatmap(input_gene_list = genes, mode = "jaccard", legend.position = "top", cluster = TRUE, use_viridis = FALSE)
    testthat::expect_true("ggplot" %in% class(p))
    
    p <- SCpubr::do_CorrelationHeatmap(input_gene_list = genes, mode = "jaccard", legend.position = "top", remove.diagonal = TRUE)
    testthat::expect_true("ggplot" %in% class(p))
    
    p <- SCpubr::do_CorrelationHeatmap(input_gene_list = genes, mode = "jaccard", legend.position = "top", remove.diagonal = FALSE)
    testthat::expect_true("ggplot" %in% class(p))
    
    p <- SCpubr::do_CorrelationHeatmap(input_gene_list = genes, mode = "jaccard", legend.position = "right")
    testthat::expect_true("ggplot" %in% class(p))
  })

  testthat::test_that("do_CorrelationHeatmap: PASS - group.by", {
    testthat::skip_on_cran()

    genes <- list("A" = Seurat::VariableFeatures(sample)[1:5],
                  "B" = Seurat::VariableFeatures(sample)[6:10],
                  "C" = Seurat::VariableFeatures(sample)[11:15])

    p <- SCpubr::do_CorrelationHeatmap(sample = sample,
                                    group.by = "seurat_clusters")
    testthat::expect_true("ggplot" %in% class(p))
  })

  testthat::test_that("do_CorrelationHeatmap: PASS - group.by - rotate axis labels", {
    testthat::skip_on_cran()

    genes <- list("A" = Seurat::VariableFeatures(sample)[1:5],
                  "B" = Seurat::VariableFeatures(sample)[6:10],
                  "C" = Seurat::VariableFeatures(sample)[11:15])

    p <- SCpubr::do_CorrelationHeatmap(sample = sample,
                                    group.by = "seurat_clusters",
                                    axis.text.x.angle = 0)
    testthat::expect_true("ggplot" %in% class(p))
  })
}