library(dplyr)
data(iowa)
iowa <- redist_map(iowa, existing_plan = cd_2010, pop_tol = 0.05, total_pop = pop)
plans <- redist_smc(iowa, nsims = 100, silent = TRUE)
plans %>%
mutate(comp = distr_compactness(iowa)) %>%
group_by(draw) %>%
summarize(pop_dev = max(abs(total_pop/mean(total_pop) - 1)),
comp = comp[1]) %>%
redist.plot.scatter(pop_dev, comp)
Run the code above in your browser using DataLab