tests/testthat/testdeltaa.R

In paleomorph: Geometric Morphometric Tools for Paleobiology

context('deltaa')

test_that('deltaa works', {

  # delta should be 0
  olda <- array( rep(1:6, 3, each = 3), dim = c(2, 3, 3))
  newa <- array( rep(1:6, 3, each = 3), dim = c(2, 3, 3))

  m <- 3
  n <- 2

  expect_equal(0, deltaa(olda, newa, m, n, scaleDelta = FALSE))


  # small example.
  # each coord in newa is +1 coord in olda.
  # Each cooord has distance of 1.73205080756887719318.
  # Total is 1.732051 * 2 * 2

  olda <- array( rep(1:4, 3), dim = c(2, 3, 2))
  newa <- olda + 1

  m <- 2
  n <- 2

  expect_equal(1.73205080756887719318 * 2 * 2, deltaa(olda, newa, m, n, scaleDelta = FALSE))

})



test_that('deltaa works for negative coords', {

  # small example.
  # each coord in newa is +1 coord in olda.
  # Each cooord has distance of 1.73205080756887719318.
  # Total is 1.732051 * 2 * 2

  olda <- array( rep(0, 3), dim = c(2, 3, 2))
  newa <- array( rep(1, 3), dim = c(2, 3, 2))

  m <- 2
  n <- 2

  positive <- deltaa(olda, newa, m, n, scaleDelta = FALSE)

  newa <- array( rep(-1, 3), dim = c(2, 3, 2))
  negative <- deltaa(olda, newa, m, n, scaleDelta = FALSE) 
 
  expect_equal(positive, negative)


})





test_that('deltaa works with scaleDelta = TRUE', {

  # small example.
  # each coord in newa is +1 coord in olda.
  # Each cooord has distance of 1.73205080756887719318.
  # Total is 1.732051 * 2 * 2

  olda <- array( rep(0, 3), dim = c(2, 3, 2))
  newa <- array( rep(1, 3), dim = c(2, 3, 2))

  m <- 2
  n <- 2

  notScaled <- deltaa(olda, newa, m, n, scaleDelta = FALSE)

  scaled <- deltaa(olda, newa, m, n, scaleDelta = TRUE)

  expect_true(scaled < notScaled)
  expect_true(scaled * dim(olda)[3] * dim(olda)[1] == notScaled)


})