X-Git-Url: https://git.auder.net/doc/html/img/rock_paper_scissors_lizard_spock.gif?a=blobdiff_plain;f=pkg%2FR%2FcomputeGridLambda.R;h=e2b6303b69ebbe6c140c1aa9a7daeffc28082477;hb=19041906b8d80eb9a7dac7bffebf3992bcec6ccf;hp=a051441201ef7793d337dea0811be4d9a913c888;hpb=0eb161e3f3d018bce7d98fc85622d14910f89d43;p=valse.git
diff --git a/pkg/R/computeGridLambda.R b/pkg/R/computeGridLambda.R
index a051441..e2b6303 100644
--- a/pkg/R/computeGridLambda.R
+++ b/pkg/R/computeGridLambda.R
@@ -24,9 +24,8 @@ computeGridLambda = function(phiInit, rhoInit, piInit, gamInit, X, Y,
m = dim(phiInit)[2]
k = dim(phiInit)[3]
- # 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)
+ list_EMG = EMGLLF(phiInit, rhoInit, piInit, gamInit, mini, maxi,
+ gamma, lambda=0, X, Y, tau)
grid = array(0, dim=c(p,m,k))
for (i in 1:p)
{
@@ -34,6 +33,5 @@ computeGridLambda = function(phiInit, rhoInit, piInit, gamInit, X, Y,
grid[i,j,] = abs(list_EMG$S[i,j,]) / (n*list_EMG$pi^gamma)
}
grid = unique(grid)
- grid = grid[grid <= 1]
grid
}