It is creating `NA` values due to using `st_within`, migh be better just to check if centroids are inside the cluster regions.
It is creating
NAvalues due to usingst_within, migh be better just to check if centroids are inside the cluster regions.