tests/testthat/test-ew_colless.R

In treestats: Phylogenetic Tree Statistics

context("ew_colless")

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

    colless <- treestats::ew_colless(focal_tree)
    colless_check <- treebalance::ewCollessI(focal_tree)
    testthat::expect_equal(colless, colless_check)

    focal_ltab <- treestats::phylo_to_l(focal_tree)

    colless <- treestats::ew_colless(focal_ltab)
    testthat::expect_equal(colless, colless_check)

    ## with extinct lineages:
    focal_tree <- ape::rphylo(n = 100, birth = 1, death = 0.2, fossils = TRUE)

    colless <- treestats::ew_colless(focal_tree)
    colless_check <- treebalance::ewCollessI(focal_tree)
    testthat::expect_equal(colless, colless_check)

    focal_ltab <- treestats::phylo_to_l(focal_tree)

    colless <- treestats::ew_colless(focal_ltab)
    testthat::expect_equal(colless, colless_check)
  }
})

test_that("wrong_object", {
  testthat::expect_error(
    treestats::ew_colless(10),
    "input object has to be phylo or ltable"
  )

  testthat::expect_error(
    treestats::ew_colless(list()),
    "input object has to be phylo or ltable"
  )
})