tests/testthat/test_initialize_ploidy.R

In SCOPE: A normalization and copy number estimation method for single-cell DNA sequencing

context("PreEst_Ploidy")
library(SCOPE)

test_that("Ploidy Initialization works", {
  Gini <- get_gini(Y_sim)
  normObj.sim <- normalize_codex2_ns_noK(Y_qc = Y_sim,
                                         gc_qc = ref_sim$gc,
                                         norm_index = which(Gini<=0.12))

  Yhat.noK.sim <- normObj.sim$Yhat
  ploidy.sim <- initialize_ploidy(Y = Y_sim, Yhat = Yhat.noK.sim, ref = ref_sim)
  expect_equal(sum(Gini<=0.12), sum(ploidy.sim==2))
})