REVIEW 3 major objections 5 minor 39 references
Relative Explanations for Contextual Problems with Endogenous Uncertainty: An Application to Competitive Facility Location
T0 review · 3 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read A Wasserstein-regularized mixed-integer bilinear program computes relative counterfactual explanations for binary-decision contextual stochastic programs with endogenous uncertainty, demonstrated on competitive facility location.
desk verdict Solid formulation contribution, but an impossible runtime row and a missing advertised case study undercut the empirical headline; worth a serious referee after rework. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the relative counterfactual explanation: a new context $x$ for which some feasible decision $z\in\mathcal{D}$ earns a factor $\alpha$ of the factual expected reward. The argument runs through the standard mixed-integer linear reformulation of multinomial logit choice probabilities, extended so that the transformed attractiveness values $\phi_d(x)=\exp(\bar{v}_d(x_d))$ are decision variables rather than fixed parameters; this produces bilinear constraints linking the opening decisions, the attractiveness values, and the choice probabilities. The Wasserstein regularization enters as a finite optimal transport problem: with $\pi^n$ coupling the factual and counterfactual choice distributions, the objective contains $\lambda\sum_n\sum_{c,c'}\pi^n_{cc'}\delta(c,c')^2$. A model-free version that optimizes directly over probabilities gives a lower bound, and a greedy feasibility heuristic supplies the warm start that makes the bilinear program solvable within the hour limit on the test instances.
What would settle it
Run the instance families of Table 1 with identical warm starts and a one-hour limit for $\lambda=0$ and $\lambda=0.1$: if the unregularized model reaches optimality in less time on a majority of the larger instances, the claimed computational benefit of Wasserstein regularization is not general. A sharper test is to take instances where the factual choice distribution is nearly uniform, where the Wasserstein penalty should force the largest distortions and the unregularized objective should look relatively cheap; if $\lambda=0.1$ still dominates there, the speedup is robust to distribution shape.
Extended reading notes
Core claim
The paper's central claim is that including the squared 2-Wasserstein distance between the factual and counterfactual choice distributions in the objective is not just a way to make counterfactuals realistic; it also accelerates solving the explanation problem and produces sparser solutions. In the multinomial logit facility location model, the explanation problem with a desired facility set is formulated as a single mixed-integer bilinear program over attractiveness variables $\phi_d(x)$, opening variables $z_d$, and optimal transport plans $\pi^n$, where the Wasserstein term is the transport cost between the factual and counterfactual choice distributions. A model-free relaxation over probability distributions provides a lower bound, and a constructive greedy heuristic provides a feasible warm start. On randomly generated instances, the regularized version with $\lambda=0.1$ reaches the one-hour time limit less often and yields sparser solutions than the unregularized version, while keeping the counterfactual choice probabilities close to the factual ones.
Load-bearing premise
The framework assumes the conditional distribution $P(y\mid z,x)$ is known, and every explanation is relative to that fixed model; if the model is misspecified, the computed counterfactual changes are not reliable evidence about what would happen in the world.
Editorial extensions
If this is right
- For any desired facility configuration, planners get a concrete, minimal context shift that makes that configuration acceptable, expressed in the same attractiveness or environmental variables the demand model uses.
- The same mixed-integer bilinear template transfers to other binary-decision problems with endogenous uncertainty—network design, assortment planning, school siting—as long as the outcome distribution is known and depends on covariates through a monotone transformation.
- The model-free lower bound gives a certificate of how much structural realism costs: if the bound is far below the model-based value, no context can realize the cheap distributional shift the bound finds.
- Regularized explanations are sparser and produce smoother demand transitions, so they are more likely to be perceived as plausible by stakeholders than explanations that change many small features.
- Because the desired space $\mathcal{D}$ can be any subset of configurations, the same code can answer what would make any one of a family of alternatives acceptable, not just a single expert-chosen target.
Reading between the lines
- I infer that the Wasserstein penalty will matter most when the factual choice probabilities are concentrated: with sharply peaked distributions the transport cost forces probability mass to move along short paths, whereas the unregularized $\ell_1$ objective can jump to an arbitrary distribution at the same covariate cost; the paper's experiments do not vary the concentration of the factual distri
- A natural extension the paper leaves open is to estimate $P(y\mid z,x)$ from data and propagate estimation error into the explanations; I would expect the slack in the demand constraint to be a useful diagnostic for when a found counterfactual is an artifact of a misspecified model.
- I infer from the reported instances that the speedup from regularization is most reliable for small budgets: several $r=8$ instances still reach the time limit even with $\lambda=0.1$, so the computational benefit may weaken as the budget grows.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies relative counterfactual explanations for contextual stochastic optimization problems in which the outcome distribution depends on both the context and the decisions. It formulates the explanation task as a non-convex optimization problem (2), introduces a cost that combines ℓ1 dissimilarity in a transformed feature space with a 2-Wasserstein regularization term (4), and proposes a model-free lower bound plus a greedy warm start. The method is specialized to a choice-based competitive facility location problem under a multinomial logit model; after a Haase-type reformulation, explanations are computed by solving the mixed-integer bilinear program (9). Section 5 reports experiments on a small illustrative instance and on synthetic instances of varying size, claiming that moderate Wasserstein regularization reduces computation time and yields sparser counterfactuals. The supplied abstract also claims validation on a real-world electric vehicle charging station planning case study in Montreal, but this case study is not present in the manuscript body.
Significance. If the computational claims held, this would be a useful contribution: it extends relative counterfactual explanations from continuous expert-provided target solutions to problems with binary decisions and endogenous uncertainty, and it allows a desired feasible set rather than a single target solution. The mathematical derivation from the MNL model to formulation (9) follows Haase (2009) and is internally consistent, and the authors make code and instances publicly available, which supports reproducibility. The lower-bound and warm-start scheme is a sensible practical addition. The significance of the empirical claims is, however, currently undermined by an arithmetically impossible row in Table 1 and by an abstract that advertises a real-world validation absent from the body; the modeling contribution itself appears sound.
major comments (3)
- [Section 5, Table 1] The row for N=100, D=10, r=8, λ=0.1 reports TL=5 and average time 1,094.136 s. Since five runs hit the 3,600 s time limit, any valid average over the ten instances must be at least 1,800 s, so the reported average is arithmetically impossible. This row directly supports the paper's headline claim that λ=0.1 reduces computation time, and the table as reported therefore does not support that claim. Please correct the data, recalculate the summary statistics, or rerun the affected instances before the efficiency claim can be accepted.
- [Abstract; Section 5] The abstract supplied with the paper states that the framework is validated on a real-world case study of electric vehicle charging station planning in Montreal. The full text contains no such case study: Section 5 contains only the illustrative example and the synthetic instances of Tables 1–2, and Section 6 does not mention Montreal. Either the case study and its results must be added to the manuscript, or the abstract's validation claim must be removed or corrected.
- [Section 5, Table 1] The statement that Wasserstein regularization yields sparser counterfactuals is stronger than what the table shows. For N=200, D=20, r=8, sparsity increases from 0.075 at λ=0 to 0.105 at λ=0.1 and 0.205 at λ=1, and several λ=1 rows, such as N=100, D=10, r=8, are less sparse than the corresponding λ=0 row. The qualitative claim in the text should be qualified as instance-dependent, or the analysis should explain when and why the sparsity improvement holds.
minor comments (5)
- [Section 4.1] The symbol D is used both for the set of candidate locations and for the desired feasible space; constraint (8g), written as "zd∈D", appears to be a typo for "z∈D" (or "z∈\mathcal{D}"), and the notation should be disambiguated.
- [Section 4.2] In the warm-start description, Step 2 says to calculate the factual captured demand "for each facility d," but the displayed quantity Q_factual = ∑_n q_n ∑_d P0_n(d) is a single total over all facilities; please clarify the wording.
- [Section 5] The phrase "the Wasserstein term is scaled by its bound" is not defined; please specify the exact scaling factor used in the objective before reporting the experiments.
- [Equations (9g)–(9i)] The constraint ranges "d∈D e∈E n∈N" appear to include e∈E unintentionally in (9g), which does not depend on e; this should be corrected for readability.
- [Figures 1 and 2] The captions contain typos: "solid dimonds" should be "solid diamonds" in Figure 1, and "a a solution" should be "a solution" in Figure 2.
Circularity Check
No significant circularity: the counterfactual formulation and bilinear reformulation are derived from explicit definitions and an external reformulation, with self-citations confined to background and notation.
full rationale
The derivation chain is self-contained. Definition 1.1 fixes the notion of a relative explanation; Problem (2) is the literal optimization formulation of that definition, and (3) merely restricts the distribution P to a family H. In the CFLP application, (8) is obtained by substituting the MNL choice probabilities (7) and the cost L = l1 + lambda*W2 into (2)-(3), and (9) is an algebraic reformulation of (8) using the external Haase (2009) linearization plus the added constraint (9i). None of these steps fits a parameter to the quantity being predicted. The Wasserstein term is an explicit design choice governed by the hyperparameter lambda, and the efficiency/sparsity claims are empirical comparisons of lambda = 0 versus lambda > 0 on the same instances; observing a lower Wasserstein distance under regularization is the intended effect of the objective, and the paper labels it 'as expected' rather than presenting it as an independent prediction. The self-citations (Bogetoft et al. 2024; Carrizosa et al. 2024; Legault and Frejinger 2025; Sadana et al. 2025) appear only in related-work discussion, as notation, or as a survey pointer, and no load-bearing uniqueness theorem or ansatz is imported from the authors' own prior work. The known-distribution assumption is explicitly scoped in Section 1 and is a modeling limitation, not a circular input. Separate empirical-validity concerns exist: Table 1's N=100, D=10, r=8, lambda=0.1 row reports TL=5 with a 1094-second average, which is arithmetically impossible under a 3600-second time limit, and the abstract's advertised Montreal case study does not appear in the body. These concerns affect the factual support for the runtime and validation claims, but they are not circularity in the derivation chain.
Assumptions & free parameters
free parameters (5)
- λ (Wasserstein regularization weight) =
0, 0.1, 1 (experiments)
- α (allowed performance loss factor) =
1 in all experiments
- ϵ (positive lower bound in lower-bound formulation) =
not reported
- Wasserstein scaling factor =
not reported
- greedy warm-start attractiveness constant c =
increased until demand constraint met
assumptions (5)
- domain assumption P(y|z,x) is known, with MNL parameters β given
- domain assumption Contextual covariates are continuous and unbounded above (X ⊆ R_+^{dx})
- domain assumption Utility \bar{v}_d is linear in features for the post-processing recovery of x
- domain assumption MNL proportional substitution pattern is a valid model of customer choice
- domain assumption Squared 2-Wasserstein with Euclidean ground distance on facility locations is the appropriate dissimilarity for distributions
Cite this review
Pith. "Pith review of Relative Explanations for Contextual Problems with Endogenous Uncertainty: An Application to Competitive Facility Location." pith.science (2026). https://pith.science/paper/X2QOQCJ6
@misc{pith2026250619155,
author = {Pith},
title = {Pith review of: Relative Explanations for Contextual Problems with Endogenous Uncertainty: An Application to Competitive Facility Location},
year = {2026},
howpublished = {\url{https://pith.science/paper/X2QOQCJ6}},
note = {Machine review of arXiv:2506.19155}
}
read the original abstract
In this paper, we consider a contextual stochastic optimization problem in which unknown parameters follow distributions that depend on contextual covariates and decisions. The problem is motivated by transportation infrastructure decisions such as facility location or network design. In such high-stakes settings, decisions must often be communicated, justified, and reconsidered under alternative stakeholder requirements. To this end, we propose a framework for computing relative counterfactual explanations. These explanations identify the smallest changes in the covariates required for a solution to satisfy prescribed constraints while limiting the performance loss to a controlled level. Whereas relative explanations have been introduced in prior literature, to the best of our knowledge, this is the first work focusing on problems with binary decision variables and endogenous uncertainty. We propose a methodology that uses the Wasserstein distance as a regularization term in the objective. Beyond improving tractability, this regularization yields explanations with desirable structural properties: it produces sparser counterfactuals, induces smoother transitions in the underlying choice distributions, and keeps the counterfactual behavior close to realistic demand patterns. We illustrate the method using a choice-based competitive facility location problem and present numerical experiments that demonstrate its ability to efficiently compute sparse, plausible, and interpretable explanations. We further validate the framework on a real-world case study of electric vehicle charging station planning in Montreal, where the explanations reveal the minimal capacity investments and environmental conditions required to justify including a candidate location in the charging network.
Figures
Reference graph
Works this paper leans on
-
[1]
Aigner, K.-M., Goerigk, M., Hartisch, M., Liers, F., and Miehlich, A. (2024). A framework for data-driven explainability in mathematical optimization. In Proceedings of the AAAI Conference on Artificial Intelligence , volume 38, pages 20912--20920
work page 2024
-
[2]
Aros-Vera, F., Marianov, V., and Mitchell, J. E. (2013). p-hub approach for the optimal park-and-ride facility location problem. European Journal of Operational Research , 226(2):277--285
work page 2013
-
[3]
Azizi, M., Vayanos, P., Wilder, B., Rice, E., and Tambe, M. (2018). Designing fair, efficient, and interpretable policies for prioritizing homeless youth for housing resources. In 15th International Conference Integration of Constraint Programming, Artificial Intelligence, and Operations Research, CPAIOR , pages 35--51. Springer
work page 2018
-
[4]
Benati, S. and Hansen, P. (2002). The maximum capture problem with random utilities: Problem formulation and algorithms. European Journal of Operational Research , 143(3):518--530
work page 2002
-
[5]
Bertsimas, D., Delarue, A., Eger, W., Hanlon, J., and Martin, S. (2020). Bus routing optimization helps B oston public schools design better policies. INFORMS Journal on Applied Analytics , 50(1):37--49
work page 2020
-
[6]
Bertsimas, D. and Mundru, N. (2023). Optimization-based scenario reduction for data-driven two-stage stochastic optimization. Operations Research , 71(4):1343--1361
work page 2023
-
[7]
Blanchet, J., Kang, Y., and Murthy, K. (2019). Robust W asserstein profile inference and applications to machine learning. Journal of Applied Probability , 56(3):830--857
2019
-
[8]
Bogetoft, P., Ram \' rez-Ayerbe, J., and Morales, D. R. (2024). Counterfactual analysis and target setting in benchmarking. European Journal of Operational Research , 315(3):1083--1095
work page 2024
Show all 39 references
-
[9]
Carrizosa, E., Ramírez-Ayerbe, J., and Romero Morales , D. (2024). Mathematical optimization modelling for group counterfactual explanations. European Journal of Operational Research , 319(2):399--412
2024
-
[10]
Eubanks, V. (2018). Automating inequality: How high-tech tools profile, police, and punish the poor . St. Martin's Press
2018
-
[11]
White Paper on Artificial Intelligence: A European approach to excellence and trust
European Commission (2020). White Paper on Artificial Intelligence: A European approach to excellence and trust . https://ec.europa.eu/info/publications/white-paper-artificial-intelligence-european-approach-excellence-and-trust_en. Accessed on June 22, 2025
2020
-
[12]
Forel, A., Parmentier, A., and Vidal, T. (2023). Explainable data-driven optimization: from context to decision and back again. In International Conference on Machine Learning , pages 10170--10187. PMLR
2023
-
[13]
Gao, R., Chen, X., and Kleywegt, A. J. (2024). Wasserstein distributionally robust optimization and variation regularization. Operations Research , 72(3):1177--1191
2024
-
[14]
and Hartisch, M
Goerigk, M. and Hartisch, M. (2023). A framework for inherently interpretable optimization models. European Journal of Operational Research , 310(3):1312--1324
2023
-
[15]
Guidotti, R. (2024). Counterfactual explanations and how to find them: literature review and benchmarking. Data Mining and Knowledge Discovery , 38:2770--2824
2024
-
[16]
Gurobi Optimizer Reference Manual
Gurobi Optimization, LLC (2024). Gurobi Optimizer Reference Manual
2024
-
[17]
Haase, K. (2009). Discrete location planning. Technical Report WP-09-07, Institute for Transport and Logistics Studies, University of Sydney
2009
-
[18]
and Müller, S
Haase, K. and Müller, S. (2014). A comparison of linear reformulations for multinomial logit choice probabilities in facility location models. European Journal of Operational Research , 232(3):689--691
2014
-
[19]
Voluntary code of conduct on the responsible development and management of advanced generative ai systems
Innovation, Science and Economic Development Canada (2023). Voluntary code of conduct on the responsible development and management of advanced generative ai systems. Government of Canada
2023
-
[20]
Karimi, A.-H., Barthe, G., Sch \"o lkopf, B., and Valera, I. (2022). A survey of algorithmic recourse: C ontrastive explanations and consequential recommendations. ACM Computing Surveys , 55(5):1--29
2022
-
[21]
and Beck, J
Korikov, A. and Beck, J. C. (2021). Counterfactual explanations via inverse constraint programming. In 27th International Conference on Principles and Practice of Constraint Programming (CP 2021) , pages 35--1. Schloss Dagstuhl--Leibniz-Zentrum f \"u r Informatik
2021
-
[22]
and Beck, J
Korikov, A. and Beck, J. C. (2023). Objective-based counterfactual explanations for linear discrete optimization. In International Conference on Integration of Constraint Programming, Artificial Intelligence, and Operations Research , pages 18--34. Springer
2023
-
[23]
Korikov, A., Shleyfman, A., and Beck, C. (2021). Counterfactual explanations for optimization-based decisions in the context of the GDPR . In ICAPS 2021 workshop on explainable AI planning
2021
-
[24]
\.I ., and Hertog, D
Kurtz, J., Birbil, S . \.I ., and Hertog, D. d. (2024). Counterfactual explanations for linear optimization. arXiv preprint arXiv:2405.15431
2024 arXiv
-
[25]
and Schmidt, M
Lefebvre, H. and Schmidt, M. (2025). Computing counterfactual explanations for linear optimization: A new class of bilevel models and a tailored penalty alternating direction method. https://optimization-online.org/?p=28803. Working paper, Optimization Online
2025
-
[26]
and Frejinger, E
Legault, R. and Frejinger, E. (2024). A model-free approach for solving choice-based competitive facility location problems using simulation and submodularity. INFORMS Journal on Computing
2024
-
[27]
and Provost, F
Martens, D. and Provost, F. (2014). Explaining data-driven document classifications. MIS Quarterly , 38(1):73--99
2014
-
[28]
Miller, T. (2019). Explanation in artificial intelligence: Insights from the social sciences. Artificial Intelligence , 267:1--38
2019
-
[29]
Mi s i \'c , V. V. and Perakis, G. (2020). Data analytics in operations management: A review. Manufacturing & Service Operations Management , 22(1):158--169
2020
-
[30]
Panaretos, V. M. and Zemel, Y. (2019). Statistical aspects of wasserstein distances. Annual Review of Statistics and Its Application , 6(Volume 6, 2019):405--431
2019
-
[31]
Peyr \'e , G., Cuturi, M., et al. (2019). Computational optimal transport: With applications to data science. Foundations and Trends in Machine Learning , 11(5-6):355--607
2019
-
[32]
Sadana, U., Chenreddy, A., Delage, E., Forel, A., Frejinger, E., and Vidal, T. (2025). A survey of contextual optimization methods for decision-making under uncertainty. European Journal of Operational Research , 320(2):271--289
2025
-
[33]
Solomon, J. (2018). Optimal transport on discrete domains. AMS Short Course on Discrete Differential Geometry
2018
-
[34]
Villani, C. et al. (2009). Optimal transport: old and new , volume 338. Springer
2009
-
[35]
Vivier-Ardisson, G., Forel, A., Parmentier, A., and Vidal, T. (2024). Cf-opt: Counterfactual explanations for structured prediction. In Forty-first International Conference on Machine Learning
2024
-
[36]
Wachter, S., Mittelstadt, B., and Russell, C. (2017). Counterfactual explanations without opening the black box: Automated decisions and the GDPR . Harvard Journal of Law & Technology , 31:841--887
2017
-
[37]
Blueprint for an ai bill of rights: Making automated systems work for the american people
White House Office of Science and Technology Policy (2022). Blueprint for an ai bill of rights: Making automated systems work for the american people. https://bidenwhitehouse.archives.gov/ostp/ai-bill-of-rights/ Retrieved from the White House archives. Accessed on June 22, 2025
2022
-
[38]
Xie, W. (2021). On distributionally robust chance constrained programs with wasserstein distance. Mathematical Programming , 186(1):115--155
2021
-
[39]
You, L., Cao, L., Nilsson, M., Zhao, B., and Lei, L. (2024). Distributional counterfactual explanation with optimal transport. arXiv preprint arXiv:2401.13112
2024 arXiv
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.