'update'
[valse.git] / R / suppressionmodelesegaux.R
CommitLineData
39046da6
BA
1suppressionmodelesegaux = 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}