X-Git-Url: https://git.auder.net/?a=blobdiff_plain;f=TODO;h=c5f4081dd34ea2e61f684bfe3639e4a8c9e63617;hb=e06420d26899e68416f782e83ed586e9e8f0b6e6;hp=53b4c97bc958969fdbfb4593dab3f13a388aadde;hpb=a174b8ea1f322992068ab42810df017a2b9620ee;p=epclust.git diff --git a/TODO b/TODO index 53b4c97..c5f4081 100644 --- a/TODO +++ b/TODO @@ -67,3 +67,38 @@ cwt : trim R part // : clever by rows retenir cwt... Stockage matrices : en colonnes systématiquement ? + +TODO: revoir les arguments, simplifier (dans les clustering...), + permettre algos de clustering quelconques, args: medoids (curves puis dists), K + +Plutôt ça que hack auto notebook generation: + https://mail.scipy.org/pipermail/ipython-dev/2014-March/013550.html + +#si pas lissage --> on doit trouver 1 : TOCHECK +#regarder biwavelets smooth.wavelet : c'est comme ça qu'il faut coder +#Options : smooth_lvl + filtre gaussien (smooth_method ?) +(this file : +https://github.com/tgouhier/biwavelet/blob/master/R/smooth.wavelet.R ) + +===== + +Voici le code : + +library(epclust) + +n <- 5 +N <- 128 +M <- matrix(runif(n * N), nrow = n) #séries en lignes +M <- t(apply(M, 1, cumsum)) / sqrt(n * N) +matplot(t(M), type = 'l', lty = 1) + +dists3a <- computeWerDists(1:n, function(inds) M[,inds], n) + +# je passe de distance Wer à corrélations : moralement tout est à 0 +summary((1 - dists3a[lower.tri(dists3a)]^2) / N / n) + +# or, les courbes sont très similaires, nous devrions être proche de 1 ! +corM <- cor(t(M)) +summary(corM[lower.tri(corM)]) + +---------> à checker