On the Ross recovery under the single-factor spot rate model
|
|
- Owen Mitchell
- 6 years ago
- Views:
Transcription
1 .... On the Ross recovery under the single-factor spot rate model M. Kijima Tokyo Metropolitan University 11/08/2016 Kijima (TMU) Ross Recovery August 11, / 35
2 Plan of My Talk..1 Introduction: Motivation, Literature Review..2 The Ross Recovery; Review..3 The Setup One-Dimensional Diffusion Finite Discrete Approximation; Skip-Free Random Walk..4 Numerical Method for the Ross Recovery..5 Numerical Example: Vasicek Model..6 Extensions Quadratic Gaussian Model One-Dimensional Diffusion with One-Side Jumps Two-Dimension Models..7 Conclusion and Future Research Kijima (TMU) Ross Recovery August 11, / 35
3 Motivation Risk management of interest-rate sensitive products has become more important than ever due to the drastic change of the monetary policies, in particular after the credit crunch. For example, many countries have adopted the so-called zero interest-rate policy (ZIRP) and now started introducing negative interest rates as in Japan. As a result, it becomes harder and harder to predict the future dynamics of the interest rates because of the untraditional policies. Scenarios based on historical data often become useless to predict such drastic changes (e.g., Swiss Franc shock) Practitioners want to have forward-looking scenarios for the risk management rather than the traditional backward-looking ones. Kijima (TMU) Ross Recovery August 11, / 35
4 Motivation, Continued For the risk evaluation purpose, it is important to distinguish the risk-neutral measure Q from the physical (observed) measure P. While Q is needed for the pricing of financial assets, P is used to generate future interest-rate scenarios. In practice, the common way to construct two models under Q and P simultaneously is first to assume a financial model under either P or Q and then identify the other model under Q or P, respectively, by estimating the market price of risk from market data. E.g., Kijima et al. (2014) assume that the spot-rate process under P follows a QG model whose parameters are estimated by GMM, and assuming the market price of risk is a piecewise linear function of time, they calibrate the parameters of the model under Q. Kijima (TMU) Ross Recovery August 11, / 35
5 Motivation, Continued On the other hand, as in Dai and Singleton (2000), it is possible to start from a (affine) model under Q and, by assuming that the market price of risk follows a (affine) model, a (affine) model under P can be estimated from the term structure observed in the market. Note that, in either models, estimation of the parameters typically requires heavy computation, because the optimization problem involved is usually non-linear and difficult to converge because of the existence of many local-optima. Also, in these approaches, the estimation relies on the historical data and so, the resulting risk evaluation is a backward-looking. Practitioners require a simple and fast method to evaluate, e.g., P/L distributions of derivatives that are consistent to the current derivative prices based on forward-looking risk scenarios. For this purpose, they typically assume P = Q or a zero drift model. Kijima (TMU) Ross Recovery August 11, / 35
6 The Ross Recovery Ross (2015) proposes the idea to recover the model under P from the model under Q; i.e., it can avoid the estimation of the market price of risk, which is presumably the most difficult and contentious. Because this method uses the market data of derivatives only and derivative prices are determined by future dynamics of the underlying assets, the resulting risk evaluation is considered to be a forward-looking, in contrast to the existing methods. Qin and Linetsky (2016) state that such an identification would be of great interest to finance researchers and market participants, as it would open avenues for extracting market s assessment of physical probabilities that could be incorporated in investment decisions and supply scenarios for risk management. Kijima (TMU) Ross Recovery August 11, / 35
7 Literature Review Following Ross (2015), there appear a vast literature that discuss the Ross recovery theorem. The first set of papers such as Audrino et al. (2015) and Kiriu and Hibiki (2015) use the market data of option prices to estimate the transition law under Q by using the MP method and apply the theorem to translate them to that under P. However, the resulting MP problem is rather huge and ill-posed, which makes the problem very difficult to solve. The second set of papers such as Audrino et al. (2015) and Boroviǒka, et al. (2014) provide empirical analyses of the theorem. They test whether or not the recovery yields predictive information beyond what can be gleaned from risk-neutral densities. The results often appear very opposite; some are affirmative and some are very negative. Kijima (TMU) Ross Recovery August 11, / 35
8 Literature Review, Continued The third set of papers such as Park (2015) and Qin and Linetsky (2016) intend to extend the recovery theorem to a continuous-time setting. The continuous-time model fails to recover a physical measure from a risk-neutral measure in general. If X t is recurrent under P, then recovery is possible. Park (2015) investigates what information is sufficient to recover when X t is transient. Qin and Linetsky (2016) extend the theorem to recurrent Borel right processes (time-homogeneous). For the one-dimensional diffusion including the CIR and Vasicek models (with positive exponential jumps), they provide a complete answer to the problem. Kijima (TMU) Ross Recovery August 11, / 35
9 In This Paper We study the Ross recovery when the spot-rate process follows a one-dimensional diffusion (possibly inhomogeneous). In particular, we study the (inhomogeneous) Hull White model, which is the most common term-structure model for practitioners. We construct a discrete-time skip-free random walk under Q from the market data of option prices. The transition laws of the random walk under P can be calculated by the recovery theorem by applying a novel bisection search. A numerical example is given to support the usefulness of our method for the risk evaluation purpose. Our method can be extended to the case that the spot-rate process follows a one-dimensional diffusion with one-side jumps. Kijima (TMU) Ross Recovery August 11, / 35
10 The Ross Recovery; Setup Consider a discrete-time economy with time index {0, 1, 2,..., T }. Let r n be the riskfree spot rate at time n and denote B n = e n 1 u=0 r u t, where t > 0 is the time interval. For any asset price S t, we have [ ] S n = E P ηn+1 n S n+1 = E Q n η n [ ] Bn S n+1 B n+1 where η n denotes the state price density at time n. Let X n be a finite Markov chain (possibly inhomogeneous) with transition probability matrix Q n = (q ij (n)) under Q. Its transition matrix under P is denoted by P n = (p ij (n)). Uncertainty is driven by X n only and we shall denote S n = s n (X n ), r n = r n (X n ), η n = η n (X n ) Note that these functions are time-dependent. Kijima (TMU) Ross Recovery August 11, / 35
11 The Ross Recovery; Review Suppose X n = i, so that s n (i) = j η n+1 (j) η n (i) p ij(n)s n+1 (j) = 1 R i (n) q ij (n)s n+1 (j) j where R i (n) = e r n(i) t, which implies η n+1 (j) η n (i) p ij(n) = 1 R i (n) q ij(n) = q ij (n), i, j Ross (2015) basically assumes there are δ(n) and ξ i (n) such that η n+1 (j) η n (i) = δ(n) ξ j(n) ξ i (n) δ(n)ξ j(n) ξ i (n) p ij(n) = q ij (n), i, j In other words, there is δ(n) such that η n+1 (j) = δ(n)η n (j) for all j; i.e., δ(n) is independent of j, called transition independent. Kijima (TMU) Ross Recovery August 11, / 35
12 The Ross Recovery Theorem Let D = diag(r i ) and E = diag(ξ i ) be the diagonal matrices with diagonals R i and ξ i, respectively (suppress (n) for simplicity). The above equations can be formally written in matrix form as δe 1 P E = D 1 Q P = δ 1 E(D 1 Q)E 1 If D 1 Q is primitive, then there exist λ and z such that (D 1 Q)z = λz, λ > 0, z > 0 due to the Perron-Frobenius (PF) theorem. Further suppose that P is stochastic; we then have P e = δ 1 E(D 1 Q)y = e, y = E 1 e, where e denotes the column vector with all entities being unity. By the uniqueness, we obtain δ = λ and y = z. It follows that ξ i = 1/z i and we thus have p ij = z j λr i z i q ij Kijima (TMU) Ross Recovery August 11, / 35
13 Some Remarks..1 Suppose that R i = R, i.e., the spot rate does not depend on the underlying Markov chain X n. Then, since the transition matrix Q is stochastic (if so is P ), we have z = e, λ = 1/R, which implies that p ij = q ij...2 It may be common to assume that η t = η(x t ) for some time-independent function η(x). But, this case yields δ = 1 so that we have p ij = q ij again...3 If we assume S n = s(x n ) for some time-independent function s(x), we have arbitrage opportunities. Hence, the asset price S n should be dependent on time such as derivatives with finite maturities. Kijima (TMU) Ross Recovery August 11, / 35
14 The Setup Consider the one-dimensional (possibly inhomogeneous) recurrent diffusion X t under Q with appropriate boundary conditions: dx t = µ(t, X t )dt + σ(t, X t )dz t Such a diffusion can be approximated by a finite random walk. Let t > 0 and define x = ξ t for some ξ > 0. We discretize the model by t n = n t and x i = i x. The transition probabilities q ij (n) are determined by µ(t n, x i ) = 1 t [q i,i+1(n) x q i,i 1 (n) x] and σ 2 (t n, x i ) = 1 t [q i,i+1(n)( x) 2 + q i,i 1 (n)( x) 2 ] µ 2 (t n, x i ) t Kijima (TMU) Ross Recovery August 11, / 35
15 The Ross Recovery Consider the tri-diagonal matrix Q n = (q ij (n)) under Q. The transition probability at time T from i to j of the random walk under Q can be obtained by (Q 0 Q 1 Q T 1 ) ij. Let R i (n) = e r n(x i ) t and D n = diag(r i (n)). From the recovery theorem, we then have, for each n, p ij (n) = z j (n) λ(n)r i (n)z i (n) q ij(n), i, j where λ(n) and z(n) = (z i (n)) are the PF eigenvalue and the associated eigenvector of D 1 n Q n. The transition probability under P is obtained by (P 0 P 1 P T 1 ) ij. Kijima (TMU) Ross Recovery August 11, / 35
16 The PF Solution In order to solve the PF equation, consider in general λ z 1 z 2. z N 1 z N m 1 u d 2 m 2 u 2 0 = d N 1 m N 1 u N d N m N z 1 z 2. z N 1 z N Namely, we solve the following system of linear equations: λz 1 = m 1 z 1 + u 1 z 2 λz 2 = d 2 z 1 + m 2 z 2 + u 2 z 3. λz N 1 = d N 1z N 2 + m N 1z N 1 + u N 1z N λz N = d N z N 1 + m N z N Kijima (TMU) Ross Recovery August 11, / 35
17 The PF Solution, Continued Defining χ i = z i+1 /z i, the above equations become λ = m 1 + u 1χ 1 λ = d 2 χ m 2 + u 2 χ 2. λ = d N 1χ 1 N 2 + mn 1 + un 1χN 1 λ = d N χ 1 N 1 + mn Choose λ > 0 arbitrary, and solve the equations from above as χ 1 (λ) = [λ m 1 ]/u 1 χ 2 (λ) = [λ d 2 χ 1 1 (λ) m 2 ]/u 2. χ N 1 (λ) = [λ d N 1 χ 1 N 2 (λ) m N 1]/u N 1 χ N (λ) = λ d N χ 1 N 1 (λ) + mn Kijima (TMU) Ross Recovery August 11, / 35
18 Algorithm to Find the PF Solution If χ i (λ) > 0 and χ N (λ) = 0, then this λ is the PF eigenvalue. By induction, it is shown that the functions χ i (λ) are continuous and strictly increasing in λ for all i. Moreover, if χ i (λ) > 0, we have the following relationship: χ N (x) < 0 x < λ, χ N (x) = 0 x = λ, χ N (x) > 0 x > λ Therefore, by using the standard search (e.g., bisection method), we can easily find the PF eigenvalue λ that satisfies χ N (λ) = 0, and at the same time χ i = χ i (λ) numerically. Kijima (TMU) Ross Recovery August 11, / 35
19 The Hull White Model Suppose that the spot rate r t follows dr t = (ϕ(t) ar t )dt + σ(t)dz t where ϕ(t) and σ(t) are deterministic functions of time t, a is a positive constant, and z t denotes a standard BM under Q. ϕ(t) is chosen to be consistent with the current term structure. It is well known that r t = θ(t) + x t, dx t = ax t dt + σ(t)dz t The shift function θ(t) is given by ( t ) θ(t) = e at e as ϕ(s)ds + r 0 Note that r 0 = θ(0) implies x 0 = 0. 0 Kijima (TMU) Ross Recovery August 11, / 35
20 Discretization For the time interval t, according to Hull and White (1994), we discretize the state as x = σ(t n ) 3 t, i.e., ξ = σ(t n ) 3 Construct the nodes (n, i) when time is n and state is i so that the realization of x t is given by x(n, i) = x i i x. The state space is S = { I, I + 1,..., I 1, I} The transition probabilities under Q are defined as (i ±I q i,i+1 (n) = a2 i 2 ( t) 2 ai t 2 q i,i (n) = 2 3 a2 i 2 ( t) 2 q i,i 1 (n) = a2 i 2 ( t) 2 + ai t 2 In order for these probabilities to be positive, we must have 6 6 3a t < I i I < 3a t q ij (n) = q ij on S are determined as far as a is estimated. Kijima (TMU) Ross Recovery August 11, / 35
21 Ross Recovery Suppose that the risk-neutral transition probabilities Q = (q ij ) are estimated by some means (equivalently, a is given). It follows that R i (n) = e rn(i) t = e (θn+xi) t = e θn t R x i (n), Rx i (n) = exi t Note that the constant term in spot rates does not affect the estimate of the transition probabilities p ij (n) under P. Denoting Dn x = diag(rx i (n)), we have the PF equation ((D x n ) 1 Q)z n = λ n z n, λ n > 0, z n = (z i (n)) > 0 Defining χ i (n) = z i+1 (n)/z i (n) for the solution z n, we have p i,i+1 (n) = χ i(n) λ n R x i (n)q i,i+1, p i,i 1 (n) = 1 λ n R x i (n)χ i 1(n) q i,i 1 Kijima (TMU) Ross Recovery August 11, / 35
22 Remarks..1 The transition probabilities q ij depend on a, but not on the volatility σ(t). That is, once the speed of mean reversion a is estimated by some means, the parameters t and a alone determine the transition probabilities under Q...2 The volatility σ(t) affects the discretization of the state space...3 Since Ri x(n) depends on x which depends on σ(t n), the solution of the PF equation also depends on the volatility σ(t)...4 Hence, the transition probabilities p ij (n) depend on the volatility σ(t), and the recovery is inhomogeneous in time...5 Also, note that Ri x < 1 for i < 0. Kijima (TMU) Ross Recovery August 11, / 35
23 Comparison with Qin and Linetsky (2016) When θ(t) = θ and σ(t) = σ, the Hull White is reduced to the Vasicek: dr t = a(θ r t )dt + σdz Q t They show that, under P, the spot rate process follows ) dr t = a (θ σ2 a r 2 t dt + σdzt P for a choice of eigenvalue λ = θ σ 2 /(2a 2 ) In the Hull White model, their results are not directly applicable, because it is not time-homogeneous. Note. Bloomberg recommends to fix a = 3% and calibrate σ(t) as a piecewise constant. The volatility function is usually calibrated by swaption prices in the market. Kijima (TMU) Ross Recovery August 11, / 35
24 Figure: Future Distribution in Vasicek Date: Sep. 30, 2015 (Source: Reuters) Estimated Result: a = 3%, σ = 1.02% (σ/a) 2 = 11.56% Y P Q Kijima (TMU) Ross Recovery August 11, / 35
25 Extension 1: Shifted QG Model Suppose that the spot rate r t follows the SDE under Q r t = (x t + α(t)) 2 + θ(t), dx t = ax t dt + σdz t Here, θ(t) is the shift function (used to fit the initial term structure) and, when θ(t) < 0, the spot rates can be negative. When θ(t) = 0, the model is the QG model of Pelsser (1997). Kijima et al. (2014) derive the closed form solution for discount bond as well as cap/floor prices when α(t) is a piecewise linear function. Suppose that Q n = (q ij (n)) is estimated by some means. Then, R i (n) = e rn(x i) t = e θn t R x i (n), Rx i (n) = e(x i+α n) 2 t Hence, the Ross Recovery can be applied by the same way as in the Hull White model. Kijima (TMU) Ross Recovery August 11, / 35
26 Extension 2: With One-Side Jumps Consider the one-dimensional recurrent diffusion with jumps under Q with appropriate boundary conditions: dx t = µ(t, X t )dt + σ(t, X t )dz t + dj t where J t is a compound Poisson with negative (or positive) jumps. J t and z t are mutually independent. Such a jump-diffusion process can be approximated by a skip-free random walk (positive or negative, respectively) The associated transition probabilities can be calculated easily because of the independence. Kijima (TMU) Ross Recovery August 11, / 35
27 Extension 2: With One-Side Jumps, Continued The skip-free positive Markov chain has the transition matrix of the form m 1 u d 2 m 2 u 2 0 Q = d N 1,1 d N 1 m N 1 u N 1 d N,1 d N,N 2 d N,N 1 m N Kijima (1993) develops a bisection algorithm to solve the PF equation, when the Markov chain is skip-free to the one direction. Using the solution of the PF equation, we can apply the recovery theorem to obtain the transition matrix under P. Kijima (TMU) Ross Recovery August 11, / 35
28 Extension 3: Two-Dim. Model Consider a two-dimensional recurrent Markov chain (X t, Y t ) with transition matrix P = (p (i,j),(k,l) ) and Q = (q (i,j),(k,l) ) under P and Q, respectively (assume time-homogeneity for simplicity). Assume that S t = s t (X t, Y t ), r t = r(x t ), η t = η t (X t, Y t ) so, X t drives the interest-rate dynamics. For X t = i, Y t = j, we have R i = e r(i) t and s t (i, j) = k,l η T (k, l) η t (i, j) p (i,j),(k,l)s T (k, l) = It follows that, under the Ross assumption, 1 q (i,j),(k,l) s T (k, l) R i δ η k,l η i,j p (i,j),(k,l) = 1 R i q (i,j),(k,l), i, j, k, l k,l Kijima (TMU) Ross Recovery August 11, / 35
29 The Setup for the Two-Dim. Model We label the states (i, j) in a lexicographical order. Suppose we have N states for X t and M states for Y t, and define the subspace S i = {(i, 1), (i, 2),..., (i, M 1), (i, M)} The state space is S = N i=1 S i, and the total number is N M. The transition matrix under Q has the following form: Q 11 Q 12 Q 1N Q 21 Q 22 Q 2N Q =......, Q ik = ( ) M q (i,j),(k,l) j,l=1 Q N1 Q N2 Q NN Q ik is the transition matrix from S i to S k. Kijima (TMU) Ross Recovery August 11, / 35
30 The Recovery Theorem: Two-Dim. Model Define the diagonal matrix D = diag(r i e) and consider the PF equation ( D 1 Q)z = λz, λ > 0, z > 0 The recovery is obtained by p (i,j),(k,l) = z k,l λr i z i,j q (i,j),(k,l) Here, λ is the PF eigenvalue and z = (z i ) is the associated eigenvector with z i = (z i,j ) being the sub-vector on the sub-space S i. Kijima (TMU) Ross Recovery August 11, / 35
31 Equity Option under Stochastic Interest Rates Consider the model ds t = r t dt + σ S dw t, S t dr t = a(θ r t )dt + σ r dz t, where dw t dz t = ρdt. The option pricing formula is available for the model, so assume that the parameters are all estimated under Q. Define so that dy t = y t = 1 σ S log S t ρ σ r r t (( 1 + aρ ) r t aρ θ σ ) S dt + dw t ρdz t σ S σ r σ r 2 Kijima (TMU) Ross Recovery August 11, / 35
32 Equity Option under SIR, Continued Also, define ΣdW t = dw t ρdz t, Σ = 1 ρ 2 so that dw t dz t = 0 We then have mutually independent processes dx t = ax t dt + σ r dz t (( 1 dy t = + aρ σ S σ r ) (θ + x t ) aρ θ σ S σ r 2 ) dt + ΣdW t Using these processes, we have r t = θ + x t { S t = exp σ S y t + ρσ } S (θ + x t ) σ r Kijima (TMU) Ross Recovery August 11, / 35
33 Two-Dim. Lattice: Equity Option under SIR We construct a two-dim. tri-nominal model for (x t, y t ) as follows. Take t > 0 and define x = σ r 3 t, y = Σ 3 t Transition probabilities of X n are given as before, because x t is the OU (r t is the Vasicek). Transition probabilities of Y n are obtained by µ y (y j ) = 1 ] [q yj,j+1 t y qyj,j 1 y [ ] σ 2 y (y j) = 3Σ 2 q y j,j+1 + qy j,j 1 µ 2 y (y j) t Kijima (TMU) Ross Recovery August 11, / 35
34 Two-Dim. Lattice, Continued Solving these equations, we get q y j,j+1 = m2 i ( t)2 + m i t 2 q y j,j 1 = m2 i ( t)2 m i t 2 Because X n and Y n are mutually independent, we finally have q (i,j),(k,l) = q x ik qy jl The resulting transition matrix under Q becomes a block tri-diagonal of the form M I U I O O D I+1 M I+1 U I+1 O D 1 Q = O D I 1 M I 1 U I 1 O O D I M I Kijima (TMU) Ross Recovery August 11, / 35
35 Conclusion and Future Research In this paper, we develop a numerical method for the Ross recovery when the spot rate follows a one-dimensional diffusion (possibly inhomogeneous and with one-side jumps). We approximate such a process by a discrete-time, finite random walk with skip-free nature. The PF equation can be solved easily by applying a novel bisection search. Our method can be applied to equity options under stochastic interest-rate economy. There remain many problems; in particular, lots of numerical examples to test the usefulness of the Ross recovery in risk management (and also investment decisions) Kijima (TMU) Ross Recovery August 11, / 35
36 Thank You for Your Attention Kijima (TMU) Ross Recovery August 11, / 35
The 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 informationSteven Heston: Recovering the Variance Premium. Discussion by Jaroslav Borovička November 2017
Steven Heston: Recovering the Variance Premium Discussion by Jaroslav Borovička November 2017 WHAT IS THE RECOVERY PROBLEM? Using observed cross-section(s) of prices (of Arrow Debreu securities), infer
More informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationMulti-dimensional Term Structure Models
Multi-dimensional Term Structure Models We will focus on the affine class. But first some motivation. A generic one-dimensional model for zero-coupon yields, y(t; τ), looks like this dy(t; τ) =... dt +
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationValuation of derivative assets Lecture 8
Valuation of derivative assets Lecture 8 Magnus Wiktorsson September 27, 2018 Magnus Wiktorsson L8 September 27, 2018 1 / 14 The risk neutral valuation formula Let X be contingent claim with maturity T.
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More informationθ(t ) = T f(0, T ) + σ2 T
1 Derivatives Pricing and Financial Modelling Andrew Cairns: room M3.08 E-mail: A.Cairns@ma.hw.ac.uk Tutorial 10 1. (Ho-Lee) Let X(T ) = T 0 W t dt. (a) What is the distribution of X(T )? (b) Find E[exp(
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 informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
More informationMultiname and Multiscale Default Modeling
Multiname and Multiscale Default Modeling Jean-Pierre Fouque University of California Santa Barbara Joint work with R. Sircar (Princeton) and K. Sølna (UC Irvine) Special Semester on Stochastics with Emphasis
More informationLecture 17. The model is parametrized by the time period, δt, and three fixed constant parameters, v, σ and the riskless rate r.
Lecture 7 Overture to continuous models Before rigorously deriving the acclaimed Black-Scholes pricing formula for the value of a European option, we developed a substantial body of material, in continuous
More informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More informationLecture 5: Review of interest rate models
Lecture 5: Review of interest rate models Xiaoguang Wang STAT 598W January 30th, 2014 (STAT 598W) Lecture 5 1 / 46 Outline 1 Bonds and Interest Rates 2 Short Rate Models 3 Forward Rate Models 4 LIBOR and
More informationImplementing an Agent-Based General Equilibrium Model
Implementing an Agent-Based General Equilibrium Model 1 2 3 Pure Exchange General Equilibrium We shall take N dividend processes δ n (t) as exogenous with a distribution which is known to all agents There
More informationModelling Default Correlations in a Two-Firm Model by Dynamic Leverage Ratios Following Jump Diffusion Processes
Modelling Default Correlations in a Two-Firm Model by Dynamic Leverage Ratios Following Jump Diffusion Processes Presented by: Ming Xi (Nicole) Huang Co-author: Carl Chiarella University of Technology,
More informationModelling Credit Spread Behaviour. FIRST Credit, Insurance and Risk. Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent
Modelling Credit Spread Behaviour Insurance and Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent ICBI Counterparty & Default Forum 29 September 1999, Paris Overview Part I Need for Credit Models Part II
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 informationCrashcourse Interest Rate Models
Crashcourse Interest Rate Models Stefan Gerhold August 30, 2006 Interest Rate Models Model the evolution of the yield curve Can be used for forecasting the future yield curve or for pricing interest rate
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
More informationEquilibrium Term Structure Models. c 2008 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 854
Equilibrium Term Structure Models c 2008 Prof. Yuh-Dauh Lyuu, National Taiwan University Page 854 8. What s your problem? Any moron can understand bond pricing models. Top Ten Lies Finance Professors Tell
More informationHedging Credit Derivatives in Intensity Based Models
Hedging Credit Derivatives in Intensity Based Models PETER CARR Head of Quantitative Financial Research, Bloomberg LP, New York Director of the Masters Program in Math Finance, Courant Institute, NYU Stanford
More informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
More 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 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 informationThe Term Structure of Interest Rates under Regime Shifts and Jumps
The Term Structure of Interest Rates under Regime Shifts and Jumps Shu Wu and Yong Zeng September 2005 Abstract This paper develops a tractable dynamic term structure models under jump-diffusion and regime
More informationInterest rate models in continuous time
slides for the course Interest rate theory, University of Ljubljana, 2012-13/I, part IV József Gáll University of Debrecen Nov. 2012 Jan. 2013, Ljubljana Continuous time markets General assumptions, notations
More informationOrder book resilience, price manipulations, and the positive portfolio problem
Order book resilience, price manipulations, and the positive portfolio problem Alexander Schied Mannheim University PRisMa Workshop Vienna, September 28, 2009 Joint work with Aurélien Alfonsi and Alla
More informationInterest Rate Volatility
Interest Rate Volatility III. Working with SABR Andrew Lesniewski Baruch College and Posnania Inc First Baruch Volatility Workshop New York June 16-18, 2015 Outline Arbitrage free SABR 1 Arbitrage free
More informationTerm Structure Lattice Models
IEOR E4706: Foundations of Financial Engineering c 2016 by Martin Haugh Term Structure Lattice Models These lecture notes introduce fixed income derivative securities and the modeling philosophy used to
More informationApplication of Stochastic Calculus to Price a Quanto Spread
Application of Stochastic Calculus to Price a Quanto Spread Christopher Ting http://www.mysmu.edu/faculty/christophert/ Algorithmic Quantitative Finance July 15, 2017 Christopher Ting July 15, 2017 1/33
More informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
More 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 informationRECURSIVE VALUATION AND SENTIMENTS
1 / 32 RECURSIVE VALUATION AND SENTIMENTS Lars Peter Hansen Bendheim Lectures, Princeton University 2 / 32 RECURSIVE VALUATION AND SENTIMENTS ABSTRACT Expectations and uncertainty about growth rates that
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 informationRisk Neutral Measures
CHPTER 4 Risk Neutral Measures Our aim in this section is to show how risk neutral measures can be used to price derivative securities. The key advantage is that under a risk neutral measure the discounted
More informationAD in Monte Carlo for finance
AD in Monte Carlo for finance Mike Giles giles@comlab.ox.ac.uk Oxford University Computing Laboratory AD & Monte Carlo p. 1/30 Overview overview of computational finance stochastic o.d.e. s Monte Carlo
More 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 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 informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationFunctional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs.
Functional vs Banach space stochastic calculus & strong-viscosity solutions to semilinear parabolic path-dependent PDEs Andrea Cosso LPMA, Université Paris Diderot joint work with Francesco Russo ENSTA,
More informationOne-Factor Models { 1 Key features of one-factor (equilibrium) models: { All bond prices are a function of a single state variable, the short rate. {
Fixed Income Analysis Term-Structure Models in Continuous Time Multi-factor equilibrium models (general theory) The Brennan and Schwartz model Exponential-ane models Jesper Lund April 14, 1998 1 Outline
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 informationRoss Recovery theorem and its extension
Ross Recovery theorem and its extension Ho Man Tsui Kellogg College University of Oxford A thesis submitted in partial fulfillment of the MSc in Mathematical Finance April 22, 2013 Acknowledgements I am
More informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
More informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More informationM.I.T Fall Practice Problems
M.I.T. 15.450-Fall 2010 Sloan School of Management Professor Leonid Kogan Practice Problems 1. Consider a 3-period model with t = 0, 1, 2, 3. There are a stock and a risk-free asset. The initial stock
More informationFinite dimensional realizations of HJM models
Finite dimensional realizations of HJM models Tomas Björk Stockholm School of Economics Camilla Landén KTH, Stockholm Lars Svensson KTH, Stockholm UTS, December 2008, 1 Definitions: p t (x) : Price, at
More informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
More informationSensitivity Analysis on Long-term Cash flows
Sensitivity Analysis on Long-term Cash flows Hyungbin Park Worcester Polytechnic Institute 19 March 2016 Eastern Conference on Mathematical Finance Worcester Polytechnic Institute, Worceseter, MA 1 / 49
More informationRisk, Return, and Ross Recovery
Risk, Return, and Ross Recovery Peter Carr and Jiming Yu Courant Institute, New York University September 13, 2012 Carr/Yu (NYU Courant) Risk, Return, and Ross Recovery September 13, 2012 1 / 30 P, Q,
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 3, Mikhail Zaslavsky 3 University of Michigan, Ann
More informationImplementing the HJM model by Monte Carlo Simulation
Implementing the HJM model by Monte Carlo Simulation A CQF Project - 2010 June Cohort Bob Flagg Email: bob@calcworks.net January 14, 2011 Abstract We discuss an implementation of the Heath-Jarrow-Morton
More informationThe Binomial Model. The analytical framework can be nicely illustrated with the binomial model.
The Binomial Model The analytical framework can be nicely illustrated with the binomial model. Suppose the bond price P can move with probability q to P u and probability 1 q to P d, where u > d: P 1 q
More informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
More 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 informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationStochastic volatility modeling in energy markets
Stochastic volatility modeling in energy markets Fred Espen Benth Centre of Mathematics for Applications (CMA) University of Oslo, Norway Joint work with Linda Vos, CMA Energy Finance Seminar, Essen 18
More informationA Recovery That We Can Trust? Deducing and Testing the Restrictions of the Recovery Theorem
A Recovery That We Can Trust? Deducing and Testing the Restrictions ohe Recovery Theorem Gurdip Bakshi Fousseni Chabi-Yo Xiaohui Gao University of Houston December 4, 2015 Bakshi & Chabi-Yo & Gao Test
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 informationTopic 2 Implied binomial trees and calibration of interest rate trees. 2.1 Implied binomial trees of fitting market data of option prices
MAFS5250 Computational Methods for Pricing Structured Products Topic 2 Implied binomial trees and calibration of interest rate trees 2.1 Implied binomial trees of fitting market data of option prices Arrow-Debreu
More informationAn Introduction to Point Processes. from a. Martingale Point of View
An Introduction to Point Processes from a Martingale Point of View Tomas Björk KTH, 211 Preliminary, incomplete, and probably with lots of typos 2 Contents I The Mathematics of Counting Processes 5 1 Counting
More informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
More informationA Hybrid Commodity and Interest Rate Market Model
A Hybrid Commodity and Interest Rate Market Model University of Technology, Sydney June 1 Literature A Hybrid Market Model Recall: The basic LIBOR Market Model The cross currency LIBOR Market Model LIBOR
More informationA model reduction approach to numerical inversion for parabolic partial differential equations
A model reduction approach to numerical inversion for parabolic partial differential equations Liliana Borcea Alexander V. Mamonov 2, Vladimir Druskin 2, Mikhail Zaslavsky 2 University of Michigan, Ann
More informationGeneralized Recovery
Generalized Recovery Christian Skov Jensen Copenhagen Business School David Lando Copenhagen Business School and CEPR Lasse Heje Pedersen AQR Capital Management, Copenhagen Business School, NYU, CEPR December,
More informationSmoking Adjoints: fast evaluation of Greeks in Monte Carlo calculations
Report no. 05/15 Smoking Adjoints: fast evaluation of Greeks in Monte Carlo calculations Michael Giles Oxford University Computing Laboratory, Parks Road, Oxford, U.K. Paul Glasserman Columbia Business
More informationCredit-Equity Modeling under a Latent Lévy Firm Process
.... Credit-Equity Modeling under a Latent Lévy Firm Process Masaaki Kijima a Chi Chung Siu b a Graduate School of Social Sciences, Tokyo Metropolitan University b University of Technology, Sydney September
More informationLinear-Rational Term-Structure Models
Linear-Rational Term-Structure Models Anders Trolle (joint with Damir Filipović and Martin Larsson) Ecole Polytechnique Fédérale de Lausanne Swiss Finance Institute AMaMeF and Swissquote Conference, September
More informationCalibration of Ornstein-Uhlenbeck Mean Reverting Process
Calibration of Ornstein-Uhlenbeck Mean Reverting Process Description The model is used for calibrating an Ornstein-Uhlenbeck (OU) process with mean reverting drift. The process can be considered to be
More informationLikelihood Estimation of Jump-Diffusions
Likelihood Estimation of Jump-Diffusions Extensions from Diffusions to Jump-Diffusions, Implementation with Automatic Differentiation, and Applications Berent Ånund Strømnes Lunde DEPARTMENT OF MATHEMATICS
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 information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
More informationState processes and their role in design and implementation of financial models
State processes and their role in design and implementation of financial models Dmitry Kramkov Carnegie Mellon University, Pittsburgh, USA Implementing Derivative Valuation Models, FORC, Warwick, February
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 informationMODELING INVESTMENT RETURNS WITH A MULTIVARIATE ORNSTEIN-UHLENBECK PROCESS
MODELING INVESTMENT RETURNS WITH A MULTIVARIATE ORNSTEIN-UHLENBECK PROCESS by Zhong Wan B.Econ., Nankai University, 27 a Project submitted in partial fulfillment of the requirements for the degree of Master
More informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
More informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationStochastic Local Volatility: Excursions in Finite Differences
Stochastic Local Volatility: Excursions in Finite Differences ICBI Global Derivatives Paris April 0 Jesper Andreasen Danske Markets, Copenhagen kwant.daddy@danskebank.dk Outline Motivation: Part A & B.
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 informationLog-Robust Portfolio Management
Log-Robust Portfolio Management Dr. Aurélie Thiele Lehigh University Joint work with Elcin Cetinkaya and Ban Kawas Research partially supported by the National Science Foundation Grant CMMI-0757983 Dr.
More informationParameter estimation in SDE:s
Lund University Faculty of Engineering Statistics in Finance Centre for Mathematical Sciences, Mathematical Statistics HT 2011 Parameter estimation in SDE:s This computer exercise concerns some estimation
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationMATH 3200 Exam 3 Dr. Syring
. Suppose n eligible voters are polled (randomly sampled) from a population of size N. The poll asks voters whether they support or do not support increasing local taxes to fund public parks. Let M be
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 information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More informationConstructing Markov models for barrier options
Constructing Markov models for barrier options Gerard Brunick joint work with Steven Shreve Department of Mathematics University of Texas at Austin Nov. 14 th, 2009 3 rd Western Conference on Mathematical
More informationHigh Frequency Trading in a Regime-switching Model. Yoontae Jeon
High Frequency Trading in a Regime-switching Model by Yoontae Jeon A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Mathematics University
More informationSTOCHASTIC INTEGRALS
Stat 391/FinMath 346 Lecture 8 STOCHASTIC INTEGRALS X t = CONTINUOUS PROCESS θ t = PORTFOLIO: #X t HELD AT t { St : STOCK PRICE M t : MG W t : BROWNIAN MOTION DISCRETE TIME: = t < t 1
More informationApproximating a multifactor di usion on a tree.
Approximating a multifactor di usion on a tree. September 2004 Abstract A new method of approximating a multifactor Brownian di usion on a tree is presented. The method is based on local coupling of the
More informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
More informationdt+ ρσ 2 1 ρ2 σ 2 κ i and that A is a rather lengthy expression that we may or may not need. (Brigo & Mercurio Lemma Thm , p. 135.
A 2D Gaussian model (akin to Brigo & Mercurio Section 4.2) Suppose where ( κ1 0 dx(t) = 0 κ 2 r(t) = δ 0 +X 1 (t)+x 2 (t) )( X1 (t) X 2 (t) ) ( σ1 0 dt+ ρσ 2 1 ρ2 σ 2 )( dw Q 1 (t) dw Q 2 (t) ) In this
More information