X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=SelectModel%2Fselectionmodele.R;fp=SelectModel%2Fselectionmodele.R;h=b92c4e4bf00aef0871accb7a30f5c3d4e4d17b9d;hb=35b42a4bd37b162a3d579693b2b5fa4913a52ed5;hp=0000000000000000000000000000000000000000;hpb=1d3c1faaef57d906a7f12490040398b252ff049d;p=valse.git

diff --git a/SelectModel/selectionmodele.R b/SelectModel/selectionmodele.R
new file mode 100644
index 0000000..b92c4e4
--- /dev/null
+++ b/SelectModel/selectionmodele.R
@@ -0,0 +1,38 @@
+vec_bin = function(X,r){
+  Z = c()
+  indice = c()
+  j=1
+  for(i in 1:length(X)){
+    if(X[i] == r){
+      Z[i] = 1
+      indice[j] = i
+      j=j+1
+    }
+    else{
+      Z[i] = 0
+    }
+  }
+  return(list(Z,indice))
+}
+
+selectionmodele = function(vraisemblance){
+  D = vraimsemblance[,2]
+  D1 = unique(D)
+  
+  indice = rep(1, length(D1))
+  
+  #select argmax MLE
+  if(length(D1)>2){
+    for(i in 1:length(D1)){
+      A = c()
+      for(j in 1:length(D)){
+        if(D[[j]]==D1[[i]]){
+          a = c(a, vraimsemblance[j,1])
+        }
+      }
+      b = max(a)
+      indice[i] = which.max(vec_bin(vraimsemblance,b)[[1]]) #retourne le premier indice du vecteur binaire où u_i ==1
+    }
+  }
+  return(list(indice,D1))
+}
\ No newline at end of file