inst/tinytest/test_hadamard.R

In KlausVigo/phangorn: Phylogenetic Reconstruction and Analysis

v <- 1:8
data(yeast)
dm <- dist.hamming(yeast)
# RY-coding
yeast_ry <- acgt2ry(yeast)


# test Hadamard conjugation
# fast fft like multiplication
expect_true(inherits(H <- hadamard(3), "matrix"))
expect_equal(as.vector(H %*% v), fhm(v))
expect_true(inherits(spl_ry <- h2st(yeast_ry), "splits"))
expect_true(inherits(spl_dm <- distanceHadamard(dm), "splits"))
expect_true(inherits(fit4 <- h4st(yeast)[[1]], "splits"))