Inputparametters.R + SelecModel

[valse.git] / SelectModel / selectionindice.R

diff --git a/SelectModel/selectionindice.R b/SelectModel/selectionindice.R
new file mode 100644 (file)
index 0000000..6782b2e
--- /dev/null
+++ b/SelectModel/selectionindice.R
@@ -0,0 +1,24 @@
+selectionindice = function(phi, seuil){
+  dim_phi = dim(phi)
+  pp = dim_phi[1]
+  m = dim_phi[2]
+  
+  A = matrix(0, pp, m)
+  B = matrix(0, pp, m)
+  
+  for(j in 1:pp){
+    cpt1 = 0
+    cpt2 = 0
+    for(mm in 1:m){
+      if(max(phi[j,mm,])> seuil){
+        cpt1 = cpt1 + 1
+        A[j,cpt] = mm
+      }
+      else{
+        cpt2 = cpt2+1
+        B[j, cpt2] = mm
+      }
+    }
+  }
+  return(list(A,B))
+}
\ No newline at end of file