| 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 | } |