tests/testthat/test-RunHdpxParallel-multitypeF-fast96.R

In steverozen/mSigHdp: Mutational signature discovery using HDP (hierarchical Dirichlet process)

test_that("RunHdpxParallel-fast96", {

  input.catalog <-
    ICAMS::ReadCatalog("SBS96.ground.truth/ground.truth.syn.catalog.csv")

  reg <- new.env()
  load("RunhdpInternal.testdata/NewRunHdpParallel-fast96-2-cores.Rdata",
       envir = reg)

  retvalx <- RunHdpxParallel(
    input.catalog     = input.catalog[1:10,1:15],
    CPU.cores         = 2,
    seedNumber        = 44,
    K.guess           = 5,
    multi.types       = FALSE,
    verbose           = FALSE,
    num.child.process =  2,
    burnin            = 50, # Super low for fast testing
    post.space        = 5,  # Low for fast testing
    post.cpiter       = 1,  # Low for fast testing
    overwrite         = TRUE,
    checkpoint        = FALSE,
    out.dir           = tempfile()
  )

  if (FALSE) { # To regenerate test data
    save(retvalx,
         file = "RunhdpInternal.testdata/NewRunHdpParallel-fast96-2-cores.Rdata")
  }
  expect_equal(retvalx, reg$retvalx)

})