projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
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
})
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)
sumLogLLH <- sumLogLLH + log(norm_fact) - log((2 * base::pi)^(m/2))
}
llhLambda <- c(sumLogLLH/n, (dimension + m + 1) * k - 1)