tests/testthat/test_filterAge.R

In rmsharp/nprcmanager: Genetic Tools for Colony Management

#' Copyright(c) 2017-2020 R. Mark Sharp
#' This file is part of nprcgenekeepr
context("filterAge")
library(testthat)

ped <- nprcgenekeepr::qcPed
ped$gen <- findGeneration(ped$id, ped$sire, ped$dam)
kmat <- kinship(ped$id, ped$sire, ped$dam, ped$gen,
                sparse = FALSE)
kin <- kinMatrix2LongForm(kmat)
filteredKin1 <- filterAge(kin, ped, minAge = 1)
filteredKin5 <- filterAge(kin, ped, minAge = 5)
lessThan5 <- c("C1ICXL", "2KULR3", "RI0O7F", "7M51X5", "170ZTZ", "Y7PPEZ",
               "CFPEEU", "ZC5SCR", "218FOV", "2IXJ2N", "CAST4W", "JGPN6K",
               "HOYW0S", "DD1U77", "0DAV0I", "HLI95R", "TZ5NUB", "DR5GXB",
               "EUG3WE", "FHV13N", "OUM6QF", "6Z7MD9", "309VM2", "8KM1MP",
               "I9TQ0T", "INGWI7")
test_that("filterAge removes the correct rows", {
  expect_false(any("Y7PPEZ" %in% filteredKin1$id1))
  expect_true(any("C1ICXL" %in% filteredKin1$id1))
  expect_false(any("C1ICXL" %in% filteredKin5$id1))
}
)