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

diff --git a/SelectModel/suppressionmodelesegaux2.R b/SelectModel/suppressionmodelesegaux2.R
new file mode 100644
index 0000000..4574afc
--- /dev/null
+++ b/SelectModel/suppressionmodelesegaux2.R
@@ -0,0 +1,26 @@
+suppressionmodelesegaux2 = function(B1,rho,pi){
+  ind = c()
+  dim_B1 = dim(B1)
+  B2 = array(0,dim=c(dim_B1[1],dim_B1[2],dim_B1[3]))
+  nombreLambda=dim_B1[[2]]
+  glambda = rep(0,nombreLambda)
+  
+  #for(j in 1:nombreLambda){
+  #  for(ll in 1:(l-1)){
+  #    if(B1[,,l] == B1[,,ll]){
+  #      ind = c(ind, l)
+  #    }
+  #  }
+  #}
+  #ind = unique(ind)
+  #B1 = B1[,,-ind]
+  #rho = rho[,,,-ind] 
+  #pi = pi[,-ind]
+  
+  suppressmodel = suppressionmodelesegaux(B1,B2,glambda,rho,pi)
+  B1 = suppressmodel[[1]]
+  ind = suppressmodel[[4]]
+  rho = suppressmodel[[5]]
+  pi = suppressmodel[[6]]
+  return(list(B1,ind,rho,pi))
+}
\ No newline at end of file