#' @param k number of components
#' @param X matrix of covariates (of size n*p)
#' @param Y matrix of responses (of size n*m)
+#' @param fast boolean to enable or not the C function call
#'
#' @return a list with phiInit, rhoInit, piInit, gamInit
#'
-#' @importFrom methods new
#' @importFrom stats cutree dist hclust runif
#' @export
initSmallEM <- function(k, X, Y, fast)
piInit <- piInit1[b, ]
gamInit <- gamInit1[, , b]
- return(list(phiInit = phiInit, rhoInit = rhoInit, piInit = piInit, gamInit = gamInit))
+ list(phiInit = phiInit, rhoInit = rhoInit, piInit = piInit, gamInit = gamInit)
}