tests/testthat/test_match_two_sig_sets.R

In mSigTools: Mutational Signature Analysis Tools

test_that("test match_two_sig_sets", {
  ex.sigs <- matrix(c(0.2, 0.8, 0.3, 0.7), nrow = 2)
  colnames(ex.sigs) <- c("ex1", "ex2")
  ref.sigs <- matrix(c(0.21, 0.79, 0.19, 0.81), nrow = 2)
  colnames(ref.sigs) <- c("ref1", "ref2")
  rr <- match_two_sig_sets(ex.sigs, ref.sigs, cutoff = .985)
  expect_equal(
    rr,
    list(
      table = structure(
        list(
          x1 = c("ex1", "ex2"), x2 = c("ref2", "ref1"),
          dist = c(0.000106231316866934, 0.0105052831822057)
        ),
        row.names = 1:2, class = "data.frame"
      ),
      orig.matrix = structure(c(
        0.000110047078473174, 0.0105052831822057,
        0.000106231316866934, 0.0151846812080705
      ),
      .Dim = c(2L, 2L),
      .Dimnames = list(c("ex1", "ex2"), c("ref1", "ref2"))
      ),
      modified.matrix =
        structure(c(
          0.000110047078473174,
          0.0105052831822057,
          0.000106231316866934, 9e+99
        ),
        .Dim = c(2L, 2L),
        .Dimnames = list(c("ex1", "ex2"), c("ref1", "ref2"))
        )
    )
  )
  ex.sigs2 <- cbind(ex.sigs, ex3 = c(0.18, 0.82))
  rr <- match_two_sig_sets(ex.sigs2, ref.sigs, cutoff = .985)
  testthat::expect_equal(
    rr,
    list(
      table = structure(
        list(
          x1 = c("ex1", "ex3"), x2 = c("ref1", "ref2"),
          dist = c(
            0.000110047078473174,
            0.000102493133434356
          )
        ),
        row.names = 1:2,
        class = "data.frame"
      ),
      orig.matrix =
        structure(c(
          0.000110047078473174,
          0.0105052831822057,
          0.000955977750583648,
          0.000106231316866934,
          0.0151846812080705,
          0.000102493133434356
        ),
        .Dim = 3:2,
        .Dimnames = list(c("ex1", "ex2", "ex3"), c("ref1", "ref2"))
        ),
      modified.matrix =
        structure(c(
          0.000110047078473174, 0.0105052831822057,
          0.000955977750583648, 0.000106231316866934,
          9e+99, 0.000102493133434356
        ),
        .Dim = 3:2,
        .Dimnames =
          list(c("ex1", "ex2", "ex3"), c("ref1", "ref2"))
        )
    )
  )
})