tests/testthat/test_nipals.R

In pimentel/fscca: Fast SCCA

context("NIPALS CCA")

n <- 100
p <- 15
q <- 10
X <- matrix(rnorm(n*p), ncol = p)
Y <- matrix(rnorm(n*q), ncol = q)

test_that("input dimensions are ok",
    {
        expect_that(nipals(X[1:10, ], Y[1:3,]), throws_error())
    })

test_that("out dimensions are correct",
    {
        res <- nipals(X, Y)
        expect_equal(nrow(res$a), p)
        expect_equal(nrow(res$b), q)
    })


context("Sparse NIPALS")
test_that("lambdas are >= 0",
    {
        expect_that(sparse_nipals(X, Y, -2.0, 0.0), throws_error())
        expect_that(sparse_nipals(X, Y, 0.0, -0.1), throws_error())
        expect_that(sparse_nipals(X, Y, -10.0, -0.1), throws_error())
    })