gtclust
gtclust.Rmd
data("modesshare.idf")
res=gtclust_poly(modesshare.idf,"ward",scaling="zscore")
#> Warning: Some features were not numeric and have been removed from the
#> clustering.
plot(res)
idf_agg=geocutree(res,k=80) |> rowwise(cl) |> mutate(total = sum(c_across(nodep:tcom)))
idf = modesshare.idf |> rowwise(CODE_IRIS) |> mutate(total = sum(c_across(nodep:tcom)))
#> Warning in st_centroid.sf(idf_agg): st_centroid assumes attributes are constant
#> over geometries of x
#res=geohclust_poly(modesshare,method = "ward")
#fr_agg2k=geocutree(res,2000)
#fr_agg10k=geocutree(res,1000)