'update'
[epclust.git] / TODO
CommitLineData
572d139a 1simulateur : complètement aléatoire ? from Irish dataset ?
62deb424 2wmtsa : wavBootstrap
572d139a
BA
3wavelets methods in statistics with R - p180
4
572d139a
BA
5Essayer distance wdist du package biwavelet ?
6
7geometric structure of high dim data and dim reduction 2011
1568db29 8
5edda192 9Piste à explorer pour les comparaisons: H20
e2a7484d 10
8702eb86
BA
11#https://stat.ethz.ch/pipermail/r-help/2011-June/280133.html
12#randCov = function(d)
13#{
14# x <- matrix(rnorm(d*d), nrow=d)
15# x <- x / sqrt(rowSums(x^2))
16# x %*% t(x)
17#}
492cd9e7 18
492cd9e7 19#TODO: map-reduce more appropriate R/clustering.R ligne 88
2c14dbea
BA
20
21#TODO: use dbs(),
e906736e
BA
22#https://www.r-bloggers.com/debugging-parallel-code-with-dbs/
23#http://gforge.se/2015/02/how-to-go-parallel-in-r-basics-tips/
2c14dbea
BA
24
25PLOT:
26plot manifold 2D distances WER /
27fenetre tempo forme des courbes /
28medoids /
29gain en prevision: clust puis full --> enercast
6ad3f3fd 30
e906736e 31-------
e06420d2
BA
32
33Voici le code :
34
35library(epclust)
36
37n <- 5
38N <- 128
39M <- matrix(runif(n * N), nrow = n) #séries en lignes
40M <- t(apply(M, 1, cumsum)) / sqrt(n * N)
41matplot(t(M), type = 'l', lty = 1)
42
43dists3a <- computeWerDists(1:n, function(inds) M[,inds], n)
44
45# je passe de distance Wer à corrélations : moralement tout est à 0
46summary((1 - dists3a[lower.tri(dists3a)]^2) / N / n)
47
48# or, les courbes sont très similaires, nous devrions être proche de 1 !
49corM <- cor(t(M))
50summary(corM[lower.tri(corM)])
51
b7cd987d
BA
52#si pas lissage --> on doit trouver 1 : TOCHECK
53#regarder biwavelets smooth.wavelet : s'en inspirer