Multiname and Multiscale Default Modeling
|
|
- Jared Ross
- 5 years ago
- Views:
Transcription
1 Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis on Finance Kick-off-Workshop September 8-12, 2008 RICAM, Linz, Austria 1
2 Basics of Collaterized Debt Obligations (CDOs) e rt k f l (T k 1 )α l (T k T k 1 ) = k k e rt k (f l (T k 1 ) f l (T k )) (1 R) Predetermined payment times: T k, k {1, 2,,K} Constant risk-free short rate r Constant recovery R Tranches: {0-3, 3-7, 7-10, 10-15, 15-30}% (CDX tranches) α l is the yield associated with tranche l: the rate at which the insurance buyer pays for protection of tranche l f l (T k ) is the risk-neutral expected fraction of tranche l left at time T k (defaults between T k 1 and T k being accounted at T k ) 2
3 Doubly Stochastic Modeling We consider N obligors or underlying names. A particular obligor i defaults at the first arrival τ i of a Cox process with stochastic intensity or hazard rate X (i). Conditioned on the paths of the hazard rates, the default times τ i of the firms are independent, and the probability that obligor i has survived till time T, is given by exp( T 0 X(i) s ds). The unconditional survival probability is { IP {τ i > T } = IE e T } 0 X(i) s. Consider a subset {i 1,,i n } of obligors. The probability of the joint survival of this set till time T is then IP {τ i1 > T,,τ in > T } = IE { e n j=1 T } 0 X(i j ) s ds. 3
4 Vasicek Intensities and Survival Probabilities Negative intensities? Duffie-Singleton, Credit Risk 2003: the computational advantage with explicit solutions may be worth the approximation error associated with this Gaussian formulation. where the (W (i) t ( ) dx (i) t = κ i θ i X (i) t dt + σ i dw (i) t, ) are correlated Brownian motions with d W (i), W (j) = c ij dt. t We denote the survival probability for name i by { S i (T; x i ) = IP (τ i > T X (i) 0 = x i ) = IE e T 0 X(i) s ds X (i) 0 = x i }. We also denote the joint survival probability of all N names by { S(T;x, N) = IE e T N } 0 X(i) s ds X (1) 0 = x 1,,X (N) 0 = x N, 4
5 PDE for the Survival Probabilities From the Feynman-Kac formula it follows that the joint survival probability from time t till time T { u(t,x) = IE e T N } t X(i) s ds X t = x, solves the partial differential equation u t i,j=1 (σ i σ j c ij ) 2 u x i x j + κ i (θ i x i ) u x i ( N x i )u = 0, with terminal condition u(t,x) = 1. 5
6 Independent Case Assume first that the covariance matrix c is the identity matrix, corresponding to the components of the intensity process X being independent. Then, as is well known, or can be readily checked, the solution is given by the product-affine formula where B i (s) = s 0 A i (s) = e θ i u(t,x) = N A i (T t)e B i(t t)x i, e κiξ dξ = 1 e κ iτ, κ i s and θ d i = θ i σ2 i 2κ 2 i 0 κ ib i (ξ) dξ+ 1 2 σ2 i. s 0 B2 i (ξ) dξ = e ( θ d i (s B i(s))+ σ2 i 4κ i B 2 i (s) ), 6
7 Correlated Case In the general correlated case, we find with N u(t,x) = A c (T t) A i (T t)e B i(t t)x i, A c (s) = e 1 2 N N j (σ iσ j c ij ) s 0 B i(ξ)b j (ξ)dξ. The last integral is given explicitly by s 0 B i (ξ)b j (ξ) dξ = s κ i κ j B i (s) κ i (κ i + κ j ) B j (s) κ j (κ i + κ j ) B i(s)b j (s) (κ i + κ j ). 7
8 Symmetric Name Case The dynamics and the starting points of the intensities are the same for all the names. This is convenient to understand the effects of the correlation and the size of the portfolio. Specifically, we have ( = κ dx (i) t θ X (i) t ) dt + σ dw (i) t, X (i) 0 = x, with the parameters κ, θ σ and x assumed constant and positive. Moreover, we assume that the correlation matrix is defined by c ij = ρ X, for i j, with ρ X 0, and ones on the diagonal. Remark that such a correlation structure can be obtained by letting W (i) t = 1 ρ X W (i) t + ρ X W (0) t, where W (i), i = 0, 1,, N, are independent standard Brownian motions. 8
9 Explicit Survival Probabilities The joint survival probability for n given names, say the first n names, is S(T;x, n) = IE { e T with 0 (X(1) s } + +X(n) s ) ds X (1) 0 = x,,x (n) 0 = x = e n [θ (T B(T))+[1+(n 1)ρ X ]σ 2 B 2 (T)/(4κ)+xB(T)], B(T) = 1 e κt κ θ = θ [1 + (n 1)ρ X ] σ2 2κ 2. This expression shows explicitly how the joint survival probability depends on the correlation ρ X and the basket size n. Note in particular how the basket size enhances the correlation effect., 9
10 The Loss Distribution The loss distribution at time T of a basket of size N is given by its mass function p n = IP {(#names defaulted at time T) = n}, n = 0, 1,,N, and, in the symmetric correlated case, is explicitly p n = N n n S N+j n ( 1) j, n j j=0 using the short hand notation S n = S(T; (x,, x), n) for the joint survival probability of n names. This gives rise to an O(N 2 ) procedure for calculating the loss distribution. However, a direct implementation of this formula is not numerically stable due to catastrophic cancellation errors in finite precision arithmetics. 10
11 Alternative Implementation Note first that, from the explicit formula for the joint survival probabilities, we can write with S n = e n [θ (T B(T))+[1+(n 1)ρ X ]σ 2 B 2 (T)/(4κ)+xB(T)] = e d 1n+d 2 n 2, d 1 = d 1 (T, x) = θ(t B(T)) + xb(t) 1 2 σ2 (1 ρ X )B (2) (T), d 2 = d 2 (T) = 1 2 σ2 ρ X B (2) (T), B (2) (T) = T 0 B 2 (s) ds = (T B(T)) κ 2 B(T)2 2κ, where we assume that the model parameters are chosen so that d 1 > 0. 11
12 Alternative Implementation (continued) In the independent case ρ X = 0, we get the binomial distribution: p n = N n (1 e d 1 ) n e (N n)d 1 =: p n (d 1 ). In the general case, we can write { S n = IE e d 1n n } 2d 2 Z = IE { e n (d 1 + 2d 2 Z) }, for Z (or Z) a zero mean unit variance Gaussian random variable. Therefore, in the general case we find { p n = IE p n (d 1 + } 2d 2 Z). Thus, we get the loss distribution stably and fast by integrating (non-negative) binomial distributions with respect to the Gaussian density. 12
13 Alternative Implementation (final comments) We remark that p n = IE { p n (d 1 + } 2d 2 Z) essentially corresponds to conditioning with respect to the correlating Brownian motion W (0), in W (i) t = 1 ρ X W (i) t + ρ X W (0) t. The argument d 1 + 2d 2 Z will be negative for Z negative and with large magnitude. This reflects the fact that we are using a Vasicek model where the intensity may be negative. In the following numerical illustrations, we condition the Gaussian density to Z > d 1 / 2d 2, and choose parameters such that the complementary event has probability less than
14 Example with Constant Parameters and Strong Correlations 0.12 LOSS DISTRIBUTION INDEPENDENT STRONG CORRELATIONS TRANCHE PREMIA INDEPENDENT STRONG CORRELATIONS INDEPENDENT STRONG CORRELATIONS θ =.02, κ =.5, σ =.015, x =.02, T = 5, N = 125, r = 3%, R = 40%, ρ X = 0, or.75 14
15 Multiscale Stochastic Volatility Modeling Under the risk-neutral probability measure we assume the model dx (i) t = κ i (θ i X (i) t )dt + σ (i) t dw (i) t, for 1 i N where the W (j) s are correlated Brownian motions as before. The volatilities are stochastic and depend on a fast evolving factor Y and slowly evolving factor Z: σ (i) t = σ i (Y t, Z t ), where the functions σ i (y, z) are positive, bounded and bounded away from zero, and smooth in the second variable. The fast process is modeled for instance by dy t = 1 ε (m Y t)dt + ν 2 dw (y) ε t, ε << 1, d W (i), W (y) = ρ Y dt, for 1 i N. t 15
16 Multiscale Stochastic Volatility Modeling (continued) The slow factor evolves as dz t = δc(z)dt + δg(z) dw (z) t, 1/δ >> 1, d W (i), W (z) = ρ Z dt. t The functions and c and g are assumed to be smooth. We denote by ρ Y Z the correlation coefficient defined by d W (y), W (z) = ρ Y Z dt t We assume the coefficients ρ Y, ρ Z, ρ Y Z and the matrix (c ij ) are such that the joint covariance matrix of the Brownian motions W (i) (i = 1,,n), W (y) and W (z) is non-negative definite. 16
17 S(T;x, y, z, N) = IE { e N Joint Survival Probabilities T 0 X(i) s } ds X 0 = x, Y 0 = y, Z 0 = z. In this case the joint survival probability from time t { u ε,δ (t,x; y, z, N) = IE e T N } t X(i) s ds X t = x, Y t = y, Z t = z, solves the partial differential equation with the notation L ε,δ u ε,δ = 0, u ε,δ (T,x, y, z) = 1 L ε,δ = t + L (x,y,z) ( N ) x i where L (x,y,z) denotes the infinitesimal generator of the Markov process (X t, Y t, Z t ) under the risk-neutral measure., 17
18 Operator Notations L ε,δ = 1 ε L ε L 1 + L 2 + δm 1 + δm 2 + δ ε M 3, L 0 = ν 2 2 L 1 = 2νρ Y L 2 = t M 1 = ρ Z g(z) y 2 + (m y) y, i,j=1 2 σ i (y, z) x i y, 2 c ij σ i (y, z)σ j (y, z) + x i x j σ i (y,z) 2 x i z, M 2 = 1 2 g2 (z) 2 z 2 + c(z) z, M 3 = 2νg(z)ρ Y Z 2 y z. ( κ i (θ i x i ) N x i x i ), 18
19 Expansions Starting with u ε,δ = u ε 0 + δ u ε 1 + δ u ε 2 =, we get ( 1 ε L ) L 1 + L 2 u ε ε 0 = 0, u ε 0(T) = 1, ( 1 ε L ) L 1 + L 2 u ε ε 1 = (M 1 + ε 1 ) M 3 u ε 0, u ε 1(T) = 0. Then we expand u ε 0 = u 0 + ε u 1,0 + εu 2,0 + ε ε u 3,0 +. and get L 0 u 0 = 0, (choose u 0 independent of y), L 0 u 1,0 + L 1 u 0 =0 = 0, (choose u 1,0 independent of y), L 0 u 2,0 + L 1 u 1,0 =0 + L 2 u 0 = 0, (Poisson equation in u 2,0 ), Centering condition: L 2 u 0 = 0, where the triangular brackets denote an average with respect to the invariant distribution of the fast process Y. 19
20 Recall that L 2 = t +1 2 i,j=1 Leading Order Term 2 c ij σ i (y, z)σ j (y, z) + x i x j ( κ i (θ i x i ) N ) x i x i The leading order term u 0 of the joint survival probability is obtained by solving the constant volatility problem with the effective diffusion [d(z)] ij = d ij (z) := c ij σ i (y, z)σ j (y, z)φ(y) dy, with Φ being the invariant distribution for the Y process. Since the process Y evolves on the fast scale its leading order effect is obtained by averaging with respect to Φ. The process Z evolves on a relatively slow scale and at this level of approximation its effect corresponds to just evaluating this process at its current frozen level z. 20
21 Formula for the Leading Order Term The leading order term u 0 is the survival probability which solves where L e (d(z)) = t +1 2 L e (d(z))u 0 = 0, u 0 (T,x; z) = 1, i,j=1 2 d ij (z) + x i x j ( κ i (θ i x i ) N x i x i ) As previously in the constant volatility case, the solution is given by with u 0 (t,x; z) = A c (T t) N A i (T t)e B i(t t)x i, B i (s) = 1 e κ iτ, κ i A c (s) = e 1 2 N N j (d ij(z)) s 0 B i(ξ)b j (ξ) dξ. 21
22 Expansions (continued) Coming back to the expansion u ε 0 = u 0 + ε u 1,0 + εu 2,0 + ε ε u 3,0 +, we get u 2,0 = L 1 0 (L 2 L 2 ) u 0, L 0 u 3,0 + L 1 u 2,0 + L 2 u 1,0 = 0, (Poisson equation in u 3,0 ), so that the centering condition becomes L 1 u 2,0 + L 2 u 1,0 = 0, L 2 u 1,0 = L 1 u 2,0 = L 1 L 1 0 (L 2 L 2 ) u 0 22
23 Correction u 1,0 due to the fast volatility factor Y We introduce the operator A 1,0 = L 1 L 1 0 (L 2 L 2 ), which will be given explicitly below. The function u 1,0 (t,x, z) solves the inhomogeneous problem L e (d(z))u 1,0 = A 1,0 u 0, u 1,0 (T,x; z) = 0. Thus, u 1,0 solves a linear equation with the effective operator L e (d(z)), but now the problem involves a source term A 1,0 u 0, defined in terms of the leading order survival probability u 0, and with a zero terminal condition. 23
24 Expansions (continued) Coming back to the term δ u ε 1, we have ( 1 ε L ) L 1 + L 2 u ε ε 1 = (M 1 + ε 1 ) M 3 u ε 0, u ε 1(T) = 0. We expand u ε 1 = u 0,1 + ε u 1,1 + εu 2,1 +, and keep only the terms of order one, to obtain as previously: L 2 u 0,1 = M 1 u 0 where we have used that M 3 takes derivatives with respect to y. 24
25 Correction u 0,1 due to the slow volatility factor Z In this case we introduce the operator and obtain: A 0,1 = M 1, The function u 0,1 (t,x, z) is the solution of the problem L e (d(z))u 0,1 = A 0,1 u 0, u 0,1 (T,x; z) = 0, which is again a source problem with respect to the operator L e (d(z)) and with a zero terminal condition. 25
26 Explicit Formulas and Group Parameters We introduce the symmetric matrix Ψ(y, z) satisfying the Poisson equations L 0 Ψ i1,i 2 = c i1,i 2 σ i1 (y, z)σ i2 (y, z) d i1,i 2 (z), and the coefficients V ε 3 (z, i 1, i 2, i 3 ) = ε ρ Y ν 2 σ i3 Ψ i1 i 2 y One then obtains that the scaled operator εa 1,0 can be written εa1,0 = ε L 1 L 1 0 (L 2 L 2 ) = i 1,i 2,i 3 =1. V3 ε 3 (z, i 1, i 2, i 3 ). x i1 x i2 x i3 26
27 Explicit Formula for the Correction ε u 1,0 Setting ε u 1,0 (t,x; z) = D(T t)u 0 (t,x; z), leads to the following ODE for D: D = i 1,i 2,i 3 =1 V ε 3 (z, i 1, i 2, i 3 )B i1 B i2 B i3, D(0) = 0, which gives εu1,0 = i 1,i 2,i 3 =1 T V3 ε (z, i 1, i 2, i 3 ) 0 B i1 (s)b i2 (s)b i3 (s) ds u 0. It depends on the underlying model structure in a complicated way, but only the effective market group parameters V ε 3 ( ) are needed to compute the fast time scale correction u 1,0. 27
28 Explicit Formula for the Correction δ u 0,1 Using M 1 = ρ Z g(z) N σ i(y, z) 2 x i z, we find that δa0,1 = δ M 1 = where we introduced the coefficients V δ V δ 1 (z, i) = δ g(z)ρ Z σ i. 2 1 (z, i), z x i As for the fast correction, it follows that δ u0,1 = ( 1 2 i 1 =1 V δ 1 (z, i 1 ) i 2 =1 i 3 =1 d T dz (d i 1 i 2 ) 0 B i1 (v) v 0 B i1 (s)b i2 (s) ds dv ) u 0. 28
29 Accuracy Result For any fixed t < T, x IR N and y, z IR, ( u ε,δ (t,x, y, z) u 0 (t,x, z) + ε u 1,0 (t,x, z) + δ u 0,1 (t,x, z)) = O(ε+δ), where u ε,δ is the solution of the original problem L ε,δ u ε,δ = 0, u ε,δ (T,x, y, z) = 1, and u 0, ε u 1,0 and δ u 0,1 are given by the formulas above. One of the difficulties is that the potential ( x i i) in ( N L ε,δ = t + L (x,y,z) x i ) is unbounded from above since the X (i) s are unbounded from below., 29
30 The transformation Accuracy Result (sketch of proof) u ε,δ (t,x, y, z) = M ε,δ (t, y, z) N e B i(t t)x i reduces to a Feynman Kac equation for M with a bounded potential in (y, z), bounded time-dependent coefficients, and smooth terminal condition M ε,δ (T, y, z) = 1. The proof given in the one-dimensional equity case by Fouque Papanicolaou Sircar Sølna is then easily extended. 30
31 Stochastic Volatility Effects in the Symmetric Case In the symmetric case with σ (i) t σ(y t, Z t ), we obtain S(T;x, y, z, n) = u ε,δ (0; (x,,x), y, z, n) ũ ε,δ (0; (x,,x), z, n) = [ ] (1 + D ε (T; z, n) + D δ (T; z, n))e n θ (z)(t B(T))+[1+(n 1)ρ X ] σ 2 (z) B2 (T) 4κ +xb(t) where σ 2 (z) = σ(, z) 2, D ε (T; z, n) = v 3 (z)n 2 (1 + (n 1)ρ X )B (3) (T), B (3) (T) = D δ (T; z,n) = v 1 (z)n 2 (1 + (n 1)ρ X ) B (3) (T), B(3) (T) = v 1 (z) = δ 2 g(z)ρ Z σ(, z) z L 0 Ψ = σ 2 (y, z) σ 2 (, z), T 0 T σ 2 (, z), v 3 (z) = ε 2 ρ Y ν 0 θ (z) = θ [1 + (n 1)ρ X ] σ2 (z) 2κ 2. B 3 (s) ds, B(s)B (2) (s) ds, Ψ(, z) σ(,z) y 31
32 Stochastic Volatility Effects in the Symmetric Case (continued) Using the constant volatility notation, we have with S(T; x, z, n) ( )) (1 + n 3 ρ X v 3 (z)b (3) (T) + v 1 (z) B (3) (T) e nd 1(T,x,z)+n 2 d 2 (T,z), d 1 (T,x,z) = θ(t B(T)) + xb(t) 1 2 (1 ρ X) σ(z) 2 B (2) (T), d 2 (T, z) = 1 2 ρ X σ(z) 2 B (2) (T), d 2 (T,z) = d 2 (T, z) + (1 ρ X ) ( ) v 3 (z)b (3) (T) + v 1 (z) B (3) (T), 32
33 Stochastic Volatility Effects on the Loss Distribution Recall the case with constant volatility and ρ X = 0: p n (x ) = N n n e x (N+j n) ( 1) j n j j=0 = N n (1 e x ) n e (N n)x. In the general case with ρ X 0, we introduce as before Z N(0, 1), we replace d 2 by d 2, to obtain the loss distribution: { p n = IE p n (d d ( ) 2 Z) + ρ X (v 3 B (3) + v B(3) 1 ) p n (d 1 + } 2 d 2 Z), Note that n=0 p n (x) = d3 dx 3 n=0 p n (x) = 0 p n = 1 (but p n may be < 0). n=0 33
34 Effect of Combined Name Correlation and Stochastic Volatility LOSS DISTRIBUTION CONSTANT PARAMETERS STOCHASTIC VOLATILITY TRANCHE PREMIA CONSTANT PARAMETERS STOCHASTIC VOLATILITY CONSTANT PARAMETERS STOCHASTIC VOLATILITY θ =.03, κ =.5, σ =.02, x =.03, T = 5, N = 125, r = 3%, R = 40%, ρ X =.01, and v 3 = , v 1 =
35 Summary Models with Vasicek intensities are computationally tractable. Stochastic volatility can be implemented via perturbation techniques (regular and singular). Strong effect on loss distribution and on senior tranches. Generalizations: Fluctuations in the hazard rate level (θ stochastic). Short rate term structure effect (r stochastic). Name heterogeneity via grouping. 35
Multiscale Stochastic Volatility Models
Multiscale Stochastic Volatility Models Jean-Pierre Fouque University of California Santa Barbara 6th World Congress of the Bachelier Finance Society Toronto, June 25, 2010 Multiscale Stochastic Volatility
More informationMultiscale Stochastic Volatility Models Heston 1.5
Multiscale Stochastic Volatility Models Heston 1.5 Jean-Pierre Fouque Department of Statistics & Applied Probability University of California Santa Barbara Modeling and Managing Financial Risks Paris,
More informationApplications to Fixed Income and Credit Markets
Applications to Fixed Income and Credit Markets Jean-Pierre Fouque University of California Santa Barbara 28 Daiwa Lecture Series July 29 - August 1, 28 Kyoto University, Kyoto 1 Fixed Income Perturbations
More informationAsian Options under Multiscale Stochastic Volatility
Contemporary Mathematics Asian Options under Multiscale Stochastic Volatility Jean-Pierre Fouque and Chuan-Hsiang Han Abstract. We study the problem of pricing arithmetic Asian options when the underlying
More informationCalibration to Implied Volatility Data
Calibration to Implied Volatility Data Jean-Pierre Fouque University of California Santa Barbara 2008 Daiwa Lecture Series July 29 - August 1, 2008 Kyoto University, Kyoto 1 Calibration Formulas The implied
More informationVolatility Time Scales and. Perturbations
Volatility Time Scales and Perturbations Jean-Pierre Fouque NC State University, soon UC Santa Barbara Collaborators: George Papanicolaou Stanford University Ronnie Sircar Princeton University Knut Solna
More informationStochastic Volatility Effects on Defaultable Bonds
Stochastic Volatility Effects on Defaultable Bonds Jean-Pierre Fouque Ronnie Sircar Knut Sølna December 24; revised October 24, 25 Abstract We study the effect of introducing stochastic volatility in the
More informationEvaluation of compound options using perturbation approximation
Evaluation of compound options using perturbation approximation Jean-Pierre Fouque and Chuan-Hsiang Han April 11, 2004 Abstract This paper proposes a fast, efficient and robust way to compute the prices
More informationMULTISCALE STOCHASTIC VOLATILITY ASYMPTOTICS
MULTISCALE STOCHASTIC VOLATILITY ASYMPTOTICS JEAN-PIERRE FOUQUE, GEORGE PAPANICOLAOU, RONNIE SIRCAR, AND KNUT SOLNA Abstract. In this paper we propose to use a combination of regular and singular perturbations
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 informationVariance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment
Variance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment Jean-Pierre Fouque Tracey Andrew Tullie December 11, 21 Abstract We propose a variance reduction method for Monte Carlo
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
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 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 informationStochastic Volatility Modeling
Stochastic Volatility Modeling Jean-Pierre Fouque University of California Santa Barbara 28 Daiwa Lecture Series July 29 - August 1, 28 Kyoto University, Kyoto 1 References: Derivatives in Financial Markets
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 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 informationMULTISCALE STOCHASTIC VOLATILITY FOR EQUITY, INTEREST RATE, AND CREDIT DERIVATIVES
MULTISCALE STOCHASTIC VOLATILITY FOR EQUITY, INTEREST RATE, AND CREDIT DERIVATIVES Building upon the ideas introduced in their previous book, Derivatives in Financial Markets with Stochastic Volatility,
More informationCalibration of Stock Betas from Skews of Implied Volatilities
Calibration of Stock Betas from Skews of Implied Volatilities Jean-Pierre Fouque Eli Kollman January 4, 010 Abstract We develop call option price approximations for both the market index and an individual
More informationOption Pricing Under a Stressed-Beta Model
Option Pricing Under a Stressed-Beta Model Adam Tashman in collaboration with Jean-Pierre Fouque University of California, Santa Barbara Department of Statistics and Applied Probability Center for Research
More informationOption Pricing Models for European Options
Chapter 2 Option Pricing Models for European Options 2.1 Continuous-time Model: Black-Scholes Model 2.1.1 Black-Scholes Assumptions We list the assumptions that we make for most of this notes. 1. The underlying
More informationAsymptotic Pricing of Commodity Derivatives using Stochastic Volatility Spot Models
Asymptotic Pricing of Commodity Derivatives using Stochastic Volatility Spot Models Samuel Hikspoors and Sebastian Jaimungal a a Department of Statistics and Mathematical Finance Program, University of
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationTwo and Three factor models for Spread Options Pricing
Two and Three factor models for Spread Options Pricing COMMIDITIES 2007, Birkbeck College, University of London January 17-19, 2007 Sebastian Jaimungal, Associate Director, Mathematical Finance Program,
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationSingular Perturbations in Option Pricing
Singular Perturbations in Option Pricing J.-P. Fouque G. Papanicolaou R. Sircar K. Solna March 4, 2003 Abstract After the celebrated Black-Scholes formula for pricing call options under constant volatility,
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 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 informationThe Term Structure of Interest Rates under Regime Shifts and Jumps
The Term Structure of Interest Rates under Regime Shifts and Jumps Shu Wu and Yong Zeng September 2005 Abstract This paper develops a tractable dynamic term structure models under jump-diffusion and regime
More informationHedging under Arbitrage
Hedging under Arbitrage Johannes Ruf Columbia University, Department of Statistics Modeling and Managing Financial Risks January 12, 2011 Motivation Given: a frictionless market of stocks with continuous
More informationFast narrow bounds on the value of Asian options
Fast narrow bounds on the value of Asian options G. W. P. Thompson Centre for Financial Research, Judge Institute of Management, University of Cambridge Abstract We consider the problem of finding bounds
More informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
More informationRisk Neutral Measures
CHPTER 4 Risk Neutral Measures Our aim in this section is to show how risk neutral measures can be used to price derivative securities. The key advantage is that under a risk neutral measure the discounted
More informationThe Lognormal Interest Rate Model and Eurodollar Futures
GLOBAL RESEARCH ANALYTICS The Lognormal Interest Rate Model and Eurodollar Futures CITICORP SECURITIES,INC. 399 Park Avenue New York, NY 143 Keith Weintraub Director, Analytics 1-559-97 Michael Hogan Ex
More informationSaddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance
Saddlepoint Approximation Methods for Pricing Financial Options on Discrete Realized Variance Yue Kuen KWOK Department of Mathematics Hong Kong University of Science and Technology Hong Kong * This is
More informationIntroduction to Affine Processes. Applications to Mathematical Finance
and Its Applications to Mathematical Finance Department of Mathematical Science, KAIST Workshop for Young Mathematicians in Korea, 2010 Outline Motivation 1 Motivation 2 Preliminary : Stochastic Calculus
More informationMonte Carlo Simulation of a Two-Factor Stochastic Volatility Model
Monte Carlo Simulation of a Two-Factor Stochastic Volatility Model asymptotic approximation formula for the vanilla European call option price. A class of multi-factor volatility models has been introduced
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationMulti-dimensional Term Structure Models
Multi-dimensional Term Structure Models We will focus on the affine class. But first some motivation. A generic one-dimensional model for zero-coupon yields, y(t; τ), looks like this dy(t; τ) =... dt +
More informationMartingale Approach to Pricing and Hedging
Introduction and echniques Lecture 9 in Financial Mathematics UiO-SK451 Autumn 15 eacher:s. Ortiz-Latorre Martingale Approach to Pricing and Hedging 1 Risk Neutral Pricing Assume that we are in the basic
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 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 informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More informationdt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135.
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where ( κ1 0 dx(t) = 0 κ 2 r(t) = δ 0 +X 1 (t)+x 2 (t) )( X1 (t) X 2 (t) ) ( σ1 0 dt+ ρσ 2 1 ρ2 σ 2 )( dw Q 1 (t) dw Q 2 (t) ) In this
More informationSupplementary Appendix to The Risk Premia Embedded in Index Options
Supplementary Appendix to The Risk Premia Embedded in Index Options Torben G. Andersen Nicola Fusari Viktor Todorov December 214 Contents A The Non-Linear Factor Structure of Option Surfaces 2 B Additional
More information16. Inflation-Indexed Swaps
6. Inflation-Indexed Swaps Given a set of dates T,...,T M, an Inflation-Indexed Swap (IIS) is a swap where, on each payment date, Party A pays Party B the inflation rate over a predefined period, while
More informationCredit Derivatives and Risk Aversion
Credit Derivatives and Risk Aversion Tim Leung Ronnie Sircar Thaleia Zariphopoulou October 27, revised December 27 Abstract We discuss the valuation of credit derivatives in extreme regimes such as when
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 informationHedging Credit Derivatives in Intensity Based Models
Hedging Credit Derivatives in Intensity Based Models PETER CARR Head of Quantitative Financial Research, Bloomberg LP, New York Director of the Masters Program in Math Finance, Courant Institute, NYU Stanford
More informationAsymptotic Method for Singularity in Path-Dependent Option Pricing
Asymptotic Method for Singularity in Path-Dependent Option Pricing Sang-Hyeon Park, Jeong-Hoon Kim Dept. Math. Yonsei University June 2010 Singularity in Path-Dependent June 2010 Option Pricing 1 / 21
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 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 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 informationMASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.
MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS May/June 2006 Time allowed: 2 HOURS. Examiner: Dr N.P. Byott This is a CLOSED
More informationRMSC 4005 Stochastic Calculus for Finance and Risk. 1 Exercises. (c) Let X = {X n } n=0 be a {F n }-supermartingale. Show that.
1. EXERCISES RMSC 45 Stochastic Calculus for Finance and Risk Exercises 1 Exercises 1. (a) Let X = {X n } n= be a {F n }-martingale. Show that E(X n ) = E(X ) n N (b) Let X = {X n } n= be a {F n }-submartingale.
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 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 informationPricing Asian Options with Stochastic Volatility
Pricing Asian Options with Stochastic Volatility Jean-Pierre Fouque and Chuan-Hsiang Han June 5, 23 Abstract In this paper, we generalize the recently developed dimension reduction technique of Vecer for
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationPRICING TIMER OPTIONS UNDER FAST MEAN-REVERTING STOCHASTIC VOLATILITY
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 4, Winter 009 PRICING TIMER OPTIONS UNDER FAST MEAN-REVERTING STOCHASTIC VOLATILITY DAVID SAUNDERS ABSTRACT. Timer options are derivative securities
More informationIn this lecture we will solve the final-value problem derived in the previous lecture 4, V (1) + rs = rv (t < T )
MSC FINANCIAL ENGINEERING PRICING I, AUTUMN 2010-2011 LECTURE 5: THE BLACK AND SCHOLES FORMULA AND ITS GREEKS RAYMOND BRUMMELHUIS DEPARTMENT EMS BIRKBECK In this lecture we will solve the final-value problem
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 informationSensitivity Analysis on Long-term Cash flows
Sensitivity Analysis on Long-term Cash flows Hyungbin Park Worcester Polytechnic Institute 19 March 2016 Eastern Conference on Mathematical Finance Worcester Polytechnic Institute, Worceseter, MA 1 / 49
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 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 informationDr. Maddah ENMG 625 Financial Eng g II 10/16/06
Dr. Maddah ENMG 65 Financial Eng g II 10/16/06 Chapter 11 Models of Asset Dynamics () Random Walk A random process, z, is an additive process defined over times t 0, t 1,, t k, t k+1,, such that z( t )
More informationOne-Factor Models { 1 Key features of one-factor (equilibrium) models: { All bond prices are a function of a single state variable, the short rate. {
Fixed Income Analysis Term-Structure Models in Continuous Time Multi-factor equilibrium models (general theory) The Brennan and Schwartz model Exponential-ane models Jesper Lund April 14, 1998 1 Outline
More informationSTOCHASTIC INTEGRALS
Stat 391/FinMath 346 Lecture 8 STOCHASTIC INTEGRALS X t = CONTINUOUS PROCESS θ t = PORTFOLIO: #X t HELD AT t { St : STOCK PRICE M t : MG W t : BROWNIAN MOTION DISCRETE TIME: = t < t 1
More informationLocal Volatility Dynamic Models
René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Columbia November 9, 27 Contents Joint work with Sergey Nadtochyi Motivation 1 Understanding
More informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationOn the Ross recovery under the single-factor spot rate model
.... On the Ross recovery under the single-factor spot rate model M. Kijima Tokyo Metropolitan University 11/08/2016 Kijima (TMU) Ross Recovery SMU @ August 11, 2016 1 / 35 Plan of My Talk..1 Introduction:
More informationOptimal robust bounds for variance options and asymptotically extreme models
Optimal robust bounds for variance options and asymptotically extreme models Alexander Cox 1 Jiajie Wang 2 1 University of Bath 2 Università di Roma La Sapienza Advances in Financial Mathematics, 9th January,
More informationdt + ρσ 2 1 ρ2 σ 2 B i (τ) = 1 e κ iτ κ i
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where dx(t) = ( κ1 0 0 κ 2 ) ( X1 (t) X 2 (t) In this case we find (BLACKBOARD) that r(t) = δ 0 + X 1 (t) + X 2 (t) ) ( σ1 0 dt + ρσ 2
More informationSLE and CFT. Mitsuhiro QFT2005
SLE and CFT Mitsuhiro Kato @ QFT2005 1. Introduction Critical phenomena Conformal Field Theory (CFT) Algebraic approach, field theory, BPZ(1984) Stochastic Loewner Evolution (SLE) Geometrical approach,
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 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 informationOptimal trading strategies under arbitrage
Optimal trading strategies under arbitrage Johannes Ruf Columbia University, Department of Statistics The Third Western Conference in Mathematical Finance November 14, 2009 How should an investor trade
More informationPDE Approach to Credit Derivatives
PDE Approach to Credit Derivatives Marek Rutkowski School of Mathematics and Statistics University of New South Wales Joint work with T. Bielecki, M. Jeanblanc and K. Yousiph Seminar 26 September, 2007
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 informationHedging of Credit Derivatives in Models with Totally Unexpected Default
Hedging of Credit Derivatives in Models with Totally Unexpected Default T. Bielecki, M. Jeanblanc and M. Rutkowski Carnegie Mellon University Pittsburgh, 6 February 2006 1 Based on N. Vaillant (2001) A
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 informationEstimation of dynamic term structure models
Estimation of dynamic term structure models Greg Duffee Haas School of Business, UC-Berkeley Joint with Richard Stanton, Haas School Presentation at IMA Workshop, May 2004 (full paper at http://faculty.haas.berkeley.edu/duffee)
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 informationPortfolio Optimization Under a Stressed-Beta Model
Portfolio Optimization Under a Stressed-Beta Model Jean-Pierre Fouque Adam P. Tashman 2 Abstract This paper presents a closed-form solution to the portfolio optimization problem where an agent wishes to
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 informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationApproximation Methods in Derivatives Pricing
Approximation Methods in Derivatives Pricing Minqiang Li Bloomberg LP September 24, 2013 1 / 27 Outline of the talk A brief overview of approximation methods Timer option price approximation Perpetual
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 informationA NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ
A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS Burhaneddin İZGİ Department of Mathematics, Istanbul Technical University, Istanbul, Turkey
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 information13.3 A Stochastic Production Planning Model
13.3. A Stochastic Production Planning Model 347 From (13.9), we can formally write (dx t ) = f (dt) + G (dz t ) + fgdz t dt, (13.3) dx t dt = f(dt) + Gdz t dt. (13.33) The exact meaning of these expressions
More informationExam Quantitative Finance (35V5A1)
Exam Quantitative Finance (35V5A1) Part I: Discrete-time finance Exercise 1 (20 points) a. Provide the definition of the pricing kernel k q. Relate this pricing kernel to the set of discount factors D
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 informationCredit Risk Models with Filtered Market Information
Credit Risk Models with Filtered Market Information Rüdiger Frey Universität Leipzig Bressanone, July 2007 ruediger.frey@math.uni-leipzig.de www.math.uni-leipzig.de/~frey joint with Abdel Gabih and Thorsten
More informationModel-independent bounds for Asian options
Model-independent bounds for Asian options A dynamic programming approach Alexander M. G. Cox 1 Sigrid Källblad 2 1 University of Bath 2 CMAP, École Polytechnique University of Michigan, 2nd December,
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationHigh Frequency Trading in a Regime-switching Model. Yoontae Jeon
High Frequency Trading in a Regime-switching Model by Yoontae Jeon A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Mathematics University
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
More informationReplication and Absence of Arbitrage in Non-Semimartingale Models
Replication and Absence of Arbitrage in Non-Semimartingale Models Matematiikan päivät, Tampere, 4-5. January 2006 Tommi Sottinen University of Helsinki 4.1.2006 Outline 1. The classical pricing model:
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU IFS, Chengdu, China, July 30, 2018 Peter Carr (NYU) Volatility Smiles and Yield Frowns 7/30/2018 1 / 35 Interest Rates and Volatility Practitioners and
More information