Merge branch 'update_unit_tests' into dev

Author: itsrainingdata

tests/testthat/test-ccdr_call.R

@@ -0,0 +1,60 @@
+context("ccdr_call")
+
+pp <- 10L
+nn <- 5L
+X.test <- data.frame(matrix(runif(pp*pp), ncol = pp))
+cors.test <- sparsebnUtils::cor_vector(X.test)
+betas.test <- matrix(runif(pp*pp), ncol = pp)
+indexj.test <- rep(0L, pp + 1)
+nj.test <- as.integer(rep(nn, pp))
+cors.length <- pp*(pp+1)/2
+lambda.test <- pi
+gamma.test <- 2.0
+eps.test <- 0.1
+maxIters.test <- 1000L
+alpha.test <- 10
+
+### The call for ccdr_call
+# ccdr_call <- function(data,
+#                       ivn = NULL,
+#                       betas,
+#                       lambdas,
+#                       lambdas.length,
+#                       gamma,
+#                       error.tol,
+#                       rlam,
+#                       max.iters,
+#                       alpha,
+#                       verbose = FALSE)
+
+test_that("ccdr_call runs as expected", {
+    ### lambdas unspecified
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = 20, rlam = 0.01, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), NA)
+
+    ### lambdas specified
+    expect_error(ccdr_call(data = X.test, lambdas = c(10, 5, 1), gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), NA)
+})
+
+test_that("Check input: data", {
+    expect_error(ccdr_call(data = 1, lambdas = NULL, lambdas.length = 20, rlam = 0.01, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "'data' argument must be a data.frame or matrix!")
+    expect_error(ccdr_call(data = matrix(c(NA,1,1,1), nrow = 2), lambdas = NULL, lambdas.length = 20, rlam = 0.01, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "missing values detected!")
+})
+
+test_that("Check input: lambdas, lambdas.length", {
+    # lambdas.length
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = NULL, rlam = 0.01, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "Both lambdas and lambdas.length unspecified: Must specify a value for at least one of these arguments!")
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = "test", rlam = 0.01, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "lambdas.length must be numeric!")
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = -1, rlam = 0.01, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "lambdas.length must be positive!")
+
+    # lambdas
+    expect_error(ccdr_call(data = X.test, lambdas = c(10, "test", 1), gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "lambdas must be a numeric vector!")
+    expect_error(ccdr_call(data = X.test, lambdas = c(10, 5, -1), gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "lambdas must contain only nonnegative values!")
+})
+
+test_that("Check input: rlam", {
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = 20, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "rlam must be specified if lambdas is not explicitly specified.")
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = 20, rlam = "test", gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "rlam must be numeric!")
+    expect_error(ccdr_call(data = X.test, lambdas = NULL, lambdas.length = 20, rlam = -1, gamma = gamma.test, error.tol = eps.test, max.iters = maxIters.test, alpha = alpha.test), "rlam must be >= 0!")
+})
+
+

tests/testthat/test-ccdr_gridR.R

@@ -1,24 +1,32 @@
 context("ccdr_gridR")
 
-pp <- 10
-nn <- 5
-X.test <- matrix(runif(pp*pp), ncol = pp)
-X.test <- data.frame(X.test)
-cors.test <- sparsebnUtils::cor_vector_ivn(X.test, NULL)
+pp <- 10L
+nn <- 5L
+X.test <- data.frame(matrix(runif(pp*pp), ncol = pp))
+cors.test <- sparsebnUtils::cor_vector(X.test)
 betas.test <- matrix(runif(pp*pp), ncol = pp)
+indexj.test <- rep(0L, pp + 1)
+nj.test <- as.integer(rep(nn, pp))
+cors.length <- pp*(pp+1)/2
+lambda.test <- pi
+gamma.test <- 2.0
+eps.test <- 0.1
+maxIters.test <- 1000L
+alpha.test <- 10
 
-# test_that("ccdr_gridR runs as expected", {
-#     ### No error / Also tests default values for betas, nlam, gamma, eps, maxIters, alpha
-#     expect_error(ccdr_gridR(cors = cors.test), NA)
-# })
+### The call for ccdr_gridR
+# ccdr_gridR <- function(cors,
+#                        pp, nn,
+#                        nj = NULL,
+#                        indexj = NULL,
+#                        betas,
+#                        lambdas,
+#                        gamma,
+#                        eps,
+#                        maxIters,
+#                        alpha,
+#                        verbose)
 
-# test_that("Check input: rlam", {
-#     ### eps is non-numeric
-#     expect_error(ccdr.run(data = dat.test, lambdas.length = lambdas.length.test, rlam = "not a number"))
-#
-#     ### rlam is negative (note that rlam = 0 is OK)
-#     expect_error(ccdr.run(data = dat.test, lambdas.length = lambdas.length.test, rlam = -1e-5))
-#
-#     ### Need to add a test for zero case
-#
-# })
+test_that("ccdr_gridR runs as expected", {
+    expect_error(ccdr_gridR(cors = as.numeric(cors.test), pp = pp, nn = nn, betas = betas.test, lambdas = lambda.test, gamma = gamma.test, eps = eps.test, maxIters = maxIters.test, alpha = alpha.test, verbose = FALSE), NA)
+})

tests/testthat/test-ccdr_singleR.R

@@ -1,17 +1,17 @@
 context("ccdr_singleR")
 
-    pp <- 10L
-    nn <- 5L
-    indexj <- rep(0L, pp + 1)
-    nj <- as.integer(rep(nn, pp))
-    cors.length <- pp*(pp+1)/2
-    cors.test <- runif(cors.length)
-    betas.test <- matrix(runif(pp*pp), ncol = pp)
-    lambda.test <- pi
-    gamma.test <- 2.0
-    eps.test <- 0.1
-    maxIters.test <- 1000L
-    alpha.test <- 10
+pp <- 10L
+nn <- 5L
+indexj <- rep(0L, pp + 1)
+nj <- as.integer(rep(nn, pp))
+cors.length <- pp*(pp+1)/2
+cors.test <- runif(cors.length)
+betas.test <- matrix(runif(pp*pp), ncol = pp)
+lambda.test <- pi
+gamma.test <- 2.0
+eps.test <- 0.1
+maxIters.test <- 1000L
+alpha.test <- 10
 
 ### The call for ccdr_singleR
 # ccdr_singleR <- function(cors,