tests/testthat/test-distance.R

context("distance")

test_that("Euclidean distance is correct", {
          set.seed(888)
          x <- matrix(runif(50), ncol=5)
          d.ecodist <- distance(x, "euclidean")
          d.base <- dist(x)
          expect_equal(attributes(d.ecodist), attributes(d.base)[names(attributes(d.base)) != "call"])
          expect_equal(as.vector(d.ecodist), as.vector(d.base))
})

test_that("Bray-Curtis distance is correct", {
          set.seed(888)
          x <- matrix(runif(50), ncol=5)
          expect_equal(distance(x, "bray"), bcdist(x))
})

test_that("Mahalanobis icov is correct", {
          set.seed(888)
          x <- matrix(runif(50), ncol=5)
          x.md <- full(distance(x, "mahal"))
          sub.md <- full(distance(x[1:5, ], "mahal", icov=cov(x)))
          expect_equal(x.md[1:5, 1:5], sub.md)
})
phiala/ecodist documentation built on Nov. 5, 2023, 10:47 a.m.