library(cluster) library(tidyverse)

energy <- read.csv(“energy_drinks.csv”) View(energy)

en <- select(energy, D1:D5) e1 <- dist(en)

h1 <- hclust(e1, method = “average”) plot(h1, hang = -1) heatmap(as.matrix(e1), Rowv = as.dendrogram(h1), Colv = ‘Rowv’)

clustere1 <- cutree(h1, k = 3)

sil1 <- silhouette(clustere1, e1) summary(sil1)

en_clus <- cbind(energy, clustere1) en_clus <- mutate(en_clus, cluster = case_when(clustere1 == 1 ~ ‘C1’, clustere1 == 2 ~ ‘C2’, clustere1 == 3 ~ ‘C3’, ))

write_csv(en_clus, ‘en_sur_cluster.csv’)