Skip to contents
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)

dep_agg = dep %>% 
  add_column(cl= cutree(res,h=7)) %>% 
  group_by(cl) %>% 
  summarise(p_gauche=mean(p_gauche))

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)