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

diff --git a/src/test/OLD_TEST_MATLAB/testConstructionModelesLassoMLE.m b/src/test/OLD_TEST_MATLAB/testConstructionModelesLassoMLE.m
new file mode 100644
index 0000000..27c1208
--- /dev/null
+++ b/src/test/OLD_TEST_MATLAB/testConstructionModelesLassoMLE.m
@@ -0,0 +1,46 @@
+function[] = testConstructionModelesLassoMLE()
+
+	testFolder = 'data/';
+	delimiter = '\n';
+
+	%get dimensions
+	dimensions = dlmread(strcat(testFolder,'dimensions'), delimiter);
+	n = dimensions(1);
+	p = dimensions(2);
+	m = dimensions(3);
+	k = dimensions(4);
+	L = dimensions(5);
+
+	%get all input arrays
+	phiInit = reshape(dlmread(strcat(testFolder,'phiInit'), delimiter), p, m, k);
+	rhoInit = reshape(dlmread(strcat(testFolder,'rhoInit'), delimiter), m, m, k);
+	piInit = transpose(dlmread(strcat(testFolder,'piInit'), delimiter));
+	gamInit = reshape(dlmread(strcat(testFolder,'gamInit'), delimiter), n, k);
+	mini = int64(dlmread(strcat(testFolder,'mini'), delimiter));
+	maxi = int64(dlmread(strcat(testFolder,'maxi'), delimiter));
+	gamma = dlmread(strcat(testFolder,'gamma'), delimiter);
+	glambda = dlmread(strcat(testFolder,'glambda'), delimiter);
+	X = reshape(dlmread(strcat(testFolder,'X'), delimiter), n, p);
+	Y = reshape(dlmread(strcat(testFolder,'Y'), delimiter), n, m);
+	seuil = dlmread(strcat(testFolder,'seuil'), delimiter);
+	tau = dlmread(strcat(testFolder,'tau'), delimiter);
+	A1 = int64(reshape(dlmread(strcat(testFolder,'A1'), delimiter), p, m+1, L));
+	A2 = int64(reshape(dlmread(strcat(testFolder,'A2'), delimiter), p, m+1, L));
+
+	%run constructionModelesLassoMLE.m
+	[phi,rho,pi,lvraisemblance] = constructionModelesLassoMLE(phiInit,rhoInit,piInit,gamInit,mini,maxi,gamma,glambda,X,Y,seuil,tau,A1,A2);
+
+	%get all stored outputs
+	ref_phi = reshape(dlmread(strcat(testFolder,'phi'), delimiter), p, m, k, L);
+	ref_rho = reshape(dlmread(strcat(testFolder,'rho'), delimiter), m, m, k, L);
+	ref_pi = reshape(dlmread(strcat(testFolder,'pi'), delimiter), k, L);
+	ref_lvraisemblance = reshape(dlmread(strcat(testFolder,'lvraisemblance'), delimiter), L, 2);
+
+	%check that output correspond to stored output
+	tol = 1e-5;
+	checkOutput('phi',phi,ref_phi,tol);
+	checkOutput('rho',rho,ref_rho,tol);
+	checkOutput('pi',pi,ref_pi,tol);
+	checkOutput('lvraisemblance',lvraisemblance,ref_lvraisemblance,tol);
+
+end