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

pdogs <- read.csv("r:/prehist_dog.csv", header=TRUE, row.names=1)
spdogs <- scale(pdogs)
set.seed(210315)

spdogs.k2 <- kmeans(spdogs,centers=2,nstart=100)
spdogs.k2
spdogs.k2$cluster

spdogs.k3 <- kmeans(spdogs,centers=3,nstart=100)
spdogs.k3
spdogs.k3$cluster