DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Generating moment matching scenarios using optimization techniques

Abstract

An optimization based method is proposed to generate moment matching scenarios for numerical integration and its use in stochastic programming. The main advantage of the method is its flexibility: it can generate scenarios matching any prescribed set of moments of the underlying distribution rather than matching all moments up to a certain order, and the distribution can be defined over an arbitrary set. This allows for a reduction in the number of scenarios and allows the scenarios to be better tailored to the problem at hand. The method is based on a semi-infinite linear programming formulation of the problem that is shown to be solvable with polynomial iteration complexity. A practical column generation method is implemented. The column generation subproblems are polynomial optimization problems; however, they need not be solved to optimality. It is found that the columns in the column generation approach can be efficiently generated by random sampling. The number of scenarios generated matches a lower bound of Tchakaloff's. The rate of convergence of the approximation error is established for continuous integrands, and an improved bound is given for smooth integrands. Extensive numerical experiments are presented in which variants of the proposed method are compared to Monte Carlomore » and quasi-Monte Carlo methods on both numerical integration problems and stochastic optimization problems. The benefits of being able to match any prescribed set of moments, rather than all moments up to a certain order, is also demonstrated using optimization problems with 100-dimensional random vectors. Here, empirical results show that the proposed approach outperforms Monte Carlo and quasi-Monte Carlo based approaches on the tested problems.« less

Authors:
 [1];  [1]
  1. Northwestern Univ., Evanston, IL (United States)
Publication Date:
Research Org.:
Northwestern Univ., Evanston, IL (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
OSTI Identifier:
1321131
Grant/Contract Number:  
SC0005102
Resource Type:
Accepted Manuscript
Journal Name:
SIAM Journal on Optimization
Additional Journal Information:
Journal Volume: 23; Journal Issue: 2; Journal ID: ISSN 1052-6234
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; scenario generation; moment matching; cubature; column generation; convex programming; statistical bounds; semi-infinite programming

Citation Formats

Mehrotra, Sanjay, and Papp, Dávid. Generating moment matching scenarios using optimization techniques. United States: N. p., 2013. Web. doi:10.1137/110858082.
Mehrotra, Sanjay, & Papp, Dávid. Generating moment matching scenarios using optimization techniques. United States. https://doi.org/10.1137/110858082
Mehrotra, Sanjay, and Papp, Dávid. Thu . "Generating moment matching scenarios using optimization techniques". United States. https://doi.org/10.1137/110858082. https://www.osti.gov/servlets/purl/1321131.
@article{osti_1321131,
title = {Generating moment matching scenarios using optimization techniques},
author = {Mehrotra, Sanjay and Papp, Dávid},
abstractNote = {An optimization based method is proposed to generate moment matching scenarios for numerical integration and its use in stochastic programming. The main advantage of the method is its flexibility: it can generate scenarios matching any prescribed set of moments of the underlying distribution rather than matching all moments up to a certain order, and the distribution can be defined over an arbitrary set. This allows for a reduction in the number of scenarios and allows the scenarios to be better tailored to the problem at hand. The method is based on a semi-infinite linear programming formulation of the problem that is shown to be solvable with polynomial iteration complexity. A practical column generation method is implemented. The column generation subproblems are polynomial optimization problems; however, they need not be solved to optimality. It is found that the columns in the column generation approach can be efficiently generated by random sampling. The number of scenarios generated matches a lower bound of Tchakaloff's. The rate of convergence of the approximation error is established for continuous integrands, and an improved bound is given for smooth integrands. Extensive numerical experiments are presented in which variants of the proposed method are compared to Monte Carlo and quasi-Monte Carlo methods on both numerical integration problems and stochastic optimization problems. The benefits of being able to match any prescribed set of moments, rather than all moments up to a certain order, is also demonstrated using optimization problems with 100-dimensional random vectors. Here, empirical results show that the proposed approach outperforms Monte Carlo and quasi-Monte Carlo based approaches on the tested problems.},
doi = {10.1137/110858082},
journal = {SIAM Journal on Optimization},
number = 2,
volume = 23,
place = {United States},
year = {Thu May 16 00:00:00 EDT 2013},
month = {Thu May 16 00:00:00 EDT 2013}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 20 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Risk aversion in multistage stochastic programming: A modeling and algorithmic perspective
journal, February 2016

  • Homem-de-Mello, Tito; Pagnoncelli, Bernardo K.
  • European Journal of Operational Research, Vol. 249, Issue 1
  • DOI: 10.1016/j.ejor.2015.05.048

A Randomized Mirror-Prox Method for Solving Structured Large-Scale Matrix Saddle-Point Problems
journal, January 2013

  • Baes, Michel; Bürgisser, Michael; Nemirovski, Arkadi
  • SIAM Journal on Optimization, Vol. 23, Issue 2
  • DOI: 10.1137/11085801X

Artificial neural network-based methodology for short-term electric load scenario generation
conference, September 2015

  • Vagropoulos, Stylianos I.; Kardakos, Evaggelos G.; Simoglou, Christos K.
  • 2015 18th International Conference on Intelligent System Application to Power Systems (ISAP)
  • DOI: 10.1109/ISAP.2015.7325540

Monte Carlo sampling-based methods for stochastic optimization
journal, January 2014

  • Homem-de-Mello, Tito; Bayraksan, Güzin
  • Surveys in Operations Research and Management Science, Vol. 19, Issue 1
  • DOI: 10.1016/j.sorms.2014.05.001

ANN-based scenario generation methodology for stochastic variables of electric power systems
journal, May 2016

  • Vagropoulos, Stylianos I.; Kardakos, Evaggelos G.; Simoglou, Christos K.
  • Electric Power Systems Research, Vol. 134
  • DOI: 10.1016/j.epsr.2015.12.020

A Cutting Surface Algorithm for Semi-Infinite Convex Programming with an Application to Moment Robust Optimization
journal, January 2014

  • Mehrotra, Sanjay; Papp, Dávid
  • SIAM Journal on Optimization, Vol. 24, Issue 4
  • DOI: 10.1137/130925013

Containment Problems for Polytopes and Spectrahedra
journal, January 2013

  • Kellner, Kai; Theobald, Thorsten; Trabandt, Christian
  • SIAM Journal on Optimization, Vol. 23, Issue 2
  • DOI: 10.1137/120874898

An empirical analysis of scenario generation methods for stochastic optimization
journal, November 2016


Scenario generation for stochastic optimization problems via the sparse grid method
journal, April 2015

  • Chen, Michael; Mehrotra, Sanjay; Papp, Dávid
  • Computational Optimization and Applications, Vol. 62, Issue 3
  • DOI: 10.1007/s10589-015-9751-7

Multivariate simultaneous approximation
journal, December 2002


An encyclopaedia of cubature formulas
journal, June 2003


A new moment matching algorithm for sampling from partially specified symmetric distributions
journal, November 2008


The complexity of optimizing over a simplex, hypercube or sphere: a short survey
journal, December 2007


Scenario reduction in stochastic programming
journal, March 2003


Simulation and optimization approaches to scenario tree generation
journal, April 2004

  • Gülpınar, Nalan; Rustem, Berç; Settergren, Reuben
  • Journal of Economic Dynamics and Control, Vol. 28, Issue 7
  • DOI: 10.1016/S0165-1889(03)00113-1

Likelihood approximation by numerical integration on sparse grids
journal, May 2008


Scenario tree modeling for multistage stochastic programs
journal, November 2007


GloptiPoly: Global optimization over polynomials with Matlab and SeDuMi
journal, June 2003

  • Henrion, Didier; Lasserre, Jean-Bernard
  • ACM Transactions on Mathematical Software, Vol. 29, Issue 2
  • DOI: 10.1145/779359.779363

Global Optimization with Polynomials and the Problem of Moments
journal, January 2001


An Approximate Method for Sampling Correlated Random Variables from Partially-Specified Distributions
journal, February 1998


Epi-convergent discretizations of stochastic programs via integration quadratures
journal, February 2005


Discrepancy and integration of continuous functions
journal, February 1988


Scenario Generation
book, February 2011


A comparison between (quasi-)Monte Carlo and cubature rule based methods for solving high-dimensional integration problems
journal, March 2003


Sums of Squares and Semidefinite Program Relaxations for Polynomial Optimization Problems with Structured Sparsity
journal, January 2006

  • Waki, Hayato; Kim, Sunyoung; Kojima, Masakazu
  • SIAM Journal on Optimization, Vol. 17, Issue 1
  • DOI: 10.1137/050623802

Approximating the complexity measure of Vavasis-Ye algorithm is NP-hard
journal, September 1999


Works referencing / citing this record:

A copula-based scenario tree generation algorithm for multiperiod portfolio selection problems
journal, January 2019


Scenario generation for stochastic optimization problems via the sparse grid method
journal, April 2015

  • Chen, Michael; Mehrotra, Sanjay; Papp, Dávid
  • Computational Optimization and Applications, Vol. 62, Issue 3
  • DOI: 10.1007/s10589-015-9751-7

An empirical analysis of scenario generation methods for stochastic optimization
journal, November 2016


Scenario generation in stochastic programming using principal component analysis based on moment-matching approach
journal, October 2019


Performance Comparison of Scenario-Generation Methods Applied to a Stochastic Optimization Asset-Liability Management Model
journal, April 2018


ANN-based scenario generation methodology for stochastic variables of electric power systems
journal, May 2016

  • Vagropoulos, Stylianos I.; Kardakos, Evaggelos G.; Simoglou, Christos K.
  • Electric Power Systems Research, Vol. 134
  • DOI: 10.1016/j.epsr.2015.12.020

Risk aversion in multistage stochastic programming: A modeling and algorithmic perspective
journal, February 2016

  • Homem-de-Mello, Tito; Pagnoncelli, Bernardo K.
  • European Journal of Operational Research, Vol. 249, Issue 1
  • DOI: 10.1016/j.ejor.2015.05.048