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,...
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,...