tests/testthat/test_plotCounts.R

In EngeLab/CIMseq: CIMseq is a software suite for analysis of scRNAseq data that has been produced using the cellular interactions by multiplet sequencing method

context("spPlotCounts")

##run test plotCountsData
test_that("check that plotCountsData outputs the expected result", {
  #This is largley tested in the estimateCells function
  
  ###TEST1####
  nc1 <- ncol(getData(CIMseqSinglets_test, "counts"))
  nc2 <- ncol(getData(CIMseqMultiplets_test, "counts"))
  nc <- nc1 + nc2
  
  #run function
  output <- plotCountsData(CIMseqSinglets_test, CIMseqMultiplets_test)
  
  #test
  expect_is(output, c("tbl_df", "tbl", "data.frame"))
  expect_type(output, "list")
  expect_true(ncol(output) == 7)
  expect_true(nrow(output) == nc)
})

##run test plotCountsERCC
test_that("check that plotCountsERCC outputs the expected result", {
  
  #run function
  output <- plotCountsMarkers(CIMseqSinglets_test, CIMseqMultiplets_test)
  
  #test
  expect_type(output, "list")
  expect_is(output, c("gg", "ggplot"))
})

##run test plotCountsMarkers
test_that("check that plotCountsMarkers outputs the expected result", {
    
  #run function
  markers <- c("RNR2", "NME7")
  output <- plotCountsMarkers(CIMseqSinglets_test, CIMseqMultiplets_test, markers)
  
  #test
  expect_type(output, "list")
  expect_is(output, c("gg", "ggplot"))
})

##run test plotUnsupervisedClass
test_that("check that plotUnsupervisedClass outputs the expected result", {
  
  #run function
  output <- plotUnsupervisedClass(CIMseqSinglets_test, CIMseqMultiplets_test)
  
  #test
  expect_type(output, "list")
  expect_is(output, c("gg", "ggplot"))
})

##run test plotUnsupervisedMarkers
test_that("check that plotUnsupervisedMarkers outputs the expected result", {
  
  #run function
  output <- plotUnsupervisedMarkers(
    CIMseqSinglets_test, CIMseqMultiplets_test, markers = "ACTB"
  )
  
  #test
  expect_type(output, "list")
  expect_is(output, c("gg", "ggplot"))
  
  #run function
  output <- plotUnsupervisedMarkers(
    CIMseqSinglets_test, CIMseqMultiplets_test,
    markers = c("ACTB", "GAPDH")
  )
  
  #test
  expect_type(output, "list")
  expect_is(output, c("gg", "ggplot"))
})