merge with remote

[epclust.git] / pkg / tests / testthat / test-assignMedoids.R

diff --git a/pkg/tests/testthat/test-assignMedoids.R b/pkg/tests/testthat/test-assignMedoids.R
new file mode 100644 (file)
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 )
+})