code seems OK; still wavelets test to write
[epclust.git] / epclust / tests / testthat / helper.clustering.R
index 273a0b7..f39257e 100644 (file)
@@ -9,10 +9,10 @@ computeDistortion = function(series, medoids)
        if (bigmemory::is.big.matrix(medoids))
                medoids = medoids[,] #extract standard matrix
 
-       n = nrow(series) ; L = ncol(series)
+       n = ncol(series) ; L = nrow(series)
        distortion = 0.
        for (i in seq_len(n))
                distortion = distortion + min( colSums( sweep(medoids,1,series[,i],'-')^2 ) / L )
 
-       distortion / n
+       sqrt( distortion / n )
 }