X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=src%2Ftest%2FOLD_TEST_MATLAB%2FgenerateRunSaveTest_constructionModelesLassoMLE.m;fp=src%2Ftest%2FOLD_TEST_MATLAB%2FgenerateRunSaveTest_constructionModelesLassoMLE.m;h=6e48d451945a3c592bf621524d38ff632c3723c9;hb=7b2720733e9aebe177c211119a9ec160c7e7117b;hp=0000000000000000000000000000000000000000;hpb=4725af564bfdbeb40d0caf7f124ab94ac05a97eb;p=valse.git

diff --git a/src/test/OLD_TEST_MATLAB/generateRunSaveTest_constructionModelesLassoMLE.m b/src/test/OLD_TEST_MATLAB/generateRunSaveTest_constructionModelesLassoMLE.m
new file mode 100644
index 0000000..6e48d45
--- /dev/null
+++ b/src/test/OLD_TEST_MATLAB/generateRunSaveTest_constructionModelesLassoMLE.m
@@ -0,0 +1,62 @@
+function[] = generateRunSaveTest_constructionModelesLassoMLE(n, p, m, k, mini, maxi, gamma, glambda, varargin)
+	
+	%set defaults for optional inputs
+	optargs = {200 15 10 3 5 10 1.0 [0.0,0.01,0.02,0.03,0.05,0.1,0.2,0.3,0.5,0.7,0.85,0.99]};
+	%replace defaults by user parameters
+	optargs(1:length(varargin)) = varargin;
+	[n, p, m, k, mini, maxi, gamma, glambda] = optargs{:};
+	tau = 1e-6;
+	seuil = 1e-15;
+	mini = int64(mini);
+	maxi = int64(maxi);
+	L = length(glambda);
+
+	%Generate phiInit,piInit,...
+	[phiInit,rhoInit,piInit,gamInit] = basicInitParameters(n, p, m, k);
+
+	%Generate X and Y
+	[X, Y, ~] = generateIOdefault(n, p, m, k);
+
+	A2 = zeros(p,m+1,L,'int64');
+    for i=1:L
+        A2(:,1,i) = 1:p;
+        A2(1:5,2,i) = 1:5;
+    end
+	A1 = zeros(p,m+1,L,'int64');
+    for i=1:L
+        A1(:,1,i) = 1:p;
+        A1(1:5,2,i) = 1:5;
+    end
+    
+	%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
+
+	testFolder = 'data/';
+	mkdir(testFolder);
+	delimiter = '\n';
+
+	%save inputs
+	dlmwrite(strcat(testFolder,'phiInit'), reshape(phiInit,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'rhoInit'), reshape(rhoInit,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'piInit'), piInit, delimiter);
+	dlmwrite(strcat(testFolder,'gamInit'), reshape(gamInit,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'mini'), mini, delimiter);
+	dlmwrite(strcat(testFolder,'maxi'), maxi, delimiter);
+	dlmwrite(strcat(testFolder,'gamma'), gamma, delimiter);
+	dlmwrite(strcat(testFolder,'glambda'), glambda, delimiter);
+	dlmwrite(strcat(testFolder,'X'), reshape(X,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'Y'), reshape(Y,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'seuil'), seuil, delimiter);
+	dlmwrite(strcat(testFolder,'tau'), tau, delimiter);
+	dlmwrite(strcat(testFolder,'A1'), reshape(A1,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'A2'), reshape(A2,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'dimensions'), [n,p,m,k,L], delimiter);
+
+	[phi,rho,pi,lvraisemblance] = constructionModelesLassoMLE(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambda,X,Y,seuil,tau,A1,A2);
+
+	%save output
+	dlmwrite(strcat(testFolder,'phi'), reshape(phi,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'rho'), reshape(rho,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'pi'), reshape(pi,1,[]), delimiter);
+	dlmwrite(strcat(testFolder,'lvraisemblance'), reshape(lvraisemblance,1,[]), delimiter);
+
+end