Commit | Line | Data |
---|---|---|
572d139a | 1 | simulateur : complètement aléatoire ? from Irish dataset ? |
62deb424 | 2 | wmtsa : wavBootstrap |
572d139a BA |
3 | wavelets methods in statistics with R - p180 |
4 | ||
572d139a BA |
5 | Essayer distance wdist du package biwavelet ? |
6 | ||
7 | geometric structure of high dim data and dim reduction 2011 | |
1568db29 | 8 | |
5edda192 | 9 | Piste à 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 | |
25 | PLOT: | |
26 | plot manifold 2D distances WER / | |
27 | fenetre tempo forme des courbes / | |
28 | medoids / | |
29 | gain en prevision: clust puis full --> enercast | |
6ad3f3fd | 30 | |
e906736e | 31 | ------- |
e06420d2 BA |
32 | |
33 | Voici le code : | |
34 | ||
35 | library(epclust) | |
36 | ||
37 | n <- 5 | |
38 | N <- 128 | |
39 | M <- matrix(runif(n * N), nrow = n) #séries en lignes | |
40 | M <- t(apply(M, 1, cumsum)) / sqrt(n * N) | |
41 | matplot(t(M), type = 'l', lty = 1) | |
42 | ||
43 | dists3a <- computeWerDists(1:n, function(inds) M[,inds], n) | |
44 | ||
45 | # je passe de distance Wer à corrélations : moralement tout est à 0 | |
46 | summary((1 - dists3a[lower.tri(dists3a)]^2) / N / n) | |
47 | ||
48 | # or, les courbes sont très similaires, nous devrions être proche de 1 ! | |
49 | corM <- cor(t(M)) | |
50 | summary(corM[lower.tri(corM)]) | |
51 | ||
b7cd987d BA |
52 | #si pas lissage --> on doit trouver 1 : TOCHECK |
53 | #regarder biwavelets smooth.wavelet : s'en inspirer |