tests/testthat/test_fsteachsnphudson.R

context("Fst Hudson all SNPs")

test_that("calculate Fst Hudson for all SNPs",{

  data(example_SNP)

  idx1 <- which(sample_labels == 'pop1')
  idx2 <- which(sample_labels == 'pop2')
  fst.pairwise <- fst.each.snp.hudson(simsnp$snp, idx1, idx2)

  expect_length(fst.pairwise, 2000)
})
kridsadakorn/kris documentation built on April 28, 2022, 4:23 a.m.