#!/usr/bin/env Rscript


income_elec_state = data.frame(a = c(1, 3:6, 8, 10, 12, 13), b = c(1:5, seq(6, 12, by = 2)))

plot(income_elec_state)

if (!require(ggdendro)) install.packages("ggdendro", repos = "https://cran.r-project.org/")
library(ggdendro)

distance <- dist(income_elec_state, method = "euclidean")

# single - closest pair of points
# complete - farthest pair of points
# average - average distance between all pairs
h_clust <- hclust(distance, method = "singlet")
h_clust

plot(ggdendrogram(h_clust))

cutree(h_clust, k = 2)