Quasi-Convex Stochastic Dynamic Programming
|
|
- August Lewis
- 6 years ago
- Views:
Transcription
1 Quasi-Convex Stochastic Dynamic Programming John R. Birge University of Chicago Booth School of Business JRBirge SIAM FM12, MSP, 10 July
2 General Theme Many dynamic optimization problems dealing with commodity storage management have a value function with quasi-convex (or quasi- concave) properties Taking advantage of the structure u can help overcome curse of dimensionality An approximation algorithm can repeatedly identify sub-level sets to obtain convergence JRBirge SIAM FM12, MSP, 10 July
3 Outline Motivation Basic properties Level-set approximation algorithm Computational examples Conclusions and next steps JRBirge SIAM FM12, MSP, 10 July
4 Motivation Commodity storage characterized by inventory (x) and price (p): px For a one-stage objective x 0 andp 0, V(p,x)=px is quasi-concave: JRBirge SIAM FM12, MSP, 10 July
5 Challenges Quasi-convexity is not generally additive Many DP objectives are additive but others are not, e.g., risk-sensitive objective: Log[ E(exp(θ v t (x t,p t )) ] When does the property persist in the value function across time? How to take advantage of the property when present? JRBirge SIAM FM12, MSP, 10 July
6 Simple Example: Commodity Management Objective (with action u, capacity C, discount factor δ): V t (x t,p t )= max -xt u C-x p + E[V +u,p t t u δ t+1 (x t t+1 )] Quasi-concave at end of horizon (T) Induction possible? Assume: (p t+1 /p t ) ~ F (fixed stationary) ti JRBirge SIAM FM12, MSP, 10 July
7 Quasi-Concavity Result for V t Basic induction: V t is then quasi-concave if V t (x, p) min{v t (x 1,p 1 ),V t (x 2,p 2 )}, (1) for any such combination of points, (x 1, p 1 ) and (x 2, p 2 ). Without loss of generality, suppose x 1 <x 2 and let and Suppose and x = x 1 + = x 2, (2) p = Up 1 = Dp 2. (3) V t (x 1,p 1 )= u 1 p 1 + E p1 V t+1 (x 1 + u 1,p 0 ), (4) V t (x 2,p 2 )= u 2 p 2 + E p2 V t+1 (x 2 + u 2,p 0 ). (5) JRBirge SIAM FM12, MSP, 10 July
8 Quasi-concavity of Value Function Suppose policy σ 1 starting from (x, p) at stage that sets u 1 = u 1 + at stage t and then follows policy π 1 (ω) (optimal along all p t (ω) =w t (ω)p 1 at time t; σ 2 starting from (x, p) withu 2 = u 2 at t and then π 2 (ω) -optimalfor p t (ω) = v t (ω)p 2 ; σ be the policy that randomly chooses σ 1 or σ 2 from state (x, p) V t (x, p) V σ = 1 ( (u p 0 t,σ (x, p) 2 (u 1 + )Up 1 + UE p1 V t+1 (x 1 + u 1, p ) (u 2 )Dp 2 + DE p2 V t+1 (x 2 + u 2,p 0 )) = ( Up( p Dp ) + (UV t(x 1, p 1 ) + DV t t( (x 2, p 2 )) (U + D)min(V t(x 1,p 1 ),V t (x 2,p 2 )) p = p 1 + p 2 + 2p 1p 2 min(v t (x 1,p 1 ),V t (x 2,p 2 )), 4p 1 p 2 min{v t (x 1,p 1 ),V t (x 2,p 2 )} using p = Up 1 = Dp 2, U = p 1+p 2 2p 1, D = p 1+p 2 2p 2,andthat2p 1 p 2 p p 2 2. JRBirge SIAM FM12, MSP, 10 July
9 Analysis and Extensions of Result Extends to proportional transaction costs directly Key assumption: Price process increments are independent of the price level (as in a log-normal/gbm model) Extends to other price processes as well: Distributions such that V t,π1 (x 1,Up 1 )=U 0 V t (x 1,p 1 ) and V t = t,π1 2(x 2, Dp 2 ) D V t (x 1, s 1 ) where U + D 2. JRBirge SIAM FM12, MSP, 10 July
10 Extended Processes Ornstein-Uhlenbeck: d log p t = θ(μ log p t )dt + σdw t, (1) where θ, μ, andσ are scalars; W t is a standard Brownian motion. The price p t given initial value p is then p t =exp(logpe θt + μ(1 e θt )+ Z t 0 σe θ(q t) dw q ). (2) Let p = Up 1 = Dp 2, p i t be spot price at t given spot price p i at 0, i =1, 2, then p 1 U e θt = U 0 2 D e θt = D 0 t p t p t, p t p t p t, (3) the condition, U 0 + D 0 2holdsforsufficiently small θ and sufficiently large time intervals. JRBirge SIAM FM12, MSP, 10 July
11 Algorithm Motivation Assume: Quasi-concave (-convex) structure Basic idea: approximate the level sets of the value function at each stage Use the level-set approximation in stage t+1 to construct an approximation at t. Use outer linearization i i to construct each successive approximation JRBirge SIAM FM12, MSP, 10 July
12 Algorithms Quasi-concave Dynamic Program Method(QDPM) 1. Set U t (x, p) for all t. LetU T (x, p) =V T (x, p). Set t = T REPEAT: (a) Randomly select (x t,p t ); (b) Find Ût(x U = t t, p t ) max u A(xt,p t ) v t (u, x t, p t )+ P z Z t+1 (u,x t,s t ) π t+1(z, a, x t, p t )U t+1 (z); (c) Perform Bound Update; (d) Choose t. 3. UNTIL U 1 (x 1,p 1 )isstable. Bound Update Construct upper bound on the value function using outer linearization: U t (x, p) min(u t (x, p), Ût(x t,p t )) (1) ˆ for all (x, p) such that U t (x t,p t ) T (x x t,p p t )) 0. JRBirge SIAM FM12, MSP, 10 July
13 Bound Update Û t (x t,p t ) Ût(x t,p t )(x x t,p p t )) 0 U t (x, p) Û(x t,p t ) JRBirge SIAM FM12, MSP, 10 July
14 Examples Setup: Lognormal prices (log-mean: 1, logst.dev.: 0.4) Transaction cost: 5% proportional Note: bound update must achieve sufficient accuracy to also maintain quasi-concavity Use of 1000, 500, and 250 random samples per period JRBirge SIAM FM12, MSP, 10 July
15 Approximations at T-1 JRBirge SIAM FM12, MSP, 10 July
16 Upper Bounds at T-5 JRBirge SIAM FM12, MSP, 10 July
17 Observations from Example For the T-1-stage results: Close to original for N T-1 =1000 More differences for N T-1 =500 and N T-1 =250 increasingly less consistent in shape with fewer samples. For the T-5-stage results: Similar il pattern more error less consistent level sets with N T-5 =250 JRBirge SIAM FM12, MSP, 10 July \
18 Conclusions Quasi-convexity may persist in some in some commodity inventory control problems Possible resolutions with level set approximations Limitations: Accuracy needed in the approximation to preserve quasi-concavity Unclear results for higher-dimensions JRBirge SIAM FM12, MSP, 10 July
19 Future Work Find other structures (such as risk-sensitive objective) with quasi-convex structure Identify methods to maintain the property in the approximation Extend results for multiple dimensions JRBirge SIAM FM12, MSP, 10 July
20 Thank you! JRBirge SIAM FM12, MSP, 10 July
Asset-Liability Management
Asset-Liability Management John Birge University of Chicago Booth School of Business JRBirge INFORMS San Francisco, Nov. 2014 1 Overview Portfolio optimization involves: Modeling Optimization Estimation
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 informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationThe Values of Information and Solution in Stochastic Programming
The Values of Information and Solution in Stochastic Programming John R. Birge The University of Chicago Booth School of Business JRBirge ICSP, Bergamo, July 2013 1 Themes The values of information and
More informationControl Improvement for Jump-Diffusion Processes with Applications to Finance
Control Improvement for Jump-Diffusion Processes with Applications to Finance Nicole Bäuerle joint work with Ulrich Rieder Toronto, June 2010 Outline Motivation: MDPs Controlled Jump-Diffusion Processes
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationA Hybrid Commodity and Interest Rate Market Model
A Hybrid Commodity and Interest Rate Market Model University of Technology, Sydney June 1 Literature A Hybrid Market Model Recall: The basic LIBOR Market Model The cross currency LIBOR Market Model LIBOR
More informationDeterministic Income under a Stochastic Interest Rate
Deterministic Income under a Stochastic Interest Rate Julia Eisenberg, TU Vienna Scientic Day, 1 Agenda 1 Classical Problem: Maximizing Discounted Dividends in a Brownian Risk Model 2 Maximizing Discounted
More informationLecture 8: The Black-Scholes theory
Lecture 8: The Black-Scholes theory Dr. Roman V Belavkin MSO4112 Contents 1 Geometric Brownian motion 1 2 The Black-Scholes pricing 2 3 The Black-Scholes equation 3 References 5 1 Geometric Brownian motion
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationManaging Risk with Operational and Financial Instruments
Managing Risk with Operational and Financial Instruments John R. Birge The University of Chicago Booth School of Business www.chicagobooth.edu/fac/john.birge Motivation Operations (e.g., flexible production,
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 informationBACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin
BACHELIER FINANCE SOCIETY 4 th World Congress Tokyo, 26 Investments and forward utilities Thaleia Zariphopoulou The University of Texas at Austin 1 Topics Utility-based measurement of performance Utilities
More informationLecture 3. Sergei Fedotov Introduction to Financial Mathematics. Sergei Fedotov (University of Manchester) / 6
Lecture 3 Sergei Fedotov 091 - Introduction to Financial Mathematics Sergei Fedotov (University of Manchester) 091 010 1 / 6 Lecture 3 1 Distribution for lns(t) Solution to Stochastic Differential Equation
More informationStochastic Dual Dynamic Programming Algorithm for Multistage Stochastic Programming
Stochastic Dual Dynamic Programg Algorithm for Multistage Stochastic Programg Final presentation ISyE 8813 Fall 2011 Guido Lagos Wajdi Tekaya Georgia Institute of Technology November 30, 2011 Multistage
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationStochastic Finance 2010 Summer School Ulm Lecture 1: Energy Derivatives
Stochastic Finance 2010 Summer School Ulm Lecture 1: Energy Derivatives Professor Dr. Rüdiger Kiesel 21. September 2010 1 / 62 1 Energy Markets Spot Market Futures Market 2 Typical models Schwartz Model
More informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationOrnstein-Uhlenbeck Theory
Beatrice Byukusenge Department of Technomathematics Lappeenranta University of technology January 31, 2012 Definition of a stochastic process Let (Ω,F,P) be a probability space. A stochastic process is
More informationVolatility. Roberto Renò. 2 March 2010 / Scuola Normale Superiore. Dipartimento di Economia Politica Università di Siena
Dipartimento di Economia Politica Università di Siena 2 March 2010 / Scuola Normale Superiore What is? The definition of volatility may vary wildly around the idea of the standard deviation of price movements
More informationThe Self-financing Condition: Remembering the Limit Order Book
The Self-financing Condition: Remembering the Limit Order Book R. Carmona, K. Webster Bendheim Center for Finance ORFE, Princeton University November 6, 2013 Structural relationships? From LOB Models to
More informationTHE OPTIMAL ASSET ALLOCATION PROBLEMFOR AN INVESTOR THROUGH UTILITY MAXIMIZATION
THE OPTIMAL ASSET ALLOCATION PROBLEMFOR AN INVESTOR THROUGH UTILITY MAXIMIZATION SILAS A. IHEDIOHA 1, BRIGHT O. OSU 2 1 Department of Mathematics, Plateau State University, Bokkos, P. M. B. 2012, Jos,
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 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 informationModeling via Stochastic Processes in Finance
Modeling via Stochastic Processes in Finance Dimbinirina Ramarimbahoaka Department of Mathematics and Statistics University of Calgary AMAT 621 - Fall 2012 October 15, 2012 Question: What are appropriate
More informationComprehensive Exam. August 19, 2013
Comprehensive Exam August 19, 2013 You have a total of 180 minutes to complete the exam. If a question seems ambiguous, state why, sharpen it up and answer the sharpened-up question. Good luck! 1 1 Menu
More informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationVariable Annuities with Lifelong Guaranteed Withdrawal Benefits
Variable Annuities with Lifelong Guaranteed Withdrawal Benefits presented by Yue Kuen Kwok Department of Mathematics Hong Kong University of Science and Technology Hong Kong, China * This is a joint work
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay. Solutions to Final Exam.
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2011, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (32 pts) Answer briefly the following questions. 1. Suppose
More informationConditional Full Support and No Arbitrage
Gen. Math. Notes, Vol. 32, No. 2, February 216, pp.54-64 ISSN 2219-7184; Copyright c ICSRS Publication, 216 www.i-csrs.org Available free online at http://www.geman.in Conditional Full Support and No Arbitrage
More informationValuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model
Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility
More informationProblem Set 3. Thomas Philippon. April 19, Human Wealth, Financial Wealth and Consumption
Problem Set 3 Thomas Philippon April 19, 2002 1 Human Wealth, Financial Wealth and Consumption The goal of the question is to derive the formulas on p13 of Topic 2. This is a partial equilibrium analysis
More informationLecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.
Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous
More informationHedging with Life and General Insurance Products
Hedging with Life and General Insurance Products June 2016 2 Hedging with Life and General Insurance Products Jungmin Choi Department of Mathematics East Carolina University Abstract In this study, a hybrid
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
More informationSOLVING ROBUST SUPPLY CHAIN PROBLEMS
SOLVING ROBUST SUPPLY CHAIN PROBLEMS Daniel Bienstock Nuri Sercan Özbay Columbia University, New York November 13, 2005 Project with Lucent Technologies Optimize the inventory buffer levels in a complicated
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 informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
More informationConditional Density Method in the Computation of the Delta with Application to Power Market
Conditional Density Method in the Computation of the Delta with Application to Power Market Asma Khedher Centre of Mathematics for Applications Department of Mathematics University of Oslo A joint work
More informationθ(t ) = T f(0, T ) + σ2 T
1 Derivatives Pricing and Financial Modelling Andrew Cairns: room M3.08 E-mail: A.Cairns@ma.hw.ac.uk Tutorial 10 1. (Ho-Lee) Let X(T ) = T 0 W t dt. (a) What is the distribution of X(T )? (b) Find E[exp(
More informationParameters Estimation in Stochastic Process Model
Parameters Estimation in Stochastic Process Model A Quasi-Likelihood Approach Ziliang Li University of Maryland, College Park GEE RIT, Spring 28 Outline 1 Model Review The Big Model in Mind: Signal + Noise
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 informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 2-3 Haijun Li An Introduction to Stochastic Calculus Week 2-3 1 / 24 Outline
More informationValue of Flexibility in Managing R&D Projects Revisited
Value of Flexibility in Managing R&D Projects Revisited Leonardo P. Santiago & Pirooz Vakili November 2004 Abstract In this paper we consider the question of whether an increase in uncertainty increases
More informationThe Black-Scholes PDE from Scratch
The Black-Scholes PDE from Scratch chris bemis November 27, 2006 0-0 Goal: Derive the Black-Scholes PDE To do this, we will need to: Come up with some dynamics for the stock returns Discuss Brownian motion
More informationThe University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam
The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider
More informationInvestigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs. Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2
Investigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2 1 Hacettepe University Department of Actuarial Sciences 06800, TURKEY 2 Middle
More informationAll Investors are Risk-averse Expected Utility Maximizers. Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel)
All Investors are Risk-averse Expected Utility Maximizers Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel) First Name: Waterloo, April 2013. Last Name: UW ID #:
More informationOptimization Models in Financial Engineering and Modeling Challenges
Optimization Models in Financial Engineering and Modeling Challenges John Birge University of Chicago Booth School of Business JRBirge UIUC, 25 Mar 2009 1 Introduction History of financial engineering
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
More informationFinite Memory and Imperfect Monitoring
Federal Reserve Bank of Minneapolis Research Department Finite Memory and Imperfect Monitoring Harold L. Cole and Narayana Kocherlakota Working Paper 604 September 2000 Cole: U.C.L.A. and Federal Reserve
More informationBasic Concepts in Mathematical Finance
Chapter 1 Basic Concepts in Mathematical Finance In this chapter, we give an overview of basic concepts in mathematical finance theory, and then explain those concepts in very simple cases, namely in the
More informationOverview. Transformation method Rejection method. Monte Carlo vs ordinary methods. 1 Random numbers. 2 Monte Carlo integration.
Overview 1 Random numbers Transformation method Rejection method 2 Monte Carlo integration Monte Carlo vs ordinary methods 3 Summary Transformation method Suppose X has probability distribution p X (x),
More informationSTATE UNIVERSITY OF NEW YORK AT ALBANY Department of Economics. Ph. D. Preliminary Examination: Macroeconomics Fall, 2009
STATE UNIVERSITY OF NEW YORK AT ALBANY Department of Economics Ph. D. Preliminary Examination: Macroeconomics Fall, 2009 Instructions: Read the questions carefully and make sure to show your work. You
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationStochastic Dual Dynamic Programming
1 / 43 Stochastic Dual Dynamic Programming Operations Research Anthony Papavasiliou 2 / 43 Contents [ 10.4 of BL], [Pereira, 1991] 1 Recalling the Nested L-Shaped Decomposition 2 Drawbacks of Nested Decomposition
More informationCredit Risk : Firm Value Model
Credit Risk : Firm Value Model Prof. Dr. Svetlozar Rachev Institute for Statistics and Mathematical Economics University of Karlsruhe and Karlsruhe Institute of Technology (KIT) Prof. Dr. Svetlozar Rachev
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
More informationOptimal switching problems for daily power system balancing
Optimal switching problems for daily power system balancing Dávid Zoltán Szabó University of Manchester davidzoltan.szabo@postgrad.manchester.ac.uk June 13, 2016 ávid Zoltán Szabó (University of Manchester)
More informationPakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks
Pakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks Spring 2009 Main question: How much are patents worth? Answering this question is important, because it helps
More informationQuasi-Monte Carlo for Finance
Quasi-Monte Carlo for Finance Peter Kritzer Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Linz, Austria NCTS, Taipei, November 2016 Peter Kritzer
More informationM.I.T Fall Practice Problems
M.I.T. 15.450-Fall 2010 Sloan School of Management Professor Leonid Kogan Practice Problems 1. Consider a 3-period model with t = 0, 1, 2, 3. There are a stock and a risk-free asset. The initial stock
More informationImportance sampling and Monte Carlo-based calibration for time-changed Lévy processes
Importance sampling and Monte Carlo-based calibration for time-changed Lévy processes Stefan Kassberger Thomas Liebmann BFS 2010 1 Motivation 2 Time-changed Lévy-models and Esscher transforms 3 Applications
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 informationStochastic Optimal Control
Stochastic Optimal Control Lecturer: Eilyan Bitar, Cornell ECE Scribe: Kevin Kircher, Cornell MAE These notes summarize some of the material from ECE 5555 (Stochastic Systems) at Cornell in the fall of
More informationRobustness, Model Uncertainty and Pricing
Robustness, Model Uncertainty and Pricing Antoon Pelsser 1 1 Maastricht University & Netspar Email: a.pelsser@maastrichtuniversity.nl 29 October 2010 Swissquote Conference Lausanne A. Pelsser (Maastricht
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationIs the Structural Approach More Accurate than the Statistical Approach in Bankruptcy Prediction?
Is the Structural Approach More Accurate than the Statistical Approach in Bankruptcy Prediction? Hui Hao Global Risk Management, Bank of Nova Scotia April 12, 2007 Road Map Theme: Horse racing among two
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationOptimal liquidation with market parameter shift: a forward approach
Optimal liquidation with market parameter shift: a forward approach (with S. Nadtochiy and T. Zariphopoulou) Haoran Wang Ph.D. candidate University of Texas at Austin ICERM June, 2017 Problem Setup and
More informationA Stochastic Approximation Algorithm for Making Pricing Decisions in Network Revenue Management Problems
A Stochastic Approximation Algorithm for Making ricing Decisions in Network Revenue Management roblems Sumit Kunnumkal Indian School of Business, Gachibowli, Hyderabad, 500032, India sumit kunnumkal@isb.edu
More informationConstructing Markov models for barrier options
Constructing Markov models for barrier options Gerard Brunick joint work with Steven Shreve Department of Mathematics University of Texas at Austin Nov. 14 th, 2009 3 rd Western Conference on Mathematical
More informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationOptimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models
Optimally Thresholded Realized Power Variations for Lévy Jump Diffusion Models José E. Figueroa-López 1 1 Department of Statistics Purdue University University of Missouri-Kansas City Department of Mathematics
More informationTime-changed Brownian motion and option pricing
Time-changed Brownian motion and option pricing Peter Hieber Chair of Mathematical Finance, TU Munich 6th AMaMeF Warsaw, June 13th 2013 Partially joint with Marcos Escobar (RU Toronto), Matthias Scherer
More informationAn Approximation Algorithm for Capacity Allocation over a Single Flight Leg with Fare-Locking
An Approximation Algorithm for Capacity Allocation over a Single Flight Leg with Fare-Locking Mika Sumida School of Operations Research and Information Engineering, Cornell University, Ithaca, New York
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationA Binomial Model of Asset and Option Pricing with Heterogen
A Binomial Model of Asset and Option Pricing with Heterogeneous Beliefs School of Finance and Economics, UTS 15th International Conference on Computing in Economics and Finance 15-17 July 2009 Sydney Basic
More informationEvaluating Electricity Generation, Energy Options, and Complex Networks
Evaluating Electricity Generation, Energy Options, and Complex Networks John Birge The University of Chicago Graduate School of Business and Quantstar 1 Outline Derivatives Real options and electricity
More informationShort & Long Run impact of volatility on the effect monetary shocks
Short & Long Run impact of volatility on the effect monetary shocks Fernando Alvarez University of Chicago & NBER Inflation: Drivers & Dynamics Conference 218 Cleveland Fed Alvarez Volatility & Monetary
More informationValuation of performance-dependent options in a Black- Scholes framework
Valuation of performance-dependent options in a Black- Scholes framework Thomas Gerstner, Markus Holtz Institut für Numerische Simulation, Universität Bonn, Germany Ralf Korn Fachbereich Mathematik, TU
More informationS t d with probability (1 p), where
Stochastic Calculus Week 3 Topics: Towards Black-Scholes Stochastic Processes Brownian Motion Conditional Expectations Continuous-time Martingales Towards Black Scholes Suppose again that S t+δt equals
More informationInterest rate models and Solvency II
www.nr.no Outline Desired properties of interest rate models in a Solvency II setting. A review of three well-known interest rate models A real example from a Norwegian insurance company 2 Interest rate
More informationMartingales, Part II, with Exercise Due 9/21
Econ. 487a Fall 1998 C.Sims Martingales, Part II, with Exercise Due 9/21 1. Brownian Motion A process {X t } is a Brownian Motion if and only if i. it is a martingale, ii. t is a continuous time parameter
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationBasic Concepts and Examples in Finance
Basic Concepts and Examples in Finance Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 5, 2017 ICMAT / UC3M The Financial Market The Financial Market We assume there are
More information( ) since this is the benefit of buying the asset at the strike price rather
Review of some financial models for MAT 483 Parity and Other Option Relationships The basic parity relationship for European options with the same strike price and the same time to expiration is: C( KT
More information1 The EOQ and Extensions
IEOR4000: Production Management Lecture 2 Professor Guillermo Gallego September 16, 2003 Lecture Plan 1. The EOQ and Extensions 2. Multi-Item EOQ Model 1 The EOQ and Extensions We have explored some of
More informationForward Dynamic Utility
Forward Dynamic Utility El Karoui Nicole & M RAD Mohamed UnivParis VI / École Polytechnique,CMAP elkaroui@cmapx.polytechnique.fr with the financial support of the "Fondation du Risque" and the Fédération
More informationLiquidity and Risk Management
Liquidity and Risk Management By Nicolae Gârleanu and Lasse Heje Pedersen Risk management plays a central role in institutional investors allocation of capital to trading. For instance, a risk manager
More informationRegression estimation in continuous time with a view towards pricing Bermudan options
with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten
More informationPricing theory of financial derivatives
Pricing theory of financial derivatives One-period securities model S denotes the price process {S(t) : t = 0, 1}, where S(t) = (S 1 (t) S 2 (t) S M (t)). Here, M is the number of securities. At t = 1,
More informationESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY
ESTIMATION OF UTILITY FUNCTIONS: MARKET VS. REPRESENTATIVE AGENT THEORY Kai Detlefsen Wolfgang K. Härdle Rouslan A. Moro, Deutsches Institut für Wirtschaftsforschung (DIW) Center for Applied Statistics
More informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
More information