R: Cluster spatial polygons to bigger areas with attribute constraint of minimal population

Question

Given the dataset with the spatial polygons below, I want to cluster these polygons to bigger areas, based on a similarity matrix. For the computation, I would like to use R