SelectModel/selectionmodele.m

   1 function [indice,D1]=selectionmodele(vraisemblance)
   2
   3     D=vraisemblance(:,2);
   4     [D1]=unique(D);
   5     indice=ones(1,length(D1));
   6     %On ne sélectionne que celui qui maximise : l'EMV
   7     if length(D1)>2
   8         for i=1:length(D1)
   9             a=[];
  10             for j=1:length(D)
  11                 if D(j)==D1(i)
  12                     a=[a,vraisemblance(j,1)];
  13                 end
  14             end
  15             b=max(a);
  16             indice(i)=find(vraisemblance(:,1)==b,1);
  17         end
  18     end
  19
  20 end