1 function[] = generateRunSaveTest_EMGrank(n, p, m, k, mini, maxi, gamma, rank, varargin)
3 %set defaults for optional inputs
4 optargs = {200 15 10 3 5 10 1.0 1:3};
5 %replace defaults by user parameters
6 optargs(1:length(varargin)) = varargin;
7 [n, p, m, k, mini, maxi, gamma, rank] = optargs{:};
13 Pi = (1.0/k)*ones(1,k);
20 [X, Y, ~] = generateIOdefault(n, p, m, k);
22 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
29 dlmwrite(strcat(testFolder,'Rho'), reshape(Rho,1,[]), delimiter);
30 dlmwrite(strcat(testFolder,'Pi'), Pi, delimiter);
31 dlmwrite(strcat(testFolder,'mini'), mini, delimiter);
32 dlmwrite(strcat(testFolder,'maxi'), maxi, delimiter);
33 dlmwrite(strcat(testFolder,'X'), reshape(X,1,[]), delimiter);
34 dlmwrite(strcat(testFolder,'Y'), reshape(Y,1,[]), delimiter);
35 dlmwrite(strcat(testFolder,'tau'), tau, delimiter);
36 dlmwrite(strcat(testFolder,'rank'), rank, delimiter);
37 dlmwrite(strcat(testFolder,'dimensions'), [n,p,m,k], delimiter);;
39 [phi,LLF] = EMGrank(Pi,Rho,mini,maxi,X,Y,tau,rank);
42 dlmwrite(strcat(testFolder,'phi'), reshape(phi,1,[]), delimiter);
43 dlmwrite(strcat(testFolder,'LLF'), LLF, delimiter);