tests/testthat/test-phylo2l.R

In treestats: Phylogenetic Tree Statistics

context("phylo2L")

test_that("usage", {
  set.seed(42)
  focal_tree <- ape::rphylo(n = 500, birth = 1, death = 0)

  ltable_1 <- DDD::phylo2L(focal_tree)
  ltable_2 <- treestats::phylo_to_l(focal_tree)

  diff <- ltable_1 - ltable_2

  for (i in 1:4) {
    testthat::expect_lt(mean(diff[, i]), 1e-6)
  }

  # again, but with extinct lineages:
  set.seed(42)
  focal_tree <- ape::rphylo(n = 100, birth = 1, death = 0.3, fossils = TRUE)

  ltable_1 <- DDD::phylo2L(focal_tree)
  ltable_2 <- treestats::phylo_to_l(focal_tree)

  diff <- ltable_1 - ltable_2
  for (i in 1:4) {
    testthat::expect_lt(mean(diff[, i]), 1e-6)
  }
})