X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=src%2Ftest%2Fgenerate_test_data%2Fhelpers%2FconstructionModelesLassoRank.m;h=6279416b7f105a0a8d9f9d3ec0636ffdaa075e5b;hp=415ab125acce542110372f98f757c9d0eb648e1e;hb=c3bc47052f3ccb659659c59a82e9a99ea842398d;hpb=e39bc178cf5de02489ea2dce3869ba6323e18492 diff --git a/src/test/generate_test_data/helpers/constructionModelesLassoRank.m b/src/test/generate_test_data/helpers/constructionModelesLassoRank.m index 415ab12..6279416 100644 --- a/src/test/generate_test_data/helpers/constructionModelesLassoRank.m +++ b/src/test/generate_test_data/helpers/constructionModelesLassoRank.m @@ -1,4 +1,4 @@ -function[phi,lvraisemblance] = constructionModelesLassoRank(Pi,Rho,mini,maxi,X,Y,tau,A1,rangmin,rangmax) +function[phi,llh] = constructionModelesLassoRank(Pi,Rho,mini,maxi,X,Y,tau,A1,rangmin,rangmax) PI = 4.0 * atan(1.0); @@ -22,7 +22,7 @@ function[phi,lvraisemblance] = constructionModelesLassoRank(Pi,Rho,mini,maxi,X,Y %output parameters phi = zeros(p,m,k,L*Size); - lvraisemblance = zeros(L*Size,2); + llh = zeros(L*Size,2); for lambdaIndex=1:L %On ne garde que les colonnes actives %active sera l'ensemble des variables informatives @@ -31,7 +31,7 @@ function[phi,lvraisemblance] = constructionModelesLassoRank(Pi,Rho,mini,maxi,X,Y if length(active) > 0 for j=1:Size [phiLambda,LLF] = EMGrank(Pi(:,lambdaIndex),Rho(:,:,:,lambdaIndex),mini,maxi,X(:,active),Y,tau,Rank(j,:)); - lvraisemblance((lambdaIndex-1)*Size+j,:) = [LLF, sum(Rank(j,:) .* (length(active)-Rank(j,:)+m))]; + llh((lambdaIndex-1)*Size+j,:) = [LLF, sum(Rank(j,:) .* (length(active)-Rank(j,:)+m))]; phi(active,:,:,(lambdaIndex-1)*Size+j) = phiLambda; end end