ttimbers/between_ss_total_SS_ratio_kmeans.R

## between_ss_total_SS_ratio_kmeans.R
set.seed(123)
for(n in c(50, 200, 500, 1000, 2000, 10000) ) {
  x1 <- data.frame(x1 = rnorm(n, 3), x2 = rnorm(n, 3))
  x2 <- data.frame(x1 = rnorm(n, -2), x2 = rnorm(n, -2))
  x <- rbind(x1, x2)
  a <- kmeans(x, centers= 2)
  print(paste0('Sample size: ', n, ' - between_SS / total_SS: ', round(100 * a$betweenss/a$totss, 2)))
}

## between_ss_total_SS_ratio_kmeans_output.md

      
    Raw
  

              between_ss_total_SS_ratio_kmeans_output.md
            
          
    [1] "Sample size: 50 - between_SS / total_SS: 88.67"
[1] "Sample size: 200 - between_SS / total_SS: 86.34"
[1] "Sample size: 500 - between_SS / total_SS: 86.68"
[1] "Sample size: 1000 - between_SS / total_SS: 86.23"
[1] "Sample size: 2000 - between_SS / total_SS: 86.1"
[1] "Sample size: 10000 - between_SS / total_SS: 86.11"
	set.seed(123)
	for(n in c(50, 200, 500, 1000, 2000, 10000) ) {
	x1 <- data.frame(x1 = rnorm(n, 3), x2 = rnorm(n, 3))
	x2 <- data.frame(x1 = rnorm(n, -2), x2 = rnorm(n, -2))
	x <- rbind(x1, x2)
	a <- kmeans(x, centers= 2)
	print(paste0('Sample size: ', n, ' - between_SS / total_SS: ', round(100 * a$betweenss/a$totss, 2)))
	}