prepare structure for R package
[valse.git] / R / suppressionmodelesegaux2.R
CommitLineData
35b42a4b
B
1suppressionmodelesegaux2 = 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}