Testing several clustering methods in spatial_clustering_cv

make_clusters.R

# remotes::install_github('tidymodels/spatialsample')
library(spatialsample)
library(cluster)

make_clusters <- function(func) { 
  set.seed(123)
  boston_canopy |> 
    spatial_clustering_cv(cluster_function = func) |> 
    autoplot() 
}

function_list <- list(
  "kmeans",
  "hclust",
  \(dists, v, ...) fanny(dists, v, maxit = 1e4)$clustering,
  \(dists, v, ...) pam(dists, v)$clustering,
  \(dists, v, ...) cutree(agnes(dists), v),
  \(dists, v, ...) cutree(diana(dists), v)
)

lapply(function_list, make_clusters)