tests/testthat/test-fit_to_signatures_strict.R

In MutationalPatterns: Comprehensive genome-wide analysis of mutational processes

context("test-fit_to_signatures_strict")

# Get mut_mat
mut_mat <- readRDS(system.file("states/mut_mat_data.rds",
  package = "MutationalPatterns"
))

# Get signatures
signatures <- get_known_signatures()

output <- fit_to_signatures_strict(mut_mat, signatures, max_delta = 0.05)

expected <- readRDS(system.file("states/strict_snv_refit.rds",
  package = "MutationalPatterns"
))

test_that("Output has correct class", {
  expect_true(inherits(output, "list"))
  expect_true(inherits(output$fit_res, "list"))
  expect_true(inherits(output$fit_res$contribution, "matrix"))
  expect_true(inherits(output$fit_res$reconstructed, "matrix"))
  expect_true(inherits(output$sim_decay_fig, "list"))
  expect_true(inherits(output$sim_decay_fig[[1]], "gg"))
})

test_that("Output is equal to expected", {
  expect_equal(output$fit_res, expected)
})