- # TODO: explain why gamma=1 instad of just 'gamma'?
- list_EMG = EMGLLF(phiInit, rhoInit, piInit, gamInit, mini, maxi,
- gamma=1, lamba=0, X, Y, tau)
- grid = array(0, dim=c(p,m,k))
- for (i in 1:p)
- {
- for (j in 1:m)
- grid[i,j,] = abs(list_EMG$S[i,j,]) / (n*list_EMG$pi^gamma)
- }
- grid = unique(grid)
- grid = grid[grid <= 1]
- grid
+ list_EMG <- EMGLLF(phiInit, rhoInit, piInit, gamInit, mini, maxi, gamma, lambda = 0,
+ X, Y, tau, fast)
+ grid <- array(0, dim = c(p, m, k))
+ for (j in 1:p)
+ {
+ for (mm in 1:m)
+ grid[j, mm, ] <- abs(list_EMG$S[j, mm, ])/(n * list_EMG$pi^gamma)
+ }
+ sort(unique(grid))