AN EULER-TYPE METHOD FOR THE STRONG APPROXIMATION OF THE COX-INGERSOLL-ROSS PROCESS
|
|
- Roberta Hensley
- 5 years ago
- Views:
Transcription
1 AN ULR-TYP MTHOD FOR TH STRONG APPROXIMATION OF TH COX-INGRSOLL-ROSS PROCSS STFFN DRICH, ANDRAS NUNKIRCH AND LUKASZ SZPRUCH Abstract. We analyze the strong approximation of the Cox-Ingersoll-Ross CIR process in the regime where the process does not hit zero by a positivity preserving drift-implicit uler-type method. As an error criterion we use the p-th mean of the maximum distance between the CIR process and its approximation on a finite time interval. We show that under mild assumptions on the parameters of the CIR process the proposed method attains, up to a logarithmic term, the convergence of order /2. This agrees with the standard rate of the strong convergence for global approximations of stochastic differential equations SDs with Lipschitz coefficients despite the fact that the CIR process has a non-lipschitz diffusion coefficient. Keywords. Cox-Ingersoll-Ross process; strong global approximation; square-root coefficient; Lamperti transformation; drift-implicit method; one-sided Lipschitz condition Mathematics Subject Classification 2. 65C3; 6H35. Introduction and Main Result In computational finance a lot of effort has been given to the so called Cox-Ingersoll- Ross process CIR recently. The CIR process under consideration has the following form dx t = κλ X t dt + θ X t dw t, X = x, t. Here W = W t t, is a one-dimensional Brownian motion, κ, λ, θ > and x >. It is well known that equation admits a unique strong solution which is non-negative, see e.g. Chapter 5 in [24]. The CIR process was originally proposed by Cox et al. [9] in 985 as a model for short-term interest rates. Nowadays, this model is widely used in financial modeling, e.g. as volatility process in the Heston model [8]. One of the main objectives in mathematical finance is the pricing of path-dependent derivatives. If the asset prices or the interest rates dynamics are modeled by a d- dimensional SD with solution S t t [,T ], then this corresponds to the quadrature problem p = F S where F : C[, T ], R d R is the discounted payoff of the derivative. Typically, explicit formulae for such quantities are unknown and have to be approximated by Monte Carlo methods that are based on approximate solutions of the stochastic differential equation on [, T ]. Here, the knowledge of the global strong convergence rate of the approximation is important, in particular if the efficient Multi-level Monte-Carlo method [2, 3] is to be used. Date: October 27, 2.
2 2 STFFN DRICH, ANDRAS NUNKIRCH AND LUKASZ SZPRUCH The strong global approximation of equation has been studied in several articles. Strong convergence without a rate or with a logarithmic rate of several discretization schemes has been shown in [2, 7,, 7, 9]. In [2], a general framework for the analysis of strong approximation of the CIR process is presented along with extensive simulation studies. Moreover, the strong approximation of more general Ait-Sahalia-type interest rate models is analyzed in [2]. However, only in [4] non-logarithmic convergence rates are obtained. In [4] it is shown that a symmetrized uler method has strong convergence order /2 under restrictive assumptions on the parameters of the equation, see Section 2. The difficulties to obtain strong convergence rates for the approximations of the CIR process are due to its square-root coefficient. Thus the standard theory which relies on the global Lipschitz assumption does not apply [25, 26]. In this article we focus on the regime where the CIR process does not hit zero, i.e. where PX t > for t = an assumption which is often fulfilled in interest rate models. By the Feller test, it is true if and only if 2κλ, see e.g. Chapter 5 in [24]. We will not directly do a numerical analysis for the CIR process, but for a coordinate transformation thereof. We consider the process Y t = X t, which by Itô s formula satisfies dy t = α Y t dt + βy t dt + γ dw t, t, Y = x, 2 with 4κλ θ2 α =, β = κ 8 2, γ =. This transformation, known as Lamperti transformation, allows us to shift the nonlinearity from the diffusion coefficient into the drift coefficient. Note that the drift fx = α + βx, x >, x satisfies for α >, β R the one-sided Lipschitz condition x yfx fy βx y 2, x, y >. This property is crucial to control the error propagation of drift-implicit uler schemes, see [2, 2]. The drift-implicit uler method with stepsize > for equation 2 leads to the numerical scheme α y k+ = y k + + βy k+ + γ k W, k =,,... 3 y k+ with y = x and k W = W k+ W k, k =,,.... Recalling that α, γ > and β <, equation 3 has the unique positive solution y k+ = y k + γ k W 2 β + y k + γ k W 2 + α 4 β 2 β, which we call drift-implicit square-root uler method. Transforming back, i.e. x k = y 2 k, k =,,...,
3 AN ULR-TYP MTHOD FOR TH CIR PROCSS 3 gives a strictly positive approximation of the original CIR process. This scheme has already been suggested in [2], but convergence results have not been established. Using piecewise linear interpolation, i.e. x t = k + t x k + t k x k+, t [k, k + ], we obtain a global approximation x t t [,T ] of the CIR process on [, T ]. xploiting the structure of SD 2 we establish the following theorem. Theorem.. Let 2κλ >, x > and T >. Then, for all p < 2κλ there exists a constant K p > such that /p max X t x t p K p log, t [,T ] for all, /2]. Hence we obtain the optimal strong convergence rate for the approximation of SDs with Lipschitz coefficients, see [27]. The only price to be paid is the restriction on p that arises from the need to control the inverse p-th moments of the CIR process, which are infinite for p 2κλ. The remainder of this article is structured as follows: In the next section we give a short overview on discretization schemes for the CIR process, while the proof of Theorem. is given in Section Numerical Methods for the CIR Process Discretization schemes for the CIR process were proposed in numerous articles, among these are [2, 3, 4, 7,, 5, 7, 9, 23, 28]. In this short summary we mainly focus on uler-type methods. In [4] the authors study a symmetrized uler method for the approximation of equation, i.e. x k+ = x k + κλ x k + θ x k k W 4 for k =,,.... Under the assumption 2κλ > + { } κ 8 max 6p, 6p 2 θ they found that max X k x k 2p C p p, where the constant C p > depends only p, κ, λ, θ, x and T. As a consequence of Lemma 3.5 in Section 3 the piecewise linear interpolation of this scheme satisfies the same error estimate with respect to the p-th mean maximum distance as our drift-implicit squareroot uler scheme. However for the symmetrized uler scheme the assumptions on the parameters of the CIR process are clearly more restrictive. Moreover, it is shown in [7] that the symmetrized uler scheme converges weakly with rate if 2κλ 2.
4 4 STFFN DRICH, ANDRAS NUNKIRCH AND LUKASZ SZPRUCH The truncated uler scheme x k+ = x k + κλ x k + θ x + k kw, k =,,... 5 was analyzed in [], while the scheme x k+ = x k + κλ x k + θ x k k W, k =,,... 6 was studied in [9]. Both schemes do not preserve positivity, but satisfy max X k x k 2 for without further restrictions on the parameters of the equation. In particular an application of Theorem 3. in [7] yields a logarithmic convergence rate for the uler schemes 5 and 6. Moreover, if 2κλ it follows from [6] that the uler schemes 4 6 have a pathwise convergence rate of /2 ε for all ε >, i.e. we have max X /2 ε k x k P a.s. for. The asymptotic error distribution of these schemes can be deduced from [29]: For 2κλ it holds that max X L k x k max t t [,T ] Φ t db s 8 Φ s for, where B t t is a Brownian motion independent of W t t and Φ t t is given by t Φ t = exp κt θ2 ds + θ t dw s, t. 8 X s 2 Xs While no explicit solution of equation is known, the finite dimensional distributions can be characterized in terms of a non-central chi-square distribution, see e.g. [9]. Thus, equation can be simulated exactly at a finite number of time points using the Markov property, see e.g. [8, 4]. However, the algorithms for the exact simulation of the CIR process are strongly problem dependent: The number of degrees of freedom of the non-central chi-square random variable, which has to be simulated in each step, is 4κλ/. Thus the computational cost of the algorithms depends strongly on κ, λ and θ. The same problem, i.e. strong dependence of the computational cost of the algorithm on the parameters of the equation, arises also for the exact sampling algorithm introduced in [5], which can be also applied to equation, see [6]. While for the simulation of the CIR process at a single point the exact simulation methods are useful, discretization schemes remain erior if a full sample path of the CIR process has to be simulated or if the CIR process is part of a system of stochastic differential equations, see e.g. [9]. Moreover, results on strong convergence rates for approximations of equation have an interest of its own, since it is one of the most prominent examples for a stochastic differential equation, whose coefficients do not satisfy the standard global Lipschitz assumption.
5 AN ULR-TYP MTHOD FOR TH CIR PROCSS 5 3. Proof of Theorem. In the following we will denote by c constants regardless of their value. 3.. Preliminaries. For our error analysis we need to control the inverse moments of the CIR process. Since X t follows a non-central chi-square distribution, we have 2κ X p t = x exp κt p x expκt 2κλ Γ + p θ 2 F Γ 2κλ p p, 2κλ, 2κ x expκt for p > 2κλ, where F denotes the confluent hypergeometric function, and else, see e.g. Theorem 3. in [22]. Since F a, b, z = X p t = Γb Γb a z a + O z, z, see formula 3..5 on page 54 in [], it follows for p > 2κλ that t X p t is bounded on [, T ], i.e. X p t < for p > 2κλ t [,T ] θ. 7 2 Moreover, from Theorem 3. [22] or Lemma A.2 in [7] we also have T exp ε Xs ds < 2κλ if and only if ε θ2 2 8 θ. In general, one can estimate polynomial moments 2 against exponential moments, since for q, ε > there exists c > such that x q ce εx for x. Hence we arrive at Lemma 3.. Let 2κλ >, T > and q. It holds T q Xs ds <. Further we need a smoothness result for equation 2: Lemma 3.2. Let 2κλ > and T >. Then, for all q we have and s<t T, t s Yt Y s q c t s q/2, for s, t [, T ], Y t Y s q c log q/2, for, /2], Y t q <. t [,T ]
6 6 STFFN DRICH, ANDRAS NUNKIRCH AND LUKASZ SZPRUCH Proof. For s t, we have t α Y t Y s = du + s Y u By the Cauchy-Schwarz inequality one has Y t Y s γ W t W s + t s /2 T t s α 2 Yu 2 βy u du + γ W t W s. /2 T /2 du + t s /2 β 2 Yu du 2. Since Y u = X u the first assertion follows now from 7, Lemma 3., Minkowski s inequality and the smoothness of Brownian motion in the q-th mean. For the second assertion in addition we use the well known fact that the modulus of continuity of Brownian motion satisfies s<t T, t s W t W s q c log q/2 for, /2], see e.g. []. The third assertion follows from T Y t y + γ W t + t /2 α 2 /2 T du + t /2 and Y 2 u W t q <. t [,T ] 3.2. rror Bound for the implicit uler Scheme for Y. /2 β 2 Yu 2 du Proposition 3.3. Let 2κλ >. For T > and p < 2κλ, there exists c > such that /p Y k y k p c, for, /2]. Proof. Without loss of generality we assume that < T. For the error at time point k we have the recursion with e = e k = Y k y k e k+ = e k + fy k+ fy k+ + r k r k = k+ k Multiplying both sides with e k+ we obtain fy k+ fy t dt. e 2 k+ 2 e2 k + 2 e2 k+ + e k+ fyk+ fy k+ + e k+ r k. Since e k+ fyk+ fy k+ βe 2 k+
7 AN ULR-TYP MTHOD FOR TH CIR PROCSS 7 we have and it follows that n e 2 n 2 e k+ r k, n =, 2,... k= T/ e k 2 r k. 8 k= Now, it remains to analyze the local error k+ r k = fy k+ fy t dt. k Note that for a, b >. Thus fa fb = βa b + α b a ab fa fb c + a b ab and we obtain k+ r k c + k Y t Y k+ Y k+ Y t dt. An application of Hölder s and Minkowski s inequality yields rk p k+ /p pq Yk+ c + Y Y t Y k+ pq t pq for q, q > with q + q for all q >. Since k =. Now, Lemma 3.2 gives r k p /p c k+ k + Y t Y k+ pq Y t Y k+ pq Y t + 2pq Y k+, 2pq pq pq dt dt, 9 applying 7 with q > such that pq < 2κλ yields r k p /p c 3/2, for p < 2κλ and using 8 completes the proof of the Proposition.
8 8 STFFN DRICH, ANDRAS NUNKIRCH AND LUKASZ SZPRUCH 3.3. Moment Bounds for the implicit uler Scheme for Y. Now we show that all moments of the approximation scheme are uniformly bounded. Multiplying 3 with y k+ yields It follows with and we obtain by induction that y 2 k+ = α + βy 2 k+ + y k+ y k + γ k W. for all k =,,..., T/. This allows us to show: y 2 k+ 2α + γ 2 + y 2 k + M k M k = 2γy k k W + γ 2 k W 2, y 2 k+ y 2 + 2α + γ 2 k + 2 Lemma 3.4. Let 2κλ and >, T >. Then for all p we have Proof. From and we obtain that with y k p <. y k 2 c + c k m t dw t 3 m t = 2γy l + 2γ 2 W t W l, t [l, l +. Since t [, T/ ] m t p c + c the Burkholder-Davis-Gundy inequality, i.e. s p m τ dw τ c s [,t] Jensen s inequality and 3 give that l=,..., T/ t y l p m 2 τ dτ p/2 y k 2p c + c y l p. l=,..., T/, So 2 now yields y k 4 < and the assertion follows from an induction procedure in p.
9 AN ULR-TYP MTHOD FOR TH CIR PROCSS rror Bound for the drift-implicit square-root uler Scheme. Now denote by X the piecewise linear interpolation of the CIR process with stepsize >, i.e. X t = k + t t X k + k X k+, t [k, k + ]. Lemma 3.5. Let 2κλ >, T > and, /2]. Then, for all q we have Proof. Combining the equality max t [,T ] X t X t q c log q/2. X t X s = Y t + Y s Y t Y s with the Cauchy-Schwarz inequality and Lemma 3.2 we get s<t T, t s Now the assertion follows from X t X s q c log q/2. X t X t X t X s. t [,T ] s<t T, t s Due to X t x t X t X t + X k yk 2 t [,T ] t [,T ] and the above Lemma it only remains to control /p. X k yk 2 p Since X k yk 2 p Y k + y k p Y k y k, p an application of Hölder s inequality with ε > such that + εp < 2κλ and Proposition 3.3 give X k yk 2 p Y k + y k p +ε ε +ε ε p/2. It remains to apply Lemma 3.2 and Lemma 3.4 to finish the proof of Theorem.. Acknowledgements. The authors would like to thank Martin Altmayer for valuable comments on an earlier version of the manuscript.
10 STFFN DRICH, ANDRAS NUNKIRCH AND LUKASZ SZPRUCH References [] Abramowitz, M., Stegun, I.A.: Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables. U.S. Department of Commerce, Washington 964 [2] Alfonsi, A.: On the discretization schemes for the CIR and Bessel squared processes. Monte Carlo Methods Appl. 4, [3] Alfonsi, A.: High order discretization schemes for the CIR process: Application to affine term structure and Heston models. Math. Comput , [4] Berkaoui, A., Bossy, M., Diop, A.: uler scheme for SDs with non-lipschitz diffusion coefficient: strong convergence. SAIM, Probab. Stat. 2, - 28 [5] Beskos, A., Roberts, G.O.: xact simulation of diffusions. Ann. Appl. Probab. 54, [6] Beskos, A., Papaspiliopoulos, O., Roberts, G.O., Fearnhead, P.: xact and computationally efficient likelihood-based estimation for discretely observed diffusion processes With discussion. J. R. Stat. Soc., Ser. B, Stat. Methodol. 683, [7] Bossy, M., Diop, A.: An efficient discretization scheme for one dimensional SDs with a diffusion coefficient function of the form x a, a [/2,. Working paper, INRIA 27 [8] Broadie, M., Kaya, Ö.: xact simulation of stochastic volatility and other affine jump diffusion processes. Oper. Res. 54, [9] Cox, J., Ingersoll, J., Ross, S.: A theory of the term structure of interest rates. conometrica 53, [] Deelstra, G., Delbaen, F.: Convergence of discretized stochastic interest rate processes with stochastic drift term. Appl. Stochastic Models Data Anal. 4, [] Fischer, M., Nappo, G.: On the moments of the modulus of continuity of Itô processes. Stochastic Anal. Appl. 28, [2] Giles, M.: Multi-level Monte Carlo path simulation, Oper. Res. 563, [3] Giles, M.: Improved multilevel Monte Carlo convergence using the Milstein scheme. In: Keller, A. ed. et al., Monte Carlo and Quasi-Monte Carlo Methods 26, Proceedings. Berlin, Springer, [4] Glassermann, P.: Monte Carlo Methods in Financial ngineering. Springer, New York 24 [5] Günther, M., Kahl, C., Roßberg, T.: Structure preserving stochastic integration schemes in interest rate derivative modeling. Appl. Numer. Math. 583, [6] Gyöngy, I.: A note on uler s approximations. Potential Anal. 8, [7] Gyöngy, I., Rásonyi, M.: A note on uler approximations for SDs with Hölder continuous diffusion coefficients. Stoch. Proc. Appl. 2, [8] Heston, S.L.: A closed form solution for options with stochastic volatility, with applications to bonds and currency options. Rev. Financial Stud. 6, [9] Higham, D.J., Mao, X.: Convergence of Monte Carlo simulations involving the mean-reverting square root process. J. Comp. Fin. 8, [2] Higham, D.J., Mao, X., Stuart, A.M.: Strong convergence of uler-type methods for nonlinear stochastic differential equations. SIAM J. Numer. Anal. 43, [2] Higham, D.J., Mao, X., Szpruch, L., Pan. J.: Numerical simulation of a strongly nonlinear Ait- Sahalia-type interest rate model. BIT 52, [22] Hurd, T.R., Kuznetsov, A.: xplicit formulas for Laplace transforms of stochastic integrals. Markov Process. Relat. Fields 42, [23] Jäckel, P., Kahl, C.: Fast strong approximation Monte-Carlo schemes for stochastic volatility models. Quant. Finance 6, [24] Karatzas, I., Shreve, S..: Brownian Motion and Stochastic Calculus. Springer, New York, 2nd edition 99 [25] Kloeden, P.., Platen,.: Numerical Solution of Stochastic Differential quations. Springer, Berlin, 3rd edition 999 [26] Milstein, G.N.: Numerical Integration of Stochastic Differential quations. Kluwer, Doordrecht 995
11 AN ULR-TYP MTHOD FOR TH CIR PROCSS [27] Müller-Gronbach, T.: The optimal uniform approximation of systems of stochastic differential equations. Ann. Appl. Probab. 22, [28] Moro,., Schurz, H.: Boundary preserving semianalytic numerical algorithms for stochastic differential equations. SIAM J. Sci. Comput. 29, [29] Neuenkirch, A., Zähle, H.: Asymptotic error distribution of the uler method for SDs with non-lipschitz coefficients. Monte Carlo Methods Appl. 54, Steffen Dereich, Fachbereich Mathematik und Informatik, Philipps-Universität Marburg, Hans-Meerwein-Straße, D-3532 Marburg, Germany Andreas Neuenkirch, Fachbereich Mathematik, TU Kaiserslautern, Postfach 349, D Kaiserslautern, Germany Lukasz Szpruch, Mathematical Institute, University of Oxford, St Giles, Oxford OX 3LB, ngland
An 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 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 informationAnalysing multi-level Monte Carlo for options with non-globally Lipschitz payoff
Finance Stoch 2009 13: 403 413 DOI 10.1007/s00780-009-0092-1 Analysing multi-level Monte Carlo for options with non-globally Lipschitz payoff Michael B. Giles Desmond J. Higham Xuerong Mao Received: 1
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 informationAffine term structures for interest rate models
Stefan Tappe Albert Ludwig University of Freiburg, Germany UNSW-Macquarie WORKSHOP Risk: modelling, optimization and inference Sydney, December 7th, 2017 Introduction Affine processes in finance: R = a
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 informationEquivalence between Semimartingales and Itô Processes
International Journal of Mathematical Analysis Vol. 9, 215, no. 16, 787-791 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.215.411358 Equivalence between Semimartingales and Itô Processes
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 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 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 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 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 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 informationConvergence of Discretized Stochastic (Interest Rate) Processes with Stochastic Drift Term.
Convergence of Discretized Stochastic (Interest Rate) Processes with Stochastic Drift Term. G. Deelstra F. Delbaen Free University of Brussels, Department of Mathematics, Pleinlaan 2, B-15 Brussels, Belgium
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 informationA THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES
Proceedings of ALGORITMY 01 pp. 95 104 A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES BEÁTA STEHLÍKOVÁ AND ZUZANA ZÍKOVÁ Abstract. A convergence model of interest rates explains the evolution of the
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
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 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 informationarxiv: v1 [math.pr] 15 Dec 2011
Parameter Estimation of Fiber Lay down in Nonwoven Production An Occupation Time Approach Wolfgang Bock, Thomas Götz, Uditha Prabhath Liyanage arxiv:2.355v [math.pr] 5 Dec 2 Dept. of Mathematics, University
More informationSensitivity of American Option Prices with Different Strikes, Maturities and Volatilities
Applied Mathematical Sciences, Vol. 6, 2012, no. 112, 5597-5602 Sensitivity of American Option Prices with Different Strikes, Maturities and Volatilities Nasir Rehman Department of Mathematics and Statistics
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 informationIntroduction to Affine Processes. Applications to Mathematical Finance
and Its Applications to Mathematical Finance Department of Mathematical Science, KAIST Workshop for Young Mathematicians in Korea, 2010 Outline Motivation 1 Motivation 2 Preliminary : Stochastic Calculus
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 informationStochastic Calculus of Heston s Stochastic-Volatility Model
Stochastic Calculus of Heston s Stochastic-Volatility Model MTNS 21 Floyd B. Hanson Departments of Mathematics University of Illinois and University of Chicago hanson@math.uic.edu March 5, 21 Abstract
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 informationPDE Methods for the Maximum Drawdown
PDE Methods for the Maximum Drawdown Libor Pospisil, Jan Vecer Columbia University, Department of Statistics, New York, NY 127, USA April 1, 28 Abstract Maximum drawdown is a risk measure that plays an
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 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 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 informationMean-square stability properties of an adaptive time-stepping SDE solver
Mean-square stability properties of an adaptive time-stepping SDE solver H. Lamba and T. Seaman Department of Mathematical Sciences, George Mason University, MS 3F2, 4400 University Drive, Fairfax, VA
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 informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationNumerical Solution of Stochastic Differential Equations with Jumps in Finance
Numerical Solution of Stochastic Differential Equations with Jumps in Finance Eckhard Platen School of Finance and Economics and School of Mathematical Sciences University of Technology, Sydney Kloeden,
More informationMODELLING 1-MONTH EURIBOR INTEREST RATE BY USING DIFFERENTIAL EQUATIONS WITH UNCERTAINTY
Applied Mathematical and Computational Sciences Volume 7, Issue 3, 015, Pages 37-50 015 Mili Publications MODELLING 1-MONTH EURIBOR INTEREST RATE BY USING DIFFERENTIAL EQUATIONS WITH UNCERTAINTY J. C.
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 informationAmerican Foreign Exchange Options and some Continuity Estimates of the Optimal Exercise Boundary with respect to Volatility
American Foreign Exchange Options and some Continuity Estimates of the Optimal Exercise Boundary with respect to Volatility Nasir Rehman Allam Iqbal Open University Islamabad, Pakistan. Outline Mathematical
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 informationA note on the existence of unique equivalent martingale measures in a Markovian setting
Finance Stochast. 1, 251 257 1997 c Springer-Verlag 1997 A note on the existence of unique equivalent martingale measures in a Markovian setting Tina Hviid Rydberg University of Aarhus, Department of Theoretical
More informationLogarithmic derivatives of densities for jump processes
Logarithmic derivatives of densities for jump processes Atsushi AKEUCHI Osaka City University (JAPAN) June 3, 29 City University of Hong Kong Workshop on Stochastic Analysis and Finance (June 29 - July
More informationOption Pricing under Delay Geometric Brownian Motion with Regime Switching
Science Journal of Applied Mathematics and Statistics 2016; 4(6): 263-268 http://www.sciencepublishinggroup.com/j/sjams doi: 10.11648/j.sjams.20160406.13 ISSN: 2376-9491 (Print); ISSN: 2376-9513 (Online)
More informationStochastic Runge Kutta Methods with the Constant Elasticity of Variance (CEV) Diffusion Model for Pricing Option
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 18, 849-856 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4381 Stochastic Runge Kutta Methods with the Constant Elasticity of Variance
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 informationPortfolio optimization problem with default risk
Portfolio optimization problem with default risk M.Mazidi, A. Delavarkhalafi, A.Mokhtari mazidi.3635@gmail.com delavarkh@yazduni.ac.ir ahmokhtari20@gmail.com Faculty of Mathematics, Yazd University, P.O.
More informationRegression estimation in continuous time with a view towards pricing Bermudan options
with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten
More informationCalibration of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationDrawdowns Preceding Rallies in the Brownian Motion Model
Drawdowns receding Rallies in the Brownian Motion Model Olympia Hadjiliadis rinceton University Department of Electrical Engineering. Jan Večeř Columbia University Department of Statistics. This version:
More informationApplication of MCMC Algorithm in Interest Rate Modeling
Application of MCMC Algorithm in Interest Rate Modeling Xiaoxia Feng and Dejun Xie Abstract Interest rate modeling is a challenging but important problem in financial econometrics. This work is concerned
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 informationThe Azema Yor embedding in non-singular diusions
Stochastic Processes and their Applications 96 2001 305 312 www.elsevier.com/locate/spa The Azema Yor embedding in non-singular diusions J.L. Pedersen a;, G. Peskir b a Department of Mathematics, ETH-Zentrum,
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 informationWeierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions
Weierstrass Institute for Applied Analysis and Stochastics Maximum likelihood estimation for jump diffusions Hilmar Mai Mohrenstrasse 39 1117 Berlin Germany Tel. +49 3 2372 www.wias-berlin.de Haindorf
More informationPRICING TIMER OPTIONS UNDER FAST MEAN-REVERTING STOCHASTIC VOLATILITY
CANADIAN APPLIED MATHEMATICS QUARTERLY Volume 17, Number 4, Winter 009 PRICING TIMER OPTIONS UNDER FAST MEAN-REVERTING STOCHASTIC VOLATILITY DAVID SAUNDERS ABSTRACT. Timer options are derivative securities
More informationOPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF FINITE
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 WeA11.6 OPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF
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 informationGeneralized Affine Transform Formulae and Exact Simulation of the WMSV Model
On of Affine Processes on S + d Generalized Affine and Exact Simulation of the WMSV Model Department of Mathematical Science, KAIST, Republic of Korea 2012 SIAM Financial Math and Engineering joint work
More informationCONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES
CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES D. S. SILVESTROV, H. JÖNSSON, AND F. STENBERG Abstract. A general price process represented by a two-component
More informationEstimation of Value at Risk and ruin probability for diffusion processes with jumps
Estimation of Value at Risk and ruin probability for diffusion processes with jumps Begoña Fernández Universidad Nacional Autónoma de México joint work with Laurent Denis and Ana Meda PASI, May 21 Begoña
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 informationSOME APPLICATIONS OF OCCUPATION TIMES OF BROWNIAN MOTION WITH DRIFT IN MATHEMATICAL FINANCE
c Applied Mathematics & Decision Sciences, 31, 63 73 1999 Reprints Available directly from the Editor. Printed in New Zealand. SOME APPLICAIONS OF OCCUPAION IMES OF BROWNIAN MOION WIH DRIF IN MAHEMAICAL
More informationEnlargement of filtration
Enlargement of filtration Bernardo D Auria email: bernardo.dauria@uc3m.es web: www.est.uc3m.es/bdauria July 6, 2017 ICMAT / UC3M Enlargement of Filtration Enlargement of Filtration ([1] 5.9) If G is a
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationPolynomial processes in stochastic portofolio theory
Polynomial processes in stochastic portofolio theory Christa Cuchiero University of Vienna 9 th Bachelier World Congress July 15, 2016 Christa Cuchiero (University of Vienna) Polynomial processes in SPT
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 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 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 informationRichardson Extrapolation Techniques for the Pricing of American-style Options
Richardson Extrapolation Techniques for the Pricing of American-style Options June 1, 2005 Abstract Richardson Extrapolation Techniques for the Pricing of American-style Options In this paper we re-examine
More informationConditional Density Method in the Computation of the Delta with Application to Power Market
Conditional Density Method in the Computation of the Delta with Application to Power Market Asma Khedher Centre of Mathematics for Applications Department of Mathematics University of Oslo A joint work
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 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 informationarxiv: v1 [q-fin.pr] 23 Feb 2014
Time-dependent Heston model. G. S. Vasilev, Department of Physics, Sofia University, James Bourchier 5 blvd, 64 Sofia, Bulgaria CloudRisk Ltd (Dated: February 5, 04) This work presents an exact solution
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationarxiv: v1 [q-fin.cp] 1 Nov 2016
Essentially high-order compact schemes with application to stochastic volatility models on non-uniform grids arxiv:1611.00316v1 [q-fin.cp] 1 Nov 016 Bertram Düring Christof Heuer November, 016 Abstract
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 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 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 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 informationPricing Dynamic Solvency Insurance and Investment Fund Protection
Pricing Dynamic Solvency Insurance and Investment Fund Protection Hans U. Gerber and Gérard Pafumi Switzerland Abstract In the first part of the paper the surplus of a company is modelled by a Wiener process.
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 informationModule 10:Application of stochastic processes in areas like finance Lecture 36:Black-Scholes Model. Stochastic Differential Equation.
Stochastic Differential Equation Consider. Moreover partition the interval into and define, where. Now by Rieman Integral we know that, where. Moreover. Using the fundamentals mentioned above we can easily
More informationStochastic Differential Equations in Finance and Monte Carlo Simulations
Stochastic Differential Equations in Finance and Department of Statistics and Modelling Science University of Strathclyde Glasgow, G1 1XH China 2009 Outline Stochastic Modelling in Asset Prices 1 Stochastic
More informationA NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS. Burhaneddin İZGİ
A NEW NOTION OF TRANSITIVE RELATIVE RETURN RATE AND ITS APPLICATIONS USING STOCHASTIC DIFFERENTIAL EQUATIONS Burhaneddin İZGİ Department of Mathematics, Istanbul Technical University, Istanbul, Turkey
More informationSADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD
The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint
More informationForward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions
Forward Monte-Carlo Scheme for PDEs: Multi-Type Marked Branching Diffusions Pierre Henry-Labordère 1 1 Global markets Quantitative Research, SOCIÉTÉ GÉNÉRALE Outline 1 Introduction 2 Semi-linear PDEs 3
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 informationJournal of Mathematical Analysis and Applications
J Math Anal Appl 389 (01 968 978 Contents lists available at SciVerse Scienceirect Journal of Mathematical Analysis and Applications wwwelseviercom/locate/jmaa Cross a barrier to reach barrier options
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 informationarxiv: v1 [q-fin.pm] 13 Mar 2014
MERTON PORTFOLIO PROBLEM WITH ONE INDIVISIBLE ASSET JAKUB TRYBU LA arxiv:143.3223v1 [q-fin.pm] 13 Mar 214 Abstract. In this paper we consider a modification of the classical Merton portfolio optimization
More informationRough volatility models: When population processes become a new tool for trading and risk management
Rough volatility models: When population processes become a new tool for trading and risk management Omar El Euch and Mathieu Rosenbaum École Polytechnique 4 October 2017 Omar El Euch and Mathieu Rosenbaum
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 informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More informationDiffusions, Markov Processes, and Martingales
Diffusions, Markov Processes, and Martingales Volume 2: ITO 2nd Edition CALCULUS L. C. G. ROGERS School of Mathematical Sciences, University of Bath and DAVID WILLIAMS Department of Mathematics, University
More informationOn worst-case investment with applications in finance and insurance mathematics
On worst-case investment with applications in finance and insurance mathematics Ralf Korn and Olaf Menkens Fachbereich Mathematik, Universität Kaiserslautern, 67653 Kaiserslautern Summary. We review recent
More informationStochastic Modelling in Finance
in Finance Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH April 2010 Outline and Probability 1 and Probability 2 Linear modelling Nonlinear modelling 3 The Black Scholes
More informationConvexity Theory for the Term Structure Equation
Convexity Theory for the Term Structure Equation Erik Ekström Joint work with Johan Tysk Department of Mathematics, Uppsala University October 15, 2007, Paris Convexity Theory for the Black-Scholes Equation
More informationLaw of the Minimal Price
Law of the Minimal Price Eckhard Platen School of Finance and Economics and Department of Mathematical Sciences University of Technology, Sydney Lit: Platen, E. & Heath, D.: A Benchmark Approach to Quantitative
More informationHedging of Contingent Claims under Incomplete Information
Projektbereich B Discussion Paper No. B 166 Hedging of Contingent Claims under Incomplete Information by Hans Föllmer ) Martin Schweizer ) October 199 ) Financial support by Deutsche Forschungsgemeinschaft,
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 information