X-Git-Url: https://git.auder.net/?p=epclust.git;a=blobdiff_plain;f=pkg%2Ftests%2Ftestthat%2Ftest-assignMedoids.R;fp=pkg%2Ftests%2Ftestthat%2Ftest-assignMedoids.R;h=9b55d78d5029f716f1e56e38d0d7b7917ed8fe23;hp=0000000000000000000000000000000000000000;hb=e906736ea27105237e84c904dce6170353726292;hpb=57f337af19cd6251815bb1ff2d62f4c58e8b6078

diff --git a/pkg/tests/testthat/test-assignMedoids.R b/pkg/tests/testthat/test-assignMedoids.R
new file mode 100644
index 0000000..9b55d78
--- /dev/null
+++ b/pkg/tests/testthat/test-assignMedoids.R
@@ -0,0 +1,19 @@
+context("assignMedoids")
+
+test_that("assignMedoids behave as expected",
+{
+	# Generate a gaussian mixture
+	n <- 999
+	L <- 7
+	medoids <- cbind( rep(0,L), rep(-5,L), rep(5,L) )
+	# short series...
+	require("MASS", quietly=TRUE)
+	series <- t( rbind( MASS::mvrnorm(n/3, medoids[,1], diag(L)),
+		MASS::mvrnorm(n/3, medoids[,2], diag(L)),
+		MASS::mvrnorm(n/3, medoids[,3], diag(L)) ) )
+
+	# With high probability, medoids indices should resemble 1,1,1,...,2,2,2,...,3,3,3,...
+	mi <- assignMedoids(series, medoids)
+	mi_ref <- rep(1:3, each=n/3)
+	expect_lt( mean(mi != mi_ref), 0.01 )
+})