choicedesign.design

Classes for constructing efficient experimental designs.

The main user-facing class is EffDesign, which combines initial design generation, constrained random-swap optimisation, and optional blocking.

Typical workflow:

design = EffDesign(X=[alt1_A, alt2_A], ncs=18)
init   = design.gen_initdesign(cond=['alt1_A > alt2_A'], seed=42)
result = design.optimise(init, V={1: V1, 2: V2}, time_lim=1)
# result = (optimal_design, init_derr, final_derr, n_iter, utility_balance)

Classes

EffDesign(X, ncs)

Efficient design for a discrete choice experiment.

FullFactDesign(X)

Full-factorial design covering all combinations of attribute levels.
class choicedesign.design.EffDesign(X: dict, ncs: int)

Efficient design for a discrete choice experiment.

Combines initial design generation, random-swap optimisation, and optional blocking in a single object.

Parameters:

  • X (list[Attribute]) – Ordered list of Attribute objects that define the design columns and their discrete levels.

  • ncs (int) – Number of choice situations (rows) in the design matrix.

Examples

>>> from choicedesign.design import EffDesign
>>> from choicedesign.expressions import Attribute, Parameter
>>> alt1_A = Attribute('alt1_A', [1, 2, 3])
>>> alt2_A = Attribute('alt2_A', [1, 2, 3])
>>> beta_A  = Parameter('beta_A', -0.1)
>>> V = {1: beta_A * alt1_A, 2: beta_A * alt2_A}
>>> design = EffDesign(X=[alt1_A, alt2_A], ncs=18)
>>> init   = design.gen_initdesign(seed=42)
>>> result = design.optimise(init, V=V, iter_lim=500)
evaluate(design: DataFrame, V: dict, model: str = 'mnl', criterion: str = 'd', cost_param=None, wtp_params=None, bayes_draws: int = None, seed: int = None)

Evaluate design

Evaluates a design stored in a Pandas data frame

Parameters:

  • design (pd.DataFrame) – Design to evaluate

  • V (dict) – A dictionary with the utility function.

  • model (str) – The base model for the efficient design, by default ‘mnl’

  • criterion (str) – Optimality criterion: 'd' (D-error, default), 'a' (A-error), or 'c' (C-error / WTP variance). When 'c', cost_param and wtp_params are required.

  • cost_param (Parameter, optional) – The cost (denominator) parameter. Required when criterion='c'.

  • wtp_params (list[Parameter], optional) – Parameters whose WTP variances are evaluated. Required when criterion='c'.

  • bayes_draws (int, optional) – Number of Monte Carlo draws for Db-error evaluation. Only valid with criterion='d'.

  • seed (int, optional) – Random seed for Bayesian draws, by default None

Returns:

  • perf (float) – The criterion value of the design

  • ubalance_ratio (float) – Utility balance ratio

gen_blocks(design: DataFrame, n_blocks: int, n_iter: int = 1000)

Assign choice situations to blocks.

Minimises the correlation between the block assignment and all attribute columns by evaluating n_iter random permutations and keeping the best one.

Parameters:

  • design (pandas.DataFrame) – Optimised design from optimise() (must include a CS column).

  • n_blocks (int) – Number of blocks.

  • n_iter (int, optional) – Number of random permutations evaluated by the search, by default 1000.

Returns:

  • design (pandas.DataFrame) – Design with an additional Block column.

  • corr_list (list[float]) – History of best total absolute correlation found at each improvement.

gen_initdesign(cond: list = None, seed: bool = None)

Generate initial design matrix

It generates the initial design matrix. The user can define a set of conditions that must be satisfied.

Parameters:

  • cond (list[str], optional) – List of conditions that the final design must hold. Each element is a string that contains a single condition. Conditions can be of the form of binary relations (e.g., X > Y where X and Y are attributes of a specific alternative) or conditional relations (e.g., if X > a then Y < b where a and b are values). Users can specify multiple conditions when the operator if is defined, separated by the operator &, by default None

  • seed (bool, None) – Random seed, by default None

Returns:

init_design – A Pandas DataFrame with the initial design matrix.

Return type:

pandas.DataFrame

optimise(init_design: DataFrame, V: dict, model: str = 'mnl', algorithm: str = 'swap', criterion: str = 'd', cost_param=None, wtp_params=None, bayes_draws: int = None, iter_lim: int = None, noimprov_lim: int = None, time_lim: int = None, seed: int = None, verbose: bool = False)

Optimise the design using a random-search algorithm.

Starts from an initial design and iteratively improves it according to the selected criterion and stopping rules. At least one stopping criterion (iter_lim, noimprov_lim, or time_lim) must be supplied.

Parameters:

  • init_design (pandas.DataFrame) – The initial design matrix, typically from gen_initdesign().

  • V (dict) – A dictionary with the utility functions, keyed by alternative index. e.g. {1: V1, 2: V2}

  • model (str) – The base model for the efficient design, by default ‘mnl’

  • algorithm (str) – Optimisation algorithm to use. Options: 'swap' (random swapping, default), 'rsc' (random Relabelling, Swapping, Cycling), or 'federov' (Modified Federov — tries all full-factorial candidates per row per iteration; slower per iteration but more systematic).

  • criterion (str) – Optimality criterion: 'd' (D-error, default), 'a' (A-error), or 'c' (C-error / WTP variance). When 'c', cost_param and wtp_params are required.

  • cost_param (Parameter, optional) – The cost (denominator) parameter used to compute WTP ratios. Required when criterion='c'.

  • wtp_params (list[Parameter], optional) – Parameters whose WTP ratios are minimised. Required when criterion='c'.

  • bayes_draws (int, optional) – Number of Monte Carlo draws for Db-efficient (Bayesian) design. Only valid with criterion='d'. When set, the optimizer minimises the expected D-error averaged over draws from the prior distributions of parameters that have prior_std defined.

  • iter_lim (int, optional) – Number of iterations before the algorithm stops, by default None

  • noimprov_lim (int, optional) – Number of iterations without improvement before the algorithm stops, by default None

  • time_lim (int, optional) – Time (in minutes) before the algorithm stops, by default None

  • seed (int, optional) – Random seed, by default None

  • verbose (bool, optional) – Whether status messages and progress are shown, by default False

Returns:

  • optimal_design (pandas.DataFrame) – The final (optimal) design

  • init_perf (float) – Criterion value of the initial design

  • final_perf (float) – Criterion value of the final design

  • final_iter (int) – Total number of iterations

  • ubalance_ratio (float) – Utility balance ratio

class choicedesign.design.FullFactDesign(X: list)

Full-factorial design covering all combinations of attribute levels.

Parameters:

X (list[Attribute]) – List of Attribute objects.

Notes

This class is provided for completeness. Its dependency pyDOE2 is incompatible with Python 3.12+ and is not installed by default.

gen_design()

Generate full-factorial design matrix

Returns:

design – A Pandas DataFrame with all combinations of attribute levels.

Return type:

pandas.DataFrame