tests/testthat/test_optimize_exposure_QP.R

In mSigTools: Mutational Signature Analysis Tools

test_that("optimize_exposure_QP", {
  skip_if("" == system.file(package = "ICAMS"))
  skip_if("" == system.file(package = "cosmicsig"))
  spectrum <- ICAMS::ReadCatalog(file = "testdata/SBS96_spectrum.csv")
  rr <- optimize_exposure_QP(spectrum = spectrum,
                             signatures = cosmicsig::COSMIC_v3.0$signature$GRCh37$SBS96[ , 1:5])
  expect_equal(rr,
               c(SBS1 = 1383.896, SBS2 = 1684.294,
                 SBS3 = 6958.131,  SBS4 = 1196.792,
                 SBS5 = 3676.888), tolerance = 1e-5)

})