X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=pkg%2FR%2FconstructionModelesLassoMLE.R;h=0197e1a32d785a74efeda4efa7e6b666deb389e6;hb=57a8adeb0e558a54f9481e0256e258525aeb8cd1;hp=3967dfcf623892689d186b6ab7ee8979ddde35ed;hpb=ca277ac5ab51fef149014eb5e4610403fdb3227b;p=valse.git

diff --git a/pkg/R/constructionModelesLassoMLE.R b/pkg/R/constructionModelesLassoMLE.R
index 3967dfc..0197e1a 100644
--- a/pkg/R/constructionModelesLassoMLE.R
+++ b/pkg/R/constructionModelesLassoMLE.R
@@ -76,12 +76,11 @@ 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)
+      norm_fact <- sum(gam)
       sumLogLLH <- sumLogLLH + log(norm_fact) - log((2 * base::pi)^(m/2))
     }
-    llhLambda <- c(sumLogLLH/n, (dimension + m + 1) * k - 1)
+    llhLambda <- c(-sumLogLLH/n, (dimension + m + 1) * k - 1)
     # densite <- vector("double", n)
     # for (r in 1:k)
     # {