sel.lambda = selected[[lambda]]
col.sel = which(colSums(sel.lambda)!=0)
res_EM = EMGLLF(phiInit[col.sel,,],rhoInit,piInit,gamInit,mini,maxi,gamma,0.,X[,col.sel],Y,tau)
sel.lambda = selected[[lambda]]
col.sel = which(colSums(sel.lambda)!=0)
res_EM = EMGLLF(phiInit[col.sel,,],rhoInit,piInit,gamInit,mini,maxi,gamma,0.,X[,col.sel],Y,tau)