test_that("hard cases work", {
  x <- c(6L, 87L, 429L, 392L, 86L)
  g1 <- 3
  g2 <- 2
  par = c(1e-5, 1e-5, 1 - 1e-5)
  expect_no_error(
    obj_dr_pp(par = par, x = x, g1 = g1, g2 = g2)
  )
})

test_that("hard case works for tetraploids", {
  # load("./mout_Cc.rda")
  # mout_Cc$snpdf$snp[3055:3056]
  # sub <- filter_snp(mout_Cc, snp %in% c("chr4c_2876569", "chr4c_2876589"))
  # mhard_g <- multidog_to_g(mout = sub, ploidy = 4, type = "off_g")
  # mhard_gl <- multidog_to_g(mout = sub, ploidy = 4, type = "off_gl")
  # save(mhard_g, mhard_gl, file = "./mhard.RData")

  load("./mhard.RData")

  expect_no_error(
    {
      #seg_lrt(x = mhard_gl$g[1, ,], p1_ploidy = 4, p2_ploidy = 4, p1 = mhard_gl$p1[[1]], p2 = mhard_gl$p1[[1]])
      #seg_multi(g = mhard_gl$g, p1_ploidy = 4, p2_ploidy = 4, p1 = mhard_gl$p1, p2 = mhard_gl$p2)
      em_li(mhard_gl$g[1, ,])
    }
  )
})