'update'
[epclust.git] / pkg / R / simulateCurve.R
1 Ref_Centroids <<- NULL
2 Tot_Nb_Curves <<- NA
3
4 #' initCentroids
5 #'
6 #' Must be called before getDetRandCurve()
7 #'
8 #' @param N Number of curves
9 #' @param D Number of sample points
10 #' @param K Number of clusters
11 #'
12 #' @export
13 initClustersParams <- function(N, D, K)
14 {
15 # Generate K centroids
16 Ref_Centroids <<- sapply(1:K, function(k) cumsum(rnorm(D)))
17 Tot_Nb_Curves <<- N
18 }
19
20 #' getDetRandCurve
21 #'
22 #' TODO: improve this (using known clusters centers...)
23 #'
24 #' @param i Some index
25 #'
26 #' @examples
27 #' \dontrun{
28 #' initCentroids(N=250000, D=15000, K=15)
29 #' res_ascii <- claws(getDetRandCurve, K1=50, K2=15, nb_series_per_chunk=500,
30 #' nb_items_clust=100, random=FALSE, verbose=TRUE, ncores_clust=3)
31 #' }
32 #' @export
33 getDetRandCurve <- function(indices)
34 {
35 sapply(indices, function(i) {
36 if (i > Tot_Nb_Curves)
37 return (NULL)
38 set.seed(i)
39 j <- sample(ncol(Ref_Centroids), 1)
40 Ref_Centroids[,j] + rnorm(nrow(Ref_Centroids))
41 })
42 }