X-Git-Url: https://git.auder.net/variants/Cwda/current/git-favicon.png?a=blobdiff_plain;ds=sidebyside;f=src%2Fsources%2FconstructionModelesLassoMLE.c;h=34e5808e37b6dda447190399d282266957251aa3;hb=46a2e676b2d85eef5a1811a6e623b65327fc451d;hp=ad2a71821b94a23cea46dd6ea66efc9c95358dcc;hpb=ef67d338c7f28ba041abe40ca9a8ab69f8365a90;p=valse.git diff --git a/src/sources/constructionModelesLassoMLE.c b/src/sources/constructionModelesLassoMLE.c index ad2a718..34e5808 100644 --- a/src/sources/constructionModelesLassoMLE.c +++ b/src/sources/constructionModelesLassoMLE.c @@ -13,7 +13,8 @@ void constructionModelesLassoMLE_core( const Real* gamInit, // paramètre initial des probabilités a posteriori de chaque échantillon int mini,// nombre minimal d'itérations dans l'algorithme EM int maxi,// nombre maximal d'itérations dans l'algorithme EM - Real gamma,// valeur de gamma : puissance des proportions dans la pénalisation pour un Lasso adaptatif + Real gamma,// valeur de gamma : puissance des proportions dans la pénalisation + //pour un Lasso adaptatif const Real* glambda, // valeur des paramètres de régularisation du Lasso const Real* X, // régresseurs const Real* Y, // réponse @@ -33,9 +34,15 @@ void constructionModelesLassoMLE_core( int k, // nombre de composantes int L) // taille de glambda { - //preparation: phi = 0 + //preparation: phi,rho,pi = 0, llh=+Inf for (int u=0; u
0
- //~ phi(A2(j,1,lambdaIndex),b,:,lambdaIndex) = 0.0;
- //~ end
if (lengthB > 0)
{
+ //phi[A2[j,1,lambdaIndex],b,,lambdaIndex] = 0.
for (int mm=0; mm