#' @importFrom stats cutree dist hclust runif
initSmallEM <- function(k, X, Y, fast)
{
- n <- nrow(Y)
- m <- ncol(Y)
+ n <- nrow(X)
p <- ncol(X)
+ m <- ncol(Y)
nIte <- 20
Zinit1 <- array(0, dim = c(n, nIte))
betaInit1 <- array(0, dim = c(p, m, k, nIte))