tests/test_that_mfa.R

In suren-rathnayake/EMMIXmfa: Mixture Models with Component-Wise Factor Analyzers

rm(list=ls())
require(testthat)
require(EMMIXmfa)

set.seed(1)
Y <- iris[, -5]
g <- 3
q <- 2
p <- ncol(Y)
n <- nrow(Y)

context("mfa")
model <- mfa(Y, g, q, nkmeans = 2, nrandom = 2, tol = 1.e-1)

expect_that(0, equals(minmis(model$clust, model$clust)))
expect_that(1, equals(ari(model$clust, model$clust)))
expect_that(model$clust, equals(predict(model, Y)))


mod2 <- mfa (Y, g, q, init_para = model)

test_that("mfa with init_para as past model", {
  expect_that(mod2, is_a("mfa"))
})

test_that("mfa with init_para as past model", {
  expect_that(mod2, is_a("emmix"))
})

X <- Y
X[1, 4] <- NA
expect_that(mfa(X, g, q),
            throws_error("`Y' has missing value"))

X[1, 4] <- "a"
expect_that(mfa(X, g, q),
            throws_error("`Y' has a non-numeric element"))

X <- X[,1]
expect_that(mfa(X, g, q),
            throws_error("The data must have more than one variable."))

X <- Y
expect_that(mfa(X, g, ncol(X)),
            throws_error("The number of factors must be less than the number of variables."))

expect_that(mfa(X, g, ncol(X)+1),
            throws_error("The number of factors must be less than the number of variables."))

expect_that(mfa(X, g = g, q = -1),
            throws_error("q must be a positive integer."))

expect_that(mfa(X, g = -2, q = q),
            throws_error("g must be a positive integer."))

expect_that(mfa(X, g = 2.5, q = q),
            throws_error("g must be a positive integer."))

expect_that(mfa(X, g = g, q = 1.5),
            throws_error("q must be a positive integer."))

test_that("mfa with itmax NULL", {
  expect_that(mfa (Y, g, q, itmax=NULL), throws_error())
})

test_that("mfa with itmax NULL", {
  expect_that(mfa (Y, g, q, itmax=NULL), throws_error())
})

test_that("mfa with itmax neg", {
  expect_that(mfa (Y, g, q, itmax=-1), throws_error())
})

context("q = 1")
q <- 1
context("mfa")
model <- mfa(Y, g, q, nkmeans = 2, nrandom = 2, tol = 1.e-5)

X <- Y
X[1, 4] <- NA
expect_that(mfa(X, g, q),
            throws_error("`Y' has missing value"))

X[1, 4] <- "a"
expect_that(mfa(X, g, q),
            throws_error("`Y' has a non-numeric element"))

X <- X[,1]
expect_that(mfa(X, g, q),
            throws_error("The data must have more than one variable."))

X <- Y
expect_that(mfa(X, g, ncol(X)),
            throws_error("The number of factors must be less than the number of variables."))

expect_that(mfa(X, g, ncol(X)+1),
            throws_error("The number of factors must be less than the number of variables."))

expect_that(mfa(X, g = g, q = -1),
            throws_error("q must be a positive integer."))

expect_that(mfa(X, g = -2, q = q),
            throws_error("g must be a positive integer."))

expect_that(mfa(X, g = 2.5, q = q),
            throws_error("g must be a positive integer."))

expect_that(mfa(X, g = g, q = 1.5),
            throws_error("q must be a positive integer."))

test_that("mfa with itmax NULL", {
  expect_that(mfa (Y, g, q, itmax=NULL), throws_error())
})

test_that("mfa with itmax NULL", {
  expect_that(mfa (Y, g, q, itmax=NULL), throws_error())
})

test_that("mfa with itmax neg", {
  expect_that(mfa (Y, g, q, itmax=-1), throws_error())
})

test_that("mfa with init_para as past model", {
  expect_that(mfa (Y, g, q, init_para = model), is_a("mfa"))
})

dat <- rmix(100, model)

test_that("mfa with init_para as past model", {
  expect_that(dat, is_a("matrix"))
})

test_that("mfa with init_para as past model", {
  expect_that(nrow(dat), equals(100))
})

test_that("mfa with init_para as past model", {
  expect_that(ncol(dat), equals(p))
})