+++ /dev/null
-# R-equivalent of computeMedoidsIndices, requiring a matrix
-# (thus potentially breaking "fit-in-memory" hope)
-R_computeMedoidsIndices <- function(medoids, series)
-{
- nb_series = ncol(series) #series in columns
-
- mi = rep(NA,nb_series)
- for (i in 1:nb_series)
- mi[i] <- which.min( colSums( sweep(medoids, 1, series[,i], '-')^2 ) )
-
- mi
-}