[valse.git] / R / suppressionmodelesegaux2.R

suppressionmodelesegaux2 = function(B1,rho,pi){
  ind = c()
  dim_B1 = dim(B1)
  B2 = array(0,dim=c(dim_B1[1],dim_B1[2],dim_B1[3]))
  nombreLambda=dim_B1[[2]]
  glambda = rep(0,nombreLambda)
  
  #for(j in 1:nombreLambda){
  #  for(ll in 1:(l-1)){
  #    if(B1[,,l] == B1[,,ll]){
  #      ind = c(ind, l)
  #    }
  #  }
  #}
  #ind = unique(ind)
  #B1 = B1[,,-ind]
  #rho = rho[,,,-ind] 
  #pi = pi[,-ind]
  
  suppressmodel = suppressionmodelesegaux(B1,B2,glambda,rho,pi)
  B1 = suppressmodel[[1]]
  ind = suppressmodel[[4]]
  rho = suppressmodel[[5]]
  pi = suppressmodel[[6]]
  return(list(B1,ind,rho,pi))
}
Commit	Line	Data
35b42a4b B	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	}