This R package enables researchers to sample redistricting plans from a pre-specified target distribution using Sequential Monte Carlo and Markov Chain Monte Carlo algorithms. The package supports various constraints in the redistricting process, such as geographic compactness and population parity requirements. Tools for analysis, including computation of various summary statistics and plotting functionality, are also included.
redist is available on CRAN and can be installed using:
You can also install the most recent development version of
redist (which is usually quite stable) using the `remotes`` package.
A basic analysis has two steps. First, you define a redistricting plan using
redist_map. Then you simulate plans using one of the algorithm functions:
library(redist) library(dplyr) data(iowa) # set a 0.1% population constraint iowa_map = redist_map(iowa, existing_plan=cd_2010, pop_tol=0.001, total_pop = pop) # simulate 500 plans using the SMC algorithm iowa_plans = redist_smc(iowa_map, nsims=500) #> SEQUENTIAL MONTE CARLO #> Sampling 500 99-unit maps with 4 districts and population between 760827 and 762350.
After generating plans, you can use
redist’s plotting functions to study the geographic and partisan characteristics of the simulated ensemble.
iowa_plans = iowa_plans %>% mutate(Compactness = distr_compactness(iowa_map), `Population deviation` = plan_parity(iowa_map), `Democratic vote` = group_frac(iowa_map, dem_08, tot_08)) hist(iowa_plans, `Population deviation`) + hist(iowa_plans, Compactness) + plot_layout(guides="collect") + plot_annotation(title="Simulated plan characteristics")
redist.plot.scatter(iowa_plans, `Population deviation`, Compactness) + labs(title="Population deviation and compactness by plan")
plot(iowa_plans, `Democratic vote`, size=0.5, color_thresh=0.5) + scale_color_manual(values=c("black", "tomato2", "dodgerblue")) + labs(title="Democratic vote share by district")