test_that("predict returns scores with correct dimension", {

  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 = 2)

  pr <- predict(fit, X)

  expect_true(is.list(pr))
  expect_equal(length(pr), 2)

  expect_equal(ncol(pr[[1]]), 2)
  expect_equal(ncol(pr[[2]]), 2)

})