R/suppressionmodelesegaux.R

   1 suppressionmodelesegaux = function(B1,B2,glambda,rho,pi){
   2   ind = c()
   3   for(j in 1:length(glambda)){
   4     for(ll in 1:(l-1)){
   5       if(B1[,,l] == B1[,,ll]){
   6         ind = c(ind, l)
   7       }
   8     }
   9   }
  10   ind = unique(ind)
  11   B1 = B1[,,-ind]
  12   glambda = glambda[-ind]
  13   B2 = B2[,,-ind]
  14   rho = rho[,,,-ind]
  15   pi = pi[,-ind]
  16
  17   return(list(B1,B2,glambda,ind,rho,pi))
  18 }