EMG_rank = EMGrank(Pi[,lambdaIndex], Rho[,,,lambdaIndex], mini, maxi, X[, active], Y, tau, Rank[j,])
phiLambda = EMG_rank$phi
LLF = EMG_rank$LLF
EMG_rank = EMGrank(Pi[,lambdaIndex], Rho[,,,lambdaIndex], mini, maxi, X[, active], Y, tau, Rank[j,])
phiLambda = EMG_rank$phi
LLF = EMG_rank$LLF