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 )
}