X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=R%2FconstructionModelesLassoMLE.R;fp=R%2FconstructionModelesLassoMLE.R;h=346339ed836a7e9e2306737bdf427f62041f986a;hp=ed05b2a1c625a1735c5f99b031664bf5241e61ac;hb=621ec46ceb49d86ad3d525cfb263329f1d1faff8;hpb=5955cc25fade323296485500104dd90990e2cdc1

diff --git a/R/constructionModelesLassoMLE.R b/R/constructionModelesLassoMLE.R
index ed05b2a..346339e 100644
--- a/R/constructionModelesLassoMLE.R
+++ b/R/constructionModelesLassoMLE.R
@@ -22,13 +22,14 @@ constructionModelesLassoMLE = function(phiInit,rhoInit,piInit,gamInit,mini,maxi,
 
     res = EMGLLF(phiInit[a,,],rhoInit,piInit,gamInit,mini,maxi,gamma,0.,X[,a],Y,tau)
 
+		#TODO: supprimer ça et utiliser parLapply(...)
     for (j in 1:length(a))
       phi[a[j],,,lambdaIndex] = res$phi[j,,]
     rho[,,,lambdaIndex] = res$rho
     pi[,lambdaIndex] = res$pi
 
     dimension = 0
-    for (j in 1:p)
+    for (j in 1:p) #TODO: doit pouvoir être fait beaucoup plus simplement
 		{
       b = A2[j,2:dim(A2)[2],lambdaIndex]
       b = b[b!=0]
@@ -40,7 +41,7 @@ constructionModelesLassoMLE = function(phiInit,rhoInit,piInit,gamInit,mini,maxi,
 
     #on veut calculer l'EMV avec toutes nos estimations
 		densite = matrix(0, nrow=n, ncol=L)
-		for (i in 1:n)
+		for (i in 1:n) #TODO: pas besoin de cette boucle (vectorize)
 		{
 			for (r in 1:k)
 			{