data("pres2022_t1_communes")
deplist = c(56)
dep = pres2022_t1_communes |>
filter(INSEE_DEP %in% deplist) |>
select(p_gauche,p_droite,p_extreme_droite) |>
filter(!is.na(p_gauche))
res=geohclust_poly(dep)
#> Warning: Some regions were isolated. The hierarchy was automatically completed
#> to reach one cluster.
plot(res)


fr = pres2022_t1_communes |>
select(p_gauche,p_droite,p_extreme_droite,NOM) |>
filter(!is.na(p_gauche))
res=geohclust_poly(fr)
#> Warning: Some features were not numeric and have been removed from the
#> clustering.
#> Warning: Some regions were isolated. The hierarchy was automatically completed
#> to reach one cluster.
fr_agg2k=geocutree(res,2000)
fr_agg10k=geocutree(res,1000)


