tests/testthat/test-spectral_clustering.R

In TreeDist: Calculate and Map Distances Between Phylogenetic Trees

d <- dist(c(1:20, 60:41))

test_that("Spectral clustering fails gracefully", {
  expect_error(SpectralEigens(d, nEig = 0), "nEig must be.*positive")
})

test_that("Spectral clustering works", {
  allEig <- SpectralEigens(d, nEig = Inf)
  expect_equal(dim(allEig), c(40, 40))
  expect_equal(abs(SpectralEigens(d, nEig = 2)), abs(allEig[, 40:39]),
               tolerance = sqrt(.Machine[["double.eps"]]))

  expect_warning(deprecated <- SpectralClustering(d, nEig = Inf),
                 "'SpectralClustering' is deprecated.")
  # Can't check raw equality because signs may flip
  expect_true(all.equal(crossprod(deprecated), crossprod(allEig),
                        tolerance = sqrt(.Machine[["double.eps"]])))
})