Higher-Order Stochastic Dominance Constraints in Optimization
Rajmadan Lakshmanan, Alois Pichler, Miloš Kopa
TL;DR
This paper tackles optimization under higher-order stochastic dominance constraints, which are inherently infinite in number. It develops two equivalent formulations—one using expectation-based dominance $X \preccurlyeq^{(p)} Y$ and another using risk measures $\mathcal{R}_\beta$—and proves that the uncountable constraint set can be reduced to a finite set of test points. A Newton-based optimization framework is derived, employing a single aggregate constraint $\overline{g}_p(x) \le 0$ and gradient calculations via the implicit function theorem, enabling efficient handling of non-linear dominance constraints. Numerical experiments on portfolio optimization with real data (MSCI benchmark and standardized datasets) demonstrate that higher-order dominance yields higher returns and more concentrated allocations to top-performing assets, while significantly reducing constraint complexity relative to traditional SSD/TSD methods. The work provides practical tools for robust, higher-order decision-making under uncertainty and lays groundwork for broader applications and future extensions in multi-objective and multi-dominance settings.
Abstract
This contribution examines optimization problems that involve stochastic dominance constraints. These problems have uncountably many constraints. We develop methods to solve the optimization problem by reducing the constraints to a finite set of test points needed to verify stochastic dominance. This improves both theoretical understanding and computational efficiency. Our approach introduces two formulations of stochastic dominance$\unicode{x2013}$one employs expectation operators and another based on risk measures$\unicode{x2013}$allowing for efficient verification processes. Additionally, we develop an optimization framework incorporating these stochastic dominance constraints. Numerical results validate the robustness of our method, showcasing its effectiveness for solving higher-order stochastic dominance problems, with applications to fields such as portfolio optimization.
