tests/testthat/test-bootstrap.R

In ruqianl/comapr: Crossover analysis and genetic map construction

# test boostrap
data(snp_geno_gr)
data("parents_geno")
test_that("bootstrapDist works", {

  corrected_geno <- correctGT(gt_matrix = mcols(snp_geno_gr),
                              ref = parents_geno$ref,
                              alt = parents_geno$alt)
  GenomicRanges::mcols(snp_geno_gr) <- corrected_geno
  #  co_geno$interval_ID <- rownames(co_geno)
  # melt_geno <- melt(co_geno,id= 'interval_ID')
  # colnames(melt_geno) <- c("interval_ID","SampleID","Cross_over")
  marker_gr_cos <- countCOs(snp_geno_gr)

  # expect_equal(sum(marker_gr_cos[4:5,"X92"]$X92),1)
  # dist <- calGeneticDist(marker_gr_cos)
  # expect_true(abs(dist$kosambi_cm[4]-13.67603)<1e-3)
  #
  bootstrap_results <- suppressWarnings(bootstrapDist(marker_gr_cos,B=10,
                                                      group_by  = c("X1","X9")))
  expect_true(length(bootstrap_results)==10)
})

test_that("permuteDist works", {

  corrected_geno <- correctGT(gt_matrix = mcols(snp_geno_gr),
                              ref = parents_geno$ref,
                              alt = parents_geno$alt)
  GenomicRanges::mcols(snp_geno_gr) <- corrected_geno
  #  co_geno$interval_ID <- rownames(co_geno)
  # melt_geno <- melt(co_geno,id= 'interval_ID')
  # colnames(melt_geno) <- c("interval_ID","SampleID","Cross_over")
  marker_gr_cos <- countCOs(snp_geno_gr)

  # expect_equal(sum(marker_gr_cos[4:5,"X92"]$X92),1)
  # dist <- calGeneticDist(marker_gr_cos)
  # expect_true(abs(dist$kosambi_cm[4]-13.67603)<1e-3)
  #
  perms_results <- suppressWarnings(permuteDist(marker_gr_cos,B=10,
                                                      group_by  = c("X1","X9")))
  expect_true(!is.null(perms_results$permutes))
  expect_true(!is.null(perms_results$observed_diff))
  expect_true(!is.null(perms_results$nSample))

})