Scenario optimization

For other uses, see Scenario (disambiguation).

The scenario approach or scenario optimization approach is a technique for obtaining solutions to robust optimization and chance-constrained optimization problems based on a sample of the constraints. It also relates to inductive reasoning in modeling and decision-making. The technique has existed for decades as a heuristic approach and has more recently been given a systematic theoretical foundation.

In optimization, robustness features translate into constraints that are parameterized by the uncertain elements of the problem. In the scenario method, ^{[1] [2] [3]} a solution is obtained by only looking at a random sample of constraints (heuristic approach) called scenarios and a deeply-grounded theory tells the user how “robust” the corresponding solution is related to other constraints. This theory justifies the use of randomization in robust and chance-constrained optimization.

Data-driven optimization

At times, scenarios are obtained as random extractions from a model. More often, however, scenarios are instances of the uncertain constraints that are obtained as observations (data-driven science). In this latter case, no model of uncertainty is needed to generate scenarios. Moreover, most remarkably, also in this case scenario optimization comes accompanied by a full-fledged theory because all scenario optimization results are distribution-free and can therefore be applied even when a model of uncertainty is not available.

Theoretical results

For constraints that are convex (e.g. in semidefinite problems, involving LMIs (Linear Matrix Inequalities)), a deep theoretical analysis has been established which shows that the probability that a new constraint is not satisfied follows a distribution that is dominated by a Beta distribution. This result is tight since it is exact for a whole class of convex problems. ^[3] More generally, various empirical levels have been shown to follow a Dirichlet distribution, whose marginals are beta distribution. ^[4] The scenario approach with ${\displaystyle L_{1}}$ regularization has also been considered, ^[5] and handy algorithms with reduced computational complexity are available. ^[6] Extensions to more complex, non-convex, set-ups are still objects of active investigation.

Along the scenario approach, it is also possible to pursue a risk-return trade-off. ^{[7] [8]} Moreover, a full-fledged method can be used to apply this approach to control. ^[9] First ${\displaystyle N}$ constraints are sampled and then the user starts removing some of the constraints in succession. This can be done in different ways, even according to greedy algorithms. After elimination of one more constraint, the optimal solution is updated, and the corresponding optimal value is determined. As this procedure moves on, the user constructs an empirical “curve of values”, i.e. the curve representing the value achieved after the removing of an increasing number of constraints. The scenario theory provides precise evaluations of how robust the various solutions are.

A remarkable advance in the theory has been established by the recent wait-and-judge approach: ^[10] one assesses the complexity of the solution (as precisely defined in the referenced article) and from its value formulates precise evaluations on the robustness of the solution. These results shed light on deeply-grounded links between the concepts of complexity and risk. A related approach, named "Repetitive Scenario Design" aims at reducing the sample complexity of the solution by repeatedly alternating a scenario design phase (with reduced number of samples) with a randomized check of the feasibility of the ensuing solution. ^[11]

Example

Consider a function ${\displaystyle R_{\delta }(x)}$ which represents the return of an investment; it depends on our vector of investment choices ${\displaystyle x}$ and on the market state ${\displaystyle \delta }$ which will be experienced at the end of the investment period.

Given a stochastic model for the market conditions, we consider ${\displaystyle N}$ of the possible states ${\displaystyle \delta _{1},\dots ,\delta _{N}}$ (randomization of uncertainty). Alternatively, the scenarios ${\displaystyle \delta _{i}}$ can be obtained from a record of observations.

We set out to solve the scenario optimization program

${\displaystyle \max _{x}\min _{i=1,\dots ,N}R_{\delta _{i}}(x).\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ (1)}$

This corresponds to choosing a portfolio vector x so as to obtain the best possible return in the worst-case scenario. ^{[12] [13]}

After solving (1), an optimal investment strategy ${\displaystyle x^{\ast }}$ is achieved along with the corresponding optimal return ${\displaystyle R^{\ast }}$. While ${\displaystyle R^{\ast }}$ has been obtained by looking at ${\displaystyle N}$ possible market states only, the scenario theory tells us that the solution is robust up to a level ${\displaystyle \varepsilon }$, that is, the return ${\displaystyle R^{\ast }}$ will be achieved with probability ${\displaystyle 1-\varepsilon }$ for other market states.

In quantitative finance, the worst-case approach can be overconservative. One alternative is to discard some odd situations to reduce pessimism; ^[7] moreover, scenario optimization can be applied to other risk-measures including CVaR – Conditional Value at Risk – so adding to the flexibility of its use. ^[14]

Application fields

Fields of application include: prediction, systems theory, regression analysis (Interval Predictor Models in particular), Actuarial science, optimal control, financial mathematics, machine learning, decision making, supply chain, and management.

References

1. Calafiore, Giuseppe; Campi, M.C. (2005). "Uncertain convex programs: Randomized solutions and confidence levels". Mathematical Programming. 102: 25–46. doi:10.1007/s10107-003-0499-y. S2CID   1063933.
2. Calafiore, G.C.; Campi, M.C. (2006). "The Scenario Approach to Robust Control Design". IEEE Transactions on Automatic Control. 51 (5): 742–753. doi:10.1109/TAC.2006.875041. S2CID   49263.
3. Campi, M. C.; Garatti, S. (2008). "The Exact Feasibility of Randomized Solutions of Uncertain Convex Programs". SIAM Journal on Optimization. 19 (3): 1211–1230. doi:10.1137/07069821X.
4. Carè, A.; Garatti, S.; Campi, M. C. (2015). "Scenario Min-Max Optimization and the Risk of Empirical Costs". SIAM Journal on Optimization. 25 (4): 2061–2080. doi:10.1137/130928546. hdl: 11311/979283 .
5. Campi, M. C.; Carè, A. (2013). "Random Convex Programs with L₁-Regularization: Sparsity and Generalization". SIAM Journal on Control and Optimization. 51 (5): 3532–3557. doi:10.1137/110856204.
6. Carè, Algo; Garatti, Simone; Campi, Marco C. (2014). "FAST—Fast Algorithm for the Scenario Technique". Operations Research. 62 (3): 662–671. doi:10.1287/opre.2014.1257. hdl: 11311/937164 .
7. Campi, M. C.; Garatti, S. (2011). "A Sampling-and-Discarding Approach to Chance-Constrained Optimization: Feasibility and Optimality". Journal of Optimization Theory and Applications. 148 (2): 257–280. doi:10.1007/s10957-010-9754-6. S2CID   7856112.
8. Calafiore, Giuseppe Carlo (2010). "Random Convex Programs". SIAM Journal on Optimization. 20 (6): 3427–3464. doi:10.1137/090773490.
9. "Modulating robustness in control design: Principles and algorithms". IEEE Control Systems Magazine. 33 (2): 36–51. 2013. doi:10.1109/MCS.2012.2234964. S2CID   24072721.
10. Campi, M. C.; Garatti, S. (2018). "Wait-and-judge scenario optimization". Mathematical Programming. 167: 155–189. doi:10.1007/s10107-016-1056-9. hdl: 11311/1002492 . S2CID   39523265.
11. Calafiore, Giuseppe C. (2017). "Repetitive Scenario Design". IEEE Transactions on Automatic Control. 62 (3): 1125–1137. arXiv: 1602.03796 . doi:10.1109/TAC.2016.2575859. S2CID   47572451.
12. Pagnoncelli, B. K.; Reich, D.; Campi, M. C. (2012). "Risk-Return Trade-off with the Scenario Approach in Practice: A Case Study in Portfolio Selection". Journal of Optimization Theory and Applications. 155 (2): 707–722. doi:10.1007/s10957-012-0074-x. S2CID   1509645.
13. Calafiore, Giuseppe Carlo (2013). "Direct data-driven portfolio optimization with guaranteed shortfall probability". Automatica. 49 (2): 370–380. doi:10.1016/j.automatica.2012.11.012. S2CID   5762583.
14. Ramponi, Federico Alessandro; Campi, Marco C. (2018). "Expected shortfall: Heuristics and certificates". European Journal of Operational Research. 267 (3): 1003–1013. doi:10.1016/j.ejor.2017.11.022. S2CID   3553018.