essai fusion

[valse.git] / pkg / R / constructionModelesLassoMLE.R

diff --git a/pkg/R/constructionModelesLassoMLE.R b/pkg/R/constructionModelesLassoMLE.R
index 3967dfc..9743f0c 100644 (file)
--- 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)