Nothing
Q = 6
K = 4
N = 8
simu <- simulate_BBC(N = N, L = 100, epsilon = 0, theta_true = NULL)
dtm.full = simu$dtm.full
V = dim(dtm.full)[2]
Ytruth = simu$Ytruth
res <- mmpca_clust(simu$dtm.full, Q = 6, K = 4,
Yinit = Ytruth,
max.epochs = 1,
keep = 1,
verbose = 0,
nruns = 1,
mc.cores = 1)
test_that("check slots", code = {
expect_is(res, "mmpcaClust")
expect_equal(res@K, 4)
expect_equal(res@Q, 6)
expect_length(unique(res@clustering), Q)
expect_length(res@clustering, N)
expect_true(res@icl < 0)
expect_true(res@llhood < 0)
# check increasing bounds
expect_equal(res@logLikelihoods, sort(res@logLikelihoods, decreasing = F))
})
test_that("check dims", code = {
expect_equal(dim(res@beta), c(K,V))
expect_equal(dim(res@gamma), c(6, 4))
})
test_that("Check that Beta is proba matrix", code = {
expect_equal(Matrix::rowSums(res@beta), rep(1, K))
})
test_that("Check wrong Yinit error message", code = {
expect_error(mmpca_clust(simu$dtm.full, Q = 6, K = 4,
Yinit = c(1),
max.epochs = 1,
keep = 1,
verbose = 0))
expect_error(mmpca_clust(simu$dtm.full, Q = 6, K = 4,
Yinit = rep(1, N),
max.epochs = 1,
keep = 1,
verbose = 0))
expect_error(mmpca_clust(simu$dtm.full, Q = 6, K = 4,
Yinit = 'wrong!',
max.epochs = 1,
keep = 1,
verbose = 0))
})
test_that("Check wrong Yinit error message", code = {
init.beta = matrix(1, nrow = 4, data = V)
expect_error(mmpca_clust(simu$dtm.full, Q = 6, K = 4,
init.beta = init.beta,
max.epochs = 1,
keep = 1,
verbose = 0))
expect_error(mmpca_clust(simu$dtm.full, Q = 6, K = 4,
init.beta = 'wrong !',
max.epochs = 1,
keep = 1,
verbose = 0))
})
skip_on_cran()
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.