The data-driven COS method
|
|
- Eustace Washington
- 5 years ago
- Views:
Transcription
1 The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica Reading group, March 13, 2017 Reading group, March 13, /
2 1 The COS method 2 Learning densities 3 The data-driven COS (ddcos) method 4 Choice of parameters in ddcos method 5 Applications of the ddcos method 6 Conclusions Reading group, March 13, /
3 The COS method A lot of work behind: [FO08], [FO09], etc. Fourier-based method to price options. Starting point is risk-neutral valuation formula: v(x, t) = e r(t t) E [v(y, T ) x] = e r(t t) R v(y, T )f (y x)dy, where r is the risk-free rate and f (y x) is the density of the underlying process. Typically, we have: ( ) ( ) S(0) S(T ) x := log and y := log, K K f (y x) is unknown in most of cases. However, characteristic function available for many models. Exploit the relation between the density and the characteristic function (Fourier pair). Reading group, March 13, /
4 The COS method - European options f (y x) is approximated, on a finite interval [a, b], by a cosine series ( f (y x) = 1 ( A A k (x) cos kπ y a ) ), b a b a k=1 b ( A 0 = 1, A k (x) = f (y x) cos kπ y a ) dy, k = 1, 2,.... b a a Interchanging the summation and integration and introducing the definition V k := 2 b ( v(y, T ) cos kπ y a ) dy, b a a b a we find that the option value is given by v(x, t) e r(t t) A k (x)v k, k=0 where indicates that the first term is divided by two. Reading group, March 13, /
5 Pricing European options with the COS method Coefficients A k can be computed from the ChF. Coefficients V k are known analytically (for many types of options). Closed-form expressions for the option Greeks and Γ = v(x, t) S Γ = 2 v(x, t) S 2 = 1 v(x, t) S(0) x = exp( r(t t)) exp( r(t t)) k=0 A k (x) x V k S(0), k=0 ( A k(x) + 2 A k (x) x) x 2 ) ) Vk S 2 (0) Due to the rapid decay of the coefficients, v(x, t), and Γ can be approximated with high accuracy by truncating to N terms. Reading group, March 13, /
6 Learning densities Statistical learning theory: deals with the problem of finding a predictive function based on data. We follow the analysis about the problem of density estimation proposed by Vapnik in [Vap98]. Given independent and identical distributed samples X 1, X 2,..., X n. By definition, density f (x) is related to the cumulative distribution function, F (x), by means of the expression x f (y)dy = F (x). Function F (x) is approximated by the empirical approximation F n (x) = 1 n n η(x X i ), i=1 where η( ) is the step-function. Convergence O(1/ n). Reading group, March 13, /
7 Regularization approach The previous equation can be rewritten as a linear operator equation Cf = F F n, where the operator Ch := x h(z)dz. Stochastic ill-posed problem. Regularization method (Vapnik). Given a lower semi-continuous functional W (f ) such that: Solution of Cf = Fn belongs to D, the domain of definition of W (f ). The functional W (f ) takes real non-negative values in D. The set M c = {f : W (f ) c} is compact in H (the space where the solution exits and is unique). Then we can construct the functional R γn (f, F n ) = L 2 H(Cf, F n ) + γ n W (f ), where L H is a metric of the space H (loss function) and γ n is the parameter of regularization satisfying that γ n 0 as n. Under these conditions, a function f n minimizing the functional converges almost surely to the desired one. Reading group, March 13, /
8 Regularization and Fourier-based density estimators Assume f (x) belongs to the functions whose p-th derivatives belong to L 2 (0, π), the kernel K(z x) and ( 2 W (f ) = K(z x)f (x)dx) dz, R R The risk functional ( x ) 2 ( 2 R γn (f, F n ) = f (y)dy F n (x) dx+γ n K(z x)f (x)dx) dz. R 0 R R Denoting by ˆf (u), ˆF n (u) and ˆK(u) the Fourier transforms, by definition ˆF n (u) = 1 F n (x)e iux dx 2π R = 1 n η(x X j )e iux dx = 1 n exp( iux j ), 2nπ n iu R j=1 where i = 1 is the imaginary unit. Reading group, March 13, / j=1
9 Regularization and Fourier-based density estimators By employing the convolution theorem and Parseval s identity ˆf (u) 1 n n j=1 R γn (f, F n ) = exp( iux j) 2 + γ n iu ˆK(u)ˆf (u) 2. L 2 L 2 The condition to minimize R γn (f, F n ) is given by, ˆf (u) u 2 1 nu 2 which gives us, n exp( iux j ) + γ n ˆK(u) ˆK( u)ˆf (u) = 0, j=1 ( ) 1 1 ˆf n (u) = 1 + γ n u 2 ˆK(u) ˆK( u) n n exp( iux j ). j=1 Reading group, March 13, /
10 Regularization and Fourier-based density estimators K(x) = δ (p) (x), and the desired PDF, f (x) and its p-th derivative (p 0) belongs to L 2 (0, π), the risk functional becomes π ( x ) 2 π ( 2 R γn (f, F n ) = f (y)dy F n (x) dx + γ n f (x)) (p) dx. 0 0 Given orthonormal functions, ψ 1 (θ),..., ψ k (θ),... f n (θ) = 1 π + 2 à k ψ k (θ), π k=1 with à 0, à 1,..., à k,... expansion coefficients, à k =< f n, ψ k >. The coefficients à k cannot be directly computed from f n, but à k =< f n, ψ k >=< ˆf n, ˆψ k > π ( ) = γ n u 2 ˆK(u) ˆK( u) n 0 0 n exp( iuθ j ) ˆψ k (u)du. j=1 Reading group, March 13, /
11 Regularization and Fourier-based density estimators Using cosine series expansions, i.e., ψ k (θ) = cos(kθ), it is well-known that ˆψ k (u) = 1 (δ(u k) + δ(u + k)). 2 This facilitates the computation of à k avoiding the calculation of the integral. Thus, the minimum of R γn à k = 1 ( ) n 1 2n 1 + γ n ( k) 2 ˆK( k) exp(ikθ j ) ˆK(k) j=1 ( ) n γ n k 2 ˆK(k) exp( ikθ j ) ˆK( k) 1 1 = 1 + γ n k 2 ˆK(k) ˆK( k) n j=1 j=1 n 1 1 cos(kθ j ) = 1 + γ n k 2(p+1) n n cos(kθ j ), where θ j (0, π) are given samples of the unknown distribution. In the last step, ˆK(u) = (iu) p is used. j=1 Reading group, March 13, /
12 The data-driven COS method Employ the solution of the regularization problem for density estimation in the COS framework. In both, the density is assumed to be in the form of a cosine series expansion. The minimum of the functional is in terms of the expansion coefficients. Take advantage of the COS machinery: pricing options, Greeks, etc. The samples must follow risk-neutral measure (Monte Carlo paths). Reading group, March 13, /
13 The data-driven COS method Key idea: Ã k approximates A k. Risk neutral samples from an asset at time T, S 1 (t), S 2 (t),..., S n (t). With a logarithmic transformation, we have ( ) Sj (T ) Y j := log. K The regularization solution is defined in (0, π), by transformation θ j = π Y j a b a, where the boundaries a and b are defined as a := min 1 j n (Y j), b := max 1 j n (Y j). Reading group, March 13, /
14 The data-driven COS method - European options The A k coefficients are replaced by the data-driven Ãk 1 ) n n j=1 (kπ cos Y j a b a A k à k = 1 + γ n k 2(p+1). The ddcos pricing formula for European options ṽ(x, t) = e r(t t) k=0 1 n = e r(t t) à k V k. k=0 ) n j=1 (kπ cos Y j a b a 1 + γ n k 2(p+1) V k As in the original COS method, we must truncate the infinite sum to a finite number of terms N N ṽ(x, t) = e r(t t) à k V k, k=0 Reading group, March 13, /
15 The data-driven COS method - Greeks Data-driven expressions for the and Γ sensitivities. Define the corresponding sine coefficients as B k := 1 n ) n j=1 (kπ sin Y j a b a 1 + γ n k 2(p+1). Taking derivatives of the ddcos pricing formulat w.r.t the samples, Y j, the data-driven Greeks, and Γ, can be obtained by = e r(t t) N k=0 N ( Γ = e r(t t) B k k=0 ( B k kπ b a ) V k S(0), ( ) ) kπ kπ 2 b a Ãk b a V k S 2 (0). Reading group, March 13, /
16 The data-driven COS method - Variance reduction admits in the computation of Ãk. Here, antithetic variates (AV) to our method. Since the samples must be i.i.d., an immediate application of AV is not possible. Assume antithetic samples, Y, that can be computed without computational effort, a new estimator i Ā k := 1 2 (Ãk + à k where à k are antithetic coefficients, obtained from Y i. It can be proved that the use of Āk results in a variance reduction. Additional information to reduce the variance. For example, the martingale property S(T ) = S(T ) 1 n S j (T ) + E[S(T )], n = S(T ) 1 n j=1 ), n S j (T ) + S(0) exp(rt ). j=1 Reading group, March 13, /
17 Choice of parameters in ddcos method The choice of optimal values of γ n and p. There is no rule or procedure to obtain an optimal p. As a rule of thumb, p = 0 seems to be the most appropriate value. Fixing p, we rely on the computation of an optimal γ n True p = 0 p = 1 p = 2 p = Figure: Parameter p analysis. Reading group, March 13, /
18 Choice of γ n γ n impacts the efficiency of the ddcos method: it is related to the number of samples, n, and number of terms, N. For the regularization parameter γ n, a rule that ensures asymptotic convergence log log n γ n =. n In practical situations: not optimal. Exploit the relation between the empirical and real (unknown) CDFs. Reading group, March 13, /
19 Choice of γ n This relation can be modeled by statistical laws or statistics: Kolmogorov-Smirnov, Anderson-Darling, Smirnov-Cramér von Mises. Preferable: a measure of the distance between the F n (x) and F (x) follows a known distribution. We have chosen Smirnov-Cramér von Mises(SCvM): ω 2 = n (F (x) F n (x)) 2 df (x). i=1 R Assume we have an approximation, F γn (which depends on γ n ). An almost optimal γ n is computed by solving the equation n ( F γn ( X i ) i 0.5 ) 2 = m S i n 12n, where X 1, X 2,..., X n is the ordered array of samples X 1, X 2,..., X n and m S the mean of the ω 2. Reading group, March 13, /
20 Influence of γ n To assess the impact of γ n : Mean integrated Squared Error (MiSE): [ ] [ ] E f n f 2 2 = E (f n (x) f (x)) 2 dx. R A formula for the MiSE formula is derived in our context: MISE = 1 n N k=1 ( 1 1 ( 1 + γn k 2(p+1)) ) 2 A 2k A 2 k + k=n+1 Two main aspects influenced γ n : accuracy in n and stability in N. The quality of the approximated density can be also affected. A 2 k. Reading group, March 13, /
21 Influence of γ n 10 2 n "rule" 10 0 A k real n SCvM n = A k rule A k rule, n =0 A k SCvM A k SCvM, n = (a) Convergence in terms of n (b) Convergence in terms of N Figure: Influence of γ n :. Reading group, March 13, /
22 Optimal number of terms N Try to find a minimum optimal value of N. N considerably affects the performance. We wish to avoid the computation of any  k. We define a proxy for the MiSE and follow: MiSE 1 n N k=1 1 2 ( 1 + γn k 2(p+1)) n rule rule - addend 1 n rule - proxy n n SCvM SCvM - addend 1 n SCvM - proxy n Reading group, March 13, /
23 Optimal number of terms N Data: n, γ n N min = 5 N max = ɛ = 1 n MiSE prev = for N = N min : N max do MiSE N = 1 n N k=1 ɛ N = MiSE N MiSE prev MiSE N if ɛ N > ɛ then N op = N else Break 1 2 (1+γ nk 2(p+1) ) 2 N n Figure: Almost optimal N. MiSE prev = MiSE N Reading group, March 13, /
24 Applications of the ddcos method Pricing options (no better than Monte Carlo). Sensitivities or Greeks. Models without analytic characteristic function. SABR model. Risk measures: VaR and Expected shortfall. Combinations. Reading group, March 13, /
25 Applications of the ddcos method Unfortunately, the γ n based on the SCvM statistic does not provide any benefit. The use of the γ n rule entails faster ddcos estimators. ddcos converges with the expected convergence rate O(1/ n). The variance reduction techniques are successfully applied. In Greeks computation, Monte Carlo-based methods may require one or two extra simulations. In the convergence tests, the reported values are computed as the average of 50 experiments. Reading group, March 13, /
26 Applications of the ddcos - Option pricing 10 1 ddcos ddcos, AV MC 10 0 MC, AV 10 1 ddcos ddcos, AV MC 10 0 MC, AV (a) Call: Strike K = 100. (b) Put: Strike K = 100. Figure: Convergence in prices of the ddcos method: Antithetic Variates (AV); GBM, S(0) = 100, r = 0.1, σ = 0.3 and T = 2. Reading group, March 13, /
27 Applications of the ddcos - Greeks estimation 10-1 ddcos ddcos MCFD 10-2 MCFD, AV, AV 10-2 ddcos ddcos, AV MCFD MCFD, AV (a) (Call): Strike K = 100. (b) Γ: Strike K = 100. Figure: Convergence in Greeks of the ddcos method: Antithetic Variates (AV); GBM, S(0) = 100, r = 0.1, σ = 0.3 and T = 2. Reading group, March 13, /
28 Applications of the ddcos - Greeks estimation K (% of S(0)) 80% 90% 100% 110% 120% 0.1 Ref ddcos RE MCFD RE Γ Ref ddcos RE MCFD RE Table: GBM call option Greeks: S(0) = 100, r = 0.1, σ = 0.3 and T = 2. Reading group, March 13, /
29 Applications of the ddcos - Greeks estimation K (% of S(0)) 80% 90% 100% 110% 120% Ref ddcos RE MCFD RE Γ Ref ddcos RE MCFD RE Table: Merton jump-diffusion call option Greeks: S(0) = 100, r = 0.1, σ = 0.3, µ j = 0.2, σ j = 0.2 and λ = 8 and T = 2. Reading group, March 13, /
30 Applications of the ddcos - Greeks estimation K (% of S(0)) 80% 90% 100% 110% 120% Ref ddcos RE MCFD RE Table: Call option Greek under the SABR model: S(0) = 100, r = 0, σ 0 = 0.3, α = 0.4, β = 0.6, ρ = 0.25 and T = 2. Reading group, March 13, /
31 Applications of the ddcos - Greeks estimation K (% of S(0)) 80% 90% 100% 110% 120% Ref ddcos RE Hagan RE Table: under SABR model. Setting: Call, S(0) = 0.04, r = 0.0, σ 0 = 0.4, α = 0.8, β = 1.0, ρ = 0.5 and T = 2. Reading group, March 13, /
32 Applications of the ddcos - Greeks estimation 0.75 ddcos Ref ddcos (a) : Strike K = (b) Γ: Strike K = Figure: : Greeks convergence test. Reading group, March 13, /
33 Applications of the ddcos - Risk measures In the context of the Delta-Gamma approach (COS in [OGO14]). The change in a portfolio value can be generalized. L := V = V (S, t) V (S + S, t + t). The formal definition of the VaR reads P( V < VaR(q)) = 1 F L (VaR(q)) = q, with q a predefined confidence level. Given the VaR, the ES measure is computed as ES := E[ V V > VaR(q)]. Two portfolios with the same composition: one European call and half a European put on the same asset, maturity 60 days and K = 101. Different time horizons: 1 day (Portfolio 1) and 10 days (Portfolio 2). The asset follows a GBM with S(0) = 100, r = 0.1 and σ = 0.3. Reading group, March 13, /
34 Applications of the ddcos - Risk measures ddcos COS ddcos COS (a) Density Portfolio 1. (b) Density Portfolio 2. Figure: Recovered densities of L: ddcos vs. COS. Reading group, March 13, /
35 Applications of the ddcos - Risk measures 10 0 VaR ES VaR ES (a) Portfolio 1: q = 99%. (b) Portfolio 2: q = 90%. Figure: VaR and ES convergence in n. Reading group, March 13, /
36 Applications of the ddcos - Risk measures The oscillations can be removed. Two options: smoothing parameter or filters [RVO14] COS ddcos, p=1 ddcos, filter COS ddcos, p=1 ddcos, filter (a) Density Portfolio 1. (b) Density Portfolio 2. Figure: Smoothed densities of L. Reading group, March 13, /
37 Applications of the ddcos - Risk measures and SABR ddcos VaR ddcos ES ddcos f L ddcos F L ddcos f L, filter ddcos F L, filter (a) VaR and ES: q = 99%. (b) F L and f L. Figure: Delta-Gamma approach under the SABR model. Setting: S(0) = 100, K = 100, r = 0.0, σ 0 = 0.4, α = 0.8, β = 1.0, ρ = 0.5, T = 2, q = 99% and t = 1/365. Reading group, March 13, /
38 Applications of the ddcos - Risk measures and SABR q 10% 30% 50% 70% 90% VaR ES Table: VaR and ES under SABR model. Setting: S(0) = 100, K = 100, r = 0.0, σ 0 = 0.4, α = 0.8, β = 1.0, ρ = 0.5, T = 2, and t = 1/365. Reading group, March 13, /
39 Conclusions extends the COS method applicability to cases when only data samples of the underlying are available. The method exploits a closed-form solution, in terms of Fourier cosine expansions, of a regularization problem. It allows to develop a data-driven method which can be employed for option pricing and risk management. particularly results in an efficient method for the and Γ sensitivities computation, based solely on the samples. It can be employed within the Delta-Gamma approximation for calculating risk measures. A possible future extension may be the use of other basis functions. Haar wavelets are for example interesting since they provide positive densities and allow an efficient treatment of dynamic data. Reading group, March 13, /
40 References Fang Fang and Cornelis W. Oosterlee. A novel pricing method for European options based on Fourier-cosine series expansions. SIAM Journal on Scientific Computing, 31: , Fang Fang and Cornelis W. Oosterlee. Pricing early-exercise and discrete barrier options by Fourier-cosine series expansions. Numerische Mathematik, 114(1):27 62, Luis Ortiz-Gracia and Cornelis W. Oosterlee. Efficient VaR and Expected Shortfall computations for nonlinear portfolios within the delta-gamma approach. Applied Mathematics and Computation, 244:16 31, Maria J. Ruijter, Mark Versteegh, and Cornelis W. Oosterlee. On the application of spectral filters in a Fourier option pricing technique. Journal of Computational Finance, 19(1):75 106, Vladimir N. Vapnik. Statistical Learning Theory. Wiley-Interscience, Reading group, March 13, /
41 Suggestions, comments & questions Thank you for your attention Reading group, March 13, 2017 /
The data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica CMMSE 2017, July 6, 2017 Álvaro Leitao (CWI & TUDelft)
More informationA Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options
A Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options Luis Ortiz-Gracia Centre de Recerca Matemàtica (joint work with Cornelis W. Oosterlee, CWI) Models and Numerics
More informationEfficient Concentration Risk Measurement in Credit Portfolios with Haar Wavelets
Efficient Concentration Risk Measurement in Credit Portfolios with Haar Wavelets Josep J. Masdemont 1 and Luis Ortiz-Gracia 2 1 Universitat Politècnica de Catalunya 2 Centre de Recerca Matemàtica & Centrum
More informationValuation of Equity / FX Instruments
Technical Paper: Valuation of Equity / FX Instruments MathConsult GmbH Altenberger Straße 69 A-4040 Linz, Austria 14 th October, 2009 1 Vanilla Equity Option 1.1 Introduction A vanilla equity option is
More informationPricing European Options by Stable Fourier-Cosine Series Expansions
Pricing European Options by Stable Fourier-Cosine Series Expansions arxiv:171.886v2 [q-fin.cp] 8 Jan 217 Chunfa Wang Fiance School of Zhejiang University of Finance and Economics, Hangzhou, China, cfwang@zufe.edu.cn
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 informationAnalytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model
Analytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model Advances in Computational Economics and Finance Univerity of Zürich, Switzerland Matthias Thul 1 Ally Quan
More informationTime-changed Brownian motion and option pricing
Time-changed Brownian motion and option pricing Peter Hieber Chair of Mathematical Finance, TU Munich 6th AMaMeF Warsaw, June 13th 2013 Partially joint with Marcos Escobar (RU Toronto), Matthias Scherer
More informationVaR Estimation under Stochastic Volatility Models
VaR Estimation under Stochastic Volatility Models Chuan-Hsiang Han Dept. of Quantitative Finance Natl. Tsing-Hua University TMS Meeting, Chia-Yi (Joint work with Wei-Han Liu) December 5, 2009 Outline Risk
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
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 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 informationValue at Risk Ch.12. PAK Study Manual
Value at Risk Ch.12 Related Learning Objectives 3a) Apply and construct risk metrics to quantify major types of risk exposure such as market risk, credit risk, liquidity risk, regulatory risk etc., and
More informationA New Hybrid Estimation Method for the Generalized Pareto Distribution
A New Hybrid Estimation Method for the Generalized Pareto Distribution Chunlin Wang Department of Mathematics and Statistics University of Calgary May 18, 2011 A New Hybrid Estimation Method for the GPD
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 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 informationStochastic Grid Bundling Method
Stochastic Grid Bundling Method GPU Acceleration Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee London - December 17, 2015 A. Leitao &
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 informationAnalysis of truncated data with application to the operational risk estimation
Analysis of truncated data with application to the operational risk estimation Petr Volf 1 Abstract. Researchers interested in the estimation of operational risk often face problems arising from the structure
More informationFast Convergence of Regress-later Series Estimators
Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser
More information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
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 informationStrategies for Improving the Efficiency of Monte-Carlo Methods
Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful
More informationSTOCHASTIC VOLATILITY AND OPTION PRICING
STOCHASTIC VOLATILITY AND OPTION PRICING Daniel Dufresne Centre for Actuarial Studies University of Melbourne November 29 (To appear in Risks and Rewards, the Society of Actuaries Investment Section Newsletter)
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 information2 f. f t S 2. Delta measures the sensitivityof the portfolio value to changes in the price of the underlying
Sensitivity analysis Simulating the Greeks Meet the Greeks he value of a derivative on a single underlying asset depends upon the current asset price S and its volatility Σ, the risk-free interest rate
More informationEfficient Valuation of Equity-Indexed Annuities Under Lévy Processes Using Fourier-Cosine Series
Efficient Valuation of Equity-Indexed Annuities Under Lévy Processes Using Fourier-Cosine Series Geng Deng, PhD, CFA, FRM Tim Dulaney, PhD, FRM, Craig McCann, PhD, CFA Mike Yan, PhD, FRM April 1, 14 Abstract
More informationMachine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
More informationFinancial Time Series and Their Characteristics
Financial Time Series and Their Characteristics Egon Zakrajšek Division of Monetary Affairs Federal Reserve Board Summer School in Financial Mathematics Faculty of Mathematics & Physics University of Ljubljana
More informationFourier, Wavelet and Monte Carlo Methods in Computational Finance
Fourier, Wavelet and Monte Carlo Methods in Computational Finance Kees Oosterlee 1,2 1 CWI, Amsterdam 2 Delft University of Technology, the Netherlands AANMPDE-9-16, 7/7/26 Kees Oosterlee (CWI, TU Delft)
More informationLecture 1: Lévy processes
Lecture 1: Lévy processes A. E. Kyprianou Department of Mathematical Sciences, University of Bath 1/ 22 Lévy processes 2/ 22 Lévy processes A process X = {X t : t 0} defined on a probability space (Ω,
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
More informationCalculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the
VaR Pro and Contra Pro: Easy to calculate and to understand. It is a common language of communication within the organizations as well as outside (e.g. regulators, auditors, shareholders). It is not really
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 informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
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 informationPreprint núm January Robust pricing of european options with wavelets and the characteristic function. L. Ortiz-Gracia, C. W.
Preprint núm. 1144 January 213 Robust pricing of european options with wavelets and the characteristic function L. Ortiz-Gracia, C. W. Oosterlee ROBUST PRICING OF EUROPEAN OPTIONS WITH WAVELETS AND THE
More informationROBUST PRICING OF EUROPEAN OPTIONS WITH WAVELETS AND THE CHARACTERISTIC FUNCTION
ROBUST PRICING OF EUROPEAN OPTIONS WITH WAVELETS AND THE CHARACTERISTIC FUNCTION LUIS ORTIZ-GRACIA AND CORNELIS W. OOSTERLEE Abstract. We present a novel method for pricing European options based on the
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 informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationIMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS
IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
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 informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationFINITE DIFFERENCE METHODS
FINITE DIFFERENCE METHODS School of Mathematics 2013 OUTLINE Review 1 REVIEW Last time Today s Lecture OUTLINE Review 1 REVIEW Last time Today s Lecture 2 DISCRETISING THE PROBLEM Finite-difference approximations
More informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
More informationPricing Bermudan options in Lévy process models
1 1 Dept. of Industrial & Enterprise Systems Engineering University of Illinois at Urbana-Champaign Joint with Xiong Lin Bachelier Finance Society 6th Congress 6/24/2010 American options in Lévy models
More informationNormal Inverse Gaussian (NIG) Process
With Applications in Mathematical Finance The Mathematical and Computational Finance Laboratory - Lunch at the Lab March 26, 2009 1 Limitations of Gaussian Driven Processes Background and Definition IG
More informationAsian options and meropmorphic Lévy processes
Asian options and meropmorphic Lévy processes July 9 13, 2013 Outline Introduction Existing pricing approaches and the exponential functional I q Methodology overview Meromorphic Lévy processes Theoretical
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 informationFrom Characteristic Functions and Fourier Transforms to PDFs/CDFs and Option Prices
From Characteristic Functions and Fourier Transforms to PDFs/CDFs and Option Prices Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Fourier Transforms Option Pricing, Fall, 2007
More informationFast Pricing and Calculation of Sensitivities of OTM European Options Under Lévy Processes
Fast Pricing and Calculation of Sensitivities of OTM European Options Under Lévy Processes Sergei Levendorskĭi Jiayao Xie Department of Mathematics University of Leicester Toronto, June 24, 2010 Levendorskĭi
More informationA Continuity Correction under Jump-Diffusion Models with Applications in Finance
A Continuity Correction under Jump-Diffusion Models with Applications in Finance Cheng-Der Fuh 1, Sheng-Feng Luo 2 and Ju-Fang Yen 3 1 Institute of Statistical Science, Academia Sinica, and Graduate Institute
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 10-03 Acceleration of Option Pricing Technique on Graphics Processing Units Bowen Zhang Cornelis. W. Oosterlee ISSN 1389-6520 Reports of the Department of Applied
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 informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -5 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Moments of a distribubon Measures of
More informationPricing of minimum interest guarantees: Is the arbitrage free price fair?
Pricing of minimum interest guarantees: Is the arbitrage free price fair? Pål Lillevold and Dag Svege 17. 10. 2002 Pricing of minimum interest guarantees: Is the arbitrage free price fair? 1 1 Outline
More informationThe Evaluation Of Barrier Option Prices Under Stochastic Volatility
QUANTITATIVE FINANCE RESEARCH CENTRE QUANTITATIVE F INANCE RESEARCH CENTRE QUANTITATIVE FINANCE RESEARCH CENTRE Research Paper 266 January 21 The Evaluation Of Barrier Option Prices Under Stochastic Volatility
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More informationChapter 2 Uncertainty Analysis and Sampling Techniques
Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying
More 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 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 informationCredit Risk and Underlying Asset Risk *
Seoul Journal of Business Volume 4, Number (December 018) Credit Risk and Underlying Asset Risk * JONG-RYONG LEE **1) Kangwon National University Gangwondo, Korea Abstract This paper develops the credit
More informationPricing and Hedging of Commodity Derivatives using the Fast Fourier Transform
Pricing and Hedging of Commodity Derivatives using the Fast Fourier Transform Vladimir Surkov vladimir.surkov@utoronto.ca Department of Statistical and Actuarial Sciences, University of Western Ontario
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
More informationImplied Lévy Volatility
Joint work with José Manuel Corcuera, Peter Leoni and Wim Schoutens July 15, 2009 - Eurandom 1 2 The Black-Scholes model The Lévy models 3 4 5 6 7 Delta Hedging at versus at Implied Black-Scholes Volatility
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Risk Measures Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Reference: Chapter 8
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More informationOptimized Least-squares Monte Carlo (OLSM) for Measuring Counterparty Credit Exposure of American-style Options
Optimized Least-squares Monte Carlo (OLSM) for Measuring Counterparty Credit Exposure of American-style Options Kin Hung (Felix) Kan 1 Greg Frank 3 Victor Mozgin 3 Mark Reesor 2 1 Department of Applied
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 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 informationSparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models
Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Norbert Hilber Seminar of Applied Mathematics ETH Zürich Workshop on Financial Modeling with Jump Processes p. 1/18 Outline
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 informationMultilevel Monte Carlo for VaR
Multilevel Monte Carlo for VaR Mike Giles, Wenhui Gou, Abdul-Lateef Haji-Ali Mathematical Institute, University of Oxford (BNP Paribas, Hong Kong) (also discussions with Ralf Korn, Klaus Ritter) Advances
More informationOption Pricing for Discrete Hedging and Non-Gaussian Processes
Option Pricing for Discrete Hedging and Non-Gaussian Processes Kellogg College University of Oxford A thesis submitted in partial fulfillment of the requirements for the MSc in Mathematical Finance November
More informationSADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD
The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint
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 informationGRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS
GRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS Patrick GAGLIARDINI and Christian GOURIÉROUX INTRODUCTION Risk measures such as Value-at-Risk (VaR) Expected
More informationMSc in Financial Engineering
Department of Economics, Mathematics and Statistics MSc in Financial Engineering On Numerical Methods for the Pricing of Commodity Spread Options Damien Deville September 11, 2009 Supervisor: Dr. Steve
More informationINTEREST RATES AND FX MODELS
INTEREST RATES AND FX MODELS 3. The Volatility Cube Andrew Lesniewski Courant Institute of Mathematics New York University New York February 17, 2011 2 Interest Rates & FX Models Contents 1 Dynamics of
More informationA Hybrid Importance Sampling Algorithm for VaR
A Hybrid Importance Sampling Algorithm for VaR No Author Given No Institute Given Abstract. Value at Risk (VaR) provides a number that measures the risk of a financial portfolio under significant loss.
More informationSYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data
SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu September 5, 2015
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 informationTechnische Universiteit Delft Faculteit Elektrotechniek, Wiskunde en Informatica Delft Institute of Applied Mathematics
Technische Universiteit Delft Faculteit Elektrotechniek, Wiskunde en Informatica Delft Institute of Applied Mathematics Het nauwkeurig bepalen van de verlieskans van een portfolio van risicovolle leningen
More informationAsymptotic results discrete time martingales and stochastic algorithms
Asymptotic results discrete time martingales and stochastic algorithms Bernard Bercu Bordeaux University, France IFCAM Summer School Bangalore, India, July 2015 Bernard Bercu Asymptotic results for discrete
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
More informationOption Pricing and Calibration with Time-changed Lévy processes
Option Pricing and Calibration with Time-changed Lévy processes Yan Wang and Kevin Zhang Warwick Business School 12th Feb. 2013 Objectives 1. How to find a perfect model that captures essential features
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 informationIdiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective
Idiosyncratic risk, insurance, and aggregate consumption dynamics: a likelihood perspective Alisdair McKay Boston University June 2013 Microeconomic evidence on insurance - Consumption responds to idiosyncratic
More informationDistributed Computing in Finance: Case Model Calibration
Distributed Computing in Finance: Case Model Calibration Global Derivatives Trading & Risk Management 19 May 2010 Techila Technologies, Tampere University of Technology juho.kanniainen@techila.fi juho.kanniainen@tut.fi
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 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 informationOverview. Transformation method Rejection method. Monte Carlo vs ordinary methods. 1 Random numbers. 2 Monte Carlo integration.
Overview 1 Random numbers Transformation method Rejection method 2 Monte Carlo integration Monte Carlo vs ordinary methods 3 Summary Transformation method Suppose X has probability distribution p X (x),
More 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 informationFinancial Times Series. Lecture 6
Financial Times Series Lecture 6 Extensions of the GARCH There are numerous extensions of the GARCH Among the more well known are EGARCH (Nelson 1991) and GJR (Glosten et al 1993) Both models allow for
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 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 informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More information