tests/testthat/testGetSummaryPatterns.R

In AlexJanse/cMut: Identify Clustered Single Nucleotide Mutations and Link Them to Mutational Patterns

context("Check the getSummaryPatterns function")

testGroupClusters <- groupClusters(identifyClusters(testDataSet,20000,linkPatterns = T),
                                   patternIntersect = T,searchClusterPatterns = T)
testSummary <- getSummaryPatterns(testGroupClusters)

test_that("Check if the results match with the expectation",{
  expect_equal(ncol(testSummary),
               3)
  expect_equal(any(testSummary$process == "Unidentified"),
               TRUE)
  expect_equal(testSummary[testSummary$process == "AID","frequency"][[1]],
               28)
  expect_equal(testSummary[testSummary$process == "A3G","frequency"][[1]],
               2)
  expect_equal(testSummary[testSummary$process == "PolEta","frequency"][[1]],
               2)
  expect_equal(testSummary[testSummary$process == "PolZeta","frequency"][[1]],
               6)
  expect_equal(testSummary[testSummary$process == "PolZeta.endOnly","frequency"][[1]],
               12)
  expect_equal(nrow(testSummary),
               10)
})