Guaranteed Fixed-Width Confidence Intervals for Monte Carlo and Quasi-Monte Carlo Simulation
|
|
- Caren Griffin
- 5 years ago
- Views:
Transcription
1 Guaranteed Fixed-Width Confidence Intervals for Monte Carlo and Quasi-Monte Carlo Simulation Fred J. Hickernell Department of Applied Mathematics, Illinois Institute of Technology Joint work with Sou-Cheng Choi (NORC, U Chicago), Lan Jiang (IIT), Lluís Antoni Jiménez Rugama (IIT), and Art Owen (Stanford) Supported by NSF-DMS Thank you for your kind invitation to speak. It is a privilege to visit the country of my forefather, who came from Haguenau in Alsace June 1, 215 hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 1 / 15
2 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ^µ n := 1 n nÿ Y i, Y 1, Y 2,... Y Monte Carlo = sampling with computers originated at Los Alamos National Laboratory ca (Eckhardt, 1987) µ = E(Y) =?, where Y complicated» Y(4) %random # generator ans = hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 2 / 15
3 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ^µ n := 1 n nÿ Y i, Y 1, Y 2,... Y Monte Carlo = sampling with computers Given ε a, want n such that P[ µ ^µ n ď ε a ] ě 99% a fixed-width confidence interval µ = E(Y) =?, where Y complicated» Y(4) %random # generator ans = hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 2 / 15
4 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n nÿ Y i, Y 1, Y 2,... Y Monte Carlo = sampling with computers Given ε a, want n such that P[ µ ^µ n ď ε a ] ě 99% Central Limit Theorem: ^µ n «N(µ, σ 2 /n), suggesting S (2.58 ) W 2 ˆ 1.2 ˆ ^σ n = ^σ 2 ε a = sample variance option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray µ = E(Y) =?, where Y complicated» Y(4) %random # generator ans = hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 2 / 15
5 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n nÿ Y i, Y 1, Y 2,... Y Monte Carlo = sampling with computers Given ε a, want n such that P[ µ ^µ n ď ε a ] ě 99% Central Limit Theorem: ^µ n «N(µ, σ 2 /n), suggesting S (2.58 ) W 2 ˆ 1.2 ˆ ^σ n = ^σ 2 ε a = sample variance option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray µ = E(Y) =?, where Y complicated» Y(4) %random # generator ans = » meanmc_clt(y,.2) ans = 1.1 +/-.2» meanmc_clt(y,.2) ans = /-.2» meanmc_clt(y,.2) ans =.99 +/-.2 What is wrong? hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 2 / 15
6 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n nÿ Y i, Y 1, Y 2,... Y Monte Carlo = sampling with computers Given ε a, want n such that P[ µ ^µ n ď ε a ] ě 99% Central Limit Theorem: ^µ n «N(µ, σ 2 /n), suggesting S (2.58 ) W 2 ˆ 1.2 ˆ ^σ n = ε a ^σ 2 = sample variance using n σ samples option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray µ = E(Y) =?, where Y complicated» Y(4) %random # generator ans = » meanmc_clt(y,.2,.1,1) ans = 1. +/-.2» meanmc_clt(y,.2,.1,1) ans = 1. +/-.2» meanmc_clt(y,.2,.1,1) ans = 1. +/-.2 for n σ = 1, rather than 1 hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 2 / 15
7 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n nÿ Y i, Y 1, Y 2,... Y Monte Carlo = sampling with computers Given ε a, want n such that P[ µ ^µ n ď ε a ] ě 99% Central Limit Theorem: ^µ n «N(µ, σ 2 /n), suggesting S (2.58 ) W 2 ˆ 1.2 ˆ ^σ n = ^σ 2 ε a = sample variance option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray µ = E(Y) =?, where Y complicated» Y(4) %random # generator ans = PDF of Y y hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 2 / 15
8 Guaranteed Adaptive Monte Carlo Simulation Adaptive Algorithm meanmc_g. (Choi et al., ) Given i) IID random generator for Y 1, Y 2,..., ii) absolute error tolerance, ε a, and iii) n σ P N and C ą 1, Step 1. Bound the Variance. Compute the sample variance, ^σ 2, based on n σ samples of Y. We know P[C 2^σ 2 ě var(y)] ě 99.5% for kurt(y) := E[(Y µ) 4 ]/ var(y) 2 ď κ max (n σ, C) by Cantelli s inequality. Step 2. Estimate the Mean. Use a Berry-Esseen Inequality to determine the sample size, n, needed for the sample mean by solving Φ (?nε a /(C^σ) ) loooooooooomoooooooooon + n (? loooooooooooomoooooooooooon nε a /(C^σ), κ max ) ď.25. CLT part Berry-Esseen extra part Compute the sample mean, ^µ n, of an independent sample of size n. Theorem. (H. et al., 214) Assuming ď, we must have P[ µ ^µ n ď ε a ] ě 99%. The computational cost is n σ + n = O(var(Y)/ε 2 a) with high probability, even though var(y) is unknown a priori. hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 3 / 15
9 Examples for meanmc_g Asian Geometric Mean Call Option d = 1, 2, 4,..., 64 time steps Failure.8 c d ż e x 2 R d cos( x ) dx = 1 d = 1, 2, 3, 4 (Keister, 1996) Failure.8 Time (seconds) Time (seconds) Success Error 1 3 Success Error «1% success «1% success Need n σ ě 1 3 to get a reasonable κ max meanmc_g is conservative, see Bayer et al. (214) for a heuristic alternative 1 3 hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 4 / 15
10 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d.75 Given ε a, want n such that µ ^µ n ď ε a a fixed-width confidence interval We may sample the X i as IID hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
11 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
12 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
13 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
14 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
15 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
16 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
17 Stopping the Simulation When the Error is Small Enough option price probability of process failure pixel intensity ^µ n := 1 n option payoff under random scenario random failure or not, i.e., 1 or intensity of random incident light ray ż µ = E(Y) = f (x) dx =?, where Y = f (X), X U[, 1] d [,1] d nÿ Y i = 1 nÿ 1 f (X i ), n X 1, X 2,... U[, 1] d Given ε a, want n such that µ ^µ n ď ε a We may sample the X i as IID Or by highly stratified sampling, also known as quasi-monte Carlo sampling or low discrepancy sampling hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 5 / 15
18 Highly Stratified Sampling ż f (x) dx 1 nÿ f (X i ) loooooomoooooon [,1) n ď ε a, d µ=e[f (X)] highly stratified n = 2 m n =? samples needed X 1, X 2,... dependent, can be deterministic or random (Owen, 2) Sobol sequences are a popular choice (Dick and Pillichshammer, 21) Typical error analysis requires L p -norms of mixed partial derivatives of f (Niederreiter, 1992; H., 1998) hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 6 / 15
19 Highly Stratified Sampling ż f (x) dx 1 nÿ f (X i ) loooooomoooooon [,1) n ď d µ=e[f (X)] highly stratified ÿ k(κ)pdual net p fκ n = 2 m Walsh functions & coefficients f (x) = ÿ κ= ( 1) xk(κ),xy p fκ hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 6 / 15
20 Monte Carlo Simulation Adaptive Monte Carlo Adaptive Quasi-Monte Carlo Discussion References Highly Stratified Sampling n ÿ 1ÿ pfκ ď C(m)Sm (f ) ď εa, f (X i ) ď f (x) dx n d [,1) loooooomoooooon k(κ)pdual net ż µ=e[f (X)] highly stratified 1 Walsh functions & coefficients f (x) = ÿ n = 2m ( 1)xk(κ),xypfκ 1 5 Sm (f ) = data based sum of approximate p fκ w/ high κ H. and Jiménez Rugama (214) construct an algorithm that ^ n ď εa guarantees µ µ fˆκ κ= 1 1 error bound S12 (f) κ conditions: p fκ does not decay erratically as κ Ñ hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 6 / 15
21 Monte Carlo Simulation Adaptive Monte Carlo Adaptive Quasi-Monte Carlo Discussion References When Does Our Error Bound Work? May Fail f (x) f (x) Succeeds x x fˆκ 1 fˆκ error bound S12 (f) κ hickernell@iit.edu κ Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 7 / 15
22 Examples for meanmc_g and cubsobol_g Asian Geometric Mean Call Option, d = 1, 2, 4,..., 64 time steps meanmc_g Failure cubsobol_g Failure Time (seconds) Success Error Time (seconds) Success Error «1% success «1% success hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 8 / 15
23 Examples for meanmc_g and cubsobol_g c d ż e x 2 R d cos( x ) dx = 1 (Keister, 1996) d = 1,..., 4 d = 1,..., 4 meanmc_g Failure cubsobol_g Failure Time (seconds) Time (seconds) Success Error 1 3 Success Error «1% success «1% success 1 3 hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 9 / 15
24 Examples for meanmc_g and cubsobol_g c d ż e x 2 R d cos( x ) dx = 1 (Keister, 1996) d = 1,..., 4 d = 1,..., 2 meanmc_g Failure cubsobol_g Failure Time (seconds) Time (seconds) Success Error 1 3 Success Error «1% success «97% success 1 3 hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 9 / 15
25 Observations on meanmc_g and cubsobol_g Highly stratified sampling (quasi-monte Carlo) can provide great efficiency gains. cublattice_g uses a integration lattices to achieve a high degree of stratification (Jiménez Rugama and H., 214). Y or f must lie in, which describes how nasty they are allowed to be. Bigger =ñ more robust algorithms with greater cost. One may not be able to verify if Y or f lies in a priori, but the theorems explain what might go wrong and provide guidance on the choice of algorithm parameters. Focus on (instead of other shapes) of Y or f because our " * " * problems are homogeneous. error bounds positively hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 1 / 15
26 Observations on meanmc_g and cubsobol_g cont d Exist general sufficient conditions under which adaption provides no advantage (Bahadur and Savage, 1956 ; Bakhvalov, 197 ; Traub et al., 1988, Chapter 4, Theorem ; Novak, 1996). To violate those conditions we consider nonconvex of Y or f. Can accommodate µ ^µ n ď max(ε a, ε r µ ). These algorithms featured in our Guaranteed Automatic Integration Library (GAIL) code.google.com/p/gail/ (Choi et al., ). >> Y =... %define your random # generator Y >> tic, muiid = meanmc_g(y,.1,), toc muiid = 1.2 Elapsed time is seconds. >> f =... %define your integrand f >> tic, musobol = cubsobol_g(f,... ), toc musobol = 1.3 Elapsed time is.3778 seconds. hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 11 / 15
27 Monte Carlo Simulation Adaptive Monte Carlo Adaptive Quasi-Monte Carlo Discussion References Guaranteed Automatic Integration Library (GAIL) Version 2.1 code.google.com/p/gail/ GAIL 3. coming this fall. Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 12 / 15
28 Future Work In progress Control variates (subtle for quasi-monte Carlo (H. et al., 25)) Adaptive importance sampling Multilevel and related methods for -dimensional problems Porting to R and other languages Need further expertise or time Y P t, 1u, i.e., Bernoulli random variables (Jiang and H., 214) Markov Chain Monte Carlo Variety of use cases, e.g., Bayesian inference hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 13 / 15
29 References I Bahadur, R. R. and L. J. Savage The nonexistence of certain statistical procedures in nonparametric problems, Ann. Math. Stat. 27, Bakhvalov, N. S On the optimality of linear methods for operator approximation in convex classes of functions (in Russian), Zh. Vychisl. Mat. i Mat. Fiz. 1, English transl.: USSR Comput. Math. Math. Phys. 11 (1971) Bayer, C., H. Hoel, E. von Schwerin, and R. Tempone On nonasymptotic optimal stopping criteria in Monte Carlo Simulations, SIAM J. Sci. Comput. 36, A869 A885. Choi, S.-C. T., Y. Ding, F. J. H., L. Jiang, Ll. A. Jiménez Rugama, X. Tong, Y. Zhang, and X. Zhou GAIL: Guaranteed Automatic Integration Library (versions ). Dick, J. and F. Pillichshammer. 21. Digital nets and sequences: Discrepancy theory and quasi-monte Carlo integration, Cambridge University Press, Cambridge. Eckhardt, R Stan Ulam, John von Neumann, and the Monte Carlo method, Los Alamos Science, H., F. J A generalized discrepancy and quadrature error bound, Math. Comp. 67, H., F. J., L. Jiang, Y. Liu, and A. B. Owen Guaranteed conservative fixed width confidence intervals via Monte Carlo sampling, Monte Carlo and quasi-monte Carlo methods 212, pp hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 14 / 15
30 References II H., F. J. and Ll. A. Jiménez Rugama Reliable adaptive cubature using digital sequences. submitted for publication, arxiv: [math.na]. H., F. J., C. Lemieux, and A. B. Owen. 25. Control variates for quasi-monte Carlo, Statist. Sci. 2, Jiang, L. and F. J. H Guaranteed conservative confidence intervals for means of Bernoulli random variables. submitted for publication, arxiv: Jiménez Rugama, Ll. A. and F. J. H Adaptive multidimensional integration based on rank-1 lattices. submitted for publication, arxiv: Keister, B. D Multidimensional quadrature algorithms, Computers in Physics 1, Niederreiter, H Random number generation and quasi-monte Carlo methods, CBMS-NSF Regional Conference Series in Applied Mathematics, SIAM, Philadelphia. Novak, E On the power of adaption, J. Complexity 12, Owen, A. B. 2. Monte Carlo, quasi-monte Carlo, and randomized quasi-monte Carlo, Monte Carlo, quasi-monte Carlo, and randomized quasi-monte Carlo, pp Traub, J. F., G. W. Wasilkowski, and H. Woźniakowski Information-based complexity, Academic Press, Boston. hickernell@iit.edu Fixed-Width Monte Carlo Confidence Intervals Journées de Statistique, 6/1/15 15 / 15
Reliable Error Estimation for Monte Carlo and Quasi-Monte Carlo Simulation When do you have enough samples to stop?
Reliable Error Estimation for Monte Carlo and Quasi-Monte Carlo Simulation When do you have enough samples to stop? Fred J. Hickernell Department of Applied Mathematics, Illinois Institute of Technology
More informationAdaptive Quasi-Monte Carlo Where Each Integrand Value Depends on All Data Sites: the American Option
Adaptive Quasi-Monte Carlo Where Each Integrand Value Depends on All Data Sites: the American Option Lluís Antoni Jiménez Rugama Joint work with Fred J. Hickernell (IIT) and Da Li (IIT) Room 12, Bldg E1,
More informationGuaranteed Conservative Fixed Width Confidence Intervals Via Monte Carlo Sampling
Guaranteed Conservative Fixed Width Confidence Intervals Via Monte Carlo Sampling Fred J. Hickernell, Lan Jiang, Yuewei Liu 2, and Art Owen 3 Department of Applied Mathematics, Illinois Institute of Technology,
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 informationEFFICIENCY IMPROVEMENT BY LATTICE RULES FOR PRICING ASIAN OPTIONS. Christiane Lemieux Pierre L Ecuyer
Proceedings of the 1998 Winter Simulation Conference D.J. Medeiros, E.F. Watson, J.S. Carson and M.S. Manivannan, eds. EFFICIENCY IMPROVEMENT BY LATTICE RULES FOR PRICING ASIAN OPTIONS Christiane Lemieux
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationQuasi-Monte Carlo for Finance
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
More informationX i. i=1. We consider the objective of choosing M sufficiently large so that the error probability satisfies
ON NON-ASYMPTOTIC OPTIMAL STOPPING CRITERIA IN MONTE CARLO SIMULATIONS CHRISTIAN BAYER, HÅKON HOEL, ERIK VON SCHWERIN, AND RAÚL TEMPONE Abstract. We consider the setting of estimating the mean of a random
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 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 informationConfidence Intervals Introduction
Confidence Intervals Introduction A point estimate provides no information about the precision and reliability of estimation. For example, the sample mean X is a point estimate of the population mean μ
More informationContents Critique 26. portfolio optimization 32
Contents Preface vii 1 Financial problems and numerical methods 3 1.1 MATLAB environment 4 1.1.1 Why MATLAB? 5 1.2 Fixed-income securities: analysis and portfolio immunization 6 1.2.1 Basic valuation of
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 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 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 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 informationUniversity of California Berkeley
University of California Berkeley Improving the Asmussen-Kroese Type Simulation Estimators Samim Ghamami and Sheldon M. Ross May 25, 2012 Abstract Asmussen-Kroese [1] Monte Carlo estimators of P (S n >
More informationMonte Carlo Simulations in the Teaching Process
Monte Carlo Simulations in the Teaching Process Blanka Šedivá Department of Mathematics, Faculty of Applied Sciences University of West Bohemia, Plzeň, Czech Republic CADGME 2018 Conference on Digital
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 and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More 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 informationChapter 7: Point Estimation and Sampling Distributions
Chapter 7: Point Estimation and Sampling Distributions Seungchul Baek Department of Statistics, University of South Carolina STAT 509: Statistics for Engineers 1 / 20 Motivation In chapter 3, we learned
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 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 informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
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 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 information3. Monte Carlo Simulation
3. Monte Carlo Simulation 3.7 Variance Reduction Techniques Math443 W08, HM Zhu Variance Reduction Procedures (Chap 4.5., 4.5.3, Brandimarte) Usually, a very large value of M is needed to estimate V with
More informationOn the Use of Quasi-Monte Carlo Methods in Computational Finance
On the Use of Quasi-Monte Carlo Methods in Computational Finance Christiane Lemieux 1 and Pierre L Ecuyer 2 1 Department of Mathematics and Statistics, University of Calgary, 2500 University Drive N.W.,
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 informationWeek 1 Quantitative Analysis of Financial Markets Distributions B
Week 1 Quantitative Analysis of Financial Markets Distributions B Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036 October
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 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 informationParallel Multilevel Monte Carlo Simulation
Parallel Simulation Mathematisches Institut Goethe-Universität Frankfurt am Main Advances in Financial Mathematics Paris January 7-10, 2014 Simulation Outline 1 Monte Carlo 2 3 4 Algorithm Numerical Results
More informationSeminar: Efficient Monte Carlo Methods for Uncertainty Quantification
Seminar: Efficient Monte Carlo Methods for Uncertainty Quantification Elisabeth Ullmann Lehrstuhl für Numerische Mathematik (M2) TU München Elisabeth Ullmann (TU München) Efficient Monte Carlo for UQ 1
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationST440/550: Applied Bayesian Analysis. (5) Multi-parameter models - Summarizing the posterior
(5) Multi-parameter models - Summarizing the posterior Models with more than one parameter Thus far we have studied single-parameter models, but most analyses have several parameters For example, consider
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 informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
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 informationResults for option pricing
Results for option pricing [o,v,b]=optimal(rand(1,100000 Estimators = 0.4619 0.4617 0.4618 0.4613 0.4619 o = 0.46151 % best linear combination (true value=0.46150 v = 1.1183e-005 %variance per uniform
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 informationChapter 7: Estimation Sections
1 / 40 Chapter 7: Estimation Sections 7.1 Statistical Inference Bayesian Methods: Chapter 7 7.2 Prior and Posterior Distributions 7.3 Conjugate Prior Distributions 7.4 Bayes Estimators Frequentist Methods:
More informationKARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI
88 P a g e B S ( B B A ) S y l l a b u s KARACHI UNIVERSITY BUSINESS SCHOOL UNIVERSITY OF KARACHI BS (BBA) VI Course Title : STATISTICS Course Number : BA(BS) 532 Credit Hours : 03 Course 1. Statistical
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 informationMATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics
MATHICSE Mathematics Institute of Computational Science and Engineering School of Basic Sciences - Section of Mathematics MATHICSE Technical Report Nr. 07.2013 February 2013 On non-asymptotic optimal stopping
More informationInterval estimation. September 29, Outline Basic ideas Sampling variation and CLT Interval estimation using X More general problems
Interval estimation September 29, 2017 STAT 151 Class 7 Slide 1 Outline of Topics 1 Basic ideas 2 Sampling variation and CLT 3 Interval estimation using X 4 More general problems STAT 151 Class 7 Slide
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 informationA Heuristic Method for Statistical Digital Circuit Sizing
A Heuristic Method for Statistical Digital Circuit Sizing Stephen Boyd Seung-Jean Kim Dinesh Patil Mark Horowitz Microlithography 06 2/23/06 Statistical variation in digital circuits growing in importance
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 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 informationSTAT Chapter 7: Confidence Intervals
STAT 515 -- Chapter 7: Confidence Intervals With a point estimate, we used a single number to estimate a parameter. We can also use a set of numbers to serve as reasonable estimates for the parameter.
More informationChapter 5. Statistical inference for Parametric Models
Chapter 5. Statistical inference for Parametric Models Outline Overview Parameter estimation Method of moments How good are method of moments estimates? Interval estimation Statistical Inference for Parametric
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 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 informationMuch of what appears here comes from ideas presented in the book:
Chapter 11 Robust statistical methods Much of what appears here comes from ideas presented in the book: Huber, Peter J. (1981), Robust statistics, John Wiley & Sons (New York; Chichester). There are many
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 information(5) Multi-parameter models - Summarizing the posterior
(5) Multi-parameter models - Summarizing the posterior Spring, 2017 Models with more than one parameter Thus far we have studied single-parameter models, but most analyses have several parameters For example,
More informationSample Size for Assessing Agreement between Two Methods of Measurement by Bland Altman Method
Meng-Jie Lu 1 / Wei-Hua Zhong 1 / Yu-Xiu Liu 1 / Hua-Zhang Miao 1 / Yong-Chang Li 1 / Mu-Huo Ji 2 Sample Size for Assessing Agreement between Two Methods of Measurement by Bland Altman Method Abstract:
More informationOn the Scrambled Sobol sequences Lecture Notes in Computer Science 3516, , Springer 2005
On the Scrambled Sobol sequences Lecture Notes in Computer Science 3516, 775-782, Springer 2005 On the Scrambled Soboĺ Sequence Hongmei Chi 1, Peter Beerli 2, Deidre W. Evans 1, and Micheal Mascagni 2
More informationarxiv: v1 [math.st] 21 Mar 2016
Stratified Monte Carlo simulation of Markov chains Rana Fakhereddine a, Rami El Haddad a, Christian Lécot b, arxiv:1603.06386v1 [math.st] 21 Mar 2016 a Université Saint-Joseph, Faculté des Sciences, BP
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More information6 Central Limit Theorem. (Chs 6.4, 6.5)
6 Central Limit Theorem (Chs 6.4, 6.5) Motivating Example In the next few weeks, we will be focusing on making statistical inference about the true mean of a population by using sample datasets. Examples?
More informationModule 4: Point Estimation Statistics (OA3102)
Module 4: Point Estimation Statistics (OA3102) Professor Ron Fricker Naval Postgraduate School Monterey, California Reading assignment: WM&S chapter 8.1-8.4 Revision: 1-12 1 Goals for this Module Define
More informationPricing Exotic Options Under a Higher-order Hidden Markov Model
Pricing Exotic Options Under a Higher-order Hidden Markov Model Wai-Ki Ching Tak-Kuen Siu Li-min Li 26 Jan. 2007 Abstract In this paper, we consider the pricing of exotic options when the price dynamic
More informationIntroduction to Sequential Monte Carlo Methods
Introduction to Sequential Monte Carlo Methods Arnaud Doucet NCSU, October 2008 Arnaud Doucet () Introduction to SMC NCSU, October 2008 1 / 36 Preliminary Remarks Sequential Monte Carlo (SMC) are a set
More informationLecture Quantitative Finance Spring Term 2015
implied Lecture Quantitative Finance Spring Term 2015 : May 7, 2015 1 / 28 implied 1 implied 2 / 28 Motivation and setup implied the goal of this chapter is to treat the implied which requires an algorithm
More informationCHAPTER 2 LITERATURE REVIEW
CHAPTER 2 LITERATURE REVIEW This chapter includes the review of research work (research papers) in short. The papers reviewed illustrate use of Monte Carlo methods in diverse fields. Some papers are also
More informationStochastic Approximation Algorithms and Applications
Harold J. Kushner G. George Yin Stochastic Approximation Algorithms and Applications With 24 Figures Springer Contents Preface and Introduction xiii 1 Introduction: Applications and Issues 1 1.0 Outline
More informationSection 8.2: Monte Carlo Estimation
Section 8.2: Monte Carlo Estimation Discrete-Event Simulation: A First Course c 2006 Pearson Ed., Inc. 0-13-142917-5 Discrete-Event Simulation: A First Course Section 8.2: Monte Carlo Estimation 1/ 19
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 informationFREDRIK BAJERS VEJ 7 G 9220 AALBORG ØST Tlf.: URL: Fax: Monte Carlo methods
INSTITUT FOR MATEMATISKE FAG AALBORG UNIVERSITET FREDRIK BAJERS VEJ 7 G 9220 AALBORG ØST Tlf.: 96 35 88 63 URL: www.math.auc.dk Fax: 98 15 81 29 E-mail: jm@math.aau.dk Monte Carlo methods Monte Carlo methods
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 informationUsing Halton Sequences. in Random Parameters Logit Models
Journal of Statistical and Econometric Methods, vol.5, no.1, 2016, 59-86 ISSN: 1792-6602 (print), 1792-6939 (online) Scienpress Ltd, 2016 Using Halton Sequences in Random Parameters Logit Models Tong Zeng
More informationResearch Statement. Dapeng Zhan
Research Statement Dapeng Zhan The Schramm-Loewner evolution (SLE), first introduced by Oded Schramm ([12]), is a oneparameter (κ (0, )) family of random non-self-crossing curves, which has received a
More informationA NEW EFFICIENT SIMULATION STRATEGY FOR PRICING PATH-DEPENDENT OPTIONS
Proceedings of the 2006 Winter Simulation Conference L. F. Perrone, F. P. Wieland, J. Liu, B. G. Lawson, D. M. Nicol, and R. M. Fujimoto, eds. A NEW EFFICIENT SIMULATION STRATEGY FOR PRICING PATH-DEPENDENT
More informationConfidence Intervals for the Difference Between Two Means with Tolerance Probability
Chapter 47 Confidence Intervals for the Difference Between Two Means with Tolerance Probability Introduction This procedure calculates the sample size necessary to achieve a specified distance from the
More informationMATH 5510 Mathematical Models of Financial Derivatives. Topic 1 Risk neutral pricing principles under single-period securities models
MATH 5510 Mathematical Models of Financial Derivatives Topic 1 Risk neutral pricing principles under single-period securities models 1.1 Law of one price and Arrow securities 1.2 No-arbitrage theory and
More informationRandom Tree Method. Monte Carlo Methods in Financial Engineering
Random Tree Method Monte Carlo Methods in Financial Engineering What is it for? solve full optimal stopping problem & estimate value of the American option simulate paths of underlying Markov chain produces
More informationThe Assumption(s) of Normality
The Assumption(s) of Normality Copyright 2000, 2011, 2016, J. Toby Mordkoff This is very complicated, so I ll provide two versions. At a minimum, you should know the short one. It would be great if you
More informationThe Binomial Lattice Model for Stocks: Introduction to Option Pricing
1/33 The Binomial Lattice Model for Stocks: Introduction to Option Pricing Professor Karl Sigman Columbia University Dept. IEOR New York City USA 2/33 Outline The Binomial Lattice Model (BLM) as a Model
More informationSTAT 509: Statistics for Engineers Dr. Dewei Wang. Copyright 2014 John Wiley & Sons, Inc. All rights reserved.
STAT 509: Statistics for Engineers Dr. Dewei Wang Applied Statistics and Probability for Engineers Sixth Edition Douglas C. Montgomery George C. Runger 7 Point CHAPTER OUTLINE 7-1 Point Estimation 7-2
More informationB. Consider the problem of evaluating the one dimensional integral
Proceedings of the 2008 Winter Simulation Conference S. J. Mason, R. R. Hill, L. Mönch, O. Rose, T. Jefferson, J. W. Fowler eds. MONOTONICITY AND STRATIFICATION Gang Zhao Division of Systems Engineering
More informationIntroduction to Algorithmic Trading Strategies Lecture 8
Introduction to Algorithmic Trading Strategies Lecture 8 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Further Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Outline
More informationBusiness Statistics 41000: Probability 4
Business Statistics 41000: Probability 4 Drew D. Creal University of Chicago, Booth School of Business February 14 and 15, 2014 1 Class information Drew D. Creal Email: dcreal@chicagobooth.edu Office:
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 informationApproximations of Stochastic Programs. Scenario Tree Reduction and Construction
Approximations of Stochastic Programs. Scenario Tree Reduction and Construction W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany www.mathematik.hu-berlin.de/~romisch
More informationInferences on Correlation Coefficients of Bivariate Log-normal Distributions
Inferences on Correlation Coefficients of Bivariate Log-normal Distributions Guoyi Zhang 1 and Zhongxue Chen 2 Abstract This article considers inference on correlation coefficients of bivariate log-normal
More informationA Macro-Finance Model of the Term Structure: the Case for a Quadratic Yield Model
Title page Outline A Macro-Finance Model of the Term Structure: the Case for a 21, June Czech National Bank Structure of the presentation Title page Outline Structure of the presentation: Model Formulation
More informationThe value of foresight
Philip Ernst Department of Statistics, Rice University Support from NSF-DMS-1811936 (co-pi F. Viens) and ONR-N00014-18-1-2192 gratefully acknowledged. IMA Financial and Economic Applications June 11, 2018
More informationUNIT 4 MATHEMATICAL METHODS
UNIT 4 MATHEMATICAL METHODS PROBABILITY Section 1: Introductory Probability Basic Probability Facts Probabilities of Simple Events Overview of Set Language Venn Diagrams Probabilities of Compound Events
More informationOvernight Index Rate: Model, calibration and simulation
Research Article Overnight Index Rate: Model, calibration and simulation Olga Yashkir and Yuri Yashkir Cogent Economics & Finance (2014), 2: 936955 Page 1 of 11 Research Article Overnight Index Rate: Model,
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulation Efficiency and an Introduction to Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University
More informationChapter 8 Statistical Intervals for a Single Sample
Chapter 8 Statistical Intervals for a Single Sample Part 1: Confidence intervals (CI) for population mean µ Section 8-1: CI for µ when σ 2 known & drawing from normal distribution Section 8-1.2: Sample
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 informationAsian Option Pricing: Monte Carlo Control Variate. A discrete arithmetic Asian call option has the payoff. S T i N N + 1
Asian Option Pricing: Monte Carlo Control Variate A discrete arithmetic Asian call option has the payoff ( 1 N N + 1 i=0 S T i N K ) + A discrete geometric Asian call option has the payoff [ N i=0 S T
More informationPortfolio Management and Optimal Execution via Convex Optimization
Portfolio Management and Optimal Execution via Convex Optimization Enzo Busseti Stanford University April 9th, 2018 Problems portfolio management choose trades with optimization minimize risk, maximize
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 informationMath Models of OR: More on Equipment Replacement
Math Models of OR: More on Equipment Replacement John E. Mitchell Department of Mathematical Sciences RPI, Troy, NY 12180 USA December 2018 Mitchell More on Equipment Replacement 1 / 9 Equipment replacement
More information