Q7.1

<aside> 💡 637,86

</aside>

# required libraries
library(cluster)

# centers vector
X = c(-4, 0, 4)
Y = c(10, 0, 10)
centers = data.frame(X,Y)

# kMeans model
model <- kmeans(data, centers)

# cohesion
model$tot.withinss

Q7.2

<aside> 💡 9496,20

</aside>

# required libraries
library(cluster)

model$betweenss

Q7.3

<aside> 💡 0,78

</aside>

model_silhouette <- silhouette(model$cluster, dist(data))
model1_mean_silhouette = mean(model_silhouette[, 3])

Q7.4

<aside> 💡 (b)

</aside>

# new centers vector 
X = c(-2, 2, 0)
Y = c(0, 0, 10)
centers = data.frame(X,Y)

# kMeans model 2
model <- kmeans(data, centers)

# model 2 silhouette
model_silhouette <- silhouette(model$cluster, dist(data))
model2_mean_silhouette = mean(model_silhouette[, 3])

# return better model
which.max(c(model1_mean_silhouette, model2_mean_silhouette))