prepare structure for R package
[valse.git] / OLD_MATLAB / SelectModel / selectionmodele.m
diff --git a/OLD_MATLAB/SelectModel/selectionmodele.m b/OLD_MATLAB/SelectModel/selectionmodele.m
new file mode 100644 (file)
index 0000000..33828ab
--- /dev/null
@@ -0,0 +1,20 @@
+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