data(pennLC)
max.prop <- 0.15
pop <- tapply(pennLC$data$population, pennLC$data$county, sum)
centroids <- latlong2grid(pennLC$geo[, 2:3])
map <- pennLC$spatial.polygon
output <- create.geo.objects(max.prop, pop, centroids, map)
## number of single zones
nrow(output$cluster.coords)
length(output$areaz)
Run the code above in your browser using DataCamp Workspace