[epclust.git] / pkg / R / simulateCurve.R

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))
	})
}
Commit	Line	Data
	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	}