tests/testthat/test_count_vcf_individuals.R

In tidypopgen: Tidy Population Genetics

skip_if_not_installed("vcfR")

test_that("count_vcf_individuals with gzfile", {
  vcf_path <- system.file(
    "extdata/ploidy/ploidy_test.vcf.gz",
    package = "tidypopgen"
  )
  n_individuals <- count_vcf_individuals(vcf_path)
  vcfr_obj <- vcfR::read.vcfR(vcf_path, verbose = FALSE)
  expect_equal(n_individuals, ncol(vcfR::extract.gt(vcfr_obj)))
})