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