prepare structure for R package

[valse.git] / OLD_MATLAB / SelectModel / selectionmodele.m

function [indice,D1]=selectionmodele(vraisemblance)

	D=vraisemblance(:,2);
	[D1]=unique(D);
	indice=ones(1,length(D1));
	%On ne sélectionne que celui qui maximise : l'EMV
	if length(D1)>2
		for i=1:length(D1)
			a=[];
			for j=1:length(D)
				if D(j)==D1(i)
					a=[a,vraisemblance(j,1)];
				end
			end
			b=max(a);
			indice(i)=find(vraisemblance(:,1)==b,1);
		end
	end

end