scripts/generate_expectations.R

library(straf)

input_diploid <- list(datapath="./tests/testthat/exampleSTRAFdiplo.txt")
dat_diploid <- straf::createGenind(input_diploid, ploidy = 2)
freq_diploid <- straf::getFreqAllPop(dat_diploid)
indices_diploid <- straf::getIndicesAllPop(dat_diploid)


input_diploid2 <- list(datapath="./tests/testthat/exampleSTRAFdiplo_empty_lines.txt")
dat_diploid2 <- straf::createGenind(input_diploid2, ploidy = 2)
freq_diploid2 <- straf::getFreqAllPop(dat_diploid2)
indices_diploid2 <- straf::getIndicesAllPop(dat_diploid2)


input_haploid <- list(datapath="./tests/testthat/exampleSTRAFhaplo.txt")
dat_haploid <- straf::createGenind(input_haploid, ploidy = 1)
freq_haploid <- straf::getFreqAllPop(dat_haploid)
indices_haploid <- straf::getIndicesAllPop(dat_haploid)


input_haploid_point <- list(datapath="./tests/testthat/exampleSTRAFhaplo_point.txt")
dat_haploid_point <- straf::createGenind(input_haploid_point, ploidy = 1)
freq_haploid_point <- straf::getFreqAllPop(dat_haploid_point)
indices_haploid_point <- straf::getIndicesAllPop(dat_haploid_point)

expectations <- list(
  dat_diploid = dat_diploid,
  freq_diploid = freq_diploid,
  indices_diploid = indices_diploid,
  dat_diploid2 = dat_diploid2,
  freq_diploid2 = freq_diploid2,
  indices_diploid2 = indices_diploid2,
  dat_haploid = dat_haploid,
  freq_haploid = freq_haploid,
  indices_haploid = indices_haploid,
  dat_haploid_point = dat_haploid_point,
  freq_haploid_point = freq_haploid_point,
  indices_haploid_point = indices_haploid_point
)

saveRDS(expectations, file = "./tests/testthat/expectations.rds")
agouy/straf documentation built on Dec. 9, 2024, 4:35 p.m.