# initialisation for the allocations probabilities in each component
tauInit,
# values for the regularization parameter grid
- gridLambda = [];
+ gridLambda = c(),
# je ne crois pas vraiment qu'il faille les mettre en sortie, d'autant plus qu'on construit
# une matrice A1 et A2 pour chaque k, et elles sont grandes, donc ca coute un peu cher ...
A1,
Pi,
#immutable
- seuil = 1e-15;
+ seuil = 1e-15
),
methods = list(
"computation of the regularization grid"
#(according to explicit formula given by EM algorithm)
- gridLambda <<- grillelambda(sx.phiInit,sx.rhoInit,sx.piInit,sx.tauInit,sx.X,sx.Y,
+ gridLambda <<- gridLambda(sx.phiInit,sx.rhoInit,sx.piInit,sx.tauInit,sx.X,sx.Y,
sx.gamma,sx.mini,sx.maxi,sx.eps);
},