Anumericalalgorithm for general HJB equations : a jump-constrained BSDE approach
|
|
- Dinah Boyd
- 5 years ago
- Views:
Transcription
1 Anumericalalgorithm for general HJB equations : a jump-constrained BSDE approach Nicolas Langrené Univ. Paris Diderot - Sorbonne Paris Cité, LPMA, FiME Joint work with Idris Kharroubi (Paris Dauphine), Huyên Pham (Paris Diderot) Workshop on Stochastic Games, Equilibrium, and Applications to Energy & Commodities Markets - Fields Institute, August 29, 2013 Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 1/20
2 Modeling volatility of an asset price S t Constant : σ>0 Deterministic : σ (t) Local : σ (t, S t ) Stochastic : dσ t =... Uncertain : σ [σ min,σ max ] Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 2/20
3 Uncertain Volatility Model Example ds t = σs t dw t σ [σ min,σ max ] uncertain Super-replication price Payoff Φ=Φ(T, S T ) P + 0 =sup E Q [Φ (T, S T )] Q Q Q = Q σ Q ; σ min σ Q σ max Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 3/20
4 Stochastic control problem with controlled driver & drift & volatility Formulation dx α s = b (X α s, α s ) ds + σ (X α s, α s ) dw s T v (t, x) =supe t,x f (Xs α, α s ) ds + g (XT α ) α A t General HJB equation v b t +sup (x, a).d x v + 12 tr σσ (x, a) Dx 2 v + f (x, a) =0 a A v (T, x) =g (x), x R d on [0,T ) R d Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 4/20
5 STEP 1 : Randomization of controls Poisson random measure µ A (dt, da) on R + A, W associated to the marked point process I (τ i,ζ i ) i,valuedina I t = ζ i, τ i t <τ i+1 Uncontrolled randomized problem Linear FBSDE Y t = g (X T )+ dx s = b (X s, I s ) ds + σ (X s, I s ) dw s T v (t, x, a) =E t,x,a f (X s, I s ) ds + g (X T ) T Y t v (t, X t, I t ) t f (X s, I s ) ds t T t Z s dw s T t A U s (a) µ A (ds, da) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 5/20
6 STEP 2 : Constraint on jumps U t (a) =v (t, X t, a) v (t, X t, I t ) Now, how to retrieve HJB? Add the constraint U t (a) 0 (t, a)! Jump-constrained BSDE Minimal solution (Y, Z, U, K) of Y t = g (X T )+ T + K T K t t f (X s, I s, Y s, Z s ) ds T subject to U t (a) 0 (t, a) t A T t Z s dw s U s (a) µ A (ds, da), 0 t T Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 6/20
7 Link with general HJB equations (X, I ) Markov v = v (t, x, a) s.t.y t = v (t, X t, I t ) Key Lemma v = v (t, x, a) does not depend on a! v = v (t, x) Theorem v = v (t, x) is solution of the HJB equation v b(x,a).d t +sup x v + 12 tr σσ (x,a)dx 2 v +f x,a,y,σ (x,a).d x v =0 a A v (T, x) =g (x), x R d on [0,T ) R d Proofs : cf. [Kharroubi, Pham, 2012] Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 7/20
8 Numerical scheme U t (a) =v (t, X t, a) v (t, X t, I t ) 0 (t, a) v (t, X t, I t ) sup a A v (t, X t, a) Minimal solution v (t, X t, I t )=sup a A v (t, X t, a) v (t, X t, I t ) v (t, X t, I t ) K t K t Forward-Backward numerical scheme Y N = g (X N ) Z i = E i Y i+1 W i / i Y i = E i [Y i+1 + f i (X i, I i, Y i+1, Z i ) i ] Y i = sup A A i E i,a [Y i ] where E i,a [.] :=E [. X i, I i = A] Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 8/20
9 Towards an implementable scheme How to compute the conditional expectations? (quantification, Malliavin calculus, empirical regression,...) cf. comparative tests in [Bouchard, Warin, 2012] Conditional expectation approximation E [U F ti ]=arg inf V L(F ti,p) E (V U) 2 where S L (F ti, P) 1 Ê [U F ti ]=arg inf V S M M (V m U m ) 2 m=1 Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 9/20
10 Empirical regression schemes First algorithm ( Tsitsiklis - van Roy ) Ŷ N = g (X N ) Ŷ i = Ê i Ŷ i+1 + f i (X i, I i ) i Ŷ i = sup A A i E i,a Ŷi Upward biased (up to Monte Carlo error & regression bias) Second algorithm ( Longstaff - Schwartz ) ˆα i =arg sup E i,a Ŷi A A i ˆX i+1 = b( ˆX i, ˆα i ) i + σ( ˆX i, ˆα i ) W i ˆv (t 0, x 0 )= 1 M N f ( ˆX i+1, ˆα i ) i + g( ˆX N ) M m=1 i=1 Downward biased (up to Monte Carlo error) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 10 /20
11 Uncertain correlation model Model dst i = σ i StdW i t i, i =1, 2 dw 1 t, dwt 2 = ρdt 1 ρ min ρ ρ max 1 Super-replication price Payoff Φ=Φ T, ST 1, ST 2 P 0 + =sup E Q Φ T, ST 1, S 2 T Q Q Q = Q ρ Q ; ρ min ρ Q ρ max Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 11 /20
12 Call spread on spread S 1 (T ) S 2 (T ) Φ=(S 1 (T) S 2 (T) K 1 ) + (S 1 (T) S 2 (T) K 2 ) + S 1 (0) S 2 (0) σ 1 σ 2 ρ min ρ max K 1 K 2 T Regression basis φ (t, s 1, s 2,ρ)=(K 2 K 1 ) S(β 0 +β 1 s 1 +β 2 s 2 +β 3 ρ+β 4 ρs 1 +β 5 ρs 2 ) S (x) =1/ (1 + exp ( x)) Bang-bang optimal control ρ (t, s 1, s 2 )=argmax α φ (t, s 1, s 2,ρ)=ρ max if β 3 +β 4 s 1 +β 5 s 2 0 = ρ min else Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 12 /20
13 Results Price of Call Spread on S1(T) S2(T) Superhedging ρ=ρ max ρ=0 ρ=ρ min Subhedging Moneyness ( = S1(0) S2(0) ) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 13 /20
14 Impact of correlation range Correlation Ranges (+ ) Price of Call Spread on S1(T) S2(T) Moneyness ( = S1(0) S2(0) ) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 14 /
15 Call Spread (S (T ) K 1 ) + (S (T ) K 2 ) + S (0) = 100, K 1 =90,K 2 =110,T =1,uncertainσ [0.1, 0.2] Estimated superreplication price: Algorithm 1 Time Step / / / / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 15 /20
16 Call Spread (S (T ) K 1 ) + (S (T ) K 2 ) + S (0) = 100, K 1 =90,K 2 =110,T =1,uncertainσ [0.1, 0.2] Estimated superreplication price: Algorithm 2 Time Step / / / / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 16 /20
17 Outperformer Spread (S 2 (T ) K 1 S 1 (T )) + (S 2 (T ) K 2 S 1 (T )) + S i (0) = 100, K 1 =0.9,K 2 =1.1,T =1,uncertainσ i [0.1, 0.2], ρ = 0.5 Estimated superreplication price: Algorithm 1 Time Step / / / / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 17 /20
18 Outperformer Spread (S 2 (T ) K 1 S 1 (T )) + (S 2 (T ) K 2 S 1 (T )) + S i (0) = 100, K 1 =0.9,K 2 =1.1,T =1,uncertainσ i [0.1, 0.2], ρ = 0.5 Estimated superreplication price: Algorithm 2 Time Step / / / / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 18 /20
19 Summary 1 Probabilistic representation of stochastic control problem with controlled volatility : jump-constrained BSDE 2 Numerical scheme for jump-constrained BSDEs 3 Application to pricing under uncertain volatility 4 Extension : stochastic games, HJB-Isaacs equations Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 19 /20
20 ???? Thank you! Questions? Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 20 /20
21 ???? References J-P. Lemor and E. Gobet and X. Warin (2006) Rate of convergence of an empirical regression method for solving generalized backward stochastic differential equations J. Guyon and P. Henry-Labordère (2011) Uncertain volatility model : a Monte Carlo approach E. Gobet and P. Turkedjiev (2011) Approximation of discrete BSDE using least-squares regression B. Bouchard and X. Warin (2012) Monte Carlo valorisation of American options : facts and new algorithms to improve existing methods I. Kharroubi and H. Pham (2012) Feynman-Kac representation for Hamilton-Jacobi-Bellman IPDE S. Alanko and M. Avellaneda (2013) Reducing variance in the numerical solution of BSDEs Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 21 /20
22 ???? Example A linear-quadratic stochastic control problem Problem dx α s =( µ 0 X α s + µ 1 α s ) ds +(σ 0 + σ 1 α s ) dw s X0 α =0 T v (t, x) =supe λ 0 α A t (α s ) 2 ds λ 1 (X α T ) 2 Set of parameters µ 0 µ 1 σ 0 σ 1 λ 0 λ 1 T Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 22 /20
23 ???? Numerical parameters n = 52 M = 10 6 time steps Monte Carlo simulations Regression basis φ (t, x,α)=β 0 + β 1 x + β 2 α + β 3 xα + β 4 x 2 + β 5 α 2 Linear optimal control α (t, x) =argmaxφ (t, x,α)=a (t) x + B (t) α A (t) = 0.5 β 3 /β 5 B (t) = 0.5 β 2 /β 5 Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 23 /20
24 ???? Estimated optimal controls Optimal Control α (t, x) Diffusion value x Time t 0 Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 24 /20
25 ???? Control impact (1/2) : no control Uncontrolled diffusion 10% 20% 30% 40% 50% 60% 70% 80% 90% 99% Interquantile Ranges Time Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 25 /20
26 ???? Control impact (2/2) : optimal control Optimally controlled diffusion 10% 20% 30% 40% 50% 60% 70% 80% 90% 99% Interquantile Ranges Time Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 26 /20
27 ???? Accuracy Comparison of control coefficients A(t) and B(t) B(t) estimated theoretical Value Function ˆv(0,0)= v(0,0)= Relative Error : 1% 4 A(t) estimated theoretical Time t Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 27 /20
28 ???? Convergence rate 1 Localizations 2 Theoretical regressions 3 Empirical regressions Assumptions p 1, L g, L f, C f,0 > 0s.t. i =0,...,N 1 g (x) g (x ) L g x p x p f i (x,a,y,z) f i (x,a,y,z ) L f ( x p x p + a p a p + y y + z z ) f i (0, 0, 0, 0) C f,0 Bounded control domain A R d : Ā > 0s.t. a A, a Ā Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 28 /20
29 ???? Theoretical regression (1/2) Definition ˆλ i (U) :=arg inf E (λ.p (X i, I i ) U) 2 λ R B P i (U) :=ˆλ i (U).p (X i, I i ) Associated scheme Ŷ N = g (X N ) ˆλ Y i = regression coefficients at time t i Ŷ i = sup A A i ˆλ Y i.p i (X i, A) Problem : Ŷ i is not itself the projection of some random variable... Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 29 /20
30 ???? Theoretical regression (2/2) Alternative definition ˆλ i,a (U) :=arg inf E (λ.p (X i, A) U A ) 2 λ R B P i,a (U) :=ˆλ i (U).p (X i, A) Regression error Yi 2 Zi 2 max Ŷ i, i Ẑ i N 1 e C(T t i ) k=i E sup Y k,a Pk,A Y (Y k,a ) 2 A A k +C k E sup Zk,A Pk,A Z (Z k,a ) 2 A A k But its empirical version cannot be (efficiently) implemented... Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 30 /20
31 ???? Outperformer (S 1 (T ) S 2 (T )) + S i (0) = 100, T =1,uncertainσ i [0.1, 0.2], ρ =0 Estimated superreplication price: Algorithm 1 Time Step 12 1/128 1/ / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 31 /20 1/8
32 ???? Outperformer (S 1 (T ) S 2 (T )) + S i (0) = 100, T =1,uncertainσ i [0.1, 0.2], ρ =0 Estimated superreplication price: Algorithm 2 Time Step 12 1/128 1/ / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 32 /20 1/8
33 ???? Outperformer Spread (S 2 (T ) K 1 S 1 (T )) + (S 2 (T ) K 2 S 1 (T )) + S i (0) = 100, K 1 =0.9,K 2 =1.1,T =1,uncertainσ i [0.1, 0.2] & ρ [ 0.5, 0.5] Estimated superreplication price: Algorithm 1 Time Step 14 1/ / / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 33 /20 1/8
34 ???? Outperformer Spread (S 2 (T ) K 1 S 1 (T )) + (S 2 (T ) K 2 S 1 (T )) + S i (0) = 100, K 1 =0.9,K 2 =1.1,T =1,uncertainσ i [0.1, 0.2] & ρ [ 0.5, 0.5] Estimated superreplication price: Algorithm 2 Time Step 14 1/ / / / log2(m) Nicolas Langrené A numerical algorithm for general HJB equations : a jump-constrained BSDE approach 34 /20 1/8
European option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
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 informationContinuous-time Stochastic Control and Optimization with Financial Applications
Huyen Pham Continuous-time Stochastic Control and Optimization with Financial Applications 4y Springer Some elements of stochastic analysis 1 1.1 Stochastic processes 1 1.1.1 Filtration and processes 1
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 information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationAn overview of some financial models using BSDE with enlarged filtrations
An overview of some financial models using BSDE with enlarged filtrations Anne EYRAUD-LOISEL Workshop : Enlargement of Filtrations and Applications to Finance and Insurance May 31st - June 4th, 2010, Jena
More informationIlliquidity, Credit risk and Merton s model
Illiquidity, Credit risk and Merton s model (joint work with J. Dong and L. Korobenko) A. Deniz Sezer University of Calgary April 28, 2016 Merton s model of corporate debt A corporate bond is a contingent
More informationAnalytical formulas for local volatility model with stochastic. Mohammed Miri
Analytical formulas for local volatility model with stochastic rates Mohammed Miri Joint work with Eric Benhamou (Pricing Partners) and Emmanuel Gobet (Ecole Polytechnique Modeling and Managing Financial
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 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 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 informationPricing in markets modeled by general processes with independent increments
Pricing in markets modeled by general processes with independent increments Tom Hurd Financial Mathematics at McMaster www.phimac.org Thanks to Tahir Choulli and Shui Feng Financial Mathematics Seminar
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More informationFinancial Mathematics and Supercomputing
GPU acceleration in early-exercise option valuation Álvaro Leitao and Cornelis W. Oosterlee Financial Mathematics and Supercomputing A Coruña - September 26, 2018 Á. Leitao & Kees Oosterlee SGBM on GPU
More informationModelling Credit Spread Behaviour. FIRST Credit, Insurance and Risk. Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent
Modelling Credit Spread Behaviour Insurance and Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent ICBI Counterparty & Default Forum 29 September 1999, Paris Overview Part I Need for Credit Models Part II
More informationContinuous Time Mean Variance Asset Allocation: A Time-consistent Strategy
Continuous Time Mean Variance Asset Allocation: A Time-consistent Strategy J. Wang, P.A. Forsyth October 24, 2009 Abstract We develop a numerical scheme for determining the optimal asset allocation strategy
More informationSemi-Markov model for market microstructure and HFT
Semi-Markov model for market microstructure and HFT LPMA, University Paris Diderot EXQIM 6th General AMaMeF and Banach Center Conference 10-15 June 2013 Joint work with Huyên PHAM LPMA, University Paris
More informationOptimal Trade Execution: Mean Variance or Mean Quadratic Variation?
Optimal Trade Execution: Mean Variance or Mean Quadratic Variation? Peter Forsyth 1 S. Tse 2 H. Windcliff 2 S. Kennedy 2 1 Cheriton School of Computer Science University of Waterloo 2 Morgan Stanley New
More informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More informationRobust Portfolio Choice and Indifference Valuation
and Indifference Valuation Mitja Stadje Dep. of Econometrics & Operations Research Tilburg University joint work with Roger Laeven July, 2012 http://alexandria.tue.nl/repository/books/733411.pdf Setting
More informationOption pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard
Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Indifference pricing and the minimal entropy martingale measure Fred Espen Benth Centre of Mathematics for Applications
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 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 informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
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 informationForward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions
Forward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions Pierre Henry-Labordère 1 1 Global markets Quantitative Research, SOCIÉTÉ GÉNÉRALE Outline 1 Introduction 2 Semi-linear PDEs 3
More informationToward a coherent Monte Carlo simulation of CVA
Toward a coherent Monte Carlo simulation of CVA Lokman Abbas-Turki (Joint work with A. I. Bouselmi & M. A. Mikou) TU Berlin January 9, 2013 Lokman (TU Berlin) Advances in Mathematical Finance 1 / 16 Plan
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 informationTwo-dimensional COS method
Two-dimensional COS method Marjon Ruijter Winterschool Lunteren 22 January 2013 1/29 Introduction PhD student since October 2010 Prof.dr.ir. C.W. Oosterlee). CWI national research center for mathematics
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 informationNumerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps
Numerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps, Senior Quantitative Analyst Motivation: Swing Options An electricity or gas SUPPLIER needs to be capable,
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
More informationCONTINUOUS TIME PRICING AND TRADING: A REVIEW, WITH SOME EXTRA PIECES
CONTINUOUS TIME PRICING AND TRADING: A REVIEW, WITH SOME EXTRA PIECES THE SOURCE OF A PRICE IS ALWAYS A TRADING STRATEGY SPECIAL CASES WHERE TRADING STRATEGY IS INDEPENDENT OF PROBABILITY MEASURE COMPLETENESS,
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More informationAsymmetric information in trading against disorderly liquidation of a large position.
Asymmetric information in trading against disorderly liquidation of a large position. Caroline Hillairet 1 Cody Hyndman 2 Ying Jiao 3 Renjie Wang 2 1 ENSAE ParisTech Crest, France 2 Concordia University,
More informationSample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models
Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson Carnegie Mellon University scottrob@andrew.cmu.edu http://www.math.cmu.edu/users/scottrob June
More informationMulti-period mean variance asset allocation: Is it bad to win the lottery?
Multi-period mean variance asset allocation: Is it bad to win the lottery? Peter Forsyth 1 D.M. Dang 1 1 Cheriton School of Computer Science University of Waterloo Guangzhou, July 28, 2014 1 / 29 The Basic
More informationConvergence Analysis of Monte Carlo Calibration of Financial Market Models
Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration
More informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
More information2.1 Mean-variance Analysis: Single-period Model
Chapter Portfolio Selection The theory of option pricing is a theory of deterministic returns: we hedge our option with the underlying to eliminate risk, and our resulting risk-free portfolio then earns
More informationMSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013
MSc Financial Engineering 2012-13 CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL To be handed in by monday January 28, 2013 Department EMS, Birkbeck Introduction The assignment consists of Reading
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationOptimal Asset Allocation with Stochastic Interest Rates in Regime-switching Models
Optimal Asset Allocation with Stochastic Interest Rates in Regime-switching Models Ruihua Liu Department of Mathematics University of Dayton, Ohio Joint Work With Cheng Ye and Dan Ren To appear in International
More informationValuing American Options by Simulation
Valuing American Options by Simulation Hansjörg Furrer Market-consistent Actuarial Valuation ETH Zürich, Frühjahrssemester 2008 Valuing American Options Course material Slides Longstaff, F. A. and Schwartz,
More informationThe Evaluation of Swing Contracts with Regime Switching. 6th World Congress of the Bachelier Finance Society Hilton, Toronto June
The Evaluation of Swing Contracts with Regime Switching Carl Chiarella, Les Clewlow and Boda Kang School of Finance and Economics University of Technology, Sydney Lacima Group, Sydney 6th World Congress
More informationHelp Session 2. David Sovich. Washington University in St. Louis
Help Session 2 David Sovich Washington University in St. Louis TODAY S AGENDA Today we will cover the Change of Numeraire toolkit We will go over the Fundamental Theorem of Asset Pricing as well EXISTENCE
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 informationGamma. The finite-difference formula for gamma is
Gamma The finite-difference formula for gamma is [ P (S + ɛ) 2 P (S) + P (S ɛ) e rτ E ɛ 2 ]. For a correlation option with multiple underlying assets, the finite-difference formula for the cross gammas
More informationThe Uncertain Volatility Model
The Uncertain Volatility Model Claude Martini, Antoine Jacquier July 14, 008 1 Black-Scholes and realised volatility What happens when a trader uses the Black-Scholes (BS in the sequel) formula to sell
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationKØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours
This question paper consists of 3 printed pages FinKont KØBENHAVNS UNIVERSITET (Blok 2, 211/212) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours This exam paper
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 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 informationExact replication under portfolio constraints: a viability approach
Exact replication under portfolio constraints: a viability approach CEREMADE, Université Paris-Dauphine Joint work with Jean-Francois Chassagneux & Idris Kharroubi Motivation Complete market with no interest
More informationExact Sampling of Jump-Diffusion Processes
1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance
More informationDoubly reflected BSDEs with jumps and generalized Dynkin games
Doubly reflected BSDEs with jumps and generalized Dynkin games Roxana DUMITRESCU (University Paris Dauphine, Crest and INRIA) Joint works with M.C. Quenez (Univ. Paris Diderot) and Agnès Sulem (INRIA Paris-Rocquecourt)
More informationPuttable Bond and Vaulation
and Vaulation Dmitry Popov FinPricing http://www.finpricing.com Summary Puttable Bond Definition The Advantages of Puttable Bonds Puttable Bond Payoffs Valuation Model Selection Criteria LGM Model LGM
More informationDynamic Portfolio Choice II
Dynamic Portfolio Choice II Dynamic Programming Leonid Kogan MIT, Sloan 15.450, Fall 2010 c Leonid Kogan ( MIT, Sloan ) Dynamic Portfolio Choice II 15.450, Fall 2010 1 / 35 Outline 1 Introduction to Dynamic
More informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
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 informationarxiv: v1 [q-fin.pm] 13 Mar 2014
MERTON PORTFOLIO PROBLEM WITH ONE INDIVISIBLE ASSET JAKUB TRYBU LA arxiv:143.3223v1 [q-fin.pm] 13 Mar 214 Abstract. In this paper we consider a modification of the classical Merton portfolio optimization
More informationOptimal Selling Strategy With Piecewise Linear Drift Function
Optimal Selling Strategy With Piecewise Linear Drift Function Yan Jiang July 3, 2009 Abstract In this paper the optimal decision to sell a stock in a given time is investigated when the drift term in Black
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 informationLECTURE 4: BID AND ASK HEDGING
LECTURE 4: BID AND ASK HEDGING 1. Introduction One of the consequences of incompleteness is that the price of derivatives is no longer unique. Various strategies for dealing with this exist, but a useful
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationModeling the dependence between a Poisson process and a continuous semimartingale
1 / 28 Modeling the dependence between a Poisson process and a continuous semimartingale Application to electricity spot prices and wind production modeling Thomas Deschatre 1,2 1 CEREMADE, University
More informationAsset Pricing Models with Underlying Time-varying Lévy Processes
Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
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 informationExponential utility maximization under partial information
Exponential utility maximization under partial information Marina Santacroce Politecnico di Torino Joint work with M. Mania AMaMeF 5-1 May, 28 Pitesti, May 1th, 28 Outline Expected utility maximization
More informationRisk minimizing strategies for tracking a stochastic target
Risk minimizing strategies for tracking a stochastic target Andrzej Palczewski Abstract We consider a stochastic control problem of beating a stochastic benchmark. The problem is considered in an incomplete
More informationPricing and hedging in incomplete markets
Pricing and hedging in incomplete markets Chapter 10 From Chapter 9: Pricing Rules: Market complete+nonarbitrage= Asset prices The idea is based on perfect hedge: H = V 0 + T 0 φ t ds t + T 0 φ 0 t ds
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 informationValuing Early Stage Investments with Market Related Timing Risk
Valuing Early Stage Investments with Market Related Timing Risk Matt Davison and Yuri Lawryshyn February 12, 216 Abstract In this work, we build on a previous real options approach that utilizes managerial
More informationOn worst-case investment with applications in finance and insurance mathematics
On worst-case investment with applications in finance and insurance mathematics Ralf Korn and Olaf Menkens Fachbereich Mathematik, Universität Kaiserslautern, 67653 Kaiserslautern Summary. We review recent
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 informationMultilevel Monte Carlo Methods for American Options
Multilevel Monte Carlo Methods for American Options Simon Gemmrich, PhD Kellog College University of Oxford A thesis submitted in partial fulfillment of the MSc in Mathematical Finance November 19, 2012
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationCallable Bond and Vaulation
and Vaulation Dmitry Popov FinPricing http://www.finpricing.com Summary Callable Bond Definition The Advantages of Callable Bonds Callable Bond Payoffs Valuation Model Selection Criteria LGM Model LGM
More informationMonte Carlo Based Numerical Pricing of Multiple Strike-Reset Options
Monte Carlo Based Numerical Pricing of Multiple Strike-Reset Options Stavros Christodoulou Linacre College University of Oxford MSc Thesis Trinity 2011 Contents List of figures ii Introduction 2 1 Strike
More informationApplication of Stochastic Calculus to Price a Quanto Spread
Application of Stochastic Calculus to Price a Quanto Spread Christopher Ting http://www.mysmu.edu/faculty/christophert/ Algorithmic Quantitative Finance July 15, 2017 Christopher Ting July 15, 2017 1/33
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationEvaluating the Longstaff-Schwartz method for pricing of American options
U.U.D.M. Project Report 2015:13 Evaluating the Longstaff-Schwartz method for pricing of American options William Gustafsson Examensarbete i matematik, 15 hp Handledare: Josef Höök, Institutionen för informationsteknologi
More informationOptimal Securitization via Impulse Control
Optimal Securitization via Impulse Control Rüdiger Frey (joint work with Roland C. Seydel) Mathematisches Institut Universität Leipzig and MPI MIS Leipzig Bachelier Finance Society, June 21 (1) Optimal
More information************* with µ, σ, and r all constant. We are also interested in more sophisticated models, such as:
Continuous Time Finance Notes, Spring 2004 Section 1. 1/21/04 Notes by Robert V. Kohn, Courant Institute of Mathematical Sciences. For use in connection with the NYU course Continuous Time Finance. This
More informationRough volatility models: When population processes become a new tool for trading and risk management
Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum
More informationOptimal decumulation into annuity after retirement: a stochastic control approach
Optimal decumulation into annuity after retirement: a stochastic control approach Nicolas Langrené, Thomas Sneddon, Geo rey Lee, Zili Zhu 2 nd Congress on Actuarial Science and Quantitative Finance, Cartagena,
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 informationOptimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error
Optimum Thresholding for Semimartingales with Lévy Jumps under the mean-square error José E. Figueroa-López Department of Mathematics Washington University in St. Louis Spring Central Sectional Meeting
More informationDynamic pricing with diffusion models
Dynamic pricing with diffusion models INFORMS revenue management & pricing conference 2017, Amsterdam Asbjørn Nilsen Riseth Supervisors: Jeff Dewynne, Chris Farmer June 29, 2017 OCIAM, University of Oxford
More informationInterest Rate Bermudan Swaption Valuation and Risk
Interest Rate Bermudan Swaption Valuation and Risk Dmitry Popov FinPricing http://www.finpricing.com Summary Bermudan Swaption Definition Bermudan Swaption Payoffs Valuation Model Selection Criteria LGM
More informationNotes. Cases on Static Optimization. Chapter 6 Algorithms Comparison: The Swing Case
Notes Chapter 2 Optimization Methods 1. Stationary points are those points where the partial derivatives of are zero. Chapter 3 Cases on Static Optimization 1. For the interested reader, we used a multivariate
More informationINSTITUTE AND FACULTY OF ACTUARIES. Curriculum 2019 SPECIMEN SOLUTIONS
INSTITUTE AND FACULTY OF ACTUARIES Curriculum 2019 SPECIMEN SOLUTIONS Subject CM1A Actuarial Mathematics Institute and Faculty of Actuaries 1 ( 91 ( 91 365 1 0.08 1 i = + 365 ( 91 365 0.980055 = 1+ i 1+
More informationAnurag Sodhi University of North Carolina at Charlotte
American Put Option pricing using Least squares Monte Carlo method under Bakshi, Cao and Chen Model Framework (1997) and comparison to alternative regression techniques in Monte Carlo Anurag Sodhi University
More informationPricing Early-exercise options
Pricing Early-exercise options GPU Acceleration of SGBM method Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee Lausanne - December 4, 2016
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 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 informationRobust Portfolio Decisions for Financial Institutions
Robust Portfolio Decisions for Financial Institutions Ioannis Baltas 1,3, Athanasios N. Yannacopoulos 2,3 & Anastasios Xepapadeas 4 1 Department of Financial and Management Engineering University of the
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 information