tests/testthat/test-updateAF.R

In GregoryMatesi/RHiddenAncestries: Hidden Ancestries Group R Package

test_that("updateAF works on ancestryData", {
  expect <- c(0.173275495, 0.001237745, 0.168316089, 0.428212624, 0.204214851, 0.004950604)
  data("ancestryData")
  testval <- updateAF(data   = ancestryData,
                      reference   = c("ref_AF_eur_1000G", "gnomAD_AF_afr"),
                      observed    = "gnomAD_AF_afr",
                      pi_target   = c(1,0),
                      pi_observed = c(.15,.18))
  expect_equal(head(testval$adjustedAF), expected = expect)
})