X-Git-Url: http://iramuteq.org/git?p=iramuteq;a=blobdiff_plain;f=Rscripts%2Fdistance-labbe.R;h=eccb3ae8105bbfc3ab9ace9672f4f9239964470f;hp=3ca73fe7643abdebaad200d3db4fb73f8af7327a;hb=26dbd550139e3c5c10d966cf10fe39f1d7d32654;hpb=45360b4516897d0714d5a2935deb97b4710f8643 diff --git a/Rscripts/distance-labbe.R b/Rscripts/distance-labbe.R index 3ca73fe..eccb3ae 100644 --- a/Rscripts/distance-labbe.R +++ b/Rscripts/distance-labbe.R @@ -23,11 +23,14 @@ compute.labbe <- function(x, y, tab) { if (plus.grand == 1) { U <- N2/N1 mini.tab[,1] <- mini.tab[,1] * U + col.plusgrand <- mini.tab[,1] + cs.plus.grand <- sum(col.plusgrand[col.plusgrand>=1]) } else { U <- N1/N2 mini.tab[,2] <- mini.tab[,2] * U + col.plusgrand <- mini.tab[,2] + cs.plus.grand <- sum(col.plusgrand[col.plusgrand>1]) } - commun <- which((mini.tab[,1] > 0) & (mini.tab[,2] > 0)) deA <- which((mini.tab[,plus.petit] > 0) & (mini.tab[,plus.grand] == 0)) deB <- which((mini.tab[,plus.petit] == 0) & (mini.tab[,plus.grand] >= 1)) @@ -37,7 +40,7 @@ compute.labbe <- function(x, y, tab) { dist.deB <- abs(mini.tab[deB, plus.petit] - mini.tab[deB, plus.grand]) dist.labbe <- sum(dist.commun) + sum(dist.deA) + sum(dist.deB) - indice.labbe <- dist.labbe/(cs[plus.petit] + sum(mini.tab[,plus.grand])) + indice.labbe <- dist.labbe/(cs[plus.petit] + cs.plus.grand) indice.labbe }