# lambda == 0 because we compute the EMV: no penalization here
res <- EMGLLF(phiInit[col.sel, , ], rhoInit, piInit, gamInit, mini, maxi,
gamma, 0, X[, col.sel], Y, eps, fast)
# lambda == 0 because we compute the EMV: no penalization here
res <- EMGLLF(phiInit[col.sel, , ], rhoInit, piInit, gamInit, mini, maxi,
gamma, 0, X[, col.sel], Y, eps, fast)