- res <- EMGLLF(phiInit[col.sel, , ], rhoInit, piInit, gamInit, mini, maxi,
- gamma, 0, X[, col.sel], Y, eps, fast)
-
+ res <- EMGLLF(array(phiInit,dim=c(p,m,k))[col.sel, , ], rhoInit, piInit, gamInit,
+ mini, maxi, gamma, 0, as.matrix(X[, col.sel]), Y, eps, fast)
+