Running ‘testthat.R’ [3s/4s]
Running the tests in ‘tests/testthat.R’ failed.
Complete output:
  > library(testthat)
  > library(kernelPSI)
  > 
  > test_check("kernelPSI")
   ----------- FAILURE REPORT -------------- 
   --- failure: the condition has length > 1 ---
   --- srcref --- 
  : 
   --- package (from environment) --- 
  kernelPSI
   --- call from context --- 
  pcaLR(K)
   --- call from argument --- 
  if (class(Ksum) != "kernelMatrix") {
      Ksum <- kernlab::as.kernelMatrix(Ksum, center = FALSE)
  }
   --- R stacktrace ---
  where 1: pcaLR(K)
  where 2: eval_bare(expr, quo_get_env(quo))
  where 3: quasi_label(enquo(object), arg = "object")
  where 4 at testthat/test_statistic.R#10: expect_type(pcaLR(K), "closure")
  where 5: eval(code, test_env)
  where 6: eval(code, test_env)
  where 7: withCallingHandlers({
      eval(code, test_env)
      if (!handled && !is.null(test)) {
          skip_empty()
      }
  }, expectation = handle_expectation, skip = handle_skip, warning = handle_warning, 
      message = handle_message, error = handle_error)
  where 8: doTryCatch(return(expr), name, parentenv, handler)
  where 9: tryCatchOne(expr, names, parentenv, handlers[[1L]])
  where 10: tryCatchList(expr, names[-nh], parentenv, handlers[-nh])
  where 11: doTryCatch(return(expr), name, parentenv, handler)
  where 12: tryCatchOne(tryCatchList(expr, names[-nh], parentenv, handlers[-nh]), 
      names[nh], parentenv, handlers[[nh]])
  where 13: tryCatchList(expr, classes, parentenv, handlers)
  where 14: tryCatch(withCallingHandlers({
      eval(code, test_env)
      if (!handled && !is.null(test)) {
          skip_empty()
      }
  }, expectation = handle_expectation, skip = handle_skip, warning = handle_warning, 
      message = handle_message, error = handle_error), error = handle_fatal, 
      skip = function(e) {
      })
  where 15: test_code(desc, code, env = parent.frame())
  where 16 at testthat/test_statistic.R#3: test_that("ridgeLR and maxLR return a closure", {
      n <- 30
      p <- 20
      K <- replicate(5, matrix(rnorm(n * p), nrow = n, ncol = p), 
          simplify = FALSE)
      K <- sapply(K, function(X) return(X %*% t(X)/dim(X)[2]), 
          simplify = FALSE)
      expect_type(ridgeLR(K, mu = 0, sigma = 1, lambda = 0.1), 
          "closure")
      expect_type(pcaLR(K), "closure")
  })
  where 17: eval(code, test_env)
  where 18: eval(code, test_env)
  where 19: withCallingHandlers({
      eval(code, test_env)
      if (!handled && !is.null(test)) {
          skip_empty()
      }
  }, expectation = handle_expectation, skip = handle_skip, warning = handle_warning, 
      message = handle_message, error = handle_error)
  where 20: doTryCatch(return(expr), name, parentenv, handler)
  where 21: tryCatchOne(expr, names, parentenv, handlers[[1L]])
  where 22: tryCatchList(expr, names[-nh], parentenv, handlers[-nh])
  where 23: doTryCatch(return(expr), name, parentenv, handler)
  where 24: tryCatchOne(tryCatchList(expr, names[-nh], parentenv, handlers[-nh]), 
      names[nh], parentenv, handlers[[nh]])
  where 25: tryCatchList(expr, classes, parentenv, handlers)
  where 26: tryCatch(withCallingHandlers({
      eval(code, test_env)
      if (!handled && !is.null(test)) {
          skip_empty()
      }
  }, expectation = handle_expectation, skip = handle_skip, warning = handle_warning, 
      message = handle_message, error = handle_error), error = handle_fatal, 
      skip = function(e) {
      })
  where 27: test_code(NULL, exprs, env)
  where 28: source_file(path, new.env(parent = env), chdir = TRUE, wrap = wrap)
  where 29: force(code)
  where 30: doWithOneRestart(return(expr), restart)
  where 31: withOneRestart(expr, restarts[[1L]])
  where 32: withRestarts(testthat_abort_reporter = function() NULL, force(code))
  where 33: with_reporter(reporter = reporter, start_end_reporter = start_end_reporter, 
      {
          reporter$start_file(basename(path))
          lister$start_file(basename(path))
          source_file(path, new.env(parent = env), chdir = TRUE, 
              wrap = wrap)
          reporter$.end_context()
          reporter$end_file()
      })
  where 34: FUN(X[[i]], ...)
  where 35: lapply(paths, test_file, env = env, reporter = current_reporter, 
      start_end_reporter = FALSE, load_helpers = FALSE, wrap = wrap)
  where 36: force(code)
  where 37: doWithOneRestart(return(expr), restart)
  where 38: withOneRestart(expr, restarts[[1L]])
  where 39: withRestarts(testthat_abort_reporter = function() NULL, force(code))
  where 40: with_reporter(reporter = current_reporter, results <- lapply(paths, 
      test_file, env = env, reporter = current_reporter, start_end_reporter = FALSE, 
      load_helpers = FALSE, wrap = wrap))
  where 41: test_files(paths, reporter = reporter, env = env, stop_on_failure = stop_on_failure, 
      stop_on_warning = stop_on_warning, wrap = wrap)
  where 42: test_dir(path = test_path, reporter = reporter, env = env, filter = filter, 
      ..., stop_on_failure = stop_on_failure, stop_on_warning = stop_on_warning, 
      wrap = wrap)
  where 43: test_package_dir(package = package, test_path = test_path, filter = filter, 
      reporter = reporter, ..., stop_on_failure = stop_on_failure, 
      stop_on_warning = stop_on_warning, wrap = wrap)
  where 44: test_check("kernelPSI")
  
   --- value of length: 2 type: logical ---
  [1] TRUE TRUE
   --- function from context --- 
  function (K, mu = 0, sigma = 1) 
  {
      if (is.list(K)) {
          Ksum <- Reduce(`+`, K)
      }
      else {
          Ksum <- K
      }
      if (class(Ksum) != "kernelMatrix") {
          Ksum <- kernlab::as.kernelMatrix(Ksum, center = FALSE)
      }
      Krot <- kernlab::kpca(Ksum)@rotated
      H <- Krot %*% pracma::pinv(Krot)
      M <- sum(eigen(H, symmetric = TRUE, only.values = TRUE)[["values"]] > 
          0.5)
      statistic <- function(Z) {
          roots <- Re(polyroot(c(-M + drop(t(Z - mu) %*% H %*% 
              Z)/sigma^2, -drop(t(2 * Z - mu) %*% H %*% Z)/sigma^2, 
              drop(t(Z) %*% H %*% Z)/sigma^2)))
          roots <- roots[roots < 1, drop = FALSE]
          R <- max(sapply(roots, function(r) 2 * M * log(1 - r) + 
              2 * r * drop(t(Z - mu) %*% H %*% Z)/sigma^2 - r^2 * 
              drop(t(Z) %*% H %*% Z)/sigma^2))
          return(R)
      }
  }
  <bytecode: 0x5636ba9a12e0>
  <environment: namespace:kernelPSI>
   --- function search by body ---
  Function pcaLR in namespace kernelPSI has this body.
   ----------- END OF FAILURE REPORT -------------- 
  Fatal error: the condition has length > 1