adaptive sparse grids and quasi monte carlo for option
play

Adaptive sparse grids and quasi Monte Carlo for option pricing under - PowerPoint PPT Presentation

Adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model Chiheb Ben Hammouda Christian Bayer Ra ul Tempone 3rd International Conference on Computational Finance (ICCF2019), A Coru na 8-12 July, 2019


  1. Adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model Chiheb Ben Hammouda Christian Bayer Ra´ ul Tempone 3rd International Conference on Computational Finance (ICCF2019), A Coru˜ na 8-12 July, 2019 0

  2. 1 Option Pricing under the Rough Bergomi Model: Motivation & Challenges 2 Our Hierarchical Deterministic Quadrature Methods 3 Numerical Experiments and Results 4 Conclusions 0

  3. Rough volatility 1 1 Jim Gatheral, Thibault Jaisson, and Mathieu Rosenbaum. “Volatility is rough”. In: Quantitative Finance 18.6 (2018), pp. 933–949 1

  4. The rough Bergomi model 2 This model, under a pricing measure, is given by ⎧ = √ v t S t dZ t , ⎪ ⎪ ⎪ ⎪ dS t = ξ 0 ( t ) exp ( η ̃ ⎨ t − 1 2 η 2 t 2 H ) , ⎪ W H √ (1) v t ⎪ t ≡ ρW 1 + ⎪ ⎪ ∶= ρW 1 t + ¯ 1 − ρ 2 W ⊥ , ⎩ ρW ⊥ Z t ( W 1 ,W ⊥ ) : two independent standard Brownian motions ̃ W H is Riemann-Liouville process, defined by ̃ = ∫ 0 K H ( t − s ) dW 1 t t ≥ 0 , W H s , t √ K H ( t − s ) = 2 H ( t − s ) H − 1 / 2 , ∀ 0 ≤ s ≤ t. H ∈ ( 0 , 1 / 2 ] controls the roughness of paths, ρ ∈ [ − 1 , 1 ] and η > 0. t ↦ ξ 0 ( t ) : forward variance curve, known at time 0. 2 Christian Bayer, Peter Friz, and Jim Gatheral. “Pricing under rough volatility”. In: Quantitative Finance 16.6 (2016), pp. 887–904 2

  5. Model challenges Numerically: ▸ The model is non-Markovian and non-affine ⇒ Standard numerical methods (PDEs, characteristic functions) seem inapplicable. ▸ The only prevalent pricing method for mere vanilla options is Monte Carlo (MC) (Bayer, Friz, and Gatheral 2016; McCrickerd and Pakkanen 2018) � still computationally expensive task. ▸ Discretization methods have a poor behavior of the strong error (strong convergence rate of order H ∈ [ 0 , 1 / 2 ] ) (Neuenkirch and Shalaiko 2016) ⇒ Variance reduction methods, such as multilevel Monte Carlo (MLMC), are inefficient for very small values of H . Theoretically: ▸ No proper weak error analysis done in the rough volatility context. 3

  6. Option pricing challenges The integration problem is challenging Issue 1: Time-discretization of the rough Bergomi process (large N (number of time steps)) ⇒ S takes values in a high-dimensional space ⇒ � Curse of dimensionality when using numerical integration methods. Issue 2: The payoff function g is typically not smooth ⇒ low regularity ⇒ � slow convergence of deterministic quadrature methods. � Curse of dimensionality: An exponential growth of the work (number of function evaluations) in terms of the dimension of the integration problem. 4

  7. Methodology 3 We design efficient hierarchical pricing methods based on 1 Analytic smoothing to uncover available regularity (inspired by (Romano and Touzi 1997) in the context of stochastic volatility models). 2 Approximating the option price using deterministic quadrature methods ▸ Adaptive sparse grids quadrature (ASGQ) . ▸ Quasi Monte Carlo (QMC) . 3 Coupling our methods with hierarchical representations ▸ Brownian bridges as a Wiener path generation method ⇒ ↘ the effective dimension of the problem. ▸ Richardson Extrapolation (Condition: weak error of order 1) ⇒ Faster convergence of the weak error ⇒ ↘ number of time steps (smaller dimension). 3 Christian Bayer, Chiheb Ben Hammouda, and Raul Tempone. “Hierarchical adaptive sparse grids and quasi Monte Carlo for option pricing under the rough Bergomi model”. In: arXiv preprint arXiv:1812.08533 (2018) 5

  8. Simulation of the rough Bergomi dynamics t , ̃ Goal: Simulate jointly ( W 1 ∶ 0 ≤ t ≤ T ) , resulting in W 1 W H t 1 ,...,W t N and ̃ t 1 ,..., ̃ t t N along a given grid t 1 < ⋅⋅⋅ < t N W H W H 1 Covariance based approach (Bayer, Friz, and Gatheral 2016) ▸ Based on Cholesky decomposition of the covariance matrix of the (2 N )-dimensional Gaussian random vector t N , ̃ t 1 ,..., ̃ W 1 t 1 ,...,W 1 W H W t N . ▸ Exact method but slow ▸ At least O ( N 2 ) . 2 The hybrid scheme (Bennedsen, Lunde, and Pakkanen 2017) ▸ Based on Euler discretization but crucially improved by moment matching for the singular term in the left point rule. ▸ Accurate scheme that is much faster than the Covariance based approach. ▸ O ( N ) up to logarithmic factors that depend on the desired error. 6

  9. On the choice of the simulation scheme Figure 1.1: The convergence of the weak error E B , using MC with 6 × 10 6 samples, for example parameters: H = 0 . 07 , K = 1 ,S 0 = 1 , T = 1 , ρ = − 0 . 9 , η = 1 . 9 , ξ 0 = 0 . 0552. The upper and lower bounds are 95% confidence intervals. a) With the hybrid scheme b) With the exact scheme. weak_error 10 −1 Lb Ub ∣ E [ g ( X Δ t ) − g ( X )] ∣ ∣ E [ g ( X Δ t ) − g ( X )] ∣ 10 −1 rate= 0.76 rate= 1.00 weak_error 10 −2 Lb 10 −2 Ub rate= 1.02 rate= 1.00 10 −3 10 −2 10 −1 10 −2 10 −1 Δ t Δ t (a) (b) 7

  10. 1 Option Pricing under the Rough Bergomi Model: Motivation & Challenges 2 Our Hierarchical Deterministic Quadrature Methods 3 Numerical Experiments and Results 4 Conclusions 7

  11. Conditional expectation for analytic smoothing C RB ( T,K ) = E [( S T − K ) + ] = E [ E [( S T − K ) + ∣ σ ( W 1 ( t ) ,t ≤ T )]] √ v t dW 1 2 ρ 2 ∫ = E [ C BS ( S 0 = exp ( ρ ∫ t − 1 v t dt ) , T T 0 0 k = K, σ 2 = ( 1 − ρ 2 ) ∫ v t dt )] T 0 ≈ ∫ R 2 N C BS ( G ( w ( 1 ) , w ( 2 ) )) ρ N ( w ( 1 ) ) ρ N ( w ( 2 ) ) d w ( 1 ) d w ( 2 ) = C N RB . (2) C BS ( S 0 ,k,σ 2 ) : the Black-Scholes call price, for initial spot price S 0 , strike price k , and volatility σ 2 . G maps 2 N independent standard Gaussian random inputs to the parameters fed to Black-Scholes formula. ρ N : the multivariate Gaussian density, N : number of time steps. 8

  12. Sparse grids I Notation: Given F ∶ R d → R and a multi-index β ∈ N d + . F β ∶ = Q m ( β ) [ F ] a quadrature operator based on a Cartesian quadrature grid ( m ( β n ) points along y n ). � Approximating E [ F ] with F β is not an appropriate option due to the well-known curse of dimensionality. The first-order difference operators ∆ i F β { F β − F β − e i , if β i > 1 if β i = 1 F β where e i denotes the i th d -dimensional unit vector The mixed (first-order tensor) difference operators ∆ [ F β ] = ⊗ d i = 1 ∆ i F β Idea: A quadrature estimate of E [ F ] is M I ℓ [ F ] = ∑ ∆ [ F β ] , (3) β ∈I ℓ 9

  13. Sparse grids II E [ F ] ≈ M I ℓ [ F ] = ∑ ∆ [ F β ] , β ∈I ℓ Product approach: I ℓ = {∣ β ∣ ∞ ≤ ℓ ; β ∈ N d + } Regular sparse grids 4 : I ℓ = {∣ β ∣ 1 ≤ ℓ + d − 1; β ∈ N d + } Adaptive sparse grids quadrature (ASGQ): I ℓ = I ASGQ (Next slides). Figure 2.1: Left are product grids ∆ β 1 ⊗ ∆ β 2 for 1 ≤ β 1 ,β 2 ≤ 3. Right is the corresponding SG construction. 4 Hans-Joachim Bungartz and Michael Griebel. “Sparse grids”. In: Acta numerica 13 (2004), pp. 147–269 10

  14. ASGQ in practice The construction of I ASGQ is done by profit thresholding I ASGQ = { β ∈ N d + ∶ P β ≥ T } . Profit of a hierarchical surplus P β = ∣ ∆ E β ∣ ∆ W β . Error contribution : ∆ E β = ∣M I∪{ β } − M I ∣ . Work contribution : ∆ W β = Work [M I∪{ β } ] − Work [M I ] . Figure 2.2: A posteriori, adaptive construction as in (Haji-Ali et al. 2016): Given an index set I k , compute the profits of the neighbor indices and select the most profitable one 11

  15. ASGQ in practice The construction of I ASGQ is done by profit thresholding I ASGQ = { β ∈ N d + ∶ P β ≥ T } . Profit of a hierarchical surplus P β = ∣ ∆ E β ∣ ∆ W β . Error contribution : ∆ E β = ∣M I∪{ β } − M I ∣ . Work contribution : ∆ W β = Work [M I∪{ β } ] − Work [M I ] . Figure 2.3: A posteriori, adaptive construction as in (Haji-Ali et al. 2016): Given an index set I k , compute the profits of the neighbor indices and select the most profitable one 11

  16. ASGQ in practice The construction of I ASGQ is done by profit thresholding I ASGQ = { β ∈ N d + ∶ P β ≥ T } . Profit of a hierarchical surplus P β = ∣ ∆ E β ∣ ∆ W β . Error contribution : ∆ E β = ∣M I∪{ β } − M I ∣ . Work contribution : ∆ W β = Work [M I∪{ β } ] − Work [M I ] . Figure 2.4: A posteriori, adaptive construction as in (Haji-Ali et al. 2016): Given an index set I k , compute the profits of the neighbor indices and select the most profitable one 11

  17. ASGQ in practice The construction of I ASGQ is done by profit thresholding I ASGQ = { β ∈ N d + ∶ P β ≥ T } . Profit of a hierarchical surplus P β = ∣ ∆ E β ∣ ∆ W β . Error contribution : ∆ E β = ∣M I∪{ β } − M I ∣ . Work contribution : ∆ W β = Work [M I∪{ β } ] − Work [M I ] . Figure 2.5: A posteriori, adaptive construction as in (Haji-Ali et al. 2016): Given an index set I k , compute the profits of the neighbor indices and select the most profitable one 11

  18. ASGQ in practice The construction of I ASGQ is done by profit thresholding I ASGQ = { β ∈ N d + ∶ P β ≥ T } . Profit of a hierarchical surplus P β = ∣ ∆ E β ∣ ∆ W β . Error contribution : ∆ E β = ∣M I∪{ β } − M I ∣ . Work contribution : ∆ W β = Work [M I∪{ β } ] − Work [M I ] . Figure 2.6: A posteriori, adaptive construction as in (Haji-Ali et al. 2016): Given an index set I k , compute the profits of the neighbor indices and select the most profitable one 11

Recommend


More recommend