# R program for additional multivariate exercise 5 (VHM 802)
# only includes code for K-means clustering

euroemp <- read.csv("r:/euroemp.csv", header=TRUE, row.names=1)
seuroemp <- euroemp
seuroemp[,2:10] <- scale(seuroemp[,2:10])
set.seed(210315)

seuroemp.k2 <- kmeans(seuroemp[,2:10],centers=2,nstart=1000)
seuroemp.k2

seuroemp.k3 <- kmeans(seuroemp[,2:10],centers=3,nstart=1000)
seuroemp.k3

seuroemp.k4 <- kmeans(seuroemp[,2:10],centers=4,nstart=1000)
seuroemp.k4

seuroemp.k5 <- kmeans(seuroemp[,2:10],centers=5,nstart=1000)
seuroemp.k5

seuroemp.k6 <- kmeans(seuroemp[,2:10],centers=6,nstart=1000)
seuroemp.k6