X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2FinitSmallEM.R;h=937ea734aa427710540b69f75585f75eb6968795;hp=7e9cce57b6c2cad973db5dc1f794ac764857b7d4;hb=6775f5b98ffc7eae7ce9d4081b23b39ce66d3c0b;hpb=0ba1b11c49d7b2a0cae493200793c1ba3fb8b8e7 diff --git a/pkg/R/initSmallEM.R b/pkg/R/initSmallEM.R index 7e9cce5..937ea73 100644 --- a/pkg/R/initSmallEM.R +++ b/pkg/R/initSmallEM.R @@ -34,7 +34,7 @@ initSmallEM <- function(k, X, Y, fast) for (r in 1:k) { Z <- Zinit1[, repet] - Z_indice <- seq_len(n)[Z == r] #renvoit les indices où Z==r + Z_indice <- seq_len(n)[Z == r] #renvoit les indices ou Z==r if (length(Z_indice) == 1) { betaInit1[, , r, repet] <- MASS::ginv(crossprod(t(X[Z_indice, ]))) %*% crossprod(t(X[Z_indice, ]), Y[Z_indice, ])