+#' 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)
{