tests/testthat/test-tune.splsda.R

In mixOmics: Omics Data Integration Project

context("tune.splsda")

test_that("tune.splsda works", {
    data(breast.tumors)
    X = breast.tumors$gene.exp
    Y = as.factor(breast.tumors$sample$treatment)
    RNGversion(.mixo_rng()) ## in case RNG changes!

    set.seed(100)
    tune = tune.splsda(X, Y, ncomp = 2, nrepeat = 3, logratio = "none",
                       test.keepX = c(5, 10, 15), folds = 3, dist = "max.dist")
    
    expect_equal(tune$choice.ncomp$ncomp, 1L)
    expect_equal(tune$choice.keepX, c(comp1 = 5, comp2 = 15))
    

    set.seed(100)
    tune2 = tune.splsda(X, Y, ncomp = 2, nrepeat = 3, logratio = "none",
                       test.keepX = c(5, 10, 15), folds = 3, dist = "max.dist", cpus = 2)
    .almost_identical(tune, tune2)
    
})