Stability-based generation of scenario trees for multistage stochastic programs
|
|
- Cordelia Fitzgerald
- 5 years ago
- Views:
Transcription
1 Stability-based generation of scenario trees for multistage stochastic programs H. Heitsch and W. Römisch Humboldt-University Berlin Institute of Mathematics Berlin, Germany Page 1 of 22 ISMP 2006, Rio de Janeiro, July 30 August 5, 2006
2 Multistage stochastic programs Let ξ = {ξ t } T t=1 be an IR d -valued discrete-time stochastic process defined on some probability space (Ω, F, IP ) and with ξ 1 deterministic. The stochastic decision x t at period t is assumed to be measurable with respect to the σ-field F t (ξ) := σ(ξ 1,..., ξ t ) (nonanticipativity). Multistage stochastic program: T min IE[ x t X t, b t (ξ t ), x t ] x t is F t (ξ) measurable, t = 1,..., T, t=1 A t,0 x t + A t,1 (ξ t )x t 1 = h t (ξ t ), t = 2,..., T where X t are nonempty and polyhedral sets, A t,0 are fixed recourse matrices and b t ( ), h t ( ) and A t,1 ( ) are affine functions depending on ξ t, where ξ varies in a polyhedral set Ξ. If the process {ξ t } T t=1 has a finite number of scenarios, they exhibit a scenario tree structure. Page 2 of 22
3 To have the multistage stochastic program well defined, we assume x t L r (Ω, F, IP ; IR m t ) and ξ t L r (Ω, F, IP ; IR d ), where r 1 and r := r r 1, if costs are random r, if only right-hand sides are random, if all technology matrices are random and r = T. Then nonanticipativity may be expressed as x N r (ξ) N r (ξ) = {x T t=1l r (Ω, F, IP ; IR m t ) : x t = IE[x t F t (ξ)], t}, i.e., as a subspace constraint, by using the conditional expectations IE[ F t (ξ)]. Page 3 of 22 For T = 2 we have N r (ξ) = IR m 1 L r (Ω, F, P ; IR m 2 ). infinite-dimensional optimization problem
4 Data process approximation by scenario trees The process {ξ t } T t=1 is approximated by a process forming a scenario tree being based on a finite set N IN of nodes. n = 1 t = 1 ξ n n n t = 2 N+ (n) N T t(n) T Page 4 of 22 Scenario tree with T = 5, N = 22 and 11 leaves n = 1 root node, n unique predecessor of node n, path(n) = {1,..., n, n}, t(n) := path(n), N + (n) set of successors to n, N T := {n N : N + (n) = } set of leaves, path(n), n N T, scenario with (given) probability π n, π n := ν N + (n) πν probability of node n, ξ n realization of ξ t(n).
5 Tree representation of the optimization model { } min π n b t(n) (ξ n ), x n xn X t(n), n N, A 1,0 x 1 = h 1 (ξ 1 ) A t(n),0 x n + A t(n),1 x n =h t(n) (ξ n ), n N n N How to solve the optimization model? - Standard software (e.g., CPLEX) - Decomposition methods for (very) large scale models (Ruszczynski/Shapiro (Eds.): Stochastic Programming, Handbook, 2003) Page 5 of 22 Questions: Under which conditions and in which sense do multistage models behave stable with respect to perturbations of ξ? Can such stability results be used to generate (multivariate) scenario trees?
6 Dynamic programming Theorem: (Evstigneev 76, Rockafellar/Wets 76) Under weak assumptions the multistage stochastic program is equivalent to the (first-stage) convex minimization problem { min f(x 1, ξ)p (dξ) : x 1 X 1 }, Ξ where f is an integrand on IR m 1 Ξ given by f(x 1, ξ):= b 1 (ξ 1 ), x 1 + Φ 2 (x 1, ξ 2 ), Φ t (x 1,..., x t 1, ξ t ):=inf{ b t (ξ t ), x t +IE [ Φ t+1 (x 1,..., x t, ξ t+1 ) F t ] : x t X t, A t,0 x t + A t,1 (ξ t )x t 1 = h t (ξ t )} for t = 2,..., T, where Φ T +1 (x 1,..., x T, ξ T +1 ) := 0. The integrand f depends on the probability measure IP and, thus, also on the probability distribution P = IP ξ 1 of ξ in a nonlinear way! Hence, earlier approaches to stability fail! Page 6 of 22
7 Quantitative Stability Let us introduce some notations. Let F denote the objective function defined on L r (Ω, F, IP ; IR s ) L r (Ω, F, IP ; IR m ) IR by F (ξ, x) := IE[ T t=1 b t(ξ t ), x t ], let X t (x t 1 ; ξ t ) := {x t X t A t,0 x t + A t,1 (ξ t )x t 1 = h t (ξ t )} denote the t-th feasibility set for every t = 2,..., T and X (ξ) := {x L r (Ω, F, IP ; IR m ) x 1 X 1, x t X t (x t 1 ; ξ t )} the set of feasible elements with input ξ. Then the multistage stochastic program may be rewritten as min{f (ξ, x) : x X (ξ) N r (ξ)}. Let v(ξ) denote its optimal value and, for any α 0, l α (F (ξ, )) := {x X (ξ) N r (ξ) : F (ξ, x) v(ξ) + α} S(ξ) := l 0 (F (ξ, )) denote the α-level set and the solution set of the stochastic program with input ξ. Page 7 of 22
8 The following conditions are imposed: (A1) ξ L r (Ω, F, IP ; IR s ) for some r 1. (A2) There exists a δ > 0 such that for any ξ L r (Ω, F, IP ; IR s ) with ξ ξ r δ, any t = 2,..., T and any x 1 X 1, x τ X τ (x τ 1 ; ξ τ ), τ = 2,..., t 1, the set X t (x t 1 ; ξ t ) is nonempty (relatively complete recourse locally around ξ). (A3) The optimal values v( ξ) are finite for all ξ L r (Ω, F, IP ; IR s ) with ξ ξ r δ and the objective function F is level-bounded locally uniformly at ξ, i.e., for some α > 0 there exists a δ > 0 and a bounded subset B of L r (Ω, F, IP ; IR m ) such that l α (F ( ξ, )) is nonempty and contained in B for all ξ Lr (Ω, F, IP ; IR s ) with ξ ξ r δ. Page 8 of 22 Norm in L r : ξ r := ( T IE[ ξ t r ]) 1 r t=1
9 Theorem: (Heitsch/Römisch/Strugarek, SIAM J. Opt. 2006) Let (A1), (A2) and (A3) be satisfied, r > 1 and X 1 be bounded. Then there exist positive constants L and δ such that v(ξ) v( ξ) L( ξ ξ r + D f (ξ, ξ)) holds for all ξ L r (Ω, F, IP ; IR s ) with ξ ξ r δ. Assume that technology matrices are non-random and the solution x of the original problem is unique. If (ξ (n) ) is a sequence in T t=1l r (Ω, F t (ξ), IP ; IR s ) such that ξ (n) ξ r and D f (ξ (n), ξ) converge to 0 and if (x (n) ) is a sequence of solutions of the approximate problems, then the sequence (x (n) ) converges to x with respect to the weak topology in L r. Page 9 of 22 Here, D f (ξ, ξ) denotes the filtration distance of ξ and ξ defined by D f (ξ, ξ)= T 1 inf max{ x t IE[x t F t ( ξ)] r, x t IE[ x t F t (ξ)] r }. x S(ξ) x S( ξ) t=2
10 Remark: The continuity property of infima in the Theorem can be supplemented by a quantitative stability property of the set S 1 (ξ) of first stage solutions. Namely, there exists a constant ˆL > 0 such that sup x S 1 ( ξ) d(x, S 1 (ξ)) Ψ 1 ξ ( ˆL( ξ ξ r + D f (ξ, ξ))), where Ψ ξ (τ) := inf {IE[f(x 1, ξ)] v(ξ) : d(x 1, S 1 (ξ)) τ, x 1 X 1 } with Ψ 1 ξ (α) := sup{τ IR + : Ψ ξ (τ) α} (α IR + ) is the growth function of the original problem near its solution set S 1 (ξ). Page 10 of 22 Remark: Simple examples show that the filtration distance is indispensable for the stability result to hold.
11 Generation of scenario trees (i) In most practical situations scenarios ξ i with known probabilities p i, i = 1,..., N, can be generated, e.g., simulation scenarios from (parametric or nonparametric) statistical models of ξ or (nearly) optimal quantizations of the probability distribution of ξ. (ii) Construction of a scenario tree out of the scenarios ξ i with probabilities p i, i = 1,..., N,. Page 11 of 22
12 Approaches for (ii): (1) Bound-based approximation methods, (Frauendorfer 96, Kuhn 05, Edirisinghe 99, Casey/Sen 05). (2) Monte Carlo-based schemes (inside or outside decomposition methods) (e.g. Shapiro 03, 06, Higle/Rayco/Sen 01, Chiralaksanakul/Morton 04). (3) the use of Quasi Monte Carlo integration quadratures (Pennanen 05, 06). (4) EVPI-based sampling schemes (inside decomposition schemes) (Corvera Poire 95, Dempster 04). (5) Moment-matching principle (Høyland/Wallace 01, Høyland/Kaut/Wallace 03). (6) (Nearly) best approximations based on probability metrics (Pflug 01, Hochreiter/Pflug 02, Mirkov/Pflug 06; Gröwe-Kuska/Heitsch/Römisch 01, 03, Heitsch/Römisch 05). Page 12 of 22 Survey: Dupačová/Consigli/Wallace 00
13 Constructing scenario trees Let ξ be the original stochastic process on some probability space (Ω, F, IP ) with parameter set {1,..., T } and state space IR d. We aim at generating a scenario tree ξ tr such that and, thus, ξ ξ tr r and D f (ξ, ξ tr ) v(ξ) v(ξ tr ) are small. To determine such a scenario tree, we start with a discrete approximation ξ f consisting of scenarios ξ i = (ξ1, i..., ξt i ) with probabilities p i, i = 1,..., N. ξ f is a fan of individual scenarios. Page 13 of 22 t = 1 t = 2 t = 3 t = 4 t = 5
14 The fan ξ f is chosen such that it is adapted to the filtration (F t (ξ)) T t=1 and ξ ξ f r ε appr. Algorithms are developed that generate a scenario tree ξ tr by deleting and bundling scenarios of ξ f (that are similar at t) such that it is also adapted to the filtration (F t (ξ)) T t=1 and satisfies Since it holds (1) ξ f ξ tr r ε r T 1 (2) inf x t IE[x t F t (ξ tr )] r ε f. x S(ξ f ) t=2 D f (ξ, ξ tr ) ε appr + T 1 inf x t IE[x t F t (ξ tr )] r, x S(ξ f ) t=2 if ξ f is sufficiently close to ξ, we obtain in case ε appr + ε r δ that v(ξ) v(ξ tr ) L(2ε appr + ε r + ε f ). Page 14 of 22
15 (1) Forward tree generation Let scenarios ξ i with probabilities p i, i = 1,..., N, fixed root ξ1 IR d, r 1, and tolerances ε r, ε t, t = 2,..., T, be given such that T ε t ε r. t=2 Step 1: Set ˆξ 1 := ξ f and C 1 = {I = {1,..., N}}. Step t: Let C t 1 = {Ct 1, 1..., C K t 1 t 1 }. Determine disjoint index sets It k and Jt k of remaining and deleted scenarios such that It k Jt k = Ct 1, k a mapping α t : I I { i k α t (j) = t (j), j Jt k, k = 1,..., K t 1, j, otherwise, where i k t (j) I k t such that i k t (j) arg min i I k t ˆξ t 1,i ˆξ t 1,j t, Page 15 of 22
16 a stochastic process ˆξ t ˆξ t,i τ = { ξ α τ (i) τ, τ t,, otherwise, ξ i τ such that ˆξ t ˆξ t 1 r,t ε t. Set I t := K t 1 k=1 Ik t and C t := {α 1 t (i) : i I k t, k = 1,..., K t 1 }. Step T+1: Let C T = {CT 1,..., CK T T }. Construct a stochastic process ξ tr having K T scenarios ξtr k such that ξtr,t k := ξ α t(i) t with probabilities πt i = p j if i CT k, k = 1,..., K T, t = 2,..., T. Proposition: j C k T ξ f ξ tr r T ε t ε r. t=2 Page 16 of 22
17 t = 1 t = 2 t = 3 t = 4 t = 5 t = 1 t = 2 t = 3 t = 4 t = 5 t = 1 t = 2 t = 3 t = 4 t = 5 Page 17 of 22 t = 1 t = 2 t = 3 t = 4 t = 5 t = 1 t = 2 t = 3 t = 4 t = 5 t = 1 t = 2 t = 3 t = 4 t = 5 Illustration of the forward tree construction for an example including T=5 time periods starting with a scenario fan containing N=58 scenarios <Start Animation>
18 (2) Bounding approximate filtration distances T 1 Aim: (ξ f, ξ tr ) := inf x t IE[x t F t (ξ tr )] r ε f x S(ξ f ) t=2 Two possibilities: (i) Estimates in terms of some solutions with input ξ f, which would require to solve a two-stage model. (ii) Estimates in terms of the input ξ f. Proposition: Let (A2) and (A3) be satisfied and X 1 be bounded. Assume that the technology matrices A t,1 are non-random, 1 r < and ξ f is sufficiently close to ξ. Then there exists a constant ˆL > 0 such that Condition: (ξ f, ξ tr ) ˆL ( i I 2 i I 2 j I 2,i p j ξ j ξ i r j I 2,i p j ξ j ξ i r ε r f ) 1 r Page 18 of 22
19 Numerical experience We consider the electricity portfolio management of a municipal power company. Data was available on the electrical load demand and on electricity prices at the market place EEX. A multivariate statistical model is developed for the yearly demandprice process ξ that allowed to generate yearly demand-price scenarios ξ i, with probabilities p i = 1 N, i = 1,..., N. These scenarios are assumed to form the process ξ f. Branching in ξ tr was allowed at most monthly. The tolerances ε t at branching points were chosen such that Page 19 of 22 ε t = ε T [1 + q(1 2 t )], t = 2,..., T, T where the parameter q [0, 1] affects the branching structure of the constructed trees. For the test runs we used q = 0.6. The test runs were performed on a PC with a 3 GHz Intel Pentium CPU and 1 GByte main memory.
20 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec a) Forward tree construction with relative filtration tolerance ε rel,f = 0.35 Page 20 of 22 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec b) Forward tree construction with relative filtration tolerance ε rel,f = 0.45 Yearly demand-price scenario trees with relative tolerance ε rel,r = 0.25
21 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec a) Forward tree construction with relative filtration tolerance ε rel,f = 0.6 Page 21 of 22 Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec b) Forward tree construction with relative filtration tolerance ε rel,f = 0.7 Yearly demand-price scenario trees with relative tolerance ε rel,r = 0.5
22 ε rel,r ε rel,f Scenarios Nodes Stages Time (sec) Numerical results for yearly demand-price scenario trees Page 22 of 22
Stability of Stochastic Programming Problems
Stability of Stochastic Programming Problems W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch Page 1 of 35 Spring School Stochastic
More informationStochastic Programming: From statistical data to optimal decisions
Stochastic Programming: From statistical data to optimal decisions W. Römisch Humboldt-University Berlin Department of Mathematics (K. Emich, H. Heitsch, A. Möller) Page 1 of 24 6th International Conference
More informationJitka Dupačová and scenario reduction
Jitka Dupačová and scenario reduction W. Römisch Humboldt-University Berlin Institute of Mathematics http://www.math.hu-berlin.de/~romisch Session in honor of Jitka Dupačová ICSP 2016, Buzios (Brazil),
More informationA note on scenario reduction for two-stage stochastic programs
A note on scenario reduction for two-stage stochastic programs Holger Heitsch a and Werner Römisch a a Humboldt-University Berlin, Institute of Mathematics, 199 Berlin, Germany We extend earlier work on
More informationAirline Network Revenue Management by Multistage Stochastic Programming
Airline Network Revenue Management by Multistage Stochastic Programming K. Emich, H. Heitsch, A. Möller, W. Römisch Humboldt-University Berlin, Department of Mathematics Page 1 of 18 GOR-Arbeitsgruppe
More informationScenario generation in stochastic programming
Scenario generation in stochastic programming Werner Römisch Humboldt-University Berlin, Institute of Mathematics, 10099 Berlin, Germany, romisch@math.hu-berlin.de Abstract Stability-based methods for
More informationScenario tree generation for multi-stage stochastic programs
Scenario tree generation for multi-stage stochastic programs Holger Heitsch 1 and Werner Römisch 2 1 Humboldt University Berlin, Institute of Mathematics, 10099 Berlin, Germany, heitsch@math.hu-berlin.de
More informationScenario Reduction: Old and New Results
Scenario Reduction: Old and New Results W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch (H. Heitsch, R. Henrion, C. Küchler) Workshop
More informationScenario generation in stochastic programming with application to energy systems
Scenario generation in stochastic programming with application to energy systems W. Römisch Humboldt-University Berlin Institute of Mathematics www.math.hu-berlin.de/~romisch SESO 2017 International Thematic
More informationGAMINGRE 8/1/ of 7
FYE 09/30/92 JULY 92 0.00 254,550.00 0.00 0 0 0 0 0 0 0 0 0 254,550.00 0.00 0.00 0.00 0.00 254,550.00 AUG 10,616,710.31 5,299.95 845,656.83 84,565.68 61,084.86 23,480.82 339,734.73 135,893.89 67,946.95
More informationStability of optimization problems with stochastic dominance constraints
Stability of optimization problems with stochastic dominance constraints D. Dentcheva and W. Römisch Stevens Institute of Technology, Hoboken Humboldt-University Berlin www.math.hu-berlin.de/~romisch SIAM
More informationScenario Generation in Stochastic Programming with Application to Optimizing Electricity Portfolios under Uncertainty
Scenario Generation in Stochastic Programming with Application to Optimizing Electricity Portfolios under Uncertainty W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch
More informationScenario generation Contents
Scenario generation W. Römisch Humboldt-University Berlin Department of Mathematics www.math.hu-berlin.de/~romisch Page 1 of 39 Tutorial, ICSP12, Halifax, August 14, 2010 (1) Aproximation issues (2) Generation
More informationOperations Research Letters. On a time consistency concept in risk averse multistage stochastic programming
Operations Research Letters 37 2009 143 147 Contents lists available at ScienceDirect Operations Research Letters journal homepage: www.elsevier.com/locate/orl On a time consistency concept in risk averse
More informationConditioning of linear-quadratic two-stage stochastic programming problems
Conditioning of linear-quadratic two-stage stochastic programming problems W. Römisch Humboldt-University Berlin Institute of Mathematics http://www.math.hu-berlin.de/~romisch (K. Emich, R. Henrion) (WIAS
More informationStochastic optimization, multivariate numerical integration and Quasi-Monte Carlo methods. W. Römisch
Stochastic optimization, multivariate numerical integration and Quasi-Monte Carlo methods W. Römisch Humboldt-University Berlin Institute of Mathematics www.math.hu-berlin.de/~romisch Chemnitzer Mathematisches
More informationComplexity of two and multi-stage stochastic programming problems
Complexity of two and multi-stage stochastic programming problems A. Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA The concept
More informationarxiv: v2 [math.oc] 18 Nov 2017
DASC: A DECOMPOSITION ALGORITHM FOR MULTISTAGE STOCHASTIC PROGRAMS WITH STRONGLY CONVEX COST FUNCTIONS arxiv:1711.04650v2 [math.oc] 18 Nov 2017 Vincent Guigues School of Applied Mathematics, FGV Praia
More informationOn the Scenario-Tree Optimal-Value Error for Stochastic Programming Problems
On the Scenario-Tree Optimal-Value Error for Stochastic Programming Problems Julien Keutchayan David Munger Michel Gendreau January 2017 CIRRELT-2017-05 On the Scenario-Tree Optimal-Value Error for Stochastic
More informationIntroduction to Stochastic Optimization Part 4: Multi-stage decision
Introduction to Stochastic Optimization Part 4: Multi-stage decision problems April 23, 29 The problem ξ = (ξ,..., ξ T ) a multivariate time series process (e.g. future interest rates, future asset prices,
More informationStochastic Programming: Tutorial Part II
Stochastic Programming: Tutorial Part II W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch Page 1 of 77 Optimization: Computation,
More informationApproximation of stochastic optimization problems and scenario generation
Approximation of stochastic optimization problems and scenario generation W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany http://www.math.hu-berlin.de/~romisch Page
More informationLecture 1. Stochastic Optimization: Introduction. January 8, 2018
Lecture 1 Stochastic Optimization: Introduction January 8, 2018 Optimization Concerned with mininmization/maximization of mathematical functions Often subject to constraints Euler (1707-1783): Nothing
More informationCONVERGENCE ANALYSIS OF SAMPLING-BASED DECOMPOSITION METHODS FOR RISK-AVERSE MULTISTAGE STOCHASTIC CONVEX PROGRAMS
CONVERGENCE ANALYSIS OF SAMPLING-BASED DECOMPOSITION METHODS FOR RISK-AVERSE MULTISTAGE STOCHASTIC CONVEX PROGRAMS VINCENT GUIGUES Abstract. We consider a class of sampling-based decomposition methods
More informationScenario estimation and generation
October 10, 2004 The one-period case Distances of Probability Measures Tensor products of trees Tree reduction A decision problem is subject to uncertainty Uncertainty is represented by probability To
More informationIn Centre, Online Classroom Live and Online Classroom Programme Prices
In Centre, and Online Classroom Programme Prices In Centre Online Classroom Foundation Certificate Bookkeeping Transactions 430 325 300 Bookkeeping Controls 320 245 225 Elements of Costing 320 245 225
More informationOptimal scenario generation and reduction in stochastic programming
Noname manuscript No. (will be inserted by the editor) Optimal scenario generation and reduction in stochastic programming R. Henrion W. Römisch the date of receipt and acceptance should be inserted later
More informationBirgit Rudloff Operations Research and Financial Engineering, Princeton University
TIME CONSISTENT RISK AVERSE DYNAMIC DECISION MODELS: AN ECONOMIC INTERPRETATION Birgit Rudloff Operations Research and Financial Engineering, Princeton University brudloff@princeton.edu Alexandre Street
More informationMonthly Trading Report July 2018
Monthly Trading Report July 218 Figure 1: July 218 (% change over previous month) % Major Market Indicators 2 2 4 USEP Forecasted Demand CCGT/Cogen/Trigen Supply ST Supply Figure 2: Summary of Trading
More informationGENERALIZED second-order cone complementarity
Stochastic Generalized Complementarity Problems in Second-Order Cone: Box-Constrained Minimization Reformulation and Solving Methods Mei-Ju Luo and Yan Zhang Abstract In this paper, we reformulate the
More informationSYSTEM BRIEF DAILY SUMMARY
SYSTEM BRIEF DAILY SUMMARY * ANNUAL MaxTemp NEL (MWH) Hr Ending Hr Ending LOAD (PEAK HOURS 7:00 AM TO 10:00 PM MON-SAT) ENERGY (MWH) INCREMENTAL COST DAY DATE Civic TOTAL MAXIMUM @Max MINIMUM @Min FACTOR
More informationAnnual Average NYMEX Strip Comparison 7/03/2017
Annual Average NYMEX Strip Comparison 7/03/2017 To Year to Year Oil Price Deck ($/bbl) change Year change 7/3/2017 6/1/2017 5/1/2017 4/3/2017 3/1/2017 2/1/2017-2.7% 2017 Average -10.4% 47.52 48.84 49.58
More informationANOVA decomposition of convex piecewise linear functions
ANOVA decomposition of convex piecewise linear functions W. Römisch Abstract Piecewise linear convex functions arise as integrands in stochastic programs. They are Lipschitz continuous on their domain,
More informationAsymptotics of minimax stochastic programs
Asymptotics of minimax stochastic programs Alexander Shapiro Abstract. We discuss in this paper asymptotics of the sample average approximation (SAA) of the optimal value of a minimax stochastic programming
More informationMIDAS: A Mixed Integer Dynamic Approximation Scheme
MIDAS: A Mixed Integer Dynamic Approximation Scheme Andy Philpott, Faisal Wahid, Frédéric Bonnans May 7, 2016 Abstract Mixed Integer Dynamic Approximation Scheme (MIDAS) is a new sampling-based algorithm
More informationMonthly Trading Report Trading Date: Dec Monthly Trading Report December 2017
Trading Date: Dec 7 Monthly Trading Report December 7 Trading Date: Dec 7 Figure : December 7 (% change over previous month) % Major Market Indicators 5 4 Figure : Summary of Trading Data USEP () Daily
More informationTime Series Analysis
Time Series Analysis A time series is a sequence of observations made: 1) over a continuous time interval, 2) of successive measurements across that interval, 3) using equal spacing between consecutive
More informationSYSTEM BRIEF DAILY SUMMARY
SYSTEM BRIEF DAILY SUMMARY * ANNUAL MaxTemp NEL (MWH) Hr Ending Hr Ending LOAD (PEAK HOURS 7:00 AM TO 10:00 PM MON-SAT) ENERGY (MWH) INCREMENTAL COST DAY DATE Civic TOTAL MAXIMUM @Max MINIMUM @Min FACTOR
More informationA look into the factor model black box Publication lags and the role of hard and soft data in forecasting GDP
A look into the factor model black box Publication lags and the role of hard and soft data in forecasting GDP Marta Bańbura and Gerhard Rünstler Directorate General Research European Central Bank November
More informationWeierstraß-Institut. für Angewandte Analysis und Stochastik. Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN
Weierstraß-Institut für Angewandte Analysis und Stochastik Leibniz-Institut im Forschungsverbund Berlin e. V. Preprint ISSN 2198-5855 Are Quasi-Monte Carlo algorithms efficient for two-stage stochastic
More informationRegularized optimization techniques for multistage stochastic programming
Regularized optimization techniques for multistage stochastic programming Felipe Beltrán 1, Welington de Oliveira 2, Guilherme Fredo 1, Erlon Finardi 1 1 UFSC/LabPlan Universidade Federal de Santa Catarina
More informationUpper bound for optimal value of risk averse multistage problems
Upper bound for optimal value of risk averse multistage problems Lingquan Ding School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, GA 30332-0205 Alexander Shapiro School
More informationDynamic Generation of Scenario Trees
Dynamic Generation of Scenario Trees Georg Ch. Pflug and Alois Pichler June 16, 015 Abstract We present new algorithms for the dynamic generation of scenario trees for multistage stochastic optimization.
More informationDEPARTMENT OF THE ARMY MILITARY SURFACE DEPLOYMENT AND DISTRIBUTION COMMAND (SDDC) 1 SOLDIER WAY SCOTT AFB, IL 62225
DEPARTMENT OF THE ARMY MILITARY SURFACE DEPLOYMENT AND DISTRIBUTION COMMAND (SDDC) 1 SOLDIER WAY SCOTT AFB, IL 62225 SDDC Operations Special Requirements Branch 1 Soldier Way Scott AFB, IL 62225 April
More informationStochastic Integer Programming
IE 495 Lecture 20 Stochastic Integer Programming Prof. Jeff Linderoth April 14, 2003 April 14, 2002 Stochastic Programming Lecture 20 Slide 1 Outline Stochastic Integer Programming Integer LShaped Method
More informationA. Shapiro Introduction
ESAIM: PROCEEDINGS, December 2003, Vol. 13, 65 73 J.P. Penot, Editor DOI: 10.1051/proc:2003003 MONTE CARLO SAMPLING APPROACH TO STOCHASTIC PROGRAMMING A. Shapiro 1 Abstract. Various stochastic programming
More informationAverage 175, , , , , , ,046 YTD Total 1,098,649 1,509,593 1,868,795 1,418, ,169 1,977,225 2,065,321
AGRICULTURE 01-Agriculture JUL 2,944-4,465 1,783-146 102 AUG 2,753 6,497 5,321 1,233 1,678 744 1,469 SEP - 4,274 4,183 1,596 - - 238 OCT 2,694 - - 1,032 340-276 NOV 1,979-5,822 637 3,221 1,923 1,532 DEC
More informationAverage 175, , , , , , ,940 YTD Total 944,460 1,284,944 1,635,177 1,183, ,954 1,744,134 1,565,640
AGRICULTURE 01-Agriculture JUL 2,944-4,465 1,783-146 102 AUG 2,753 6,497 5,321 1,233 1,678 744 1,469 SEP - 4,274 4,183 1,596 - - 238 OCT 2,694 - - 1,032 340-276 NOV 1,979-5,822 637 3,221 1,923 1,532 DEC
More informationOn a class of minimax stochastic programs
On a class of minimax stochastic programs Alexander Shapiro and Shabbir Ahmed School of Industrial & Systems Engineering Georgia Institute of Technology 765 Ferst Drive, Atlanta, GA 30332. August 29, 2003
More informationScenario Grouping and Decomposition Algorithms for Chance-constrained Programs
Scenario Grouping and Decomposition Algorithms for Chance-constrained Programs Siqian Shen Dept. of Industrial and Operations Engineering University of Michigan Joint work with Yan Deng (UMich, Google)
More informationQuantitative Stability of Two Stage Linear Second-order Conic Stochastic Programs with Full Random Recourse
Quantitative Stability of Two Stage Linear Second-order Conic Stochastic Programs with Full Random Recourse Qingsong Duan, Mengwei Xu, Shaoyan Guo, and Liwei Zhang October 7, 2017 Abstract In this paper,
More informationA Tighter Variant of Jensen s Lower Bound for Stochastic Programs and Separable Approximations to Recourse Functions
A Tighter Variant of Jensen s Lower Bound for Stochastic Programs and Separable Approximations to Recourse Functions Huseyin Topaloglu School of Operations Research and Information Engineering, Cornell
More informationConvex Optimization on Large-Scale Domains Given by Linear Minimization Oracles
Convex Optimization on Large-Scale Domains Given by Linear Minimization Oracles Arkadi Nemirovski H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute of Technology Joint research
More informationTechnical note on seasonal adjustment for Capital goods imports
Technical note on seasonal adjustment for Capital goods imports July 1, 2013 Contents 1 Capital goods imports 2 1.1 Additive versus multiplicative seasonality..................... 2 2 Steps in the seasonal
More informationLong-term Water Quality Monitoring in Estero Bay
Long-term Water Quality Monitoring in Estero Bay Keith Kibbey Laboratory Director Lee County Environmental Laboratory Division of Natural Resource Management Estero Bay Monitoring Programs Three significant
More informationApproximating Fixed-Horizon Forecasts Using Fixed-Event Forecasts
Approximating Fixed-Horizon Forecasts Using Fixed-Event Forecasts Malte Knüppel and Andreea L. Vladu Deutsche Bundesbank 9th ECB Workshop on Forecasting Techniques 4 June 216 This work represents the authors
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Data Extension & Forecasting Moving
More informationMinimax and risk averse multistage stochastic programming
Minimax and risk averse multistage stochastic programming Alexander Shapiro School of Industrial & Systems Engineering, Georgia Institute of Technology, 765 Ferst Drive, Atlanta, GA 30332. Abstract. In
More informationJayalath Ekanayake Jonas Tappolet Harald Gall Abraham Bernstein. Time variance and defect prediction in software projects: additional figures
Jayalath Ekanayake Jonas Tappolet Harald Gall Abraham Bernstein TECHNICAL REPORT No. IFI-2.4 Time variance and defect prediction in software projects: additional figures 2 University of Zurich Department
More informationSOME HISTORY OF STOCHASTIC PROGRAMMING
SOME HISTORY OF STOCHASTIC PROGRAMMING Early 1950 s: in applications of Linear Programming unknown values of coefficients: demands, technological coefficients, yields, etc. QUOTATION Dantzig, Interfaces
More information2016 Year-End Benchmark Oil and Gas Prices (Average of Previous 12 months First-Day-of-the Month [FDOM] Prices)
Oil and Gas Benchmark Prices to Estimate Year-End Petroleum Reserves and Values Using U.S. Securities and Exchange Commission Guidelines from the Modernization of Oil and Gas Reporting Effective January
More informationDecomposability and time consistency of risk averse multistage programs
Decomposability and time consistency of risk averse multistage programs arxiv:1806.01497v1 [math.oc] 5 Jun 2018 A. Shapiro School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta,
More informationIntroduction to Course
.. Introduction to Course Oran Kittithreerapronchai 1 1 Department of Industrial Engineering, Chulalongkorn University Bangkok 10330 THAILAND last updated: September 17, 2016 COMP METH v2.00: intro 1/
More informationValid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints
Valid Inequalities and Restrictions for Stochastic Programming Problems with First Order Stochastic Dominance Constraints Nilay Noyan Andrzej Ruszczyński March 21, 2006 Abstract Stochastic dominance relations
More informationSplitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches
Splitting Techniques in the Face of Huge Problem Sizes: Block-Coordinate and Block-Iterative Approaches Patrick L. Combettes joint work with J.-C. Pesquet) Laboratoire Jacques-Louis Lions Faculté de Mathématiques
More informationReformulation of chance constrained problems using penalty functions
Reformulation of chance constrained problems using penalty functions Martin Branda Charles University in Prague Faculty of Mathematics and Physics EURO XXIV July 11-14, 2010, Lisbon Martin Branda (MFF
More informationOn the maximum and minimum response to an impulse in Svars
On the maximum and minimum response to an impulse in Svars Bulat Gafarov (PSU), Matthias Meier (UBonn), and José-Luis Montiel-Olea (NYU) September 23, 2015 1 / 58 Introduction Structural VAR: Theoretical
More informationRisk neutral and risk averse approaches to multistage stochastic programming.
Risk neutral and risk averse approaches to multistage stochastic programming. A. Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA
More informationTight and Compact MILP Formulation for the Thermal Unit Commitment Problem
Online Companion for Tight and Compact MILP Formulation for the Thermal Unit Commitment Problem Germán Morales-España, Jesus M. Latorre, and Andres Ramos Universidad Pontificia Comillas, Spain Institute
More informationOptimization Tools in an Uncertain Environment
Optimization Tools in an Uncertain Environment Michael C. Ferris University of Wisconsin, Madison Uncertainty Workshop, Chicago: July 21, 2008 Michael Ferris (University of Wisconsin) Stochastic optimization
More informationRobust Process Control by Dynamic Stochastic Programming
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany MARC C. STEINBACH Robust Process Control by Dynamic Stochastic Programming ZIB-Report 04-20 (May 2004) ROBUST
More informationEmpirical likelihood and self-weighting approach for hypothesis testing of infinite variance processes and its applications
Empirical likelihood and self-weighting approach for hypothesis testing of infinite variance processes and its applications Fumiya Akashi Research Associate Department of Applied Mathematics Waseda University
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationAn Effective Heuristic for Multistage Stochastic Linear Programming
An Effective Heuristic for Multistage Stochastic Linear Programming C. Beltran-Royo L. F. Escudero J. F. Monge R. E. Rodriguez-Ravines May 28, 2013 Abstract The Multistage Stochastic Linear Programming
More informationOptimization of Dispersed Energy Supply - Stochastic Programming with Recombining Scenario Trees
1 Optimization of Dispersed Energy Supply - Stochastic Programming with Recombining Scenario Trees Alexa Epe 1, Christian Küchler 2, Werner Römisch 2, Stefan Vigerske 2, Hermann-Josef Wagner 1, Christoph
More informationOPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS
APPLICATIONES MATHEMATICAE 29,4 (22), pp. 387 398 Mariusz Michta (Zielona Góra) OPTIMAL SOLUTIONS TO STOCHASTIC DIFFERENTIAL INCLUSIONS Abstract. A martingale problem approach is used first to analyze
More informationFROM EMPIRICAL OBSERVATIONS TO TREE MODELS FOR STOCHASTIC OPTIMIZATION: CONVERGENCE PROPERTIES
SIAM J. OPTIM. Vol. 6, No. 3, pp. 1715 1740 c 016 Society for Industrial and Applied Mathematics FROM EMPIRICAL OBSERVATIONS TO TREE MODELS FOR STOCHASTIC OPTIMIZATION: CONVERGENCE PROPERTIES GEORG CH.
More informationComplexity of stochastic branch and bound methods for belief tree search in Bayesian reinforcement learning
Complexity of stochastic branch and bound methods for belief tree search in Bayesian reinforcement learning Christos Dimitrakakis Informatics Institute, University of Amsterdam, Amsterdam, The Netherlands
More informationStochastic Programming Models in Design OUTLINE
Stochastic Programming Models in Design John R. Birge University of Michigan OUTLINE Models General - Farming Structural design Design portfolio General Approximations Solutions Revisions Decision: European
More information3. Problem definition
3. Problem definition In this section, we first define a multi-dimension transaction database MD, to differentiate it between the transaction databases D in traditional association rules. Then, we discuss
More informationHow are adding integers and subtracting integers related? Work with a partner. Use integer counters to find 4 2. Remove 2 positive counters.
. How are adding integers and subtracting integers related? ACTIVITY: Work with a partner. Use integer counters to find 4. Start with 4 positive counters. Remove positive counters. What is the total number
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationForecasting the Canadian Dollar Exchange Rate Wissam Saleh & Pablo Navarro
Forecasting the Canadian Dollar Exchange Rate Wissam Saleh & Pablo Navarro Research Question: What variables effect the Canadian/US exchange rate? Do energy prices have an effect on the Canadian/US exchange
More informationRisk-Averse Dynamic Optimization. Andrzej Ruszczyński. Research supported by the NSF award CMMI
Research supported by the NSF award CMMI-0965689 Outline 1 Risk-Averse Preferences 2 Coherent Risk Measures 3 Dynamic Risk Measurement 4 Markov Risk Measures 5 Risk-Averse Control Problems 6 Solution Methods
More informationA Concise Course on Stochastic Partial Differential Equations
A Concise Course on Stochastic Partial Differential Equations Michael Röckner Reference: C. Prevot, M. Röckner: Springer LN in Math. 1905, Berlin (2007) And see the references therein for the original
More informationWHEN IS IT EVER GOING TO RAIN? Table of Average Annual Rainfall and Rainfall For Selected Arizona Cities
WHEN IS IT EVER GOING TO RAIN? Table of Average Annual Rainfall and 2001-2002 Rainfall For Selected Arizona Cities Phoenix Tucson Flagstaff Avg. 2001-2002 Avg. 2001-2002 Avg. 2001-2002 October 0.7 0.0
More informationCourse Handouts: Pages 1-20 ASSET PRICE BUBBLES AND SPECULATION. Jan Werner
Course Handouts: Pages 1-20 ASSET PRICE BUBBLES AND SPECULATION Jan Werner European University Institute May 2010 1 I. Price Bubbles: An Example Example I.1 Time is infinite; so dates are t = 0,1,2,...,.
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 011 MODULE 3 : Stochastic processes and time series Time allowed: Three Hours Candidates should answer FIVE questions. All questions carry
More informationStochastic Unit Commitment with Topology Control Recourse for Renewables Integration
1 Stochastic Unit Commitment with Topology Control Recourse for Renewables Integration Jiaying Shi and Shmuel Oren University of California, Berkeley IPAM, January 2016 33% RPS - Cumulative expected VERs
More informationStochastic Dual Dynamic Integer Programming
Stochastic Dual Dynamic Integer Programming Jikai Zou Shabbir Ahmed Xu Andy Sun December 26, 2017 Abstract Multistage stochastic integer programming (MSIP) combines the difficulty of uncertainty, dynamics,
More informationMAS3706 Topology. Revision Lectures, May I do not answer enquiries as to what material will be in the exam.
MAS3706 Topology Revision Lectures, May 208 Z.A.Lykova It is essential that you read and try to understand the lecture notes from the beginning to the end. Many questions from the exam paper will be similar
More informationProbability metrics and the stability of stochastic programs with recourse
Probability metrics and the stability of stochastic programs with recourse Michal Houda Charles University, Faculty of Mathematics and Physics, Prague, Czech Republic. Abstract In stochastic optimization
More informationGTR # VLTs GTR/VLT/Day %Δ:
MARYLAND CASINOS: MONTHLY REVENUES TOTAL REVENUE, GROSS TERMINAL REVENUE, WIN/UNIT/DAY, TABLE DATA, AND MARKET SHARE CENTER FOR GAMING RESEARCH, DECEMBER 2017 Executive Summary Since its 2010 casino debut,
More informationMr. XYZ. Stock Market Trading and Investment Astrology Report. Report Duration: 12 months. Type: Both Stocks and Option. Date: Apr 12, 2011
Mr. XYZ Stock Market Trading and Investment Astrology Report Report Duration: 12 months Type: Both Stocks and Option Date: Apr 12, 2011 KT Astrologer Website: http://www.softwareandfinance.com/magazine/astrology/kt_astrologer.php
More informationRisk averse stochastic programming: time consistency and optimal stopping
Risk averse stochastic programming: time consistency and optimal stopping A. Pichler Fakultät für Mathematik Technische Universität Chemnitz D 09111 Chemnitz, Germany A. Shapiro School of Industrial and
More informationStochastic Programming: Basic Theory
Stochastic Programming: Basic Theory John R. Birge John.Birge@ChicagoBooth.edu www.chicagobooth.edu/fac/john.birge JRBirge ICSP, Bergamo, July 2013 1 Goals Provide background on the basics of stochastic
More informationScenario grouping and decomposition algorithms for chance-constrained programs
Scenario grouping and decomposition algorithms for chance-constrained programs Yan Deng Shabbir Ahmed Jon Lee Siqian Shen Abstract A lower bound for a finite-scenario chance-constrained problem is given
More informationRevisiting some results on the complexity of multistage stochastic programs and some extensions
Revisiting some results on the complexity of multistage stochastic programs and some extensions M.M.C.R. Reaiche IMPA, Rio de Janeiro, RJ, Brazil October 30, 2015 Abstract In this work we present explicit
More informationTotally Unimodular Stochastic Programs
Totally Unimodular Stochastic Programs Nan Kong Weldon School of Biomedical Engineering, Purdue University 206 S. Martin Jischke Dr., West Lafayette, IN 47907 nkong@purdue.edu Andrew J. Schaefer Department
More informationStochastic Mathematical Programs with Equilibrium Constraints, Modeling and Sample Average Approximation
Stochastic Mathematical Programs with Equilibrium Constraints, Modeling and Sample Average Approximation Alexander Shapiro and Huifu Xu Revised June 24, 2005 Abstract. In this paper, we discuss the sample
More information