X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=pkg%2FR%2FconstructionModelesLassoMLE.R;h=c55035e8540da188b8d19147a493fedb99529af3;hp=3967dfcf623892689d186b6ab7ee8979ddde35ed;hb=465c0e07fbf8363a2625f26681fa7ba31bad82b7;hpb=ca277ac5ab51fef149014eb5e4610403fdb3227b diff --git a/pkg/R/constructionModelesLassoMLE.R b/pkg/R/constructionModelesLassoMLE.R index 3967dfc..c55035e 100644 --- a/pkg/R/constructionModelesLassoMLE.R +++ b/pkg/R/constructionModelesLassoMLE.R @@ -76,9 +76,8 @@ 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)