epclust/tests/testthat/test-assignMedoids.R

   1 context("assignMedoids")
   2
   3 test_that("assignMedoids behave as expected",
   4 {
   5     # Generate a gaussian mixture
   6     n <- 999
   7     L <- 7
   8     medoids <- cbind( rep(0,L), rep(-5,L), rep(5,L) )
   9     # short series...
  10     require("MASS", quietly=TRUE)
  11     series <- t( rbind( MASS::mvrnorm(n/3, medoids[,1], diag(L)),
  12         MASS::mvrnorm(n/3, medoids[,2], diag(L)),
  13         MASS::mvrnorm(n/3, medoids[,3], diag(L)) ) )
  14
  15     # With high probability, medoids indices should resemble 1,1,1,...,2,2,2,...,3,3,3,...
  16     mi <- assignMedoids(series, medoids)
  17     mi_ref <- rep(1:3, each=n/3)
  18     expect_lt( mean(mi != mi_ref), 0.01 )
  19 })