prepare structure for R package
[valse.git] / OLD_MATLAB / InputParameters / grillelambda.m
CommitLineData
1d3c1faa
BA
1function[gridLambda] = grillelambda(phiInit,rhoInit,piInit,gamInit,X,Y,gamma,mini,maxi,tau)
2
3 n = size(X, 1);
4 [p,m,k] = size(phiInit);
5 [phi,rho,pi,~,S] = EMGLLF(phiInit,rhoInit,piInit,gamInit,mini,maxi,1,0,X,Y,tau);
6
7 gridLambda = zeros(p,m,k);
8 for j=1:p
9 for mm=1:m
10 gridLambda(j,mm,:) = abs(reshape(S(j,mm,:),[k,1])) ./ (n*pi.^gamma);
11 end
12 end
13
14 gridLambda = unique(gridLambda);
15 gridLambda(gridLambda()>1) = [];
16
17end