Stochastic Optimization with cvxpy EE364b Project Final Report
|
|
- Melvyn Stone
- 5 years ago
- Views:
Transcription
1 Stochastic Optimization with cvpy EE364b Project Final Report Alnur Ali June 5, Introduction A stochastic program is a conve optimization problem that includes random variables, and has the form minimize E f 0 (, ω) (1) subject to E f i (, ω) 0, i = 1,..., m, where R n is the optimization variable, ω R p is a random variable, and the f i : R n R p R, i = 0,..., m are conve functions of for each value of ω. In recent work, we implemented into cvpy a (basic) framework for specifying and solving stochastic programs; we refer to this framework as cvstoc. cvstoc works by replacing the epectations in (1), which we assume are difficult to compute in general, with their sample average approimations (SAAs), i.e., N E f i (, ω) (1/N) f i (, ω j ), (2) where the ω j are (Monte Carlo) samples of ω (N is fied by the user). In this project, we investigate the quality of the SAAs induced by (2), for both unconstrained and constrained problems (by adding new diagnostic functionality to cvstoc): we implement (into cvstoc 1 ) simple estimators of the optimality gap between a candidate solution (e.g., of a SAA) and the true solution of (1) these estimators, in turn, rely on estimators (which we also implement) of (a) the true objective value at a candidate and (b) the true optimal value we evaluate all our functionality on (unconstrained, constrained, and chance constrained) stochastic programs we document cvstoc 2. j=1 1 See for code. 2 See ali supp.pdf enclosed with the final project submission. 1
2 2 Lower bounds for unconstrained problems First, we consider the unconstrained problem minimize E f 0 (, ω), (3) the optimal value of which we denote as p. Computing a lower bound on p can be useful, e.g., when checking if the objective value f 0 () of a candidate solution (e.g., of a SAA) is near p. 2.1 Mean field approimation One way to compute a lower bound on p is as follows. Assume f 0 is conve in ω for each value of ; then, by Jensen s inequality, we have [Boy15, slide 5] E f 0 (, ω) f 0 (, E ω) = inf E f 0(, ω) inf f 0(, E ω) = p ˆp mf := inf f 0(, E ω). ˆp mf is referred to as the mean field approimation. In words, replace ω in (1) with its epectation, and then minimize w.r.t. ; the optimal value is a lower bound on p. 2.2 Averaging sample average approimations We can compute another lower bound on p as follows. Let minimize (1/N) N f 0(, ω i ), (4) be a SAA for (3), where the ω i are samples, and ˆp N denotes the optimal value of (4); then, by the discussion below (5.22) in [SDR14], we have that p E ˆp N. Unfortunately, the epectation on the r.h.s. here can be hard to calculate, so we obtain an approimate lower bound by averaging M times (M is fied by the user) instead. In symbols, let ˆp (1) N,..., ˆp (M) N be the optimal values of M (different) SAAs; then, we compute ˆp M,N := (1/M) M ˆp (i) N p. 2.3 A simple estimator of the optimality gap We estimate the optimality gap := f 0 () p at a candidate solution by where ˆ M,N := ˆf 0,N () ˆp M,N, ˆf 0,N () := (1/N) N f 0 (, ω i ) (i.e., average f 0 N times). ( ˆ M,N and ˆp M,N are biased estimators of and p, respectively [SDR14, chap. 5].) 2
3 2.4 Results In this section, we evaluate the estimators ˆp mf, ˆp M,N, and ˆ M,N on a stochastic least squares problem minimize E A b 2 2, (5) where the entries of A R m n are N (1, 1), and the entries of b R m are N ( 1, 1); in this case, (5) has the analytical solution T E A T A 2 E b T A + E b T b, assuming we know the second moments of (A, b). Fig. 1 presents the results 3. In the upper left, we plot the lower bounds ˆp mf (green) and ˆp M,N (blue) on the true optimal value of (5), p (red), vs. the number of Monte Carlo samples N (the number of averaged SAAs M is fied to 10 here): ˆp mf is quite loose, but ˆp M,N gets near p after M 10, N 100. In the upper right, we plot the estimator ˆf 0,N () (blue) of the true objective value E f 0 () (red) at a (arbitrarily chosen) candidate solution (of a SAA) vs. the number of Monte Carlo samples N: ˆf0,N () converges rather rapidly to f 0 () (as epected). In the lower left, we plot the estimator ˆ M,N (blue) of the true optimality gap (red) vs. the number of averaged SAAs M (the number of Monte Carlo samples N is fied to 100 here): ˆ M,N gets near after M 10, N 100 (consistent with the upper left plot). In the lower right, we plot the error of the gap estimator ˆ M,N vs. the number of Monte Carlo samples N and the number of averaged SAAs M: we see that the number of samples affects the error more than the number of SAAs (from the tilt of the plane), but costs more (in terms of runtime). Confidence intervals In Fig. 1, we plot 95% (t-based) confidence intervals 4 in light blue: in the upper left plot, the intervals do not contain the population parameter (i.e., p ) until N 100 samples, although in the upper right plot, the intervals are generally reliable. 3 Lower bounds for constrained problems Net, we consider the constrained problem (1); by weak duality, we have that, for any λ 0, { } { m M m } p inf E f 0 (, ω) + λ i E f i (, ω) (1/M) inf ˆf (j) 0,N () + (j) λ i ˆf i,n () := ˆp M,N. In words, form the Lagrangian for (1), replace epectations with their SAAs, fi λ 0, and then minimize w.r.t. ; averaging the infimums yields a (approimate) lower bound on p. 5 We can then use this ˆp M,N just as we did in the unconstrained setting. 3 We also recreated all the plots in Fig. 1 by varying the number of Monte Carlo samples (while holding the number of SAAs fied): the results are similar, and therefore omitted due to space constraints. 4 We also implemented (percentile) bootstrap confidence intervals [Was11, page 111] for the midway report, but omit those here to keep the presentation uncluttered. 5 It is also possible to upper bound p using a saddle point interpretation of the Lagrangian [SDR14, chap. 5], but these bounds are often trivial (i.e., + ), so we do not pursue them here. 3 j=1
4 Figure 1: Upper left: various lower bounds on the true optimal value (p, red) vs. the number of Monte Carlo samples. Upper right: the estimator (blue) of the true objective value (red) vs. the number of Monte Carlo samples. Lower left: the estimator (blue) of the true optimality gap (red) vs. the number of averaged SAAs. Lower right: the error of the gap estimator (see tet for definition) vs. the number of Monte Carlo samples and the number of averaged SAAs. 95% (t-based) confidence intervals in light blue. 3.1 Quadratically constrained quadratic programs In this section, we evaluate our estimators on a stochastic quadratically constrained quadratic program minimize E T P 0 + q0 T + r 0 subject to E T P 1 + q1 T + r 1 0, where P i := A T i A i with the entries of A i N (0, 1), and q i and r i are fied, for i = 1, 2. The top row of Fig. 2 presents the results. In the upper left, we plot ˆp M,N (blue) vs. the number of Monte Carlo samples: ˆp M,N gets near and its confidence interval contains 4
5 p after N 100 samples. The upper center and right plots are just as in Fig. 1, and have similar interpretations, e.g., our estimators get near their population parameters after M 100, N 100. Additional eperiments (omitted due to space constraints) indicate that varying λ does not matter much. Figure 2: Upper left: ˆp M,N (blue) vs. the number of Monte Carlo samples. Upper center: ˆ M,N (blue) vs. the number of Monte Carlo samples. Upper right: the error of the gap estimator vs. the number of Monte Carlo samples and the number of averaged SAAs. Lower left (see Sec. 3.2): the wealth obtained (higher is better) by solving a portfolio optimization problem with various analytic and approimate risk constraints vs. the number of Monte Carlo samples. Lower right: the value of the approimate risk constraint vs. the number of Monte Carlo samples. 95% (t-based) confidence intervals in light blue. 3.2 Chance constraints A chance constraint is a constraint of the form Prob (f(, ω) 0) η, where f is conve in for each value of ω, and η is typically a large probability (e.g., 0.95). Chance constraints are typically nonconve; thus, the conve approimation E (f(, ω)/α + 1) +, (6) 5
6 where α R is a new variable, is often employed instead (e.g., by cvstoc). It is difficult in general to construct estimators of the optimality gap for chance constrained problems [SDR14, chap. 5]; thus, in this section, we directly investigate the quality of a SAA to the conve approimation (6) by studying the value-at-risk (VaR) constrained portfolio optimization problem minimize E p T subject to Prob(p T 0) 1 η, 1 T = 1, 0, (7) where the vector of returns p N ( p, Σ). In this case, the VaR constraint in (7) has the analytical form p T Φ 1 (η) Σ 1/2 2, (8) where Φ 1 ( ) is the inverse standard normal cdf, while its conve approimation induced by (6), which is equivalent to a conditional value-at-risk (CVaR) constraint, can be epressed as p T ep ( (Φ 1 (η)) 2 /2 ) ( ) / 2π(1 η) Σ 1/2 2. (9) The bottom row of Fig. 2 presents the results. In the lower left, we plot the optimal value of (7) with the analytic VaR constraint (i.e., (8), denoted as p (VaR), and shown in green), (7) with the analytic CVaR constraint (i.e., (9), denoted p (CVaR), in red), and (7) with a SAA to (6) (denoted ˆp N, in blue) vs. the number of Monte Carlo samples: ˆp N gets near p (CVaR) after N 100, and is (indeed) more conservative than p (VaR). In the lower right, we plot the approimate probability Prob( T p 0) of the constraint failing: the probability approaches our specified tolerance of 1 η = 0.05 after N API overview Lastly, to situate our code, we describe the (updated) cvpy package hierarchy below. Essentially, a user may now create either an UnconstrDiagnostics or ConstrDiagnostics object, pass in a cvpy Problem, N, and M, and receive back ˆp mf, ˆp M,N, and/or ˆ M,N : cvpy atoms epressions constraints. stochastic (i.e., cvstoc) RandomVariable RandomVariableFactory epectation() prob partial optimize() {Unconstr, Constr}Diagnostics est optval est cost est optgap 6
7 References [Boy15] Stephen Boyd. Chance constrained optimization. ee364a/lectures/stoch_prog.pdf, January [SDR14] Aleander Shapiro, Darinka Dentcheva, and Andrzej Ruszczyński. Lectures on stochastic programming: modeling and theory, volume 16. SIAM, [Was11] Larry Wasserman. All of statistics. Springer Science & Business Media,
Decomposition Methods
Decomposition Methods separable problems, complicating variables primal decomposition dual decomposition complicating constraints general decomposition structures Prof. S. Boyd, EE364b, Stanford University
More informationOn Complexity of Multistage Stochastic Programs
On Complexity of Multistage Stochastic Programs Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0205, USA e-mail: ashapiro@isye.gatech.edu
More informationSTRESS-STRENGTH RELIABILITY ESTIMATION
CHAPTER 5 STRESS-STRENGTH RELIABILITY ESTIMATION 5. Introduction There are appliances (every physical component possess an inherent strength) which survive due to their strength. These appliances receive
More informationSolving real-life portfolio problem using stochastic programming and Monte-Carlo techniques
Solving real-life portfolio problem using stochastic programming and Monte-Carlo techniques 1 Introduction Martin Branda 1 Abstract. We deal with real-life portfolio problem with Value at Risk, transaction
More informationClass 13. Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science. Marquette University MATH 1700
Class 13 Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science Copyright 017 by D.B. Rowe 1 Agenda: Recap Chapter 6.3 6.5 Lecture Chapter 7.1 7. Review Chapter 5 for Eam 3.
More informationPortfolio Management and Optimal Execution via Convex Optimization
Portfolio Management and Optimal Execution via Convex Optimization Enzo Busseti Stanford University April 9th, 2018 Problems portfolio management choose trades with optimization minimize risk, maximize
More informationinduced by the Solvency II project
Asset Les normes allocation IFRS : new en constraints assurance induced by the Solvency II project 36 th International ASTIN Colloquium Zürich September 005 Frédéric PLANCHET Pierre THÉROND ISFA Université
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationECE 295: Lecture 03 Estimation and Confidence Interval
ECE 295: Lecture 03 Estimation and Confidence Interval Spring 2018 Prof Stanley Chan School of Electrical and Computer Engineering Purdue University 1 / 23 Theme of this Lecture What is Estimation? You
More informationChapter 2 Uncertainty Analysis and Sampling Techniques
Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying
More informationMean Variance Portfolio Theory
Chapter 1 Mean Variance Portfolio Theory This book is about portfolio construction and risk analysis in the real-world context where optimization is done with constraints and penalties specified by the
More informationClass 16. Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science. Marquette University MATH 1700
Class 16 Daniel B. Rowe, Ph.D. Department of Mathematics, Statistics, and Computer Science Copyright 013 by D.B. Rowe 1 Agenda: Recap Chapter 7. - 7.3 Lecture Chapter 8.1-8. Review Chapter 6. Problem Solving
More informationApproximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications
Approximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications Anna Timonina University of Vienna, Abraham Wald PhD Program in Statistics and Operations
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x
More informationA Monte Carlo Based Analysis of Optimal Design Criteria
A Monte Carlo Based Analysis of Optimal Design Criteria H. T. Banks, Kathleen J. Holm and Franz Kappel Center for Quantitative Sciences in Biomedicine Center for Research in Scientific Computation North
More informationYao s Minimax Principle
Complexity of algorithms The complexity of an algorithm is usually measured with respect to the size of the input, where size may for example refer to the length of a binary word describing the input,
More informationEE266 Homework 5 Solutions
EE, Spring 15-1 Professor S. Lall EE Homework 5 Solutions 1. A refined inventory model. In this problem we consider an inventory model that is more refined than the one you ve seen in the lectures. The
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Risk Measures Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Reference: Chapter 8
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}
More informationQuantitative Risk Management
Quantitative Risk Management Asset Allocation and Risk Management Martin B. Haugh Department of Industrial Engineering and Operations Research Columbia University Outline Review of Mean-Variance Analysis
More informationTesting Out-of-Sample Portfolio Performance
Testing Out-of-Sample Portfolio Performance Ekaterina Kazak 1 Winfried Pohlmeier 2 1 University of Konstanz, GSDS 2 University of Konstanz, CoFE, RCEA Econometric Research in Finance Workshop 2017 SGH
More informationSlides for Risk Management
Slides for Risk Management Introduction to the modeling of assets Groll Seminar für Finanzökonometrie Prof. Mittnik, PhD Groll (Seminar für Finanzökonometrie) Slides for Risk Management Prof. Mittnik,
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationMATH 5510 Mathematical Models of Financial Derivatives. Topic 1 Risk neutral pricing principles under single-period securities models
MATH 5510 Mathematical Models of Financial Derivatives Topic 1 Risk neutral pricing principles under single-period securities models 1.1 Law of one price and Arrow securities 1.2 No-arbitrage theory and
More informationMultilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More information3.2 No-arbitrage theory and risk neutral probability measure
Mathematical Models in Economics and Finance Topic 3 Fundamental theorem of asset pricing 3.1 Law of one price and Arrow securities 3.2 No-arbitrage theory and risk neutral probability measure 3.3 Valuation
More informationMTH6154 Financial Mathematics I Stochastic Interest Rates
MTH6154 Financial Mathematics I Stochastic Interest Rates Contents 4 Stochastic Interest Rates 45 4.1 Fixed Interest Rate Model............................ 45 4.2 Varying Interest Rate Model...........................
More informationAssessing Policy Quality in Multi-stage Stochastic Programming
Assessing Policy Quality in Multi-stage Stochastic Programming Anukal Chiralaksanakul and David P. Morton Graduate Program in Operations Research The University of Texas at Austin Austin, TX 78712 January
More informationOptimal Static Hedging of Currency Risk Using FX Forwards
GJMS Special Issue for Recent Advances in Mathematical Sciences and Applications-13 GJMS Vol 2. 2 11 18 Optimal Static Hedging of Currency Risk Using FX Forwards Anil Bhatia 1, Sanjay Bhat 2, and Vijaysekhar
More informationModelling the Sharpe ratio for investment strategies
Modelling the Sharpe ratio for investment strategies Group 6 Sako Arts 0776148 Rik Coenders 0777004 Stefan Luijten 0783116 Ivo van Heck 0775551 Rik Hagelaars 0789883 Stephan van Driel 0858182 Ellen Cardinaels
More informationFinal exam solutions
EE365 Stochastic Control / MS&E251 Stochastic Decision Models Profs. S. Lall, S. Boyd June 5 6 or June 6 7, 2013 Final exam solutions This is a 24 hour take-home final. Please turn it in to one of the
More information1 Consumption and saving under uncertainty
1 Consumption and saving under uncertainty 1.1 Modelling uncertainty As in the deterministic case, we keep assuming that agents live for two periods. The novelty here is that their earnings in the second
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationHedging Commodity Processes: Problems at the Intersection of Control, Operations, and Finance
Hedging Commodity Processes: Problems at the Intersection of Control, Operations, and Finance Jeffrey Kantor, Fanhui Fan, and Fernando Garcia Department of Chemical and Biomolecular Engineering University
More informationA Monte Carlo Based Analysis of Optimal Design Criteria
A Monte Carlo Based Analysis of Optimal Design Criteria H. T. Banks, Kathleen J. Holm and Franz Kappel Center for Quantitative Sciences in Biomedicine Center for Research in Scientific Computation North
More informationRisk minimization and portfolio diversification
Risk minimization and portfolio diversification Farzad Pourbabaee Minsuk Kwak raian A. Pirvu December 16, 2014 arxiv:1411.6657v2 [q-fin.pm] 15 Dec 2014 Abstract We consider the problem of minimizing capital
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Further Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Outline
More informationPricing Volatility Derivatives with General Risk Functions. Alejandro Balbás University Carlos III of Madrid
Pricing Volatility Derivatives with General Risk Functions Alejandro Balbás University Carlos III of Madrid alejandro.balbas@uc3m.es Content Introduction. Describing volatility derivatives. Pricing and
More informationUQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.
UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.
More informationCopyright 2011 Pearson Education, Inc. Publishing as Addison-Wesley.
Appendix: Statistics in Action Part I Financial Time Series 1. These data show the effects of stock splits. If you investigate further, you ll find that most of these splits (such as in May 1970) are 3-for-1
More informationCash flow matching with risks controlled by buffered probability of exceedance and conditional value-at-risk
DOI 10.1007/s10479-016-2354-6 ADVANCES OF OR IN COMMODITIES AND FINANCIAL MODELLING Cash flow matching with risks controlled by buffered probability of exceedance and conditional value-at-risk Danjue Shang
More informationLast Time. Martingale inequalities Martingale convergence theorem Uniformly integrable martingales. Today s lecture: Sections 4.4.1, 5.
MATH136/STAT219 Lecture 21, November 12, 2008 p. 1/11 Last Time Martingale inequalities Martingale convergence theorem Uniformly integrable martingales Today s lecture: Sections 4.4.1, 5.3 MATH136/STAT219
More informationCSCI 1951-G Optimization Methods in Finance Part 00: Course Logistics Introduction to Finance Optimization Problems
CSCI 1951-G Optimization Methods in Finance Part 00: Course Logistics Introduction to Finance Optimization Problems January 26, 2018 1 / 24 Basic information All information is available in the syllabus
More informationStochastic Approximation Algorithms and Applications
Harold J. Kushner G. George Yin Stochastic Approximation Algorithms and Applications With 24 Figures Springer Contents Preface and Introduction xiii 1 Introduction: Applications and Issues 1 1.0 Outline
More informationA class of coherent risk measures based on one-sided moments
A class of coherent risk measures based on one-sided moments T. Fischer Darmstadt University of Technology November 11, 2003 Abstract This brief paper explains how to obtain upper boundaries of shortfall
More informationThe Optimization Process: An example of portfolio optimization
ISyE 6669: Deterministic Optimization The Optimization Process: An example of portfolio optimization Shabbir Ahmed Fall 2002 1 Introduction Optimization can be roughly defined as a quantitative approach
More informationA Hybrid Importance Sampling Algorithm for VaR
A Hybrid Importance Sampling Algorithm for VaR No Author Given No Institute Given Abstract. Value at Risk (VaR) provides a number that measures the risk of a financial portfolio under significant loss.
More informationECON 214 Elements of Statistics for Economists 2016/2017
ECON 214 Elements of Statistics for Economists 2016/2017 Topic The Normal Distribution Lecturer: Dr. Bernardin Senadza, Dept. of Economics bsenadza@ug.edu.gh College of Education School of Continuing and
More informationDetermining Sample Size. Slide 1 ˆ ˆ. p q n E = z α / 2. (solve for n by algebra) n = E 2
Determining Sample Size Slide 1 E = z α / 2 ˆ ˆ p q n (solve for n by algebra) n = ( zα α / 2) 2 p ˆ qˆ E 2 Sample Size for Estimating Proportion p When an estimate of ˆp is known: Slide 2 n = ˆ ˆ ( )
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationMartingales. by D. Cox December 2, 2009
Martingales by D. Cox December 2, 2009 1 Stochastic Processes. Definition 1.1 Let T be an arbitrary index set. A stochastic process indexed by T is a family of random variables (X t : t T) defined on a
More informationStatistical Methods in Financial Risk Management
Statistical Methods in Financial Risk Management Lecture 1: Mapping Risks to Risk Factors Alexander J. McNeil Maxwell Institute of Mathematical Sciences Heriot-Watt University Edinburgh 2nd Workshop on
More informationOptimal retention for a stop-loss reinsurance with incomplete information
Optimal retention for a stop-loss reinsurance with incomplete information Xiang Hu 1 Hailiang Yang 2 Lianzeng Zhang 3 1,3 Department of Risk Management and Insurance, Nankai University Weijin Road, Tianjin,
More informationOptimal Dam Management
Optimal Dam Management Michel De Lara et Vincent Leclère July 3, 2012 Contents 1 Problem statement 1 1.1 Dam dynamics.................................. 2 1.2 Intertemporal payoff criterion..........................
More informationStrategy Lines and Optimal Mixed Strategy for R
Strategy Lines and Optimal Mixed Strategy for R Best counterstrategy for C for given mixed strategy by R In the previous lecture we saw that if R plays a particular mixed strategy, [p, p, and shows no
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 2 Random number generation January 18, 2018
More informationLog-Robust Portfolio Management
Log-Robust Portfolio Management Dr. Aurélie Thiele Lehigh University Joint work with Elcin Cetinkaya and Ban Kawas Research partially supported by the National Science Foundation Grant CMMI-0757983 Dr.
More informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationWeb Appendix. Are the effects of monetary policy shocks big or small? Olivier Coibion
Web Appendix Are the effects of monetary policy shocks big or small? Olivier Coibion Appendix 1: Description of the Model-Averaging Procedure This section describes the model-averaging procedure used in
More informationA Robust Option Pricing Problem
IMA 2003 Workshop, March 12-19, 2003 A Robust Option Pricing Problem Laurent El Ghaoui Department of EECS, UC Berkeley 3 Robust optimization standard form: min x sup u U f 0 (x, u) : u U, f i (x, u) 0,
More informationAppendix to: Long-Run Asset Pricing Implications of Housing Collateral Constraints
Appendix to: Long-Run Asset Pricing Implications of Housing Collateral Constraints Hanno Lustig UCLA and NBER Stijn Van Nieuwerburgh June 27, 2006 Additional Figures and Tables Calibration of Expenditure
More informationEARLY EXERCISE OPTIONS: UPPER BOUNDS
EARLY EXERCISE OPTIONS: UPPER BOUNDS LEIF B.G. ANDERSEN AND MARK BROADIE Abstract. In this article, we discuss how to generate upper bounds for American or Bermudan securities by Monte Carlo methods. These
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationOptimal Dividend Policy of A Large Insurance Company with Solvency Constraints. Zongxia Liang
Optimal Dividend Policy of A Large Insurance Company with Solvency Constraints Zongxia Liang Department of Mathematical Sciences Tsinghua University, Beijing 100084, China zliang@math.tsinghua.edu.cn Joint
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMS091) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 3 Importance sampling January 27, 2015 M. Wiktorsson
More informationDynamic Asset and Liability Management Models for Pension Systems
Dynamic Asset and Liability Management Models for Pension Systems The Comparison between Multi-period Stochastic Programming Model and Stochastic Control Model Muneki Kawaguchi and Norio Hibiki June 1,
More informationDynamic Portfolio Execution Detailed Proofs
Dynamic Portfolio Execution Detailed Proofs Gerry Tsoukalas, Jiang Wang, Kay Giesecke March 16, 2014 1 Proofs Lemma 1 (Temporary Price Impact) A buy order of size x being executed against i s ask-side
More informationVaR vs CVaR in Risk Management and Optimization
VaR vs CVaR in Risk Management and Optimization Stan Uryasev Joint presentation with Sergey Sarykalin, Gaia Serraino and Konstantin Kalinchenko Risk Management and Financial Engineering Lab, University
More informationMVE051/MSG Lecture 7
MVE051/MSG810 2017 Lecture 7 Petter Mostad Chalmers November 20, 2017 The purpose of collecting and analyzing data Purpose: To build and select models for parts of the real world (which can be used for
More informationIntroduction to Computational Finance and Financial Econometrics Introduction to Portfolio Theory
You can t see this text! Introduction to Computational Finance and Financial Econometrics Introduction to Portfolio Theory Eric Zivot Spring 2015 Eric Zivot (Copyright 2015) Introduction to Portfolio Theory
More informationProgressive Hedging for Multi-stage Stochastic Optimization Problems
Progressive Hedging for Multi-stage Stochastic Optimization Problems David L. Woodruff Jean-Paul Watson Graduate School of Management University of California, Davis Davis, CA 95616, USA dlwoodruff@ucdavis.edu
More informationCIVL Learning Objectives. Definitions. Discrete Distributions
CIVL 3103 Discrete Distributions Learning Objectives Define discrete distributions, and identify common distributions applicable to engineering problems. Identify the appropriate distribution (i.e. binomial,
More informationOvernight Index Rate: Model, calibration and simulation
Research Article Overnight Index Rate: Model, calibration and simulation Olga Yashkir and Yuri Yashkir Cogent Economics & Finance (2014), 2: 936955 Page 1 of 11 Research Article Overnight Index Rate: Model,
More informationSolving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function?
DOI 0.007/s064-006-9073-z ORIGINAL PAPER Solving dynamic portfolio choice problems by recursing on optimized portfolio weights or on the value function? Jules H. van Binsbergen Michael W. Brandt Received:
More informationBivariate Birnbaum-Saunders Distribution
Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators
More informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationStatistical Tables Compiled by Alan J. Terry
Statistical Tables Compiled by Alan J. Terry School of Science and Sport University of the West of Scotland Paisley, Scotland Contents Table 1: Cumulative binomial probabilities Page 1 Table 2: Cumulative
More informationPortfolio Optimization Using Conditional Value-At-Risk and Conditional Drawdown-At-Risk
Portfolio Optimization Using Conditional Value-At-Risk and Conditional Drawdown-At-Risk Enn Kuutan A thesis submitted in partial fulfillment of the degree of BACHELOR OF APPLIED SCIENCE Supervisor: Dr.
More informationECO220Y Estimation: Confidence Interval Estimator for Sample Proportions Readings: Chapter 11 (skip 11.5)
ECO220Y Estimation: Confidence Interval Estimator for Sample Proportions Readings: Chapter 11 (skip 11.5) Fall 2011 Lecture 10 (Fall 2011) Estimation Lecture 10 1 / 23 Review: Sampling Distributions Sample
More informationLecture 22. Survey Sampling: an Overview
Math 408 - Mathematical Statistics Lecture 22. Survey Sampling: an Overview March 25, 2013 Konstantin Zuev (USC) Math 408, Lecture 22 March 25, 2013 1 / 16 Survey Sampling: What and Why In surveys sampling
More information1 Appendix A: Definition of equilibrium
Online Appendix to Partnerships versus Corporations: Moral Hazard, Sorting and Ownership Structure Ayca Kaya and Galina Vereshchagina Appendix A formally defines an equilibrium in our model, Appendix B
More informationAnalysis of Variance and Design of Experiments-II
Analysis of Variance and Design of Experiments-II MODULE I LECTURE - 8 INCOMPLETE BLOCK DESIGNS Dr Shalabh Department of Mathematics & Statistics Indian Institute of Technology Kanpur Generally, we are
More informationRisk aversion in multi-stage stochastic programming: a modeling and algorithmic perspective
Risk aversion in multi-stage stochastic programming: a modeling and algorithmic perspective Tito Homem-de-Mello School of Business Universidad Adolfo Ibañez, Santiago, Chile Joint work with Bernardo Pagnoncelli
More informationDebt Sustainability Risk Analysis with Analytica c
1 Debt Sustainability Risk Analysis with Analytica c Eduardo Ley & Ngoc-Bich Tran We present a user-friendly toolkit for Debt-Sustainability Risk Analysis (DSRA) which provides useful indicators to identify
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More information1 Overview. 2 The Gradient Descent Algorithm. AM 221: Advanced Optimization Spring 2016
AM 22: Advanced Optimization Spring 206 Prof. Yaron Singer Lecture 9 February 24th Overview In the previous lecture we reviewed results from multivariate calculus in preparation for our journey into convex
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More informationLarge-Scale SVM Optimization: Taking a Machine Learning Perspective
Large-Scale SVM Optimization: Taking a Machine Learning Perspective Shai Shalev-Shwartz Toyota Technological Institute at Chicago Joint work with Nati Srebro Talk at NEC Labs, Princeton, August, 2008 Shai
More informationArbitrage Theory without a Reference Probability: challenges of the model independent approach
Arbitrage Theory without a Reference Probability: challenges of the model independent approach Matteo Burzoni Marco Frittelli Marco Maggis June 30, 2015 Abstract In a model independent discrete time financial
More informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More informationChapter 7. Sampling Distributions
Chapter 7 Sampling Distributions Section 7.1 Sampling Distributions and the Central Limit Theorem Sampling Distributions Sampling distribution The probability distribution of a sample statistic. Formed
More informationModule 4: Point Estimation Statistics (OA3102)
Module 4: Point Estimation Statistics (OA3102) Professor Ron Fricker Naval Postgraduate School Monterey, California Reading assignment: WM&S chapter 8.1-8.4 Revision: 1-12 1 Goals for this Module Define
More informationUsing Monte Carlo Integration and Control Variates to Estimate π
Using Monte Carlo Integration and Control Variates to Estimate π N. Cannady, P. Faciane, D. Miksa LSU July 9, 2009 Abstract We will demonstrate the utility of Monte Carlo integration by using this algorithm
More informationEE365: Risk Averse Control
EE365: Risk Averse Control Risk averse optimization Exponential risk aversion Risk averse control 1 Outline Risk averse optimization Exponential risk aversion Risk averse control Risk averse optimization
More informationModelling of Long-Term Risk
Modelling of Long-Term Risk Roger Kaufmann Swiss Life roger.kaufmann@swisslife.ch 15th International AFIR Colloquium 6-9 September 2005, Zurich c 2005 (R. Kaufmann, Swiss Life) Contents A. Basel II B.
More informationOmitted Variables Bias in Regime-Switching Models with Slope-Constrained Estimators: Evidence from Monte Carlo Simulations
Journal of Statistical and Econometric Methods, vol. 2, no.3, 2013, 49-55 ISSN: 2051-5057 (print version), 2051-5065(online) Scienpress Ltd, 2013 Omitted Variables Bias in Regime-Switching Models with
More informationTraditional Optimization is Not Optimal for Leverage-Averse Investors
Posted SSRN 10/1/2013 Traditional Optimization is Not Optimal for Leverage-Averse Investors Bruce I. Jacobs and Kenneth N. Levy forthcoming The Journal of Portfolio Management, Winter 2014 Bruce I. Jacobs
More informationEffectiveness of CPPI Strategies under Discrete Time Trading
Effectiveness of CPPI Strategies under Discrete Time Trading S. Balder, M. Brandl 1, Antje Mahayni 2 1 Department of Banking and Finance, University of Bonn 2 Department of Accounting and Finance, Mercator
More informationWhat can we do with numerical optimization?
Optimization motivation and background Eddie Wadbro Introduction to PDE Constrained Optimization, 2016 February 15 16, 2016 Eddie Wadbro, Introduction to PDE Constrained Optimization, February 15 16, 2016
More information**BEGINNING OF EXAMINATION** A random sample of five observations from a population is:
**BEGINNING OF EXAMINATION** 1. You are given: (i) A random sample of five observations from a population is: 0.2 0.7 0.9 1.1 1.3 (ii) You use the Kolmogorov-Smirnov test for testing the null hypothesis,
More information