#' @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
-#' @export
-#' @importFrom methods new
+#'
#' @importFrom stats cutree dist hclust runif
+#' @export
initSmallEM <- function(k, X, Y, fast)
{
n <- nrow(X)
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)
}