# NOT RUN {
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