tests/testthat/test-distoredDetect.R

In ruqianl/comapr: Crossover analysis and genetic map construction

test_that("detectDisorted markers works", {
  or_geno <- snp_geno[,grep("X",colnames(snp_geno))]
  rownames(or_geno) <- paste0(snp_geno$CHR,"_",snp_geno$POS)
  or_geno[1,] <- rep("Fail",dim(or_geno)[2])
  cr_geno <- correctGT(or_geno,ref = snp_geno$C57BL.6J,
                       alt = snp_geno$FVB.NJ..i.)
  ft_gt <- filterGT(cr_geno)
  ft_gt[1,] <- "Het"
  disor_gt <- getDistortedMarkers(ft_gt)
  expect_true(disor_gt[1,2]==22)
})