X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=R%2FconstructionModelesLassoMLE.R;h=50879c935aea7e04042024e2f51a071cb554fa7f;hp=55e9419e7d89b5dd43a8917a7cfe82f4faa8d144;hb=7411013519fc6aef07280502d7d201e5ca177e0c;hpb=c7dab9ff8b95a7630c7dafdcf40d60c659290ef2 diff --git a/R/constructionModelesLassoMLE.R b/R/constructionModelesLassoMLE.R index 55e9419..50879c9 100644 --- a/R/constructionModelesLassoMLE.R +++ b/R/constructionModelesLassoMLE.R @@ -15,13 +15,13 @@ constructionModelesLassoMLE = function(phiInit,rhoInit,piInit,gamInit,mini,maxi, k = dim(phiInit)[3] #TODO: phiInit[selected] et X[selected] sont bien sûr faux; par quoi remplacer ? - #lambda == 0 c'est normal ? + #lambda == 0 c'est normal ? -> ED : oui, ici on calcule le maximum de vraisembance, donc on ne pénalise plus res = EMGLLF(phiInit[selected],rhoInit,piInit,gamInit,mini,maxi,gamma,0.,X[selected],Y,tau) #comment évaluer la dimension à partir du résultat et de [not]selected ? #dimension = ... - #on veut calculer l'EMV avec toutes nos estimations + #on veut calculer la vraisemblance avec toutes nos estimations densite = vector("double",n) for (r in 1:k) {