tests/testthat/test-coalgenhetero.R

In xiangfstats/GenieR: Inference of demographic history from phylogeny

context("Vectorised trajectory and a n by 2 matrix sample")

test_that("a vectorised trajectory and a n by 2 matrix sample", {
  sample1<-cbind(c(9,1,2,1),c(0,.008,.03,.1))
  trajectory<-function(x) 100
  example_hetero<-coalgen_hetero(sample1, Vectorize(trajectory))
  expect_error(coalgen_hetero(sample1, trajectory))
  expect_error(coalgen_hetero(sample1[,1], Vectorize(trajectory)))
})


test_that("To ensure the correct output", {
  sample1<-cbind(c(9,1,2,1),c(0,.008,.03,.1))
  trajectory<-function(x) exp(10*x)
  example_hetero<-coalgen_hetero(sample1, trajectory)
  expect_true(class(example_hetero)=="list" && length(example_hetero)==2)
})