+#' initSmallEM
+#'
#' initialization of the EM algorithm
#'
#' @param k number of components
#' @return a list with phiInit, rhoInit, piInit, gamInit
#'
#' @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)
}