test_that("hopkins_vs_components returns valid output", {

  set.seed(1)

  n <- 10
  T_len <- 6

  X <- list(
    lapply(seq_len(n), function(i) matrix(rnorm(T_len * 3), T_len, 3)),
    lapply(seq_len(n), function(i) matrix(rnorm(T_len * 2), T_len, 2))
  )

  fit <- mcca_fit(
    method = "kernel",
    X = X,
    ncomp = 3
  )

  H <- hopkins_vs_components(
    fit,
    max_comp = 2,
    nrep = 5
  )

  expect_true(is.data.frame(H))

  expect_true(all(c("components", "hopkins", "sd") %in% names(H)))

  expect_equal(nrow(H), 2)

})