R/suppressionmodelesegaux2.R

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