#' @export
#' @importFrom methods new
#' @importFrom stats cutree dist hclust runif
-initSmallEM <- function(k, X, Y, fast = TRUE)
+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))
dotProduct <- tcrossprod(Y[i, ] %*% rhoInit1[, , r, repet]
- X[i, ] %*% phiInit1[, , r, repet])
Gam[i, r] <- piInit1[repet, r] *
- det(rhoInit1[, , r, repet]) * exp(-0.5 * dotProduct)
+ gdet(rhoInit1[, , r, repet]) * exp(-0.5 * dotProduct)
}
sumGamI <- sum(Gam[i, ])
gamInit1[i, , repet] <- Gam[i, ]/sumGamI
init_EMG <- EMGLLF(phiInit1[, , , repet], rhoInit1[, , , repet], piInit1[repet, ],
gamInit1[, , repet], miniInit, maxiInit, gamma = 1, lambda = 0, X, Y,
eps = 1e-04, fast)
- LLFEessai <- init_EMG$LLF
- LLFinit1[repet] <- LLFEessai[length(LLFEessai)]
+ LLFinit1[[repet]] <- init_EMG$llh
}
b <- which.min(LLFinit1)
phiInit <- phiInit1[, , , b]