# K mean clustering
install.packages(c("ggplot2", "factoextra", "cluster"))

library(ggplot2)
library(factoextra)
library(cluster)

data("iris")
irisdata <- scale(iris[, -5])

set.seed(123)

fviz_nbclust(irisdata, kmeans, method = "wss")

model <- kmeans(irisdata, centers = 3, nstart = 25)

iris$Cluster <- as.factor(model$cluster)

print(model$centers)

table(model$cluster)

fviz_cluster(model, data = irisdata)

sil <- silhouette(model$cluster, dist(irisdata))

fviz_silhouette(sil)