1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
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)
output_best <- fit_to_signatures_strict(mut_mat, signatures[,1:5], max_delta = 0.004, method = "best_subset")
output_single_sig = fit_to_signatures_strict(mut_mat, signatures[,1, drop = F], max_delta = 0.05)
expected <- readRDS(system.file("states/strict_snv_refit.rds",
package = "MutationalPatterns"
))
expected_best <- readRDS(system.file("states/strict_best_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"))
expect_true(inherits(output_best, "list"))
expect_true(inherits(output_best$fit_res, "list"))
expect_true(inherits(output_best$fit_res$contribution, "matrix"))
expect_true(inherits(output_best$fit_res$reconstructed, "matrix"))
expect_true(inherits(output_best$sim_decay_fig, "list"))
expect_true(inherits(output_best$sim_decay_fig[[1]], "gg"))
expect_true(inherits(output_single_sig, "list"))
expect_true(inherits(output_single_sig$fit_res, "list"))
expect_true(inherits(output_single_sig$fit_res$contribution, "matrix"))
expect_true(inherits(output_single_sig$fit_res$reconstructed, "matrix"))
expect_true(inherits(output_single_sig$sim_decay_fig, "list"))
expect_true(inherits(output_single_sig$sim_decay_fig[[1]], "gg"))
})
test_that("Output is equal to expected", {
expect_equal(output$fit_res, expected)
expect_equal(output_best$fit_res, expected_best)
})
|
