R/suppressionmodelesegaux.R

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