tests/testthat/test_viz_basis_vectors_seqlogo.R

In snikumbh/archR: Identify Different Architectures of Sequence Elements

context("Test basis vectors visualized as seqlogo")
library(ggplot2)
library(ggseqlogo)


test_that("Given object is matrix", {
  testFeaturesMat <- rnorm(10000) # err
  testPositionLabels <- seq(25)
  expect_error(
    viz_bas_vec_seqlogo(testFeaturesMat, pos_lab = testPositionLabels),
    "not of type matrix"
  )
})

test_that("Handling empty matrix", {
  testFeaturesMat <- matrix()
  testPositionLabels <- seq(25)
  expect_error(
    viz_bas_vec_seqlogo(testFeaturesMat,
      pos_lab = testPositionLabels
    ),
    "Empty"
  )
})

test_that("Position labels inadequate", {
  testFeaturesMat <- matrix(rnorm(10000), nrow = 200)
  testPositionLabels <- seq(20)
  expect_error(
    viz_bas_vec_seqlogo(testFeaturesMat, pos_lab = testPositionLabels),
    "Inadequate"
  )
})

test_that("Position labels over-abundant", {
  testFeaturesMat <- matrix(rnorm(10000), nrow = 200)
  testPositionLabels <- seq(60)
  expect_error(
    viz_bas_vec_seqlogo(testFeaturesMat, pos_lab = testPositionLabels),
    "Overabundant"
  )
})

# test_that("Seqlogo plotting works", {
#   # setting seed enables proper comparison between ggplot objects since we use
#   # rnorm
#   set.seed(11223344)
#   # test variables
#   testPositionLabels <- seq(25)
#   testFeaturesMat <- matrix(rnorm(100), ncol=1)
#   # test plot, include the function directly inside because it returns 
#   nothing.
#   vdiffr::expect_doppelganger("combined heatmap seqlogo plot example",
#                 viz_all_factors_as_seqlogo(featuresMatrix = testFeaturesMat,                                                                              pos_lab = testPositionLabels))
#
# })