This function uses redist_mergesplit()
or redist_flip()
to optimize a
redistrict plan according to a user-provided criteria. It does so by running
the Markov chain for "short bursts" of usually 10 iterations, and then
starting the chain anew from the best plan in the burst, according to the
criteria. This implements the ideas in the below-referenced paper, "Voting
Rights, Markov Chains, and Optimization by Short Bursts."
redist_shortburst(
map,
score_fn = NULL,
stop_at = NULL,
burst_size = ifelse(backend == "mergesplit", 10L, 50L),
max_bursts = 500L,
maximize = TRUE,
init_plan = NULL,
counties = NULL,
constraints = redist_constr(map),
compactness = 1,
adapt_k_thresh = 0.95,
reversible = TRUE,
fixed_k = NULL,
return_all = TRUE,
thin = 1L,
backend = "mergesplit",
flip_lambda = 0,
flip_eprob = 0.05,
verbose = TRUE
)
a redist_plans
object containing the final best plan
(or the best plans after each burst, if return_all=TRUE
.
A redist_map object.
A function which takes a matrix of plans and returns a score
(or, generally, a row vector) for each plan. Can also be a purrr-style
anonymous function. See ?scorers
for some function factories
for common scoring rules.
A threshold to stop optimization at. When score_fn
returns a
row vector per plan, maximize
can be an equal-length vector specifying a
threshold for each dimension, which must all be met for the algorithm to
stop.
The size of each burst. 10 is recommended for the
mergesplit
backend and 50 for the flip
backend. Can also provide
burst schedule function which takes the current iteration (an integer)
and returns the desired burst size. This can be a random function.
The maximum number of bursts to run before returning.
If TRUE
, try to maximize the score; otherwise, try to
minimize it. When score_fn
returns a row vector per plan, maximize
can
be an equal-length vector specifying whether each dimension should be
maximized or minimized.
The initial state of the map. If not provided, will default to
the reference map of the map
object, or if none exists, will sample
a random initial state using redist_smc()
. You can also request
a random initial state by setting init_plan="sample"
.
A vector containing county (or other administrative or
geographic unit) labels for each unit, which may be integers ranging from 1
to the number of counties, or a factor or character vector. If provided, the
algorithm will only generate maps which split up to ndists-1
counties.
If no county-split constraint is desired, this parameter should be left blank.
A redist_constr
with Gibbs constraints.
Controls the compactness of the generated districts, with
higher values preferring more compact districts. Must be non-negative. See
redist_mergesplit
for more information.
The threshold value used in the heuristic to select a
value k_i
for each splitting iteration.
If FALSE
and backend="mergesplit"
, the Markov chain
used will not be reversible. This may speed up optimization.
If not NULL
, will be used to set the k
parameter for the
mergesplit
backend. If e.g. k=1
then the best edge in each spanning
tree will be used. Lower values may speed up optimization at the
cost of the Markov chain no longer targeting a known distribution.
Recommended only in conjunction with reversible=FALSE
.
Whether to return all the burst results or just the best one (generally, the Pareto frontier). Recommended for monitoring purposes.
Save every thin
-th sample. Defaults to no thinning (1). Ignored
if return_all=TRUE
.
the MCMC algorithm to use within each burst, either "mergesplit" or "flip".
The parameter determining the number of swaps to attempt each iteration of flip mcmc. The number of swaps each iteration is equal to Pois(lambda) + 1. The default is 0.
The probability of keeping an edge connected in flip mcmc. The default is 0.05.
Whether to print out intermediate information while sampling. Recommended for monitoring purposes.
Cannon, S., Goldbloom-Helzner, A., Gupta, V., Matthews, J. N., & Suwal, B. (2020). Voting Rights, Markov Chains, and Optimization by Short Bursts. arXiv preprint arXiv:2011.02288.
# \donttest{
data(iowa)
iowa_map <- redist_map(iowa, existing_plan = cd_2010, pop_tol = 0.01)
redist_shortburst(iowa_map, scorer_frac_kept(iowa_map), max_bursts = 50)
redist_shortburst(iowa_map, ~ 1 - scorer_frac_kept(iowa_map)(.), max_bursts = 50)
# }
Run the code above in your browser using DataLab