R/suppressionmodelesegaux2.R

   1 suppressionmodelesegaux2 = function(B1,rho,pi){
   2   ind = c()
   3   dim_B1 = dim(B1)
   4   B2 = array(0,dim=c(dim_B1[1],dim_B1[2],dim_B1[3]))
   5   nombreLambda=dim_B1[[2]]
   6   glambda = rep(0,nombreLambda)
   7
   8   #for(j in 1:nombreLambda){
   9   #  for(ll in 1:(l-1)){
  10   #    if(B1[,,l] == B1[,,ll]){
  11   #      ind = c(ind, l)
  12   #    }
  13   #  }
  14   #}
  15   #ind = unique(ind)
  16   #B1 = B1[,,-ind]
  17   #rho = rho[,,,-ind]
  18   #pi = pi[,-ind]
  19
  20   suppressmodel = suppressionmodelesegaux(B1,B2,glambda,rho,pi)
  21   B1 = suppressmodel[[1]]
  22   ind = suppressmodel[[4]]
  23   rho = suppressmodel[[5]]
  24   pi = suppressmodel[[6]]
  25   return(list(B1,ind,rho,pi))
  26 }