From: Benjamin Auder Date: Thu, 23 Feb 2017 20:29:54 +0000 (+0100) Subject: add a few TODOs in constructionModelesEMGLLF.R X-Git-Url: https://git.auder.net/css/vendor/img/scripts/getGraph_%22%20%20%20this.image%20%20%20%22.php?a=commitdiff_plain;h=621ec46ceb49d86ad3d525cfb263329f1d1faff8;p=valse.git add a few TODOs in constructionModelesEMGLLF.R --- 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) {