X-Git-Url: https://git.auder.net/?p=valse.git;a=blobdiff_plain;f=R%2Fsuppressionmodelesegaux.R;h=a588062a5d461449e3b9d06c9d8ddf39abc2e77d;hp=a558efbe6d4f25ca876a29af53c8b002fc1abeaf;hb=39046da6016f15d625bd99cf0303ea8beb838c79;hpb=8266149c7d93aa0543cee2a1b22e1233e7b82617 diff --git a/R/suppressionmodelesegaux.R b/R/suppressionmodelesegaux.R index a558efb..a588062 100644 --- a/R/suppressionmodelesegaux.R +++ b/R/suppressionmodelesegaux.R @@ -1,18 +1,20 @@ -suppressionmodelesegaux = function(B1,B2,glambda,rho,pi){ - ind = c() - for(j in 1:length(glambda)){ - for(ll in 1:(l-1)){ - if(B1[,,l] == B1[,,ll]){ - ind = c(ind, l) - } - } - } - ind = unique(ind) - B1 = B1[,,-ind] - glambda = glambda[-ind] - B2 = B2[,,-ind] - rho = rho[,,,-ind] - pi = pi[,-ind] - - return(list(B1,B2,glambda,ind,rho,pi)) -} \ No newline at end of file +suppressionmodelesegaux = function(B1,B2,glambda,rho,pi) +{ + ind = c() + for (j in 1:length(glambda)) + { + for (ll in 1:(l-1)) + { + if(B1[,,l] == B1[,,ll]) + ind = c(ind, l) + } + } + ind = unique(ind) + B1 = B1[,,-ind] + glambda = glambda[-ind] + B2 = B2[,,-ind] + rho = rho[,,,-ind] + pi = pi[,-ind] + + return (list(B1=B1,B2=B2,glambda=glambda,ind=ind,rho=rho,pi=pi)) +}