Simulation of an SPDE Model for a Credit Basket
|
|
- Alison Blake
- 5 years ago
- Views:
Transcription
1 Simulation of an SPDE Model for a Credit Basket Christoph Reisinger Joint work with N. Bush, B. Hambly, H. Haworth christoph.reisinger@maths.ox.ac.uk. MCFG, Mathematical Institute, Oxford University C. Reisinger p.1
2 Outline Introduction, structural models of credit Example: Credit default swap spreads A general multi-factor model Large basket limit and CDO pricing Numerical simulation of the SPDE Calibration and pricing examples Improvements, extensions C. Reisinger p.2
3 Structural setup As in Merton (1974), and Black and Cox (1976), A t the company s asset value, governed by da t A t = µdt + σ dw t µ the mean return, σ the volatility, W a standard Brownian motion denoting the default threshold barrier by b, say constant, define the distance to default, X t, as X t = 1 σ ( loga t logb ). Merton: company defaults if X T < 0 Black/Cox: default time τ is given by the first time X t hits 0 C. Reisinger p.3
4 Default probabilities By first exit time theory, the probability of survival to T is where Q(T τ X t ) = Q(X s 0, t s T X t ) H(x, s) = Φ ( = H(X t, T t), ) ( x + ms e 2mx x + ms Φ ). s s Φ the standard Gaussian CDF, m the risk-neutral drift of X s C. Reisinger p.4
5 More companies Consider firm values, for i = 1,...,N (risk-neutral measure), da i t = (r f q i )A i t dt + σ i A i t dw i (t) where r f risk-free rate q i dividend yields σ i volatilities W i (t) Brownian motions and cov(w i (t), W j (t)) = ρ ij t. C. Reisinger p.5
6 More companies Consider firm values, for i = 1,...,N (risk-neutral measure), da i t = (r f q i )A i t dt + σ i A i t dw i (t) where r f risk-free rate q i dividend yields σ i volatilities W i (t) Brownian motions and cov(w i (t), W j (t)) = ρ ij t. We assume that each company has an exponential default barrier, b i (t) = K i e γ i(t t) for constants K i, γ i. T represents the maturity of the product. C. Reisinger p.5
7 Transformation Setting ( A Xt i i = ln t A i 0 ) e γ it = α i t + σ i W i (t) with α i = r f q i γ i 1 2 σ2 i, leads to a Brownian motion with ( ) drift and constant barrier B i = ln bi (0) 0, X0 i = 0. A i 0 C. Reisinger p.6
8 Transformation Setting ( A Xt i i = ln t A i 0 ) e γ it = α i t + σ i W i (t) with α i = r f q i γ i 1 2 σ2 i, leads to a Brownian motion with ( ) drift and constant barrier B i = ln bi (0) 0, X0 i = 0. Defining the running minimum A i 0 X i t = min 0 s t Xi s, and default time, τ i, as the first hitting time of the default barrier, τ i = inf{t : X i t = B i }, the survival probability is then Q(τ i > s) = Q(X i s B i). C. Reisinger p.6
9 2D case For two firms, this problem is analytically tractable, e. g. Q(t) = Q(X 1 t B 1, X 2 t B 2) = 2 ( ) β βt ea 1B 1 +a 2 B 2 +bt e r2 0 /2t nπθ0 sin β where g n (θ) = 0 n=1 a 1 = α 1σ 2 ρα 2 σ 1 (1 ρ 2 )σ 2 1 σ 2 ( re r2 /2t e A(θ)r rr0 ) I ( nπ β ) t dr, a 2 = α 2σ 1 ρα 1 σ 2 (1 ρ 2 )σ 1 σ ( ) nπθ sin g n (θ) dθ β b = α 1 a 1 α 2 a σ2 1a ρσ 1 σ 2 a 1 a σ2 2a ρ 2 tan β =, β [0, π] ρ etc, and I ( nπ β ) ( rr0 t ) is a modified Bessel s function. C. Reisinger p.7
10 CDS spread calculation The buyer of a k th -to-default credit default swap (CDS) on a basket of n companies pays a premium, the CDS spread, for the life of the CDS until maturity or the k th default. In the event of default by the k th underlying reference company, the buyer receives a default payment and the contract terminates. C. Reisinger p.8
11 CDS spread calculation The buyer of a k th -to-default credit default swap (CDS) on a basket of n companies pays a premium, the CDS spread, for the life of the CDS until maturity or the k th default. In the event of default by the k th underlying reference company, the buyer receives a default payment and the contract terminates. Equating discounted spread payment and discounted default payment, T DSP = ck e rfs Q(τ k > s) ds DDP = (1 R)K 0 T 0 e r fs Q(s τ k s + ds) gives the market k th -to-default CDS spread, c k, as { (1 R) 1 e rft Q(τ k > T) } T 0 r fe rfs Q(τ k > s) ds c k = T. 0 e rfs Q(τ k > s) ds C. Reisinger p.8
12 Varying T First-to-default CDS, varying T Second-to-default CDS, varying T year 2 years 3 years 4 years 5 years year 2 years 3 years 4 years 5 years Spread Correlation Spread Correlation σ 1 = σ 2 = 0.2, K 1 = 100, r f = 0.05, q 1 = q 2 = 0, γ 1 = γ 2 = 0.03, initial distance-to-default = 2, R = 0.5 C. Reisinger p.9
13 Varying σ First-to-default CDS, varying σ Spread σ i =0.1 σ i =0.15 σ i =0.2 σ i =0.25 σ i = Correlation Second-to-default CDS, varying σ Spread σ i =0.1 σ i =0.15 σ i =0.2 σ i =0.25 σ i = Correlation K 1 = 100, r f = 0.05, q 1 = q 2 = 0, R = 0.5 γ 1 = γ 2 = 0.03, initial distance to default = 2, T = 5 C. Reisinger p.10
14 Portfolio model Consider a portfolio of N different companies. Denoting the companies asset values at time t by A i t, we assume that under the risk neutral measure they follow a diffusion process given by da i t = µ(t, A i t) dt + σ(t, A i t) dw i t + m σ ij (t, A i t) dm j t, j=1 W i t and M j t are Brownian motions satisfying d W i t, M j t = 0 i, j and d W i t, W j t = d M i t, M j t = δ ij dt. C. Reisinger p.11
15 Numerical PDE solution Solve Kolmogorov equation numerically? Isotropic grid for N = 2 directions (firms): C. Reisinger p.12
16 Numerical PDE solution Solve Kolmogorov equation numerically? Isotropic grid for N = 2 directions (firms): C. Reisinger p.12
17 Numerical PDE solution Solve Kolmogorov equation numerically? Isotropic grid for N = 2 directions (firms): C. Reisinger p.12
18 Numerical PDE solution Solve Kolmogorov equation numerically? Isotropic grid for N = 2 directions (firms): Curse of dimensionality : After n refinements, unknowns h N 2 nn Accuracy ǫ has complexity C(ǫ) ǫ N/2 (for an order 2 method) exponential growth of required computational resources C. Reisinger p.12
19 What is high dimensional? Example, 30 firms: If each asset is represented by only two states, the total number of variables is already 2 30 = C. Reisinger p.13
20 What is high dimensional? Example, 30 firms: If each asset is represented by only two states, the total number of variables is already 2 30 = If we choose a reasonable number of points in each direction, say 32 = 2 5, the same total number is already obtained for dim = 6. C. Reisinger p.13
21 Sparse grids Full grid, h N points Sparse, h 1 log h N 1 But: require too much smoothness and have large constants. C. Reisinger p.14
22 Correlation data (typical) C. Reisinger 0.18 p.15 0
23 Correlation data l σ(σ) λ Spectral gap after λ Exponential decay from λ Asymptotic expansion? Not enough regularity C. Reisinger 0.18 p.16 0 k
24 A different viewpoint It is irrelevant which of the firms default. Assume interchangeable. Let ν N,t the empirical measure for the entire portfolio, ν N,t = 1 N N i=1 δ A i t Let φ Cc (R) and for measure ν t write φ, ν t = φ(x)ν t (dx). Define a family of processes F N,φ t by F N,φ t = φ, ν N,t = 1 N N φ(a i t). i=1 C. Reisinger p.17
25 Evolution of measure Applying Ito s lemma to F N,φ t F N,φ t = 1 N N i=1 t N we have φ (A i s) µ(s, A i s) dt + σ(s, A i s) dws i + N i=1 t φ (A i s) m σij 2 + σ 2 ds. j=1 m σ ij (s, A i s) dms j j=1 We now assume that σ ij = σ j for all i = 1, 2..., then t t m F N,φ t = F N,φ 0 + Aφ, ν N,s ds + σ j φ, ν N,s dms j t 0 1 N 0 j=1 N φ (A i s) σ(s, A i s) dws. i i=1 A = µ x σ2 2 x 2, σ2 = m σj 2 + σ 2, j=1 C. Reisinger p.18
26 ing equation The idiosyncratic component becomes deterministic in the infinite dimensional limit, see also [Kurtz & Xiong, 1999]. As the market factors affect each company to the same degree, we can write F N,φ t F φ t = φ, ν t as N, with df φ t = Aφ, ν t dt + m σ j φ, ν t dm j t. j=1 Alternatively, we can write this in integrated form as φ, ν t = φ, ν 0 + t 0 Aφ, ν s ds + m t j=1 0 σ j φ, ν s dm j s C. Reisinger p.19
27 Density Assume the measure ν t to be absolutely continuous with respect to the Lebesgue measure, to write ν t (dx) = v(t, x)dx In differential form, dv = µ v x dt σ2 2 v x 2 dt m j=1 x (σ j(t, x)v)dm j t. This is a stochastic PDE that describes the evolution of an infinite portfolio of assets whose dynamics were given before. For solubility see [Krylov, 1994]. Can now use this to approximate the loss distribution for a portfolio of fixed size N. C. Reisinger p.20
28 Simplification Assume asset processes are correlated via a single factor. Specifically, let the SDEs for the asset processes be given by da i t A i t = r dt + 1 ρσ dw i t + ρσ dm t, A i 0 = a i, where ρ [0, 1), d W i t, M t = 0. The distance-to-default X i t = 1 σ ( log A i t log B i t ) evolves according to dx i t = µdt + 1 ρdw i t + ρdm t, X i 0 = x i with µ = 1 σ ( r 1 2 σ2) and x i = 1 σ ( log a i log B i 0 ). C. Reisinger p.21
29 Portfoilio loss Constant coefficient SPDE for v. The portfolio loss variable L N t = (1 R) 1 N measures the losses at time t. N i=1 1 {τi t} Assume that a fraction R of losses, 0 R 1, the recovery rate, is recovered after default. Absorbing barrier condition: v(t, 0) = 0 for t 0. The loss distribution is then ( L t = (1 R) 1 0 ) v(t, x) dx. C. Reisinger p.22
30 Properties Matching the initial conditions for both measures, we need v(0, x) = 1 N N δ A i 0. Given this initial condition, i=1 L 0 = N ( 1 1 N N δ A i 0 i=1 B t ) = 0, Also, 0 L t 1, for t 0 Q(L s K) Q(L t K), for s t, which ensure that there is no arbitrage in the loss distribution. C. Reisinger p.23
31 CDOs CDOs slice up the losses into tranches, defined by their attachment and detachment points. Typical numbers are 0-3%, 3%-6%, 6%-9%, 9%-12%, 12%-22%, 22%-100%. For an attachment point a and detachment point d > a, the outstanding tranche notional and tranche loss X t = max(d L t, 0) max(a L t, 0) Y t = (d a) X t = max(l t a,0) max(l t d, 0) determine the spread and default payments for that tranche. Spreads are quoted as an annual payment, as a ratio of the a notional, but assumed to be paid quarterly. a there is a variation for the equity tranche, but we do not go into details. C. Reisinger p.24
32 Spread calculation Assume the notional to be 1, c the spread payment; n the maximum number of payments up to expiry T ; T i the payment dates for 1 i n, δ = 0.25 the interval between payments. The expected sum of discounted fee payments, referred to as the fee leg, is given by cv fee = c The protection leg is n δe rt i E Q [X Ti ], i=1 V prot = n e rt i E Q [X Ti 1 X Ti ]. i=1 The fair spread payment is then given by s = V prot V fee. C. Reisinger p.25
33 Computation Recall the SPDE dv = 1 (r 12 ) σ σ2 v x dt v xx dt ρv x dm t, t > 0, x > 0 v(0, x) = v 0 (x), x > 0 v(t, 0) = 0 t 0 Monte Carlo on top of a PDE solver computationally costly. For each simulated path of the market factor, a PDE needs to be solved, e. g. by a finite difference method. From the PDE solution, quantities like tranche losses can be computed, then averaged over the simulated paths. C. Reisinger p.26
34 Discrete defaults First make the assumption that defaults are only observed at a discrete set of times, taken quarterly to coincide with the payment dates, i.e. if a firm s value is below the default barrier on one of the observation dates T i, removed from the basket. We therefore solve the modified SPDE problem dv = 1 σ (r 12 ) σ2 v x dt v xx dt ρv x dm t, t (T k, T k+1 ), v(0, x) = v 0 (x), v(t k, x) = 0 x 0, 0 < k n C. Reisinger p.27
35 Lagrangian coordinate The boundary condition is not active in intervals (T k, T k+1 ). The Brownian driver only introduces a random shift. The solution can therefore be written as v(t, x) = 8 < : v (i) (t T k, x ρ(m t M Tk )) 0 x 0 t {T k, 1 i n} else if t (T k, T k+1 ], 0 k < n where v (k) is the solution to the (deterministic) problem v (k) t = 1 (1 ρ)v(k) xx 1 2 σ (r 12 σ2 ) v (k) x, t (0, τ) = (0, T k+1 T k ) v (k) (0, x) = v(t k, x) assuming payment dates are equally spaced, τ = T k+1 T k. C. Reisinger p.28
36 Algorithm This suggests the following inductive strategy for k = 0,...,n 1: 1. Start with v (0) (0, x) = v 0 (x). 2. Solve the PDE numerically in the interval (0, T 1 ), to obtain v (0) (T 1, x). 3. Simulate M T1, evaluate v(t 1, x). 4. For k > 0, having computed v(t k, x) in the previous step, use this as initial condition for v (k), and repeat until k = n. C. Reisinger p.29
37 Finite differences The initial distribution is assumed localised. Approximate the distribution by one with support [x min, x max ] with x min < 0 and x max > 0. Can ensure that the expected error of this approximation is much smaller than the standard error of the Monte Carlo estimates. Grid x 0 = x min, x 1 = x min + x,...,x min + j x,..., x J = x min + J x = x max, where x = (x max x min )/J, timesteps t 0 = 0, t 1 = t,...,t I = I t = τ, where t = τ/i. Define an approximation v i j to v(t i, x j ) as solution to a FD/FE scheme with θ timestepping. C. Reisinger p.30
38 Stability Standard central differencing is second order accurate in x. The backward Euler scheme θ = 1 is of first order accurate (in t) and strongly stable. The Crank-Nicolson scheme θ = 1 2 is of second order accurate, and is unconditionally stable in the l 2 -norm. We deal with initial conditions of the form v(0, x) = 1 N N δ(x x i ), i=1 Crank-Nicolson timestepping gives spurious oscillations for Dirac initial conditions, and reduces the convergence order for discontinuous intitial conditions. C. Reisinger p.31
39 Rannacher timestepping Rannacher proposed the following modification: Replace first Crank-Nicolson steps with backward Euler steps. Need to balance between accuracy and stability. Analysis by Giles and Carter of the heat equation suggests to replace the first two Crank-Nicolson steps by four backward Euler steps of half the stepsize. We do this at t = 0, and also at t = T k where the interface conditions introduce discontinuities at x = 0. This restores second order convergence in time. C. Reisinger p.32
40 Averaging The sum of δ-distributions needs approximation on the grid. Could collect the firms into symmetric intervals of width x around grid points, v 0 j = 1 N x xj + x/2 x j x/2 δ(x i x) dx, however, this reduces the overall order of the finite difference scheme to 1: The approximation cannot distinguish between initial positions x i in an interval of length x. C. Reisinger p.33
41 Projection To achieve higher (i. e. second) order, the δ s are split between adjacent grid points. The correct weighting for a single firm with distance-to-default x i in the interval [x j, x j+1 ), is v 0 k = x 2 (x j+1 x i ) k = j x 2 (x i x j ) k = j else This can be written more elegantly as L 2 -projection onto the basis of hat functions Φ k 0 k N where Φ k (x) = 1 x min (max(x x k + x), max(x k x,0)), then v 0 k = Φ k, v 0 = xmax x min Φ k (x)v 0 (x) dx. Note that x N k=0 v0 k = 1. C. Reisinger p.34
42 Interface condition At t = T k, we have to evaluate the grid function at shifted arguments that do not normally coincide with the grid. First, define a piecewise linear reconstruction from the approximation vk I obtained in the last step over the previous interval [T k 1, T k ], as N v x (T k, x) = Φ j (x M)vj. I j=0 Then, approximate the shift by setting v 0 j = xj + x/2 max(x j x/2,0) v x (T k, x) dx, with M = M Tk M Tk 1. C. Reisinger p.35
43 Interface condition Use this as initial condition for the next interval (the integral is understood to be 0 if the lower limit is larger than the upper limit). This ensures that x J vj 0 = j=0 xmax 0 v h (T k, x) dx, so the cumulative density of firms with firm values greater than 0 is preserved. Also, the solution is smoothened at x = 0 C. Reisinger p.36
44 Simulations For a given realisation of the market factor, we can approximate the loss functional L Tk at time T k by L x T k = 1 xmax 0 v x (T k, x) dx Explicitly include the dependency L x T k (Φ), where Φ i are drawn independently from a standard normal distribution, Then for N sims simulations with samples Φ l = (Φ l k ) 1 k n, 1 l N sims, E Q [X Tk ] E Q [max(d L x T k (Φ), 0) max(a L x T k (Φ), 0)] 1 N sims N sims l=1 ( max(d L x T k (Φ l ), 0) max(a L x T k (Φ l ), 0) ) C. Reisinger p.37
45 Data The simulation error has two components, the discretisation error and the variance of the Monte Carlo estimate. For the following simulations we have used these data: T = 5, r = 0.027, σ = 0.24, R = 0.7, ρ = Initial positions for individual firms, calibrated to their individual CDS spreads, were well within the range [x min, x max ] = [ 10, 20]. C. Reisinger p.38
46 Grid convergence First consider the discretisation error in t and x ǫj J ǫi I Extrapolation-based estimator for the discretisation error of L Tn for increasing J (left) and I (right) for a single realisation of the path of the market factor. We clearly see second order convergence in t and x. C. Reisinger p.39
47 MC convergence s [0, 3%] [3%, 6%] [6%, 9%] [9%, 12%] [12%, 22%] [22%, e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e e Monte Carlo estimates for expected outstanding tranche notionals for N sims = 64 4 s 1 Monte Carlo runs, s = 1,..., 10. The finite difference parameters were fixed at J = 256, I = 4. C. Reisinger p.40
48 MC convergence ], a = 0, b = 0.03 E Q [XTn log 4 N sims ], a = 0.06, b = x E Q [XTn log 4 N sims ], a = 0.12, b = 0.22 E Q [XTn x 10 7 log 4 N sims Monte Carlo estimates with standard error bars for expected losses in tranches [0, 3%], [6%, 9%], [12%, 22%] for N sims = 16 4 k 1, k = 1,...,10, J = 256, I = 4. C. Reisinger p.41
49 Index spreads The fixed coupons, traded spreads and model spreads for the itraxx Main Series 10 index. Maturity Date Fixed Coupon (bp) Traded Spread (bp) Model Spread (bp) 20/12/ /12/ /12/ February 22, Parameters used for the model spreads are r = 0.042, σ = 0.22, R = 0.4. Maturity Date Fixed Coupon (bp) Traded Spread (bp) Model Spread (bp) 20/12/ /12/ /12/ December 5, Parameters used for the model spreads are r = 0.033, σ = 0.136, R = 0.4. C. Reisinger p.42
50 Implied correlation Implied Correlation Skew for itraxx Main Series 6 Tranches, Feb 22, The implied correlation for each tranche is the value of correlation that gives a model tranche spread equal to the market tranche spread. Implied Corrleation (%) Year 7 Year 10 Year % 6% 9% 12% 22% 100% Tranche Detachment Point Model parameters are r = 0.042, σ = 0.22, R = 0.4. C. Reisinger p.43
51 5 Year Tranche Market ρ = 0.3 ρ = 0.4 ρ = 0.5 ρ = 0.6 ρ = 0.7 ρ = 0.8 ρ 0%-3% 71.5 % % 75.9 % % % % % 4 3%-6% %-9% %-12% %-22% %-100% Model tranche spreads (bp) for varying values of the correlation parameter. The equity tranches are quoted as an upfront assuming a 500bp running spread. The model is calibrated to the itraxx Main Series 10 index for Dec 5, Market levels shown are for this date; model parameters are r = 0.033, σ = 0.136, R = 0.4. C. Reisinger p.44
52 7 Year Tranche Market ρ = 0.3 ρ = 0.4 ρ = 0.5 ρ = 0.6 ρ = 0.7 ρ = 0.8 ρ 0%-3% 72.9 % % % % % 60 % % 4 3%-6% %-9% %-12% %-22% %-100% Model tranche spreads (bp) for varying values of the correlation parameter. The equity tranches are quoted as an upfront assuming a 500bp running spread. The model is calibrated to the itraxx Main Series 10 index for Dec 5, Market levels shown are for this date; model parameters are r = 0.033, σ = 0.136, R = 0.4. C. Reisinger p.45
53 10 Year Tranche Market ρ = 0.3 ρ = 0.4 ρ = 0.5 ρ = 0.6 ρ = 0.7 ρ = 0.8 0%-3% 73.8 % % % % % % % 3%-6% %-9% %-12% %-22% %-100% Model tranche spreads (bp) for varying values of the correlation parameter. The equity tranches are quoted as an upfront assuming a 500bp running spread. The model is calibrated to the itraxx Main Series 10 index for Dec 5, Market levels shown are for this date; model parameters are r = 0.033, σ = 0.136, R = 0.4. C. Reisinger p.46
54 Dynamic properties, forward starting CDOs, options on tranches,... Simulation of SPDEs driven by jumps, with Karolina Bujok. Calibration, including jumps. Analytic work by Hambly/Jin on jumps, contagion, granularity. Continuous defaults, with Mike Giles, via Multi-Level Monte Carlo. Convergence analysis through mean-square stability. Importance sampling for senior tranches, with Tom Dean. C. Reisinger p.47
55 Importance sampling Using large deviations theory, see also [Glasserman et al, 1999, 2007], preliminary observations: without importance sampling with importance sampling Probability of more than 3% defaults occuring for n Monte Carlo paths C. Reisinger p.48
56 Importance sampling Using large deviations theory, see also [Glasserman et al, 1999, 2007], preliminary observations: without importance sampling 1.8 x with importance sampling 1.8 x Probability of more than 10% defaults occuring for n Monte Carlo paths C. Reisinger p.48
57 Literature References [Black & Cox(1976)] Black, F. & Cox, J. (1976) Valuing Corporate Securities: Some Effects of Bond Indenture Provisions,J. Finance, 31, [Cameron & Martin(1947)] Cameron, R. H., & Martin, W. T. (1947) The Orthogonal Development of Non-Linear Functionals in Series of Fourier-Hermite Functionals, Ann. Math, 48, [Giles & Carter(2006)] Giles, M. & Carter, R. (2006) Convergence analysis of Crank-Nicolson and Rannacher time-marching, J. Comp. Fin., 9(4), [Hambly & Jin (2008)] Hambly, B.M. & Jin, L. (2008) SPDE approximations for large basket portfolio credit modelling, in preparation. [Haworth & Reisinger(2007)] Haworth, H., & Reisinger, C. (2007) Modelling Basket Credit Default Swaps with Default Contagion, J. Credit Risk, 3(4), C. Reisinger p.49
The Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets (Hull chapter: 12, 13, 14) Liuren Wu ( c ) The Black-Scholes Model colorhmoptions Markets 1 / 17 The Black-Scholes-Merton (BSM) model Black and Scholes
More informationThe Black-Scholes Model
The Black-Scholes Model Liuren Wu Options Markets Liuren Wu ( c ) The Black-Merton-Scholes Model colorhmoptions Markets 1 / 18 The Black-Merton-Scholes-Merton (BMS) model Black and Scholes (1973) and Merton
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 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 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 Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
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 informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
More informationCalibration Lecture 4: LSV and Model Uncertainty
Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where
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 informationCredit Risk using Time Changed Brownian Motions
Credit Risk using Time Changed Brownian Motions Tom Hurd Mathematics and Statistics McMaster University Joint work with Alexey Kuznetsov (New Brunswick) and Zhuowei Zhou (Mac) 2nd Princeton Credit Conference
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 informationMultilevel Monte Carlo Simulation
Multilevel Monte Carlo p. 1/48 Multilevel Monte Carlo Simulation Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance Workshop on Computational
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More 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 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 informationModelling Bonds & Credit Default Swaps using a Structural Model with Contagion
Modelling Bonds & Credit Default Swaps using a Structural Model with Contagion Helen Haworth, Christoph Reisinger and William Shaw June 6 Abstract This paper develops a two-dimensional structural framework
More informationBIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS
BIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS PRICING EMMS014S7 Tuesday, May 31 2011, 10:00am-13.15pm
More informationComputational Finance
Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk
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 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 informationLecture 8: The Black-Scholes theory
Lecture 8: The Black-Scholes theory Dr. Roman V Belavkin MSO4112 Contents 1 Geometric Brownian motion 1 2 The Black-Scholes pricing 2 3 The Black-Scholes equation 3 References 5 1 Geometric Brownian motion
More 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 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 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 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 informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationConstructing Markov models for barrier options
Constructing Markov models for barrier options Gerard Brunick joint work with Steven Shreve Department of Mathematics University of Texas at Austin Nov. 14 th, 2009 3 rd Western Conference on Mathematical
More informationMonte Carlo Methods. Prof. Mike Giles. Oxford University Mathematical Institute. Lecture 1 p. 1.
Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Lecture 1 p. 1 Geometric Brownian Motion In the case of Geometric Brownian Motion ds t = rs t dt+σs
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 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 informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More 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 informationA new approach for scenario generation in risk management
A new approach for scenario generation in risk management Josef Teichmann TU Wien Vienna, March 2009 Scenario generators Scenarios of risk factors are needed for the daily risk analysis (1D and 10D ahead)
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationCredit Risk. June 2014
Credit Risk Dr. Sudheer Chava Professor of Finance Director, Quantitative and Computational Finance Georgia Tech, Ernest Scheller Jr. College of Business June 2014 The views expressed in the following
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 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 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 Volatility (Working Draft I)
Stochastic Volatility (Working Draft I) Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu 1 Introduction When using the Black-Scholes-Merton model to price derivative
More informationModule 4: Monte Carlo path simulation
Module 4: Monte Carlo path simulation Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Module 4: Monte Carlo p. 1 SDE Path Simulation In Module 2, looked at the case
More informationAsset Pricing Models with Underlying Time-varying Lévy Processes
Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University
More 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 informationEuropean option pricing under parameter uncertainty
European option pricing under parameter uncertainty Martin Jönsson (joint work with Samuel Cohen) University of Oxford Workshop on BSDEs, SPDEs and their Applications July 4, 2017 Introduction 2/29 Introduction
More 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 informationDYNAMIC CORRELATION MODELS FOR CREDIT PORTFOLIOS
The 8th Tartu Conference on Multivariate Statistics DYNAMIC CORRELATION MODELS FOR CREDIT PORTFOLIOS ARTUR SEPP Merrill Lynch and University of Tartu artur sepp@ml.com June 26-29, 2007 1 Plan of the Presentation
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
More informationTheoretical Problems in Credit Portfolio Modeling 2
Theoretical Problems in Credit Portfolio Modeling 2 David X. Li Shanghai Advanced Institute of Finance (SAIF) Shanghai Jiaotong University(SJTU) November 3, 2017 Presented at the University of South California
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationPartial differential approach for continuous models. Closed form pricing formulas for discretely monitored models
Advanced Topics in Derivative Pricing Models Topic 3 - Derivatives with averaging style payoffs 3.1 Pricing models of Asian options Partial differential approach for continuous models Closed form pricing
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
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 informationAsset-based Estimates for Default Probabilities for Commercial Banks
Asset-based Estimates for Default Probabilities for Commercial Banks Statistical Laboratory, University of Cambridge September 2005 Outline Structural Models Structural Models Model Inputs and Outputs
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 13. The Black-Scholes PDE Steve Yang Stevens Institute of Technology 04/25/2013 Outline 1 The Black-Scholes PDE 2 PDEs in Asset Pricing 3 Exotic
More informationThe Black-Scholes Model
IEOR E4706: Foundations of Financial Engineering c 2016 by Martin Haugh The Black-Scholes Model In these notes we will use Itô s Lemma and a replicating argument to derive the famous Black-Scholes formula
More informationBluff Your Way Through Black-Scholes
Bluff our Way Through Black-Scholes Saurav Sen December 000 Contents What is Black-Scholes?.............................. 1 The Classical Black-Scholes Model....................... 1 Some Useful Background
More informationMultilevel Monte Carlo for Basket Options
MLMC for basket options p. 1/26 Multilevel Monte Carlo for Basket Options Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance WSC09,
More informationAN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL
AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it 1 Stylized facts Traders use the Black-Scholes formula to price plain-vanilla options. An
More informationA Two Factor Forward Curve Model with Stochastic Volatility for Commodity Prices arxiv: v2 [q-fin.pr] 8 Aug 2017
A Two Factor Forward Curve Model with Stochastic Volatility for Commodity Prices arxiv:1708.01665v2 [q-fin.pr] 8 Aug 2017 Mark Higgins, PhD - Beacon Platform Incorporated August 10, 2017 Abstract We describe
More informationMultiname and Multiscale Default Modeling
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
More informationVolatility Trading Strategies: Dynamic Hedging via A Simulation
Volatility Trading Strategies: Dynamic Hedging via A Simulation Approach Antai Collage of Economics and Management Shanghai Jiao Tong University Advisor: Professor Hai Lan June 6, 2017 Outline 1 The volatility
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
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 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 information2.3 Mathematical Finance: Option pricing
CHAPTR 2. CONTINUUM MODL 8 2.3 Mathematical Finance: Option pricing Options are some of the commonest examples of derivative securities (also termed financial derivatives or simply derivatives). A uropean
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 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 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 informationValuation of Forward Starting CDOs
Valuation of Forward Starting CDOs Ken Jackson Wanhe Zhang February 10, 2007 Abstract A forward starting CDO is a single tranche CDO with a specified premium starting at a specified future time. Pricing
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 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 informationThe Black-Scholes PDE from Scratch
The Black-Scholes PDE from Scratch chris bemis November 27, 2006 0-0 Goal: Derive the Black-Scholes PDE To do this, we will need to: Come up with some dynamics for the stock returns Discuss Brownian motion
More informationLIBOR models, multi-curve extensions, and the pricing of callable structured derivatives
Weierstrass Institute for Applied Analysis and Stochastics LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives John Schoenmakers 9th Summer School in Mathematical Finance
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 informationUnified Credit-Equity Modeling
Unified Credit-Equity Modeling Rafael Mendoza-Arriaga Based on joint research with: Vadim Linetsky and Peter Carr The University of Texas at Austin McCombs School of Business (IROM) Recent Advancements
More informationCredit Modeling and Credit Derivatives
IEOR E4706: Foundations of Financial Engineering c 2016 by Martin Haugh Credit Modeling and Credit Derivatives In these lecture notes we introduce the main approaches to credit modeling and we will largely
More informationInvestigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs. Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2
Investigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2 1 Hacettepe University Department of Actuarial Sciences 06800, TURKEY 2 Middle
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 1.1287/opre.11.864ec e-companion ONLY AVAILABLE IN ELECTRONIC FORM informs 21 INFORMS Electronic Companion Risk Analysis of Collateralized Debt Obligations by Kay Giesecke and Baeho
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 informationOn the relative pricing of long maturity S&P 500 index options and CDX tranches
On the relative pricing of long maturity S&P 5 index options and CDX tranches Pierre Collin-Dufresne Robert Goldstein Fan Yang May 21 Motivation Overview CDX Market The model Results Final Thoughts Securitized
More informationConvergence Analysis of Monte Carlo Calibration of Financial Market Models
Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration
More informationIntroduction to Financial Mathematics
Department of Mathematics University of Michigan November 7, 2008 My Information E-mail address: marymorj (at) umich.edu Financial work experience includes 2 years in public finance investment banking
More informationStochastic Computation in Finance
Stochastic Computation in Finance Chuan-Hsiang Han Dept. of Quantitative Finance, NTHU Dept of Math & CS Education TMUE November 3, 2008 Outline History of Math and Finance: Fundamental Problems in Modern
More informationA Brief Introduction to Stochastic Volatility Modeling
A Brief Introduction to Stochastic Volatility Modeling Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction When using the Black-Scholes-Merton model to
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More 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 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 informationPortability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans
Portability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans An Chen University of Ulm joint with Filip Uzelac (University of Bonn) Seminar at SWUFE,
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 informationAD in Monte Carlo for finance
AD in Monte Carlo for finance Mike Giles giles@comlab.ox.ac.uk Oxford University Computing Laboratory AD & Monte Carlo p. 1/30 Overview overview of computational finance stochastic o.d.e. s Monte Carlo
More informationlast problem outlines how the Black Scholes PDE (and its derivation) may be modified to account for the payment of stock dividends.
224 10 Arbitrage and SDEs last problem outlines how the Black Scholes PDE (and its derivation) may be modified to account for the payment of stock dividends. 10.1 (Calculation of Delta First and Finest
More information"Vibrato" Monte Carlo evaluation of Greeks
"Vibrato" Monte Carlo evaluation of Greeks (Smoking Adjoints: part 3) Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance MCQMC 2008,
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 informationForeign Exchange Derivative Pricing with Stochastic Correlation
Journal of Mathematical Finance, 06, 6, 887 899 http://www.scirp.org/journal/jmf ISSN Online: 6 44 ISSN Print: 6 434 Foreign Exchange Derivative Pricing with Stochastic Correlation Topilista Nabirye, Philip
More informationThe 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 Reading group, March 13, 2017 Reading group, March
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 investment game in incomplete markets
The investment game in incomplete markets M. R. Grasselli Mathematics and Statistics McMaster University Pisa, May 23, 2008 Strategic decision making We are interested in assigning monetary values to strategic
More information