Quasi-Monte Carlo for Finance
|
|
- Marilyn Lester
- 5 years ago
- Views:
Transcription
1 Quasi-Monte Carlo for Finance Peter Kritzer Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy of Sciences Linz, Austria NCTS, Taipei, November 2016 Peter Kritzer Quasi-Monte Carlo for Finance 1
2 1 Introduction 2 Options 3 Brownian paths 4 Non-uniform random number generation using QMC 5 Generation of Brownian paths 6 Formulation as an integration problem 7 Challenges using QMC 8 Concluding remarks Peter Kritzer Quasi-Monte Carlo for Finance 2
3 Introduction Peter Kritzer Quasi-Monte Carlo for Finance 3
4 Suppose we would like to valuate a complicated financial product whose payoff depends on a stochastic process (S t ) t 0 describing the price of an asset S at time t. (1) How can we do this? (2) Where does QMC come into play? Peter Kritzer Quasi-Monte Carlo for Finance 4
5 How? (1) Find a suitable model for the (stochastic) price process S t. (2) In many cases, this will involve a stochastic differential equation and Brownian motion. (3) Find an approximation Ŝt to S t. (4) Use a suitable method to simulate Brownian motion. (4) Use a suitable method to simulate Brownian motion. Where? Here! Peter Kritzer Quasi-Monte Carlo for Finance 5
6 Options Peter Kritzer Quasi-Monte Carlo for Finance 6
7 What is an option? Simple definition: a contract that gives one party the right to buy/sell an asset at some point of time in the future at a previously fixed price. Simplest example: European call option: Gives one party the right to buy an asset (e.g. a stock) at a future point of time T > 0 at price K. T is called the maturity of the option, K the strike price, the stock the underlying. Value of European call: max(s T K, 0). Peter Kritzer Quasi-Monte Carlo for Finance 7
8 A more sophisticated example: Suppose the stock price follows a stochastic process (S t ) t 0. Consider an option that pays ( ) 1 d max S tk K, 0. d k=1 at time T, where S tk is the price of the stock at time t k, k = 1,..., d. This option is called fixed strike Asian call option. What is the value of this option at time 0 (i.e., right now)? The answer is not straightforward, because it depends on the whole path of S t in the interval [0, T ]. We need a model to simulate the stochastic process (S t ) t 0. Peter Kritzer Quasi-Monte Carlo for Finance 8
9 Brownian paths Peter Kritzer Quasi-Monte Carlo for Finance 9
10 There are different kinds of models to describe stochastic processes such as the stock price (S t ) t 0. Depending on the context, they might be simple or more complicated. Basically all models involve Brownian motion. Many models involve stochastic differential equations (SDEs). Here we use a method stemming from the Heston model. Peter Kritzer Quasi-Monte Carlo for Finance 10
11 Brownian motion: Stochastic process in continuous time on a probability space (Ω, F, P) satisfying: (1) W 0 = 0 almost surely. (2) For any t 1, t 2 0 the random variables W t2 W t1 and W t2 t 1 have the same distribution. (3) For any t 1,..., t d with t 0 := 0 < t 1 < t 2 < < t d the random variables W t1 W t0,..., W td W td 1 are independent. (4) 1/tW t is a standard normal random variable for every t > 0. (5) Fore each ω Ω the mapping t W t (ω) is continuous. Peter Kritzer Quasi-Monte Carlo for Finance 11
12 The Heston model relies on SDEs that can be solved approximately (only a few ideas here, no details). We first need to discretize the time interval [0, T ]. For the sake of simplicity, assume equal-length time intervals: Let h = T d. Discretize [0, T ] by time nodes 0, h, 2h,..., dh = T. We could also discretize such that time intervals have different length, but for the moment use the simpler version. Peter Kritzer Quasi-Monte Carlo for Finance 12
13 Then we estimate the payoff function of the option, ( ) 1 d max S kh K, 0. d k=1 We approximate the option price by estimating the prices S kh, max ( 1 d ) d Ŝ kh K, 0. k=1 The Ŝkh depend on the increments of a Brownian motion, W k+1 = W (k+1)h W kh. Peter Kritzer Quasi-Monte Carlo for Finance 13
14 Need to simulate a discrete Brownian paths, where h = T /d. How can we do this? (W h, W 2h,..., W T ). First question: How to generate non-uniform (e.g., normally distributed) random numbers? Peter Kritzer Quasi-Monte Carlo for Finance 14
15 Non-uniform random number generation using QMC Peter Kritzer Quasi-Monte Carlo for Finance 15
16 The inversion method: Suppose that F : R [0, 1] is the cumulative distribution function (CDF) of a real random variable X, i.e., F(x) = P[X x]. Assume F is invertible, i.e., G : (0, 1) R such that G(F (x)) = x for x R and F(G(u)) = u for u (0, 1). Let U be uniformly distributed on (0, 1). Then Y = G(U) has the same distribution as X. Peter Kritzer Quasi-Monte Carlo for Finance 16
17 Instead of a random U: Use elements from a quasi-random sequence in (0, 1)! Sequences of node sets from QMC methods are usually evenly distributed quasi-random sequences in (0, 1). Examples: van der Corput sequence, Kronecker sequence, etc. Peter Kritzer Quasi-Monte Carlo for Finance 17
18 We can proceed analogously for sampling from multivariate distributions if the components are independent (e.g., standard normal vectors). Again, take node sets from multivariate QMC methods (Halton sequence, Sobol sequence, Niederreiter sequence, etc.). Slight drawback: There is no closed elementary formula for the inverse CDF of a normal distribution. However, there are suitable implementations in most software packages, and actually we only need to know how to handle the standard normal CDF. Alternative methods: importance sampling, acceptance rejection. Peter Kritzer Quasi-Monte Carlo for Finance 18
19 Generation of Brownian paths Peter Kritzer Quasi-Monte Carlo for Finance 19
20 Recall our original goal to generate discrete Brownian paths, (W h, W 2h,..., W T ). where h = T /d. More generally, how to generate a Brownian path (W t1,..., W td ), with 0 < t 1 < < t d? Easy to check: (W t1,..., W td ) is a Gaussian vector with mean zero and covariance matrix t 1 t 1 t 1... t 1 t 1 t 2 t 2... t 2 (min(t j, t k )) d j,k=1 = t 1 t 2 t 3... t t 1 t 2 t 3... t d Peter Kritzer Quasi-Monte Carlo for Finance 20
21 The forward method: Suppose we are given a standard normal vector X = (X 1,..., X d ). Calculate the discrete Brownian path (W t1,..., W td ) as follows. Set W t1 = t 1 X 1, and W tk+1 = W tk + t k+1 t k X k+1. The vector (W t1,..., W td ) has the required covariance matrix, which can be checked by using E[X j X k ] = δ j,k. Advantage: generation of a path requires generation of (X 1,..., X d ), d multiplications, and d 1 additions. Peter Kritzer Quasi-Monte Carlo for Finance 21
22 The Brownian bridge construction: Let π be a permutation of the set {1,..., d}. Assume we want to compute the W ti in the order W tπ(1),..., W tπ(d). Brownian bridge construction makes this possible. Special case: π = id: forward construction. Peter Kritzer Quasi-Monte Carlo for Finance 22
23 Two key ideas underlying BB construction: Lemma Let W be a Brownian motion and let t 1 < t 2 < t 3. Then the conditional distribution of W t2 given W t1 and W t3 is N(µ, σ 2 ) with µ = t 3 t 2 t 3 t 1 W t1 + t 2 t 1 t 3 t 1 W t3 and σ 2 = (t 3 t 2 )(t 2 t 1 ) t 3 t 1. Furthermore: For computing W tπ(j), we only need knowledge of the two elements next to the right and the left, respectively. The covariance matrix is as before. Use of QMC: Use a cleverly chosen QMC sequence to determine π (e.g, van der Corput sequence). Then we can figure out beforehand which elements are next right and left for each W π(j). This can make the computation efficient. Peter Kritzer Quasi-Monte Carlo for Finance 23
24 The PCA construction: The covariance matrix of (W t1,..., W td ) can be written in the form VD 1/2 D 1/2 V 1, where D is a diagonal matrix with positive entries, and V is an orthogonal matrix, i.e., VV = V V = I. Let X be a d-dimensional standard normal vector. Then we can construct (W t1,..., W td ) by (W t1,..., W td ) = VDX. Peter Kritzer Quasi-Monte Carlo for Finance 24
25 The PCA construction principle can be extended: Due to the properties of D and V, we have VD 1/2 D 1/2 V 1 = (VD)(VD). Let Σ denote the covariance matrix of the W tj. Then the following lemma holds. Lemma Let A be any d d matrix with AA = Σ. Let X be a standard normal vector. Then Y = AX is a normal vector with covariance matrix Σ. Peter Kritzer Quasi-Monte Carlo for Finance 25
26 Recall: In our example, the goal was to generate where h = T /d. (W h, W 2h,..., W T ), Is evenly spaced discretization an advantage? Suppose we would like to generate (W 1/d, W 2/d,..., W d/d ), i.e., discretization with step-size 1/d, where we normed (w.l.o.g.) T := 1. Peter Kritzer Quasi-Monte Carlo for Finance 26
27 In this setting, the covariance matrix is D := d d Then we can easily find the Cholesky decomposition of D, D = SS, with S := d Peter Kritzer Quasi-Monte Carlo for Finance 27
28 Becoming even more general: Lemma Let A be any d d matrix with AA = D. Then there exists an orthogonal d d matrix V such that A = SV, i.e. (SV )(SV ) = D. For a standard normal X, W = AX is a Brownian path with discretization 1/d, 2/d,..., d/d. Conversely, given any orthogonal matrix V, the matrix Z := SV satisfies ZZ = D. Hence, for evenly spaced discretization, many methods can be expressed using orthogonal matrices. Peter Kritzer Quasi-Monte Carlo for Finance 28
29 Examples: For the forward construction, the corresponding orthogonal matrix is the identity matrix. Cleverly choosing the orthogonal transform can be very helpful: e.g., we can reduce the computation cost of the PCA construction (from O(d 2 ) to O(d log d)). Peter Kritzer Quasi-Monte Carlo for Finance 29
30 One more advantage: In the Asian option example, actually only the increments of the Brownian path are relevant. Observe that if then (W 1/d, W 2/d,..., W d/d ) = SVX (W 1/d, W 2/d W 1/d,..., W d/d W (d 1)/d ) = VX. Peter Kritzer Quasi-Monte Carlo for Finance 30
31 Formulation as an integration problem Peter Kritzer Quasi-Monte Carlo for Finance 31
32 How to view all of this as an integration problem? The exemplary problem was to compute the fair price P ( expected value) of the Asian option. We approximated the price process S (and thereby the price of P) by the solution of an SDE. This involved a time discretization using d steps. This corresponds to approximating E(P) by E(P D ), where E(P D ) = p(x)φ(x) dx R d where p is some complicated function of x R d and φ is the density of the standard normal distribution. φ comes into play because it occurs in the simulation of Brownian paths. Peter Kritzer Quasi-Monte Carlo for Finance 32
33 Using, e.g., the inversion method we sampled from the standard normal distribution. By this step, we more or less transform the previous integral back to the unit cube, i.e., E(P D ) = p(φ 1 (y)) dy, [0,1] d where Φ is the CDF of the standard normal distribution. Instead of computing E(P D ) precisely, we approximate it by using QMC, E(P D ) 1 N p(φ 1 (y N k )) k=1 for some QMC node set (e.g., low discrepancy sequence). This corresponds to using low discrepancy sequences in the inversion method. Peter Kritzer Quasi-Monte Carlo for Finance 33
34 Challenges using QMC Peter Kritzer Quasi-Monte Carlo for Finance 34
35 Challenge I: High dimensions: Integration problems arising from finance are frequently very high-dimensional. Theoretical error bounds for QMC are typically bad for high dimensions d. Still, the results using QMC for such problems can be quite good. This phenomenon can be explained by the lower effective dimension of the problem (Owen et al., Griebel, Kuo, Sloan,... ), or the introduction of coordinate weights (Sloan, Woźniakowski). Peter Kritzer Quasi-Monte Carlo for Finance 35
36 Challenge II: Mapping to the unit cube: Transforming into E(P D ) = p(x)φ(x) dx R d E(P D ) = p(φ 1 (y)) dy, [0,1] d is technically possible in many cases. However, the integrand p Φ 1 may become very hard to handle, even if p is nice. Peter Kritzer Quasi-Monte Carlo for Finance 36
37 E.g., the new integrand or its derivatives might be unbounded or very large near the boundary of the unit cube. Significant parts of QMC theory do not apply to such integrands. How to overcome this problem? Find assumptions on p and Φ that are practically justifiable and guarantee nice behavior of p Φ 1. Alternative: Come up with QMC theory that works on the whole space R d. No fully satisfying answers yet but much current research (see, e.g., Kuo, Nichols, Sloan, Wasilkowski, Waterhouse,... ). Peter Kritzer Quasi-Monte Carlo for Finance 37
38 Challenge III: Non-smooth functions: An even bigger problem: Sometimes, despite all efforts, the integrand is not in one of the common QMC function spaces. Simple example: payoff of a simple European call option: Payoff S T The function does not satisfy the common smoothness assumptions necessary for QMC. Peter Kritzer Quasi-Monte Carlo for Finance 38
39 What can we do then? 1. Hope that we can by chance find a different approach that works (happens only sometimes, e.g. for European options the Black-Scholes formula). 2. Ignore the problem and apply QMC nevertheless. 3. Hope that the integrand has bounded variation, and apply the Koksma-Hlawka inequality: f (x) dx 1 N f (x k ) [0,1] N V HK (f )D N (x 1,..., x N ), s k=1 where D N (x 1,..., x N ) only depends on the choice of x 1,..., x N, and V HK (f ) is the variation of f in the sense of Hardy and Krause. If V HK (f ) <, then the integration error vanishes as N. Peter Kritzer Quasi-Monte Carlo for Finance 39
40 Further challenge popping up: for many, even nice, f, V HK (f ) =, or one cannot easily check. Recent research: Find other, easier to check kinds of variation and corresponding Koksma-Hlawka-type inequalities (recently Harman, Pausinger, Svane, Tichy,... ) Peter Kritzer Quasi-Monte Carlo for Finance 40
41 Challenge IV: Error estimates. Suppose we choose to apply QMC methods after all. How to find effective error estimates? Theoretical ones often not good enough. Randomize the QMC algorithm and apply it repeatedly. Statistics then helps us to find an estimator for the error. Peter Kritzer Quasi-Monte Carlo for Finance 41
42 Most common randomization method: Random shift: Given is a set of QMC integration nodes x 1,..., x N [0, 1) d. Choose a random vector σ [0, 1) d. Use randomized integration nodes x k + σ (mod1), 1 k d. in the QMC rule instead of x 1,..., x N. Apply the QMC rule for M independent realizations of σ, take the average of the M realizations. Peter Kritzer Quasi-Monte Carlo for Finance 42
43 Concluding remarks Peter Kritzer Quasi-Monte Carlo for Finance 43
44 Multi-asset derivatives: For simplicity, we only covered financial products depending on a single asset S. We could look at s price processes S (1),..., S (s). This typically results in a ds-dimensional problem. A second covariance matrix (describing dependencies among the assets) is involved. They can be factorized using different methods. Arranging the matrices in a clever way can help in making computations faster. Peter Kritzer Quasi-Monte Carlo for Finance 44
45 Why does QMC actually work? There are not many situations in which QMC methods are guaranteed to work well theoretically. The reason is that in practice most integrands are not smooth, while QMC theory usually needs some kind of smoothness assumptions to work well. The empirical success of QMC for large classes of problems is subject of current research (see, e.g., Kuo, Sloan and their co-authors). There is no particular QMC rule that works for all (or even most) problems. Comparison of different methods is necessary. Peter Kritzer Quasi-Monte Carlo for Finance 45
46 Thank you very much for your attention. Peter Kritzer Quasi-Monte Carlo for Finance 46
47 Useful literature: P. Glasserman. Monte Carlo Methods in Financial Engineering. Springer, New York, C. Lemieux. Monte Carlo and Quasi-Monte Carlo Sampling. Springer, New York, G. Leobacher, F. Pillichshammer. Introduction to Quasi-Monte Carlo Integration and Applications. Birkhäuser, Cham, Peter Kritzer Quasi-Monte Carlo for Finance 47
Quasi-Monte Carlo Methods in Financial Engineering: An Equivalence Principle and Dimension Reduction
Quasi-Monte Carlo Methods in Financial Engineering: An Equivalence Principle and Dimension Reduction Xiaoqun Wang,2, and Ian H. Sloan 2,3 Department of Mathematical Sciences, Tsinghua University, Beijing
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 informationAn analysis of faster convergence in certain finance applications for quasi-monte Carlo
An analysis of faster convergence in certain finance applications for quasi-monte Carlo a,b a School of Mathematics and Statistics, University of NSW, Australia b Department of Computer Science, K.U.Leuven,
More informationQuasi-Monte Carlo for Finance Applications
Quasi-Monte Carlo for Finance Applications M.B. Giles F.Y. Kuo I.H. Sloan B.J. Waterhouse October 2008 Abstract Monte Carlo methods are used extensively in computational finance to estimate the price of
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationQuasi-Monte Carlo for finance applications
ANZIAM J. 50 (CTAC2008) pp.c308 C323, 2008 C308 Quasi-Monte Carlo for finance applications M. B. Giles 1 F. Y. Kuo 2 I. H. Sloan 3 B. J. Waterhouse 4 (Received 14 August 2008; revised 24 October 2008)
More informationValuation of performance-dependent options in a Black- Scholes framework
Valuation of performance-dependent options in a Black- Scholes framework Thomas Gerstner, Markus Holtz Institut für Numerische Simulation, Universität Bonn, Germany Ralf Korn Fachbereich Mathematik, TU
More informationLecture outline. Monte Carlo Methods for Uncertainty Quantification. Importance Sampling. Importance Sampling
Lecture outline Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford KU Leuven Summer School on Uncertainty Quantification Lecture 2: Variance reduction
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 informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
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 informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationMultilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More 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 informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2019 Last Time: Markov Chains We can use Markov chains for density estimation, d p(x) = p(x 1 ) p(x }{{}
More informationCPSC 540: Machine Learning
CPSC 540: Machine Learning Monte Carlo Methods Mark Schmidt University of British Columbia Winter 2018 Last Time: Markov Chains We can use Markov chains for density estimation, p(x) = p(x 1 ) }{{} d p(x
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 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 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 informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
More informationMONTE CARLO EXTENSIONS
MONTE CARLO EXTENSIONS School of Mathematics 2013 OUTLINE 1 REVIEW OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO 3 SUMMARY MONTE CARLO SO FAR... Simple to program
More informationMonte Carlo Simulation of a Two-Factor Stochastic Volatility Model
Monte Carlo Simulation of a Two-Factor Stochastic Volatility Model asymptotic approximation formula for the vanilla European call option price. A class of multi-factor volatility models has been introduced
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More 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 informationGeometric tools for the valuation of performance-dependent options
Computational Finance and its Applications II 161 Geometric tools for the valuation of performance-dependent options T. Gerstner & M. Holtz Institut für Numerische Simulation, Universität Bonn, Germany
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 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 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 informationMath Computational Finance Option pricing using Brownian bridge and Stratified samlping
. Math 623 - Computational Finance Option pricing using Brownian bridge and Stratified samlping Pratik Mehta pbmehta@eden.rutgers.edu Masters of Science in Mathematical Finance Department of Mathematics,
More informationA Matlab Program for Testing Quasi-Monte Carlo Constructions
A Matlab Program for Testing Quasi-Monte Carlo Constructions by Lynne Serré A research paper presented to the University of Waterloo in partial fulfillment of the requirements for the degree of Master
More informationAn Efficient Quasi-Monte Carlo Simulation for Pricing Asian Options under Heston's Model
An Efficient Quasi-Monte Carlo Simulation for Pricing Asian Options under Heston's Model by Kewei Yu A thesis presented to the University of Waterloo in fulfillment of the thesis requirement for the degree
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 informationHints on Some of the Exercises
Hints on Some of the Exercises of the book R. Seydel: Tools for Computational Finance. Springer, 00/004/006/009/01. Preparatory Remarks: Some of the hints suggest ideas that may simplify solving the exercises
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 informationENHANCED QUASI-MONTE CARLO METHODS WITH DIMENSION REDUCTION
Proceedings of the 2002 Winter Simulation Conference E Yücesan, C-H Chen, J L Snowdon, J M Charnes, eds ENHANCED QUASI-MONTE CARLO METHODS WITH DIMENSION REDUCTION Junichi Imai Iwate Prefectural University,
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 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 informationEfficient Deterministic Numerical Simulation of Stochastic Asset-Liability Management Models in Life Insurance
Efficient Deterministic Numerical Simulation of Stochastic Asset-Liability Management Models in Life Insurance Thomas Gerstner, Michael Griebel, Markus Holtz Institute for Numerical Simulation, University
More informationComputational Finance. Computational Finance p. 1
Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy
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 in Finance
Monte Carlo Methods in Finance Peter Jackel JOHN WILEY & SONS, LTD Preface Acknowledgements Mathematical Notation xi xiii xv 1 Introduction 1 2 The Mathematics Behind Monte Carlo Methods 5 2.1 A Few Basic
More informationMath Computational Finance Double barrier option pricing using Quasi Monte Carlo and Brownian Bridge methods
. Math 623 - Computational Finance Double barrier option pricing using Quasi Monte Carlo and Brownian Bridge methods Pratik Mehta pbmehta@eden.rutgers.edu Masters of Science in Mathematical Finance Department
More informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 2-3 Haijun Li An Introduction to Stochastic Calculus Week 2-3 1 / 24 Outline
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 informationBarrier Option. 2 of 33 3/13/2014
FPGA-based Reconfigurable Computing for Pricing Multi-Asset Barrier Options RAHUL SRIDHARAN, GEORGE COOKE, KENNETH HILL, HERMAN LAM, ALAN GEORGE, SAAHPC '12, PROCEEDINGS OF THE 2012 SYMPOSIUM ON APPLICATION
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 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 informationIntroduction to Probability Theory and Stochastic Processes for Finance Lecture Notes
Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Fabio Trojani Department of Economics, University of St. Gallen, Switzerland Correspondence address: Fabio Trojani,
More informationFinancial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA
Financial Risk Modeling on Low-power Accelerators: Experimental Performance Evaluation of TK1 with FPGA Rajesh Bordawekar and Daniel Beece IBM T. J. Watson Research Center 3/17/2015 2014 IBM Corporation
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 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 informationMOUNTAIN RANGE OPTIONS
MOUNTAIN RANGE OPTIONS Paolo Pirruccio Copyright Arkus Financial Services - 2014 Mountain Range options Page 1 Mountain Range options Introduction Originally marketed by Société Générale in 1998. Traded
More information10. Monte Carlo Methods
10. Monte Carlo Methods 1. Introduction. Monte Carlo simulation is an important tool in computational finance. It may be used to evaluate portfolio management rules, to price options, to simulate hedging
More informationNon-semimartingales in finance
Non-semimartingales in finance Pricing and Hedging Options with Quadratic Variation Tommi Sottinen University of Vaasa 1st Northern Triangular Seminar 9-11 March 2009, Helsinki University of Technology
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 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 informationMath Option pricing using Quasi Monte Carlo simulation
. Math 623 - Option pricing using Quasi Monte Carlo simulation Pratik Mehta pbmehta@eden.rutgers.edu Masters of Science in Mathematical Finance Department of Mathematics, Rutgers University This paper
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 Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationMAS3904/MAS8904 Stochastic Financial Modelling
MAS3904/MAS8904 Stochastic Financial Modelling Dr Andrew (Andy) Golightly a.golightly@ncl.ac.uk Semester 1, 2018/19 Administrative Arrangements Lectures on Tuesdays at 14:00 (PERCY G13) and Thursdays at
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 informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationSimulating more interesting stochastic processes
Chapter 7 Simulating more interesting stochastic processes 7. Generating correlated random variables The lectures contained a lot of motivation and pictures. We'll boil everything down to pure algebra
More informationEstimating the Greeks
IEOR E4703: Monte-Carlo Simulation Columbia University Estimating the Greeks c 207 by Martin Haugh In these lecture notes we discuss the use of Monte-Carlo simulation for the estimation of sensitivities
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 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 informationProbability in Options Pricing
Probability in Options Pricing Mark Cohen and Luke Skon Kenyon College cohenmj@kenyon.edu December 14, 2012 Mark Cohen and Luke Skon (Kenyon college) Probability Presentation December 14, 2012 1 / 16 What
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More information1 Geometric Brownian motion
Copyright c 05 by Karl Sigman Geometric Brownian motion Note that since BM can take on negative values, using it directly for modeling stock prices is questionable. There are other reasons too why BM is
More informationStochastic Calculus, Application of Real Analysis in Finance
, Application of Real Analysis in Finance Workshop for Young Mathematicians in Korea Seungkyu Lee Pohang University of Science and Technology August 4th, 2010 Contents 1 BINOMIAL ASSET PRICING MODEL Contents
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 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 informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
More informationFINANCIAL OPTION ANALYSIS HANDOUTS
FINANCIAL OPTION ANALYSIS HANDOUTS 1 2 FAIR PRICING There is a market for an object called S. The prevailing price today is S 0 = 100. At this price the object S can be bought or sold by anyone for any
More informationAsymptotic methods in risk management. Advances in Financial Mathematics
Asymptotic methods in risk management Peter Tankov Based on joint work with A. Gulisashvili Advances in Financial Mathematics Paris, January 7 10, 2014 Peter Tankov (Université Paris Diderot) Asymptotic
More information[D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright
Faculty and Institute of Actuaries Claims Reserving Manual v.2 (09/1997) Section D7 [D7] PROBABILITY DISTRIBUTION OF OUTSTANDING LIABILITY FROM INDIVIDUAL PAYMENTS DATA Contributed by T S Wright 1. Introduction
More information4: SINGLE-PERIOD MARKET MODELS
4: SINGLE-PERIOD MARKET MODELS Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 4: Single-Period Market Models 1 / 87 General Single-Period
More informationMulti-Asset Options. A Numerical Study VILHELM NIKLASSON FRIDA TIVEDAL. Master s thesis in Engineering Mathematics and Computational Science
Multi-Asset Options A Numerical Study Master s thesis in Engineering Mathematics and Computational Science VILHELM NIKLASSON FRIDA TIVEDAL Department of Mathematical Sciences Chalmers University of Technology
More informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
More 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 informationNo-arbitrage theorem for multi-factor uncertain stock model with floating interest rate
Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer
More informationHomework Assignments
Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)
More informationAn Adaptive Method for Evaluating Multidimensional Contingent Claims. Part II
Dept. of Math. University of Oslo Pure Mathematics ISBN 82 553 1342 7 No. 10 ISSN 0806 2439 May 2002 An Adaptive Method for Evaluating Multidimensional Contingent Claims. Part II Lars O. Dahl 28th February
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More 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 informationNumerical Methods in Option Pricing (Part III)
Numerical Methods in Option Pricing (Part III) E. Explicit Finite Differences. Use of the Forward, Central, and Symmetric Central a. In order to obtain an explicit solution for the price of the derivative,
More informationRandomness and Fractals
Randomness and Fractals Why do so many physicists become traders? Gregory F. Lawler Department of Mathematics Department of Statistics University of Chicago September 25, 2011 1 / 24 Mathematics and the
More informationModelling the Sharpe ratio for investment strategies
Modelling the Sharpe ratio for investment strategies Group 6 Sako Arts 0776148 Rik Coenders 0777004 Stefan Luijten 0783116 Ivo van Heck 0775551 Rik Hagelaars 0789883 Stephan van Driel 0858182 Ellen Cardinaels
More informationUQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.
UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.
More informationEMPIRICAL EVIDENCE ON ARBITRAGE BY CHANGING THE STOCK EXCHANGE
Advances and Applications in Statistics Volume, Number, This paper is available online at http://www.pphmj.com 9 Pushpa Publishing House EMPIRICAL EVIDENCE ON ARBITRAGE BY CHANGING THE STOCK EXCHANGE JOSÉ
More informationSYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives
SYSM 6304: Risk and Decision Analysis Lecture 6: Pricing and Hedging Financial Derivatives M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu October
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationToward a coherent Monte Carlo simulation of CVA
Toward a coherent Monte Carlo simulation of CVA Lokman Abbas-Turki (Joint work with A. I. Bouselmi & M. A. Mikou) TU Berlin January 9, 2013 Lokman (TU Berlin) Advances in Mathematical Finance 1 / 16 Plan
More 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 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 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 informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More informationProbability. An intro for calculus students P= Figure 1: A normal integral
Probability An intro for calculus students.8.6.4.2 P=.87 2 3 4 Figure : A normal integral Suppose we flip a coin 2 times; what is the probability that we get more than 2 heads? Suppose we roll a six-sided
More informationAn Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions
Journal of Numerical Mathematics and Stochastics,1 (1) : 45-55, 2009 http://www.jnmas.org/jnmas1-5.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 An Efficient Numerical Scheme for Simulation of
More informationON AN IMPLEMENTATION OF BLACK SCHOLES MODEL FOR ESTIMATION OF CALL- AND PUT-OPTION VIA PROGRAMMING ENVIRONMENT MATHEMATICA
Доклади на Българската академия на науките Comptes rendus de l Académie bulgare des Sciences Tome 66, No 5, 2013 MATHEMATIQUES Mathématiques appliquées ON AN IMPLEMENTATION OF BLACK SCHOLES MODEL FOR ESTIMATION
More information