X-Git-Url: https://git.auder.net/?a=blobdiff_plain;ds=sidebyside;f=pkg%2FR%2FconstructionModelesLassoMLE.R;h=9743f0c9eb2872ba9654491c22e5f9239831e08d;hb=228ee602a972fcac6177db0d539bf9d0c5fa477f;hp=3967dfcf623892689d186b6ab7ee8979ddde35ed;hpb=ca277ac5ab51fef149014eb5e4610403fdb3227b;p=valse.git

diff --git a/pkg/R/constructionModelesLassoMLE.R b/pkg/R/constructionModelesLassoMLE.R
index 3967dfc..9743f0c 100644
--- a/pkg/R/constructionModelesLassoMLE.R
+++ b/pkg/R/constructionModelesLassoMLE.R
@@ -76,9 +76,9 @@ constructionModelesLassoMLE <- function(phiInit, rhoInit, piInit, gamInit, mini,
       })
       
       logGam <- logGam - max(logGam) #adjust without changing proportions
-      gam[i, ] <- exp(logGam)
-      norm_fact <- sum(gam[i, ])
-      gam[i, ] <- gam[i, ] / norm_fact
+      gam <- exp(logGam)
+      print(gam)
+      norm_fact <- sum(gam)
       sumLogLLH <- sumLogLLH + log(norm_fact) - log((2 * base::pi)^(m/2))
     }
     llhLambda <- c(sumLogLLH/n, (dimension + m + 1) * k - 1)