tests/testthat/test_dOptPenalty.R

In Theorell/DepecheR: Determination of essential phenotypic elements of clusters in high-dimensional entities

###################### dOptPenalty
context("dOptPenatly")

x <- DepecheR:::generateBimodalData(observations = 200, dataCols = 10)
x_scaled <- dScale(x$samples, center = FALSE)
dOPR <- DepecheR:::dOptPenalty(x_scaled, k = 30, maxIter = 20, 
                               minARIImprovement = 0.01, sampleSize = 100, 
                               penalties = c(0, 2, 4, 8, 16, 32, 128),
                               createOutput = FALSE, optimARI = 0.95, 
                               nCores = 2, plotDir = ".")

ARI <- dOPR[[2]][, 1]
nClust <- dOPR[[2]][, 2]
bestPenalty <- dOPR$bestPenalty
rowNamesdOpr <- row.names(dOPR[[2]])
nClust[which(rowNamesdOpr == bestPenalty)]
ARI[which(rowNamesdOpr == bestPenalty)]

test_that("dOptPenalty expected output", {
    expect_equal(nClust[which(rowNamesdOpr == bestPenalty)], 2)
    expect_equal(ARI[which(rowNamesdOpr == bestPenalty)], 1)
})