redist.crsg generates redistricting plans using a random seed a grow
algorithm. This is the compact districting algorithm described in Chen and
Rodden (2013).
redist.crsg(
adj,
total_pop,
shp,
ndists,
pop_tol,
verbose = TRUE,
maxiter = 5000
)List of length N, where N is the number of precincts. Each list element is an integer vector indicating which precincts that precinct is adjacent to. It is assumed that precinct numbers start at 0.
numeric vector of length N, where N is the number of precincts. Each element lists the population total of the corresponding precinct, and is used to enforce pop_tol constraints.
An sf dataframe to compute area and centroids with.
integer, the number of districts we want to partition the precincts into.
numeric, indicating how close district population targets have to be to the target population before algorithm converges. pop_tol=0.05 for example means that all districts must be between 0.95 and 1.05 times the size of target.pop in population size.
boolean, indicating whether the time to run the algorithm is printed.
integer, indicating maximum number of iterations to attempt before convergence to population constraint fails. If it fails once, it will use a different set of start values and try again. If it fails again, redist.rsg() returns an object of all NAs, indicating that use of more iterations may be advised. Default is 5000.
list, containing three objects containing the completed redistricting plan.
plan: A vector of length N, indicating the
district membership of each precinct.
district_list A list of length Ndistrict. Each list contains a
vector of the precincts in the respective district.
district_pop A vector of length Ndistrict, containing the
population totals of the respective districts.
Jowei Chen and Jonathan Rodden (2013) ``Unintentional Gerrymandering: Political Geography and Electoral Bias in Legislatures.'' Quarterly Journal of Political Science. 8(3): 239-269.
data("fl25")
adj <- redist.adjacency(fl25)
redist.crsg(adj = adj, total_pop = fl25$pop, shp = fl25, ndists = 2, pop_tol = .1)
#>
#> ====================
#> redist.crsg(): Automated Redistricting Starts
#>
#>
#> 2 districts built using 25 precincts in 0 seconds...
#>
#> $plan
#> [1] 2 2 2 2 1 2 1 1 2 1 1 1 1 2 1 2 2 2 1 1 2 1 2 1 1
#>
#> $district_list
#> $district_list[[1]]
#> [1] 9 12 14 7 23 10 11 24 21 4 6 19 18
#>
#> $district_list[[2]]
#> [1] 22 5 13 8 16 3 1 15 2 17 20 0
#>
#>
#> $district_pop
#> [1] 82525 92518
#>