X-Git-Url: https://git.auder.net/?a=blobdiff_plain;ds=sidebyside;f=pkg%2FR%2FsimulateCurve.R;fp=pkg%2FR%2FsimulateCurve.R;h=97aeac801ec030e9c01ccd1fc990bb689ab15a71;hb=e906736ea27105237e84c904dce6170353726292;hp=0000000000000000000000000000000000000000;hpb=57f337af19cd6251815bb1ff2d62f4c58e8b6078;p=epclust.git

diff --git a/pkg/R/simulateCurve.R b/pkg/R/simulateCurve.R
new file mode 100644
index 0000000..97aeac8
--- /dev/null
+++ b/pkg/R/simulateCurve.R
@@ -0,0 +1,42 @@
+Ref_Centroids <<- NULL
+Tot_Nb_Curves <<- NA
+
+#' initCentroids
+#'
+#' Must be called before getDetRandCurve()
+#'
+#' @param N Number of curves
+#' @param D Number of sample points
+#' @param K Number of clusters
+#'
+#' @export
+initClustersParams <- function(N, D, K)
+{
+	# Generate K centroids
+	Ref_Centroids <<- sapply(1:K, function(k) cumsum(rnorm(D)))
+	Tot_Nb_Curves <<- N
+}
+
+#' getDetRandCurve
+#'
+#' TODO: improve this (using known clusters centers...)
+#'
+#' @param i Some index
+#'
+#' @examples
+#' \dontrun{
+#' initCentroids(N=250000, D=15000, K=15)
+#' res_ascii <- claws(getDetRandCurve, K1=50, K2=15, nb_series_per_chunk=500,
+#'   nb_items_clust=100, random=FALSE, verbose=TRUE, ncores_clust=3)
+#' }
+#' @export
+getDetRandCurve <- function(indices)
+{
+	sapply(indices, function(i) {
+		if (i > Tot_Nb_Curves)
+			return (NULL)
+		set.seed(i)
+		j <- sample(ncol(Ref_Centroids), 1)
+		Ref_Centroids[,j] + rnorm(nrow(Ref_Centroids))
+	})
+}