Perturbed Gaussian Copula
|
|
- Nathan Evans
- 5 years ago
- Views:
Transcription
1 Perturbed Gaussian Copula Jean-Pierre Fouque Xianwen Zhou August 8, 6 Abstract Gaussian copula is by far the most popular copula used in the financial industry in default dependency modeling. However, it has a major drawback it does not exhibit tail dependence, a very important property for copula. he essence of tail dependence is the interdependence when extreme events occur, say, defaults of corporate bonds. In this paper we show that some tail dependence can be restored by introducing stochastic volatility on a Gaussian copula. Using perturbation methods we then derive an approximate copula called perturbed Gaussian copula in this paper. A copula is a joint distribution function of uniform random variables. Sklar s heorem states that for any multivariate distribution, the univariate marginal distributions and the dependence structure can be separated. he dependence structure is completely determined by the copula. It then implies that one can borrow the dependence structure, namely the copula, of one set of dependent random variables and exchange the marginal distributions for a totally different set of marginal distributions. An important property of copula is its invariance under monotonic transformation. More precisely, if g i is strictly increasing for each i, then (g (X ),g (X ),...,g n (X n )) have the same copula as (X,X,...,X n ). From the above discussion, it is not hard to see that copula comes in default dependency modeling very naturally. For a much detailed coverage on copula, including the precise format of Sklar s heorem, as well as modeling default dependency by way of copula, the readers are referred to Schonbucher (3) 5]. Let (Z,...,Z n ) be a normal random vector with standard normal marginals and correlation matrix R, and Φ( ) be the standard normal cumulative distribution function. hen the joint distribution function of (Φ(Z ),...,Φ(Z n )) is called the Gaussian copula with correlation matrix R. Work supported by NSF grant DMS Department of Statistics and Applied Probability, University of California, Santa Barbara, CA 936-3, fouque@pstat.ucsb.edu. Mathematics Department, North Carolina State University, Raleigh, NC 7695, xzhou@ncsu.edu.
2 Gaussian copula is by far the most popular copula used in the financial industry in default dependency modeling. his is basically because of two reasons. Firstly it is easy to simulate. Secondly it requires the right number of parameters equal to the number of correlation coefficients among the underlying names. However, Gaussian copula does not exhibit any tail dependence, a very important property for copula. he essence of tail dependence is the interdependence when extreme events occur, say, defaults of corporate bonds. In fact, this is considered as a major drawback of Gaussian copula. On the other hand, by introducing stochastic volatility into the classic Black-Scholes model, Fouque, Papanicolaou and Sircar () ], by way of singular perturbation method, gave a satisfactory answer to the smile curve problem of implied volatilities in the financial market, leading to a pricing formula which is in the form of a robust simple correction to the classic Black-Scholes constant volatility formula. Furthermore, an application of this perturbation method to defaultable bond pricing has been studied by Fouque, Sircar and Solna (5) 3]. By fitting real market data, they concluded that the method works fairly well. An extension to multi-name first passage models is proposed by Fouque, Wignall and Zhou (6) 4]. In this paper we will show the effect of stochastic volatility on a Gaussian copula. Specifically, in Section, we first set up the stochastic volatility model and state out the objective the transition density functions. hen by singular perturbation, we obtain approximate transition density functions. In order to make them true probability density functions, we introduce the transformation + tanh( ). In Section, we study this new class of approximate copula density functions, first analytically and then numerically. Section 3 concludes this paper. Asymptotics. Model Setup We start with a process (X () which follows the dynamics: t,x () t,y t ) defined on the complete probability space (Ω, F, P) and dx () t = f (Y t )dw () t, dx () t = f (Y t )dw () t, dy t = ǫ (m Y t)dt + ν ǫ dw (Y ) t, where W () t,w () t and W (Y ) t are standard Brownian motions correlated as follows d W (),W () t = ρdt, d W (),W (Y ) t = ρ Y dt, d W (),W (Y ) t = ρ Y dt, with ρ,ρ Y,ρ Y and making the correlation matrix ρ ρ Y ρ ρ Y ρ Y ρ Y
3 symmetric positive definite, ǫ and ν are positive constant numbers with ǫ being small. he f i s are real functions for i =,, and are assumed here to be bounded above and below away from. It is worth noting that f i s are not explicit functions of t. hey depend on t only through Y t. Observe that Y t is a mean-reverting process and /ǫ is the rate of mean-reversion so that Y t is fast mean-reverting. Furthermore, Y t admits the unique invariant normal distribution N(m,ν ). For a fixed time >, our objective is to find, for t <, the joint distribution { } P X () ξ,x () Xt ξ = x,y t = y and the two marginal distributions { } { } P X () Xt ξ = x,y t = y, P X () Xt ξ = x,y t = y, where X t (X () t,x () t ), x (x,x ), and ξ,ξ are two arbitrary numbers. Equivalently, we need to find the following three transition densities: { } u ǫ P X () dξ,x () Xt dξ = x,y t = y, { } v ǫ P X () Xt dξ = x,y t = y, { } v ǫ P X () Xt dξ = x,y t = y, where we show the dependence on the small parameter ǫ.. PDE Representation Let us consider u ǫ first. In terms of partial differential equation (PDE), u ǫ satisfies the following Kolmogorov backward equation L ǫ u ǫ (t,x,x,y) =, u ǫ (,x,x,y) = δ(ξ ;x )δ(ξ ;x ), where δ(ξ i ;x i ) is the Dirac delta function of ξ i with spike at ξ i = x i for i =,, and operator L ǫ has the following decomposition: with the notations: L L = (m y) y L ǫ = ǫ L + ǫ L + L, + ν y, () = ν ρ Y f (y) x y + ν ρ Y f (y) x y, () L = t + f (y) x + f (y) x 3 + ρf (y)f (y). (3) x x
4 As in Fouque, Papanicolaou and Sircar () ], we expand the solution u ǫ in powers of ǫ : u ǫ = u + ǫu + ǫu + ǫ 3/ u 3 +. In the following, we will determine the first few terms appearing on the right hand side of the above expansion. Specifically, we will retain ū u + ǫu, (4) as an approximation to u ǫ (later we will propose another approximation in order to restore positiveness.).3 Leading Order erm u Following Fouque, Papanicolaou and Sircar () ], the leading order term u, which is independent of variable y, is characterized by: L u (t,x,x ) =, (5) u (,x,x ) = δ(ξ ;x )δ(ξ ;x ), where denotes the average with respect to the invariant distribution N(m,ν ) of Y t, i.e., g for a general function g of y. } g(y) { ν π exp (y m) ν dy We define the effective volatilities σ and σ, and the effective correlation ρ by: σ f, σ f, ρ ρ f f. (6) σ σ Using the definition (3) and the notations (6), equation (5) becomes u t + u σ x + u u σ x + ρ σ σ =, x x u (,x,x ) = δ(ξ ;x )δ(ξ ;x ). It can be verified that u is the transition density of two correlated scaled Brownian motions with instantaneous correlation ρ and scale factors σ and σ, respectively. hat is, u (t,x,x ) = π σ σ ( t) (7) ρ { (ξ x ) exp ( ρ ) σ ( t) ρ(ξ x )(ξ x ) σ σ ( t) + (ξ x ) ]} σ. ( t) 4
5 .4 Correction erm ǫ u Again, similar to Fouque, Papanicolaou and Sircar () ], the correction term u, which is also independent of variable y, is characterizd by: where the operator A is defined by L u (t,x,x ) = Au, (8) u (,x,x ) =, A = L L (L L ), and the inverse L is taken on the centered quantity L L. From the definition (3) of L, it is straightforward to obtain that L L = (f (y) f ) x + (f (y) f ) x + ρ(f (y)f (y) f f ). x x Let us denote by φ (y),φ (y) and φ (y) solutions of the following Poisson equations respectively L φ (y) = f (y) f, L φ (y) = f (y) f, L φ (y) = f (y)f (y) f f. heir existence (with at most polynomial growth at infinity) is guarantied by the centering property of the right hand sides and the Fredholm alternative for the infinitesimal generator L. hey are defined up to additive constants in y which will play no role after applying the operator L which takes derivatives with respect to y. It then follows that L (L L ) = φ (y) x + φ (y) x Now by the definition () of L, we have L L (L L ) = ν ρ Y f (y) φ (y) 3 +ν ρ Y f (y) herefore the operator ǫa can be written + ρφ (y). x x x φ (y) x x + ρφ 3 ] (y) x x φ (y) 3 x x + φ 3 (y) x 3 + ρφ (y) 3 x x ]. ǫa = R 3 x R x 3 + R x x + R x x where the constant parameters R,R,R and R are defined as follows: R νρ Y ǫ f φ, 5
6 R νρ Y ǫ f φ, R νρ Y ǫ f φ + ν ǫρρ Y f φ, R νρ Y ǫ f φ + ν ǫρρ Y f φ. Note that they are all small of order ǫ. It can be checked directly that u is given explicitly by u = ( t)au, and therefore ǫu = ( t) R 3 x ] + R x 3 + R x x + R x x u. (9) Explicit formulas for the third order partial derivatives of u are given in Appendix A..5 Regularity Conditions for Density Functions Since = u (t,x,x ;ξ,ξ )dξ dξ, by Lebesgue dominated convergence theorem, we then have = k +k x k = xk for integers k,k. It follows that k +k x k xk u (t,x,x ;ξ,ξ )dξ dξ, ǫu (t,x,x ;ξ,ξ )dξ dξ =, and hence ū(t,x,x ;ξ,ξ )dξ dξ =, where ū = u + ǫ u is the approximation introduced in (4). In order to guarantee that our approximated transition density function is always non-negative, which is the other regularity condition for a density function, we seek a multiplicative perturbation of the form ũ û ( + tanh( ǫû )), 6
7 where û and û are defined such that u + ǫu = û ( + ǫû ) for any ǫ >. It can be easily seen that this is achieved with the choice: û = u, û = u /u. Now instead of using ū as our approximation for u ǫ, we use ũ = u + tanh( ǫu /u )] () = u { + tanh ( ( t) u 3 u 3 u 3 u 3 ])} u R x 3 + R x 3 + R x x + R x x. Before proving that ũ given in () is indeed a probability density function, we clarify a definition first. Definition Let g be a function of n variables (x,x,...,x n ) R n. he function g is called an n-dimensional even function if g( x, x,..., x n ) = g(x,x,...,x n ) for all (x,x,...,x n ) R n, and an n-dimensional odd function if for all (x,x,...,x n ) R n. g( x, x,..., x n ) = g(x,x,...,x n ) With this definition, we can state the following proposition. Proposition Let g(x) be a probability density function on R n for n, and ϕ(x) be an odd function. If g is an even function, then the function f defined by is also a probability density function on R n. f(x) = ( + tanh(ϕ(x))) g(x) Proof We need to prove that f is globally non-negative and its integral over R n is equal to one. Observe that tanh( ) is strictly between and, and this together with the non-negativity of g justifies that f is always non-negative. On the other hand, tanh( ) is a (-dimensional) odd function, and hence tanh(ϕ(x)) is an (n-dimensional) odd function. Now by change of variables y = x, we have I tanh(ϕ(x))g(x)dx = tanh(ϕ( y))g( y)dy R n R n = tanh(ϕ(y))g(y)dy = I, R n 7
8 which implies that I =. herefore f(x)dx = g(x)dx + I = + =. R n R n he proof is complete. Now observe that u is a probability density function with respect to the variables (ξ,ξ ), and is even on (ξ x,ξ x ). In addition, ǫu /u is an odd function on (ξ x,ξ x ). By Proposition we know that ũ given in () is indeed a probability density function. As for the approximation accuracy ũ u ǫ, we first note that tanh(x) x x3 3, when x is close to. Now for fixed (t,x,x ), when ǫ is small, we have ũ = u + tanh( ǫu /u )] ǫ u u + ( ) ] 3 ǫu u 3 u = u + ( ) ( ) u ǫu ǫ 3/ 3 u 3u = ū ǫ 3/ 3 3u. herefore ũ ū is small of order ǫ 3/, while ū u ǫ is small of order ǫ (see Fouque et al. (3) ]). hus ũ u ǫ is small of the same order of ǫ as ū u ǫ, i.e., the approximation accuracy remains unchanged when replacing ū by ũ..6 Marginal ransition Densities For the marginal transition density function { } v ǫ P X () Xt dξ = x,y t = y, the above argument goes analogously, and we obtain v ǫ v p (t,x ;,ξ σ ) ( t)r 3 x 3 p (t,x ;,ξ σ ), where p (t,x ;,ξ σ ) is the transition density of the scaled Brownian motion with scale factor σ, that is, p (t,x ;,ξ σ ) = exp { (ξ x ) } π( t) σ σ. ( t) A straightforward calculation shows that 3 p 3(ξ x ) x 3 = π σ 5 ( t) + (ξ x ) 3 ]exp { (ξ x ) } 5/ π σ 7 ( t) 7/ σ. ( t) 8
9 Note again that v (t,x ;,ξ )dξ =. o guarantee the non-negativity of the approximated density function, we, again, use instead ( ṽ p 3 )] p + tanh ( t)r p x 3 as our approximation to v ǫ. where By symmetry we have v ǫ P p (t,x ;,ξ σ ) = 3 p x 3 { } X () Xt dξ = x,y t = y v p (t,x ;,ξ σ ) ( t)r 3 x 3 ( ṽ p 3 )] p + tanh ( t)r p x 3, = and ṽ is our approximation to v ǫ. exp { (ξ x ) π( t) σ σ ( t) }, p (t,x ;,ξ σ ) 3(ξ x ) π σ 5 ( t) + (ξ x ) 3 ]exp { (ξ x ) 5/ π σ 7 ( t) 7/ σ ( t) By exactly the same argument used for ũ, one can show that ṽ and ṽ are indeed probability density functions of ξ and ξ, respectively. Furthermore, the approximation accuracies remain unchanged when switching from v to ṽ, and from v to ṽ. }, Density of the Perturbed Copula. Approximated Copula Density Now suppose that conditional on {X t = x,y t = y}, (X (),X() ) admits the copula Ψ(, ), then, by Sklar s heorem, its density function ψ(, ) can be represented as where ψ(z,z ) = z z u ǫ (t,x,x,y;,ξ,ξ ) v ǫ (t,x,y;,ξ ) v ǫ (t,x,y;,ξ ), { } = P X () Xt ξ = x,y t = y, { } = P X () Xt ξ = x,y t = y. 9
10 Observe that if the volatility terms (f ( ),f ( )) for (X () t ) were constant numbers, say, the process {Y t } t was constant or the f i s were both identically constant, then Ψ would be a Gaussian copula. where Using our approximations to u ǫ,v ǫ and vǫ, we have ψ(ζ,ζ ) ψ(ζ,ζ ) t,x () ũ(t,x,x ;,ξ,ξ ) ṽ (t,x ;,ξ ) ṽ (t,x ;,ξ ), () ζ = = ζ = = ξ ξ ξ ξ ṽ (t,x ;,ξ )dξ p (t,x ;,ξ ) + tanh ( ( t)r p (t,x ;,ξ ) ṽ (t,x ;,ξ )dξ p (t,x ;,ξ ) + tanh ( ( t)r p (t,x ;,ξ ) 3 )] p (t,x ;,ξ ) x 3 dξ, 3 )] p (t,x ;,ξ ) x 3 dξ. he function ũ is given by (), and the marginals (p,p ) and their derivatives 3 p, 3 p x 3 x 3 explicitly in the previous section.6. are given Before justifying that ψ is a probability density function defined on the unit square,], we need the following proposition. Proposition Suppose function Θ(x,x,...,x n ) is an n-dimensional probability density function on R n for n, and h (x ),h (x ),...,h n (x n ) are -dimensional strictly positive probability density functions. hen the function c defined on the unit hyper-square,] n by c(z,z,...,z n ) = Θ(x,x,...,x n ) n i= h i(x i ) with z i,] given by z i = xi h i (y i )dy i is a probability density function on,] n. Furthermore, c is a copula density function if and only if h (x ),h (x ),...,h n (x n ) are the marginal density functions of Θ(x,x,...,x n ), meaning that h i (x i ) = Θ(x,x,...,x n )dx dx dx i dx i+ dx n R n for every i =,,...,n.
11 Proof Let H i be the cumulative distribution function of h i. hen H i is strictly increasing, implying the existence of its inverse function, and z i = H i (x i ), or equivalently, x i = H i (z i ) for each i. Since Θ is non-negative, and h i s are strictly positive, the function c is non-negative. On the other hand, = = c(z,z,...,z n )dz dz dz n,] n n c(h (x ),H (x ),...,H n (x n )) h i (x i ) dx dx dx n R n i= Θ(x,x,...,x n ) dx dx dx n =. R n herefore c(z,z,...,z n ) is a probability density function on,] n. Now if the additional condition is satisfied, then we have = = c(z,z,...,z n )dz dz n,] n n c(z,h (x ),...,H n (x n )) h i (x i ) dx dx n R n i= Θ(x,x,...,x n ) dx dx n =. h (x ) R n his is to say that the marginal density function for the variable z is one, and hence the marginal distribution for the variable z is uniform. Similarly, we can show that the marginal distributions for the variables z,...,z n are also uniform. By definition of copula, we know that function c is a copula density function. he converse can be obtained by reversing the above argument. he proof is complete. Now from definition () of ψ, by combining the fact that ũ, ṽ and ṽ are all probability density functions, one can see that ψ is a density function on,] by applying Proposition. However, ψ is not a copula density function in general, because the additional condition required in Proposition is not satisfied in general in our case, and hence Ψ, the copula corresponding to density function ψ, is not an exact copula in general. with Asymptotically, when ǫ goes to, for fixed (t,x,x ), the density ψ converges to φ(z,z ) z i = ξi u (t,x,x ;,ξ,ξ ) p (t,x ;,ξ ) p (t,x ;,ξ ) ( ) ξi x i p i (t,x i ;,ξ i )dξ i = N σ i t for i =,, where N( ) denotes the univariate standard normal cumulative distribution function. One should observe that φ(, ) is the two-dimensional Gaussian copula density function with correlation parameter ρ, and that it depends only on the parameter ρ, independent of any other variables/parameters, including x,x,t,, σ, σ, etc.
12 As a consequence, Ψ converges to the Gaussian copula Φ with correlation parameter ρ. Since the method used in this paper is a perturbation method, we call Ψ a perturbed Gaussian copula.. Numerical Results In this subsection, we illustrate the effectiveness of our approximation method by showing some numerical results..9.8 Marginal Densities Gauss Perturbed Gauss Simply Pert. Gauss Perturbed Gauss Figure : Perturbed Gaussian densities In Figure, we plot p (t,x ;,ξ ), v (t,x ;,ξ ) and ṽ (t,x ;,ξ ) as functions of ξ. Note that p (t,x ;,ξ ) is a standard Gaussian density without any perturbation. he upper graph demonstrates the difference between p (t,x ;,ξ ) (standard Gaussian) and ṽ (t,x ;,ξ ) (perturbed Gaussian), and the lower one between v (t,x ;,ξ ) (simply perturbed Gaussian) and ṽ (t,x ;,ξ ) (perturbed Gaussian). It can be seen from Figure that v (t,x ;,ξ ) (simply perturbed Gaussian) takes on negative values at some places; ṽ (t,x ;,ξ ) (perturbed Gaussian), however, does not take on negative values, which is guaranteed by its formation; v (t,x ;,ξ ) and ṽ (t,x ;,ξ ) are almost globally identical, which justifies the modification of the form + tanh( );
13 ṽ (t,x ;,ξ ) is considerably different from p (t,x ;,ξ ) (standard Gaussian); specifically, it shifts to the right from p (t,x ;,ξ ); Despite the difference between ṽ (t,x ;,ξ ) and p (t,x ;,ξ ), the areas under them do seem to be of the same size, which is justified by the fact that both are probability density functions and hence the overall integrals should both be one. Gaussian Copula Density Perturbed Gaussian Copula Density Figure : Gaussian copula and perturbed Gaussian copula densities In Figure, we plot ψ(, ) in the lower graph and φ(, ) in the upper graph, the Gaussian copula density that ψ converges to when ǫ goes to. It can be seen from Figure that the standard Gaussian copula density (upper graph) and the perturbed Gaussian copula density (lower graph) both present singularities at (,) and (,) but the perturbed one has more tail dependence at (, ). Our numerous numerical experiments show that this picture is extremely sensitive to the choice of parameters and gives a lot of flexibility to the shape of the perturbed Gaussian copula 3
14 density (the Matlab code is available on demand). ail dependence is a very important property for a copula, especially when this copula is to be used in modeling default correlation. he essence of tail dependence is the interdependence when extreme events occur, say, defaults of corporate bonds. he lack of tail dependence has been for years a major criticism on standard Gaussian copula. hroughout the computation, we used the following parameters: R =., R =., R =.3, R =.3, ρ =.5, t =, σ =.5, σ =.5, x =, x =. 3 Conclusion In summary, based on a stochastic volatility model, we derived an approximate copula function by way of singular perturbation that was introduced by Fouque, Papanicolaou and Sircar () ]. During the approximation, however, in order to make the candidate probability density functions globally non-negative, instead of directly using the obtained perturbation result as in ], we introduced a multiplicative modification, namely the +tanh( ) form. It turns out that this modification is both necessary (to restore positiveness) and sufficient to guarantee the resulting functions to be density functions. Finally the resulting approximate copula the so-called perturbed Gaussian copula in this paper has a very desirable property compared to standard Gaussian copula: tail dependence at point (, ). Some numerical results were provided and they strongly supported the methods described above, both the singular perturbation and the modification. References ] J. P. Fouque, G. Papanicolaou and R. Sircar, Derivatives in Financial Markets with Stochastic Volatility, Cambridge University Press,. ] J.-P. Fouque, G. Papanicolaou, R. Sircar and K. Sølna. Singular Perturbations in Option Pricing. SIAM Journal on Applied Mathematics, 63(5), (3). 3] J.P. Fouque, R. Sircar and K. Solna, Stochastic Volatility Effects on Defaultable Bonds, working paper, to appear in Applied Mathematical Finance, 5. 4] J. P. Fouque, B. Wignall and X. Zhou, Modeling Correlated Defaults: First Passage Model under Stochastic Volatility, submitted, 6. 5] P. Schonbucher, Credit Derivative Pricing Models, Wiley, 3. 4
15 A Explicit Formulas As obtained in Section.3: u (t,x,x ) = π σ σ ( t) ρ { exp ( ρ ) (ξ x ) σ ( t) ρ(ξ x )(ξ x ) σ σ ( t) + (ξ x ) ]} σ. ( t) By a straightforward calculation and symmetry, we obtain 3 { u (ξ x ) x 3 = exp ( ρ ) σ ( t) ρ(ξ x )(ξ x ) + (ξ x ) ]} σ σ ( t) σ ( t) { (ξ x ) σ + ρ(ξ ] x ) 3 σ σ 4π σ σ 3 ( t) 3 ( ρ ) 5/ (ξ x ) σ + ρ(ξ ] } x ) 3 σ σ 6π σ σ ( t) 4 ( ρ ) 7/, 3 { u (ξ x ) x x = exp ( ρ ) σ ( t) ρ(ξ x )(ξ x ) + (ξ x ) ]} σ σ ( t) σ ( t) { ρ(ξ x ) (ξ ] x ) σ σ σ 4π σ σ 3 ( t) 3 ( ρ ) 5/ (ξ x ) σ + ρ(ξ ] x ) ρ σ σ π σ σ ( t)3 ( ρ ) 5/ (ξ x ) σ + ρ(ξ ] x ) ρ(ξ x ) (ξ ] x ) σ σ σ σ σ 3 u x 3 3 u x x { = exp { = exp 6π σ σ ( t) 4 ( ρ ) 7/ (ξ x ) ( ρ ) σ ( t) ρ(ξ x )(ξ x ) + (ξ x ) ]} σ σ ( t) σ ( t) { (ξ x ) σ + ρ(ξ ] x ) 3 σ σ 4π σ σ 3 ( t) 3 ( ρ ) 5/ (ξ x ) σ + ρ(ξ ] } x ) 3 σ σ 6π σ σ ( t) 4 ( ρ ) 7/, (ξ x ) ( ρ ) σ ( t) ρ(ξ x )(ξ x ) + (ξ x ) ]} σ σ ( t) σ ( t) { ρ(ξ x ) (ξ ] x ) σ σ σ 4π σ σ 3 ( t) 3 ( ρ ) 5/ (ξ x ) σ + ρ(ξ ] x ) ρ σ σ π σ σ ( t)3 ( ρ ) 5/ (ξ x ) σ + ρ(ξ ] x ) ρ(ξ x ) (ξ ] x ) σ σ σ σ σ } 6π σ σ ( t) 4 ( ρ ) 7/. 5 }.
Stochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationMultivariate Distribution Models
Multivariate Distribution Models Model Description While the probability distribution for an individual random variable is called marginal, the probability distribution for multiple random variables is
More informationSome Aspects of Universal Portfolio
1 Some Aspects of Universal Portfolio Tomoyuki Ichiba (UC Santa Barbara) joint work with Marcel Brod (ETH Zurich) Conference on Stochastic Asymptotics & Applications Sixth Western Conference on Mathematical
More informationVariance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models
Variance Reduction Techniques for Monte Carlo Simulations with Stochastic Volatility Models Jean-Pierre Fouque North Carolina State University SAMSI November 3, 5 1 References: Variance Reduction for Monte
More informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
More informationThe Uses of Differential Geometry in Finance
The Uses of Differential Geometry in Finance Andrew Lesniewski Bloomberg, November 21 2005 The Uses of Differential Geometry in Finance p. 1 Overview Joint with P. Hagan and D. Woodward Motivation: Varadhan
More informationAsymptotic behavior of Support Vector Machine for spiked population model
Asymptotic behavior of Support Vector Machine for spiked population model Supplementary Material This note contains technical proofs of Propositions, and 3 Proof of proposition Denote X [x,, x M ] T, y
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationKolmogorov Equations and Markov Processes
Kolmogorov Equations and Markov Processes May 3, 013 1 Transition measures and functions Consider a stochastic process {X(t)} t 0 whose state space is a product of intervals contained in R n. We define
More informationStochastic Calculus. Kevin Sinclair. August 2, 2016
Stochastic Calculus Kevin Sinclair August, 16 1 Background Suppose we have a Brownian motion W. This is a process, and the value of W at a particular time T (which we write W T ) is a normally distributed
More informationPoisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation
Poisson Jumps in Credit Risk Modeling: a Partial Integro-differential Equation Formulation Jingyi Zhu Department of Mathematics University of Utah zhu@math.utah.edu Collaborator: Marco Avellaneda (Courant
More informationarxiv: v2 [math.pr] 23 Jun 2014
COMPUTATION OF COPULAS BY FOURIER METHODS ANTONIS PAPAPANTOLEON arxiv:08.26v2 [math.pr] 23 Jun 204 Abstract. We provide an integral representation for the (implied) copulas of dependent random variables
More informationA New Generalized Gumbel Copula for Multivariate Distributions
A New Generalized Gumbel Copula for Multivariate Distributions Chandra R. Bhat* The University of Texas at Austin Department of Civil, Architectural & Environmental Engineering University Station, C76,
More informationExplicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables
Explicit Bounds for the Distribution Function of the Sum of Dependent Normally Distributed Random Variables Walter Schneider July 26, 20 Abstract In this paper an analytic expression is given for the bounds
More informationRandom Variables. Random variables. A numerically valued map X of an outcome ω from a sample space Ω to the real line R
In probabilistic models, a random variable is a variable whose possible values are numerical outcomes of a random phenomenon. As a function or a map, it maps from an element (or an outcome) of a sample
More informationGENERAL MULTIVARIATE DEPENDENCE USING ASSOCIATED COPULAS
REVSTAT Statistical Journal Volume 14, Number 1, February 2016, 1 28 GENERAL MULTIVARIATE DEPENDENCE USING ASSOCIATED COPULAS Author: Yuri Salazar Flores Centre for Financial Risk, Macquarie University,
More informationChapter 3 Second Order Linear Equations
Partial Differential Equations (Math 3303) A Ë@ Õæ Aë áöß @. X. @ 2015-2014 ú GA JË@ É Ë@ Chapter 3 Second Order Linear Equations Second-order partial differential equations for an known function u(x,
More informationReflected Brownian Motion
Chapter 6 Reflected Brownian Motion Often we encounter Diffusions in regions with boundary. If the process can reach the boundary from the interior in finite time with positive probability we need to decide
More informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationCopulas. MOU Lili. December, 2014
Copulas MOU Lili December, 2014 Outline Preliminary Introduction Formal Definition Copula Functions Estimating the Parameters Example Conclusion and Discussion Preliminary MOU Lili SEKE Team 3/30 Probability
More informationOptimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience
Optimal Trade Execution with Instantaneous Price Impact and Stochastic Resilience Ulrich Horst 1 Humboldt-Universität zu Berlin Department of Mathematics and School of Business and Economics Vienna, Nov.
More informationCONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS
EVA IV, CONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS Jose Olmo Department of Economics City University, London (joint work with Jesús Gonzalo, Universidad Carlos III de Madrid) 4th Conference
More information2.3. The heat equation. A function u(x, t) solves the heat equation if (with σ>0)
12 VOLKER BETZ 2.3. The heat equation. A function u(x, t) solves the heat equation if (with σ>) (2.1) t u 1 2 σ2 u =. In the simplest case, (2.1) is supposed to hold for all x R n, and all t>, and there
More informationHJB equations. Seminar in Stochastic Modelling in Economics and Finance January 10, 2011
Department of Probability and Mathematical Statistics Faculty of Mathematics and Physics, Charles University in Prague petrasek@karlin.mff.cuni.cz Seminar in Stochastic Modelling in Economics and Finance
More informationStability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games
Stability of Feedback Solutions for Infinite Horizon Noncooperative Differential Games Alberto Bressan ) and Khai T. Nguyen ) *) Department of Mathematics, Penn State University **) Department of Mathematics,
More informationEcon 508B: Lecture 5
Econ 508B: Lecture 5 Expectation, MGF and CGF Hongyi Liu Washington University in St. Louis July 31, 2017 Hongyi Liu (Washington University in St. Louis) Math Camp 2017 Stats July 31, 2017 1 / 23 Outline
More informationSymmetry and Exact Solutions of (2+1)-Dimensional Generalized Sasa Satsuma Equation via a Modified Direct Method
Commun. Theor. Phys. Beijing, China 51 2009 pp. 97 978 c Chinese Physical Society and IOP Publishing Ltd Vol. 51, No., June 15, 2009 Symmetry and Exact Solutions of 2+1-Dimensional Generalized Sasa Satsuma
More informationSolutions of the Financial Risk Management Examination
Solutions of the Financial Risk Management Examination Thierry Roncalli January 9 th 03 Remark The first five questions are corrected in TR-GDR and in the document of exercise solutions, which is available
More informationThe Instability of Correlations: Measurement and the Implications for Market Risk
The Instability of Correlations: Measurement and the Implications for Market Risk Prof. Massimo Guidolin 20254 Advanced Quantitative Methods for Asset Pricing and Structuring Winter/Spring 2018 Threshold
More informationPrelim Examination 2010 / 2011 (Assessing Units 1 & 2) MATHEMATICS. Advanced Higher Grade. Time allowed - 2 hours
Prelim Examination 00 / 0 (Assessing Units & ) MATHEMATICS Advanced Higher Grade Time allowed - hours Read Carefully. Calculators may be used in this paper.. Candidates should answer all questions. Full
More informationMATH4210 Financial Mathematics ( ) Tutorial 7
MATH40 Financial Mathematics (05-06) Tutorial 7 Review of some basic Probability: The triple (Ω, F, P) is called a probability space, where Ω denotes the sample space and F is the set of event (σ algebra
More informationA new approach for investment performance measurement. 3rd WCMF, Santa Barbara November 2009
A new approach for investment performance measurement 3rd WCMF, Santa Barbara November 2009 Thaleia Zariphopoulou University of Oxford, Oxford-Man Institute and The University of Texas at Austin 1 Performance
More informationSystemic Risk and Stochastic Games with Delay
Systemic Risk and Stochastic Games with Delay Jean-Pierre Fouque (with René Carmona, Mostafa Mousavi and Li-Hsien Sun) PDE and Probability Methods for Interactions Sophia Antipolis (France) - March 30-31,
More informationGoodness-of-Fit Testing with Empirical Copulas
with Empirical Copulas Sami Umut Can John Einmahl Roger Laeven Department of Econometrics and Operations Research Tilburg University January 19, 2011 with Empirical Copulas Overview of Copulas with Empirical
More informationMultiple Random Variables
Multiple Random Variables Joint Probability Density Let X and Y be two random variables. Their joint distribution function is F ( XY x, y) P X x Y y. F XY ( ) 1, < x
More informationThomas Knispel Leibniz Universität Hannover
Optimal long term investment under model ambiguity Optimal long term investment under model ambiguity homas Knispel Leibniz Universität Hannover knispel@stochastik.uni-hannover.de AnStAp0 Vienna, July
More informationMathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio ( )
Mathematical Methods for Neurosciences. ENS - Master MVA Paris 6 - Master Maths-Bio (2014-2015) Etienne Tanré - Olivier Faugeras INRIA - Team Tosca November 26th, 2014 E. Tanré (INRIA - Team Tosca) Mathematical
More informationMTH739U/P: Topics in Scientific Computing Autumn 2016 Week 6
MTH739U/P: Topics in Scientific Computing Autumn 16 Week 6 4.5 Generic algorithms for non-uniform variates We have seen that sampling from a uniform distribution in [, 1] is a relatively straightforward
More informationA MODEL FOR THE LONG-TERM OPTIMAL CAPACITY LEVEL OF AN INVESTMENT PROJECT
A MODEL FOR HE LONG-ERM OPIMAL CAPACIY LEVEL OF AN INVESMEN PROJEC ARNE LØKKA AND MIHAIL ZERVOS Abstract. We consider an investment project that produces a single commodity. he project s operation yields
More informationOrnstein-Uhlenbeck processes for geophysical data analysis
Ornstein-Uhlenbeck processes for geophysical data analysis Semere Habtemicael Department of Mathematics North Dakota State University May 19, 2014 Outline 1 Introduction 2 Model 3 Characteristic function
More informationDependence. Practitioner Course: Portfolio Optimization. John Dodson. September 10, Dependence. John Dodson. Outline.
Practitioner Course: Portfolio Optimization September 10, 2008 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y ) (x,
More informationAffine Processes. Econometric specifications. Eduardo Rossi. University of Pavia. March 17, 2009
Affine Processes Econometric specifications Eduardo Rossi University of Pavia March 17, 2009 Eduardo Rossi (University of Pavia) Affine Processes March 17, 2009 1 / 40 Outline 1 Affine Processes 2 Affine
More informationLecture 3 Partial Differential Equations
Lecture 3 Partial Differential Equations Prof. Massimo Guidolin Prep Course in Investments August-September 2016 Plan of the lecture Motivation and generalities The heat equation and its applications in
More informationModelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich
Modelling Dependence with Copulas and Applications to Risk Management Filip Lindskog, RiskLab, ETH Zürich 02-07-2000 Home page: http://www.math.ethz.ch/ lindskog E-mail: lindskog@math.ethz.ch RiskLab:
More informationFormulas for probability theory and linear models SF2941
Formulas for probability theory and linear models SF2941 These pages + Appendix 2 of Gut) are permitted as assistance at the exam. 11 maj 2008 Selected formulae of probability Bivariate probability Transforms
More informationModelling and Estimation of Stochastic Dependence
Modelling and Estimation of Stochastic Dependence Uwe Schmock Based on joint work with Dr. Barbara Dengler Financial and Actuarial Mathematics and Christian Doppler Laboratory for Portfolio Risk Management
More informationON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER
ON THE FIRST TIME THAT AN ITO PROCESS HITS A BARRIER GERARDO HERNANDEZ-DEL-VALLE arxiv:1209.2411v1 [math.pr] 10 Sep 2012 Abstract. This work deals with first hitting time densities of Ito processes whose
More informationWiener Measure and Brownian Motion
Chapter 16 Wiener Measure and Brownian Motion Diffusion of particles is a product of their apparently random motion. The density u(t, x) of diffusing particles satisfies the diffusion equation (16.1) u
More informationFolland: Real Analysis, Chapter 8 Sébastien Picard
Folland: Real Analysis, Chapter 8 Sébastien Picard Problem 8.3 Let η(t) = e /t for t >, η(t) = for t. a. For k N and t >, η (k) (t) = P k (/t)e /t where P k is a polynomial of degree 2k. b. η (k) () exists
More information04. Random Variables: Concepts
University of Rhode Island DigitalCommons@URI Nonequilibrium Statistical Physics Physics Course Materials 215 4. Random Variables: Concepts Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative
More informationDynamic Risk Measures and Nonlinear Expectations with Markov Chain noise
Dynamic Risk Measures and Nonlinear Expectations with Markov Chain noise Robert J. Elliott 1 Samuel N. Cohen 2 1 Department of Commerce, University of South Australia 2 Mathematical Insitute, University
More informationPerhaps the simplest way of modeling two (discrete) random variables is by means of a joint PMF, defined as follows.
Chapter 5 Two Random Variables In a practical engineering problem, there is almost always causal relationship between different events. Some relationships are determined by physical laws, e.g., voltage
More informationAlberto Bressan. Department of Mathematics, Penn State University
Non-cooperative Differential Games A Homotopy Approach Alberto Bressan Department of Mathematics, Penn State University 1 Differential Games d dt x(t) = G(x(t), u 1(t), u 2 (t)), x(0) = y, u i (t) U i
More informationNumerical Methods with Lévy Processes
Numerical Methods with Lévy Processes 1 Objective: i) Find models of asset returns, etc ii) Get numbers out of them. Why? VaR and risk management Valuing and hedging derivatives Why not? Usual assumption:
More informationDynamic Bertrand and Cournot Competition
Dynamic Bertrand and Cournot Competition Effects of Product Differentiation Andrew Ledvina Department of Operations Research and Financial Engineering Princeton University Joint with Ronnie Sircar Princeton-Lausanne
More information100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX
100 CHAPTER 4. SYSTEMS AND ADAPTIVE STEP SIZE METHODS APPENDIX.1 Norms If we have an approximate solution at a given point and we want to calculate the absolute error, then we simply take the magnitude
More informationModulation of symmetric densities
1 Modulation of symmetric densities 1.1 Motivation This book deals with a formulation for the construction of continuous probability distributions and connected statistical aspects. Before we begin, a
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
More informationFourier Transform & Sobolev Spaces
Fourier Transform & Sobolev Spaces Michael Reiter, Arthur Schuster Summer Term 2008 Abstract We introduce the concept of weak derivative that allows us to define new interesting Hilbert spaces the Sobolev
More informationChap. 1. Some Differential Geometric Tools
Chap. 1. Some Differential Geometric Tools 1. Manifold, Diffeomorphism 1.1. The Implicit Function Theorem ϕ : U R n R n p (0 p < n), of class C k (k 1) x 0 U such that ϕ(x 0 ) = 0 rank Dϕ(x) = n p x U
More informationDependence. MFM Practitioner Module: Risk & Asset Allocation. John Dodson. September 11, Dependence. John Dodson. Outline.
MFM Practitioner Module: Risk & Asset Allocation September 11, 2013 Before we define dependence, it is useful to define Random variables X and Y are independent iff For all x, y. In particular, F (X,Y
More information3. Probability and Statistics
FE661 - Statistical Methods for Financial Engineering 3. Probability and Statistics Jitkomut Songsiri definitions, probability measures conditional expectations correlation and covariance some important
More informationMATH 131P: PRACTICE FINAL SOLUTIONS DECEMBER 12, 2012
MATH 3P: PRACTICE FINAL SOLUTIONS DECEMBER, This is a closed ook, closed notes, no calculators/computers exam. There are 6 prolems. Write your solutions to Prolems -3 in lue ook #, and your solutions to
More informationBackward martingale representation and endogenous completeness in finance
Backward martingale representation and endogenous completeness in finance Dmitry Kramkov (with Silviu Predoiu) Carnegie Mellon University 1 / 19 Bibliography Robert M. Anderson and Roberto C. Raimondo.
More informationResearch Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance
Advances in Decision Sciences Volume, Article ID 893497, 6 pages doi:.55//893497 Research Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance Junichi Hirukawa and
More informationEQUITY MARKET STABILITY
EQUITY MARKET STABILITY Adrian Banner INTECH Investment Technologies LLC, Princeton (Joint work with E. Robert Fernholz, Ioannis Karatzas, Vassilios Papathanakos and Phillip Whitman.) Talk at WCMF6 conference,
More informationInformation and Credit Risk
Information and Credit Risk M. L. Bedini Université de Bretagne Occidentale, Brest - Friedrich Schiller Universität, Jena Jena, March 2011 M. L. Bedini (Université de Bretagne Occidentale, Brest Information
More informationMoment Properties of Distributions Used in Stochastic Financial Models
Moment Properties of Distributions Used in Stochastic Financial Models Jordan Stoyanov Newcastle University (UK) & University of Ljubljana (Slovenia) e-mail: stoyanovj@gmail.com ETH Zürich, Department
More informationOn Reparametrization and the Gibbs Sampler
On Reparametrization and the Gibbs Sampler Jorge Carlos Román Department of Mathematics Vanderbilt University James P. Hobert Department of Statistics University of Florida March 2014 Brett Presnell Department
More informationOn a class of stochastic differential equations in a financial network model
1 On a class of stochastic differential equations in a financial network model Tomoyuki Ichiba Department of Statistics & Applied Probability, Center for Financial Mathematics and Actuarial Research, University
More information1 Review of di erential calculus
Review of di erential calculus This chapter presents the main elements of di erential calculus needed in probability theory. Often, students taking a course on probability theory have problems with concepts
More informationChapter 4: First-order differential equations. Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey
Chapter 4: First-order differential equations Similarity and Transport Phenomena in Fluid Dynamics Christophe Ancey Chapter 4: First-order differential equations Phase portrait Singular point Separatrix
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationSymmetry Methods for Differential Equations and Conservation Laws. Peter J. Olver University of Minnesota
Symmetry Methods for Differential Equations and Conservation Laws Peter J. Olver University of Minnesota http://www.math.umn.edu/ olver Santiago, November, 2010 Symmetry Groups of Differential Equations
More informationOPTIMAL STOPPING OF A BROWNIAN BRIDGE
OPTIMAL STOPPING OF A BROWNIAN BRIDGE ERIK EKSTRÖM AND HENRIK WANNTORP Abstract. We study several optimal stopping problems in which the gains process is a Brownian bridge or a functional of a Brownian
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationStrong log-concavity is preserved by convolution
Strong log-concavity is preserved by convolution Jon A. Wellner Abstract. We review and formulate results concerning strong-log-concavity in both discrete and continuous settings. Although four different
More informationFinancial Econometrics and Volatility Models Copulas
Financial Econometrics and Volatility Models Copulas Eric Zivot Updated: May 10, 2010 Reading MFTS, chapter 19 FMUND, chapters 6 and 7 Introduction Capturing co-movement between financial asset returns
More informationQuestion 1. The correct answers are: (a) (2) (b) (1) (c) (2) (d) (3) (e) (2) (f) (1) (g) (2) (h) (1)
Question 1 The correct answers are: a 2 b 1 c 2 d 3 e 2 f 1 g 2 h 1 Question 2 a Any probability measure Q equivalent to P on F 2 can be described by Q[{x 1, x 2 }] := q x1 q x1,x 2, 1 where q x1, q x1,x
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationShort-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility
Short-time expansions for close-to-the-money options under a Lévy jump model with stochastic volatility José Enrique Figueroa-López 1 1 Department of Statistics Purdue University Statistics, Jump Processes,
More informationJump-type Levy Processes
Jump-type Levy Processes Ernst Eberlein Handbook of Financial Time Series Outline Table of contents Probabilistic Structure of Levy Processes Levy process Levy-Ito decomposition Jump part Probabilistic
More information= F (b) F (a) F (x i ) F (x i+1 ). a x 0 x 1 x n b i
Real Analysis Problem 1. If F : R R is a monotone function, show that F T V ([a,b]) = F (b) F (a) for any interval [a, b], and that F has bounded variation on R if and only if it is bounded. Here F T V
More information5 December 2016 MAA136 Researcher presentation. Anatoliy Malyarenko. Topics for Bachelor and Master Theses. Anatoliy Malyarenko
5 December 216 MAA136 Researcher presentation 1 schemes The main problem of financial engineering: calculate E [f(x t (x))], where {X t (x): t T } is the solution of the system of X t (x) = x + Ṽ (X s
More informationStability of optimization problems with stochastic dominance constraints
Stability of optimization problems with stochastic dominance constraints D. Dentcheva and W. Römisch Stevens Institute of Technology, Hoboken Humboldt-University Berlin www.math.hu-berlin.de/~romisch SIAM
More informationMATH 425, FINAL EXAM SOLUTIONS
MATH 425, FINAL EXAM SOLUTIONS Each exercise is worth 50 points. Exercise. a The operator L is defined on smooth functions of (x, y by: Is the operator L linear? Prove your answer. L (u := arctan(xy u
More informationDifferential Equations for Dyson Processes
Differential Equations for Dyson Processes Joint work with Harold Widom I. Overview We call Dyson process any invariant process on ensembles of matrices in which the entries undergo diffusion. Dyson Brownian
More informationIntroduction to multiscale modeling and simulation. Explicit methods for ODEs : forward Euler. y n+1 = y n + tf(y n ) dy dt = f(y), y(0) = y 0
Introduction to multiscale modeling and simulation Lecture 5 Numerical methods for ODEs, SDEs and PDEs The need for multiscale methods Two generic frameworks for multiscale computation Explicit methods
More informationStein s method and zero bias transformation: Application to CDO pricing
Stein s method and zero bias transformation: Application to CDO pricing ESILV and Ecole Polytechnique Joint work with N. El Karoui Introduction CDO a portfolio credit derivative containing 100 underlying
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationRFID & Item-Level Information Visibility
Department of Information System & Operations Management University of Florida Gainesville, Florida 32601 wzhou@ufl.edu December 13, 2008 1. Manufacturing an Engine Engine Builder Replenish Policy {A1,A2}
More informationIntroduction to Probability and Statistics (Continued)
Introduction to Probability and Statistics (Continued) Prof. icholas Zabaras Center for Informatics and Computational Science https://cics.nd.edu/ University of otre Dame otre Dame, Indiana, USA Email:
More informationSystemic Risk and Stochastic Games with Delay
Systemic Risk and Stochastic Games with Delay Jean-Pierre Fouque Collaborators: R. Carmona, M. Mousavi, L.-H. Sun, N. Ning and Z. Zhang Mean Field Games 2017 IPAM at UCLA JP Fouque (UC Santa Barbara) Stochastic
More informationMultivariate Statistics
Multivariate Statistics Chapter 2: Multivariate distributions and inference Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2016/2017 Master in Mathematical
More informationNumerical Integration of SDEs: A Short Tutorial
Numerical Integration of SDEs: A Short Tutorial Thomas Schaffter January 19, 010 1 Introduction 1.1 Itô and Stratonovich SDEs 1-dimensional stochastic differentiable equation (SDE) is given by [6, 7] dx
More informationStrictly Stationary Solutions of Autoregressive Moving Average Equations
Strictly Stationary Solutions of Autoregressive Moving Average Equations Peter J. Brockwell Alexander Lindner Abstract Necessary and sufficient conditions for the existence of a strictly stationary solution
More informationPerturbative Approaches for Robust Intertemporal Optimal Portfolio Selection
Perturbative Approaches for Robust Intertemporal Optimal Portfolio Selection F. Trojani and P. Vanini ECAS Course, Lugano, October 7-13, 2001 1 Contents Introduction Merton s Model and Perturbative Solution
More informationf (x) = k=0 f (0) = k=0 k=0 a k k(0) k 1 = a 1 a 1 = f (0). a k k(k 1)x k 2, k=2 a k k(k 1)(0) k 2 = 2a 2 a 2 = f (0) 2 a k k(k 1)(k 2)x k 3, k=3
1 M 13-Lecture Contents: 1) Taylor Polynomials 2) Taylor Series Centered at x a 3) Applications of Taylor Polynomials Taylor Series The previous section served as motivation and gave some useful expansion.
More information