tests/testthat/test_hyper3.R

In RobinHankin/hyper2: The Hyperdirichlet Distribution, Mark 2

test_that("Miscellaneous hyper3 tests", {
    h1 <- hyper3_bw(B=list("a",c("a","b"),"b"),W=list(1.2,c(1.2,1),1),powers=c(3,4,-7))
    h2 <- hyper3_nv(list(c(b=3,c=1,a=9),c(x=6,y=8)),c(5,-5))

    M <- rbind(c(1.2,0,0,0,0),c(1.2,1,0,0,0),c(9,3,1,0,0),c(0,1,0,0,0),c(0,0,0,6,8))
    colnames(M) <- c("a","b","c","x","y")
    h3 <- hyper3_m(M,c(3,4,5,-7,-5))

    expect_true(h1+h2 == h3)

    expect_error(hyper3_bw(list(c("a","b"),c("a","b","d")),list(c(1,4),c(1,-3,4)),c(1,-1),letters[1:5]))

    
} )