Worst VaR Scenarios. Paul Embrechts a Andrea Höing a, Giovanni Puccetti b. 1 Introduction

Size: px
Start display at page:

Download "Worst VaR Scenarios. Paul Embrechts a Andrea Höing a, Giovanni Puccetti b. 1 Introduction"

Transcription

1 Worst VaR Scenarios Paul Embrechts a Anrea Höing a, Giovanni Puccetti b a Department of Mathematics, ETH Zürich, CH-8092 Zürich, Switzerlan b Department of Mathematics for Decisions, University of Firenze, Firenze, Italy Abstract The worst-possible Value-at-Risk for a non-ecreasing function ψ of n epenent risks is known when n = 2 or the copula of the portfolio is boune from below. In this paper we analyze the properties of the epenence structures leaing to this solution, in particular their form an the implie functional epenence between the marginals. Furthermore we criticise the assumption of the worst-possible scenario for VaR-base risk management an we provie an alternative approach supporting comonotonicity. Key wors: Value-at-Risk, epenent risks, copulas, comonotonic risks J.E.L. Subject Classification: G10 Subj, Class.: IM01, IM12, IM52 1 Introuction Consier an insurer holing a portfolio consisting of n policies with iniviual risks X 1,...,X n over a fixe time perio. Given some measurable function ψ : R n R, a relevant task in insurance mathematics is the investigation of the risk position associate with ψ(x 1,...,X n ), when the marginal istributions of the single risks are known. Actuarial examples of the function ψ inclue n i=1 x i, simply characterizing the aggregate claim amount eriving from the policies or n i=1 h i (x i ) an h( n i=1 x i ), proviing the risk positions for a reinsurance treaty with retention functions h i, i = 1,...,n an a global reinsurance treaty with global retention function h, respectively. Corresponing author. Tel.: ; fax: aress: ahoeing@math.ethz.ch (Anrea Höing). 1 The authors woul like to thank Alessanro Juri (ETH Zürich/UBS) an Pablo Parrilo (ETH Zürich) for helpful iscussions. Preprint submitte to Elsevier Science

2 The problem of fining the best-possible lower boun on the istribution function (f) of ψ(x 1,...,X n ) has receive a consierable interest in insurance mathematics. From a financial risk management point of view, the problem is equivalent to fining the worst-possible Value-at-Risk (VaR) for the corresponing aggregate position. Here we refer to the Introuction in Embrechts an Puccetti (2004) for etails. Moelling the interepenence arising in a ranom portfolio calls for the use of copulas. If a lower boun on the copula of the vector (X 1,...,X n ) is given, the above problem is fully solve an the bouns provie in Embrechts et al. (2003) are sharp. In the no-information case the latter o hol only if n = 2. Rather than treating the technical proof of such results, for which we refer to the above cite references, in this paper we analyse in more etails the properties of their solutions. We concentrate mainly on the no-information case, when a lower boun on the copula of the portfolio is not available an a solution is known only for twoimensional portfolios. Without loss of generality, we stuy the optimizing copula for the sum of two epenent risks, which is well-known to iffer from comonotonicity. In particular we iscuss its shape, its implications in terms of epenence an we criticise it as not being a rational scenario for an insurance company. Finally, we provie an alternative optimization approach leaing to a suitable measure of risk, which supports the assumption of comonotonicity for a pruent evaluation of the VaR for the aggregate position. 2 Preliminaries an funamental results In this section we present some well-known concepts about copulas an briefly recall the funamental results about the problem of bouning the VaR for functions of epenent risks. For more etails about copulas, we refer to Nelsen (1999). 2.1 Value-at-Risk an epenence structures On some probability space (Ω, A, È), let the ranom vector X := (X 1,...,X n ) represent a portfolio of risks. Given a measurable function ψ : R n R we face the problem of fining the supremum of the VaR for the aggregate position ψ(x) over the class of possible fs for X having fixe marginals F 1,...,F n. Definition 1 Let ϕ : R R be a non-ecreasing function. Its generalize left continuous inverse is the function ϕ 1 : R R efine by ϕ 1 (y) := inf{x R ϕ(x) y}. For 0 α 1 the Value-at-Risk at probability level α for a ranom variable Y with istribution function G is its α-quantile, i.e. VaR α (Y ) := G 1 (α). 2

3 Of course, quantiles of the f of ψ(x) can be compute if the joint istribution function F(x 1,...,x n ) = È[X 1 x 1,..., X n x n ] is known. At this point, the notion of copula becomes useful. Definition 2 An n -imensional copula is an n-imensional istribution function restricte to [0, 1] n having stanar uniform marginals. We enote with C n the family of n-imensional copulas. Given a copula C C n an a set of univariate marginals F 1,...,F n, we can always efine a f F on R n having these marginals by F(x 1,...,x n ) := C(F 1 (x 1 ),...,F n (x n )). (1) Hence, given n fs F 1,..., F n, we let X C = (X 1,..., X n ) be the ranom vector on R n having a copula C satisfying (1). Observe that this copula is unique for continuous marginal fs. Conversely, Sklar s Theorem (Sklar (1973), Theorem 1) states that there always exists C C n coupling the marginals of a fixe f F trough (1). Comparing copulas pointwise an efining the riskiness of a epenence structure through this comparison, we recall that any copula C lies between the lower an upper Fréchet bouns W(u 1,...,u n ) := ( n i=1 u i n + 1) + an M(u 1,...,u n ) := min 1 i n u i, namely W C M. (2) Observe that, contrary to M, the lower Fréchet boun W is not a istribution function for n > 2. Ranom variables couple through C = M (C = W, respectively) are calle comonotonic (countermonotonic). The inepenence copula is enote by Π(u 1,..., u n ) := n i=1 u i. Remark 3 Comonotonicity characterizes the risks of the portfolio as being increasing functions of a common ranom variable. It is therefore a strong epenence an measure of epenence such as Kenall s τ or Spearman s ρ will escribe M as a perfect structure, i.e. τ(m) = ρ(m) = 1. It is precisely this representation which motivates the use of the concept of comonotonicity in financial applications. Moreover, assuming comonotonicity leas to almost all the computational benefits of inepenence, yieling, in aition, a pruent scenario in many contexts as we will emphasize in Section 4. For an in epth iscussion of comonotonicity, see Dhaene et al. (2001). 2.2 Bouns on value at risk for functions of epenent risks We now recall the two funamental results being the object of our analysis. For a proof of both theorems an further iscussions, we refer to Embrechts an Puccetti 3

4 (2004) an the references therein. For a copula C an marginals F 1,...,F n, efine σ C,ψ (F 1,...,F n )(s) := C(F 1 (x 1 ),...,F n (x n )), {ψ<s} τ C,ψ (F 1,...,F n )(s) := sup C(F 1 (x 1 ),..., F n 1 (x n 1 ), Fn (ψ x n (s))), x 1,...,x n 1 R where ψ x n (s) := sup{x n R ψ(x n, x n ) < s} for x n := (x 1,...,x n 1 ) R n 1. In the following, we refer to non-ecreasing functions ψ : R n R as being non-ecreasing in each component. Remark 4 Observe that σ C,ψ (F 1,...,F n )(s) = È[ψ(X C ) < s] for X C having marginals F 1,...,F n. In the Appenix, Proposition 23, we show that the operator τ C,ψ in (3) is actually the left-continuous version of a f, i.e. there exists a ranom variable K with È[K < s] = τ CL,ψ(F 1,...,F n )(s). This result extens a claim of Denuit et al. (1999), p. 37. As first note in Schweizer an Sklar (1974) for the sum of two risks, if C L M there oes not exist a measurable real function g such that K = g(x), with X having marginals F 1,...,F n. Theorem 5 Let X C = (X 1,..., X n ) be a ranom vector on R n (n > 1) having marginal istribution functions F 1,...,F n. Assume that there exists a copula C L such that C C L. If ψ : R n R is non-ecreasing, then for every real s we have σ C,ψ (F 1,...,F n )(s) τ CL,ψ(F 1,...,F n )(s). (3) Translate in the language of VaR, the above statement becomes for every α in the unit interval. VaR α (ψ(x 1,...,X n )) τ CL,ψ(F 1,...,F n ) 1 (α) The bouns state in Theorem 5 are pointwise best-possible an cannot be tightene if n = 2 or a lower boun C L > W on the copula of the portfolio X C is assume. Theorem 6 Further to the hypotheses of Theorem 5, we assume that ψ is also right-continuous in its last argument. Define the copula C α : [0, 1] n [0, 1] as max{α, C L (u)}, if u = (u 1,..., u n ) [α, 1] n, C α (u) := min{u 1,...,u n }, otherwise, where α = τ CL,ψ(F 1,...,F n )(s). Then this copula attains boun (3), i.e. σ Cα,ψ(F 1,..., F n )(s) = α. (4) The latter theorem motivates the investigation of the epenence structures leaing to the worst-case VaR scenario when n = 2 or C L > W. For C L = W an n > 2 the boun state in (3) is still vali but no more sharp. 4

5 3 Analysis of the worst-case portfolios The aim of the present paper is to give more insight into the shape of the copula yieling the worst-possible VaR for ψ(x C ) = ψ(x 1,...,X n ) an to unerstan the implie epenence between the marginals. Uner all possible epenence structures, the worst-case scenario for the VaR at level α is given by the copula minimizing È[ψ(X C ) < s] over s-regions epening on α. Inee, accoring to Definition 1 with m ψ (s) := inf C C n{è[ψ(xc ) < s]}, s R, (5) we have that VaR α (ψ(x C )) m 1 ψ (α), α [0, 1]. The problem at han becomes also the characterization of the copula minimizing m ψ, or equivalently maximizing m ψ (s) = 1 m ψ (s) = sup C C n {È[ψ(X C ) s]}, s R. (6) Such a copula will be referre to as a worst-case scenario for the aggregate position ψ(x C ). We use the term scenario to inicate a (possibly egenerate) set of probability measures in line with Artzner et al. (1999). Analogous to the above efinitions, in the presence of partial information, we write m CL,ψ (respectively m CL,ψ) an the infimum (supremum) is taken over all C C n satisfying the bounary conition C C L. In the next subsections, we concentrate on the sum of risks (generalizations to nonecreasing continuous functions ψ being straightforwar) an we choose C L = W. See however Section 5 for some comments on the latter choice of no epenence information. 3.1 Two-imensional portfolios If we take two risks, the boun given in Theorem 5 cannot be tightene an there always exists a two-imensional copula meeting that boun at a given point s. We restate Theorem 6 in this particular case. Theorem 7 Let X C = (X 1, X 2 ) be a ranom vector on R 2 having marginal istribution functions F 1, F 2. Define the copula C α : [0, 1] 2 [0, 1], max{α, W(u)}, if u = (u 1, u 2 ) [α, 1] 2, C α (u) := min{u 1, u 2 }, otherwise, where α = τ W,+ (F 1, F 2 )(s). Then this copula attains boun (3), i.e. σ Cα,+(F 1, F 2 )(s) = α. (7) 5

6 γ B s 0.9 γ B s A s 0.7 A s Fig. 1. Support of the copula C α, sets A s,b s an curve γ for: N(0,1)-N(1,2)-normal marginals an s = 1 (which gives α = ) (left); LN(0.4,1)-LN(0.4,1)-marginals an s = 4 (α = ) (right). Proofs of Theorem 7 can be foun in Frank et al. (1987) an Rüschenorf (1982). Our aim here is to restate the problem of maximizing (6) from a geometric point of view an illustrate the properties of the optimizing copulas leaing to the worstcase scenario for VaR. Without loss of generality, in what follows, we take continuous, increasing marginals. Let moreover G s := {(x 1, x 2 ) R 2 x 1 + x 2 s} an h : R 2 [0, 1] 2, h(x 1, x 2 ) := (F 1 (x 1 ), F 2 (x 2 )). The basic iea is to use the function h to transport the optimization problem on the unit square [0, 1] 2. In fact, U C := h(x C ) is a ranom vector, with istribution function C on [0, 1] 2. The function h is invertible, hence we have that È[X C G s ] = È[h(X C ) h(g s )] = µ C (A s ), where µ C is the measure corresponing to C on [0, 1] 2 an A s := h(g s ) = {(u 1, u 2 ) [0, 1] 2 F 1 1 (u 1 ) + F 1 2 (u 2 ) s}. The maximization function (6) can now be rewritten as m + (s) = sup C C 2 {µ C (A s )}. (8) For α = 1, (3) leas to σ C,+ (F 1, F 2 )(s) = 1 for every copula C, hence take α [0, 1). The bounary of A s is the image of the curve γ : R [0, 1] 2, γ(t) := (F 1 (t), F 2 (s t)). In Figure 1 the curve γ elimiting the set A s is rawn, with the support of the copula C α, in case of normal (N) an log-normal (LN) marginals. 6

7 The copula C α is uniformly istribute on its support, hence, efining B s := {(u 1, u 2 ) [0, 1] 2 u 1 + u 2 = 1 + α} we have µ Cα (B s ) = 1 α. As note in Nelsen (1999), p. 187, this is the crucial property leaing to the statement of Theorem 7. In fact, when 0 < α < 1, the continuity of the F i s implies that α = τ W,+ (F 1, F 2 )(s) = F 1 (x 1 ) + F 2(s x 1 ) 1 (9) for some x 1. Hence the curve γ meets the segment B s at least in one point. The technical (an for general n an C L > W rather laborious) part of the proof consists in showing that γ always lies below the segment B s, hence A s B s an µ Cα (A s ) µ Cα (B s ) = 1 α. Noting that µ Cα (A s ) 1 α, from Theorem 5 we obtain (7). For α = 0, instea, the existence of a tangent point between γ an B s is not necessary, since the copula W yiels the theorem. Analogous geometric consierations can be given for the case C L > W an for non-ecreasing continuous ψ. Remark 8 Observe that the geometric properties of the support of C α, illustrate in Figure 1, can be extene to a whole family of copulas, which implies that the epenence structure leaing to the worst-case VaR is not unique. Let Ĉ2 α an C2 α enote the family of copulas leaing to the worst possible VaR an the family of copulas sharing their support on [α, 1] 2 with C α, respectively. Formally: Ĉ 2 α : = {C C2 σ C,+ (F 1, F 2 )(s) = α}, C 2 α : = {C C2 C(u 1, u 2 ) = C α (u 1, u 2 ) for α u 1, u 2 1}. Observe that we can write Ĉ2 α = {C C2 µ C (A s ) = µ Cα (A s )}. In particular, it trivially follows that, every copula in C 2 α attains boun (3), since C 2 α Ĉ2 α. We now focus on the epenence implie by the copulas in C 2 α. The support R α := {(u 1, u 2 ) [0, α) 2 u 1 = u 2 } {(u 1, u 2 ) [α, 1] 2 u 1 + u 2 = 1 + α} of the copula C α implicitly efines the epenence of the couple ranom variables by the substitution u i = F i (x i ), i = 1, 2. In fact, if the copula C α couples X 1 an X 2 into the ranom vector X Cα an if we assume F 1, F 2 to be increasing on their omain, then we have X 2 = g(x 1 ), where the function g : R R is efine as F 2 1 (F 1 (x)), if x < F1 1 (α), g(x) := (10) F 2 1 (1 + α F 1 (x)), otherwise. 7

8 Analogously, every other copula in C 2 α efines a functional epenence ientical to that of g for x F1 1 (α). For example, the copula Cα 1 given by Cα 1 (u max{c L (u 1, u 2 ), α}, when (u 1, u 2 ) [α, 1] 2, 1, u 2 ) := otherwise, u 1 u 2 α, couples two marginals, which are inepenent if the first lies below the threshol F1 1 (α) an behaves like C α otherwise. Figure 2 compares R α with the support of Cα Fig. 2. Supports of the copulas C 0.5 (left) an C (right). Merging the two marginals by C α is therefore equivalent to letting X 2 = g(x 1 ). Hence, the two risks are mutually completely epenent. Moreover, the copula C α is a so-calle shuffle-of-m an hence implies that X 1 an X 2 are strongly piecewise strictly monotone functions of each other, in the sense efine in Mikusiński et al. (1991). Nevertheless, measures of epenence such as Kenall s τ or Spearman s ρ escribe C α as a non-perfect structure when 0 α < 1, i.e. τ(c α ), ρ(c α ) < 1. This is ue to the fact that this copula only represents piecewise comonotonicity. Mathematically, the epenence structure inuce by C α is, however, as strong a the one inuce by M, since the two variables couple by C α are in a oneto-one corresponence. Finally, note that every f on R 2 efine by applying a Ĉ 2 α-copula to the given set of marginals has a singular component, i.e. is mixe with a continuous istribution having zero erivative except for a set of Lebesgue measure zero. For instance, C α is singular on its whole omain, whereas Cα 1 only on [α, 1] 2. For more etails about singular istribution functions see Billingsley (1995), Section 31 an Nelsen (1999), p. 23. At this point, it is relevant to note that, in general, M / Ĉ2 α when 0 α < 1, the case α = 1 being the trivial one in which Ĉ2 α = C2. This provies a further geometric proof that comonotonicity oes not lea to the worst possible scenario for VaR an emphasizes the non-coherence of VaR as state in Artzner et al. (1999). Suppose that X 1 an X 2 are ientically istribute with unboune, absolutely continuous f having positive ensity f. If f is eventually ecreasing, it is easy to 8

9 show that for s large enough we have that α = 2F(s/2) 1, while σ M,+ (F, F) = F(s/2) > α. (11) A necessary conition for M to be in Ĉ2 α is that the point (α, α) lies in [0, 1]2 \ A s. Equation (11) implies that this conition is not satisfie for s large enough. Finally, M Ĉ2 0 if an only if A s = [0, 1] 2, i.e. the sum X 1 + X 2 is È-a.s. boune from below by the threshol s. In this case the problem of bouning the VaR for the sum oes not arise. We conclue that, apart from pathological cases of no actuarial importance, we have that σ M,+ (F 1, F 2 )(s) > m + (s) when 0 α < 1. This equation shows that the assumption of comonotonicity among the risks of the portfolio may lea to a angerous uner-valuation of the VaR for the aggregate position. At first, the worst epenence scenario coul seem to be the one implie by M, since uner comonotonicity it is inee known that every ranom variable is a non-ecreasing function of the other, so that high values for the first imply high values for the secon. Theorem 6 provies a eeper view on this issue, stating, instea, that for every threshol s such that α < 1, there exists a copula C α yieling a value for the VaR which is higher than that of comonotonicity. The following example further stresses the fact that M oes not belong, in general, to Ĉ2 α. Example 9 Let X 1 be normally istribute N(0, 1) with f Φ an put X 2 = X 1 to obtain È[X 1 + X 2 = 0] = 1. The copula escribing this epenence is the countermonotonic copula W, uner which X 2 is a non-increasing function of X 1. Accoring to Theorem 6, m + (0) = 0. In this set-up, the maximizing solution of (8) is then the structure of epenence which is opposite to comonotonicity (note that happens whenever α = 0), for which we have instea: σ M,+ (Φ, Φ)(0) = 1/2. Figure 3 (left) illustrates how, in the case of stanar normal marginals, for every positive s R, there exists a copula C Ĉ2 α such that σ C,+ (Φ, Φ)(s) < σ M,+ (Φ, Φ)(s). In the same figure (right) we also provie the shape of the bivariate istribution obtaine by applying C α to stanar normal marginals for s = (α = ). The reaer shoul compare this figure to Figure 2 (right). 9

10 inepenence comonotonicity lower boun = worst case Fig. 3. Range for È[X 1 + X 2 < s] for a N(0,1)-N(0,1)-portfolio. Together with the inepenence an comonotonic value we plot the lower boun m + (s) (left); the ensity of the istribution of (X 1,X 2 ) obtaine by applying the copula C to a N(0,1)-N(0,1)-portfolio is given in the figure on the right. 3.2 Two-imensional uniform portfolios We now state some simple results for uniform marginals that will turn out to be useful in unerstaning the n-imensional case. Proposition 10 Let the hypotheses of Theorem 7 be satisfie with F 1, F 2 uniformly istribute on the unit interval. Then Ĉ 2 α = C2 α. (12) Proof If α = 1, Ĉ2 1 = C 2 1 = C 2. Let 0 α = s 1 < 1 an C Ĉ2 α. 1 1 H20 H 10 A s D 1 H21 H 00 E 2 u 2 u 2 H22 H 11 D 2 H23 α α E 1 0 α 1 0 u 1 α u 1 1 Fig. 4. Sets efine in Proposition 10. Observe that, for uniform marginals, the bounary of A s coincies with B s. For the region unerlying such bounary we efine E i := {(u 1, u 2 ) [0, 1] 2 u 1 + u 2 < s, u i α}, i = 1, 2 10

11 as illustrate in Figure 4 (left). By the efinition of copula an C Ĉ2 α we have that µ C (E i A s ) = µ C (E i ) + µ C (A s ) = 1 α, i = 1, 2, µ C (A s ) = 1 α, which implies µ C (E i ) = 0, i = 1, 2 an µ C ([0, α) 2 ) = α. For the upper region we introuce the following partition: i(1 α) H ni := (α α (i + 1)(1 α), α + ] 2 n 2n+1 2 n (α + (2n i)(1 α) 1 α 2 n 2, α + (2n i)(1 α) ] n+1 2 n for n 0 an i = 0,...,2 n 1. See Figure 4 (right). In particular, consier H 00 = ( 1+α 2, 1]2 an let C i : = {(u 1, u 2 ) [0, 1] 2 u 1 + u 2 > 1 + α, u i 1 + α }, i = 1, 2, 2 D 1 : = {(u 1, u 2 ) [0, 1] 2 u 1 + u 2 = 1 + α, α u α }, 2 D 2 : = {(u 1, u 2 ) [0, 1] 2 u 1 + u 2 = 1 + α, 1 + α 2 < u 1 1}. Using the properties of a copula an consiering that E 1 an E 2 have zero µ C - measure, we have that µ C (H 00 ) + µ C (C 1 ) + µ C (D 1 ) = α 2 µ C (C 1 ) + µ C (D 1 ) = 1 + α 2 = 1 α, 2 α = 1 α 2 an hence µ C (H 00 ) = 0. Analogously, applying the same arguments to the upperright triangles of the squares [α, 1 + α 2 ] [ 1 + α 2, 1] an [ 1 + α 2, 1] [α, 1 + α ], 2 respectively, we obtain that µ C (H 10 ) = µ C (H 11 ) = 0. By iteration we have that µ C (H ni ) = 0 for all n 0, i = 0,...,2 n 1 an we trivially obtain µ C ( n=0 2 n 1 i=0 H ni ) = 0. Hence the only possibility for C is to assign probability mass (1 α) to the set D 1 D 2 = B s, which implies that C C 2 α. Remark 11 With respect to (10), for 1 s 2 an X C = (X 1, X 2 ) having stanar uniform marginals an copula C = C α, X 2 = g(x 1 ), where g : [0, 1] 11

12 [0, 1] is the linear function x, if x < s 1, g(x) = s x, otherwise. The above remark, together with Lemma 10, imply that the copula C of a uniform portfolio X C = (X 1, X 2 ) belongs to Ĉ2 α if an only if È[X 1 + X 2 = s X 1 + X 2 s] = 1. (13) 3.3 Multiimensional portfolios Though the boun (3) hols in arbitrary imensions, Theorem 7 fails to be vali if we take n > 2. Proposition 12 below shows in a simple way that, if we choose uniformly istribute marginals, it is not always possible to choose a copula C so as to obtain m ψ (s) = τ C,+ (F 1,...,F n )(s) =: α. Analogously to Ĉ2 α in the previous section, we efine Ĉ n α = {C Cn σ C,+ (F 1,..., F n )(s) = α}. Proposition 12 Let X C = (X 1,...,X n ) be a ranom vector having marginal fs uniformly istribute on [0, 1]. Take n > 2 an n 1 < s < n. Then Ĉn α =. Proof Let S n := n i=1 X i an note that, for uniform marginals, we have α = s n+1. If there exists k {1,...,n 2} such that È[S n k < s k] = 1 we have È[S n s] = 0 an the statement trivially hols. Suppose then È[S n k s k] > 0 for all k {1,...,n 2}. In this case we have È[S n s] = È[S n s, S n 1 s 1] + È[S n s, S n 1 < s 1] = È[S n s S n 1 s 1] È[S n 1 s 1], since X n is uniformly istribute on [0, 1]. Proceeing by iteration we obtain È[S n s] = È[S n s S n 1 s 1]... È[S 3 s n + 3 S 2 s n + 2] È[S 2 s n + 2 X 1 s n + 1](n s). (14) Assume now that Ĉn α, i.e. there exists XC = (X 1,...,X n ) with copula C Ĉ n α. It immeiately follows that È[S n s] = È[X X n s] = n s an hence all factors in (14), apart from the last one, must be equal to one. In particular, 12

13 inepenence 0.1 comonotonicity worst case lower boun s Fig. 5. Range for È[X 1 + X 2 + X 3 < s] for a stanar uniform portfolio. Together with the inepenence an comonotonic scenario, we plot the worst-case value m + (s) which iffers from the lower boun τ CL,+(s) given by (3). this yiels that È[X 1 + X 2 s n + 2 X 1 s n + 1] = 1, (15) È[S 3 s n + 3 S 2 s n + 2] = 1. (16) Accoring to (13), (15) implies that which, together with (16), leas to È[S 2 = s n + 2 S 2 s n + 2] = 1, 1 = È[S 3 s n + 3 S 2 s n + 2] = È[X 3 1 S 2 s n + 2] = È[X 3 1, S 2 s n + 2]). È[S 2 s n + 2] The latter equation is clearly a contraiction to the fact that X 3 is uniformly istribute on [0, 1]. Remark 13 The boun given in (3) fails to be sharp when n > 2 an C L = W. This erives from the fact that W is not a copula for n > 2, i.e. for more than two ranom variables it is impossible for each of them to be almost surely a nonincreasing function of each of the remaining ones. In Rüschenorf (1982), the worst-case VaR for uniform an binomial marginals is provie. Till now, this is the only known analytical result. In fact, the optimum epenence for uniform marginals oes not solve the general problem, showing that, contrary to the twoimensional case for n > 2, the epenence structure maximizing (6) may epen upon the choice of the marginals. In Embrechts an Puccetti (2004), however, an improve boun for the VaR is provie. Figure 5 illustrates the optimum values for uniform portfolios. 13

14 4 Evaluating risk through comonotonicity In the following, we show that the assumption of comonotonicity among the X i s may lea to a pruent evaluation of the risk associate with the aggregate position ψ(x). To this purpose, we first illustrate that such kin of epenence leas to the more angerous scenario with respect to both stop-loss an supermoular orer. Then, changing the optimization approach iscusse in the previous sections, we show that comonotonicity also arises as a suitable epenence assumption for our original VaR problem. 4.1 Stochastic orers an comonotonicity In this section we provie some motivation for the assumption of comonotonicity among risks base on stochastic orers. In this framework we illustrate an important application in actuarial mathematics. We first recall some concepts about stochastic orers. Definition 14 Let X an Y be two real ranom variables. We say that X is smaller than Y in stop-loss orer an we write X sl Y if for all non-ecreasing convex functions g : R R we have provie the expectations E[g(X)], E[g(Y )] are finite. E[g(X)] E[g(Y )], (17) The stop loss-orer compares one-imensional ranom variables. A multiimensional stochastic orer implying stop-loss orer is the so-calle supermoular orer, i.e. the orer base on a comparison of integrals of supermoular functions f : R n R satisfying f(x y) + f(x y) f(x) + f(y), for all x, y R n, where x y (x y) is the componentwise maximum (minimum) of x, y. Definition 15 Let X an Y be two n-imensional ranom vectors. We say that X is smaller than Y in supermoular orer an write X sm Y if for all supermoular functions g : R n R, E[g(X)] E[g(Y )], (18) provie the expectations E[g(X)], E[g(Y )] are finite. The next theorem recalls two important results about supermoular an stop-loss orers. In particular it states that comonotonicity represents the worst possible epenence scenario with respect to both such orers. 14

15 Theorem 16 Let X C = (X 1,...,X n ) be a n-imensional ranom vector having marginal istributions F 1,...,F n an copula C. Let ψ : R n R be a nonecreasing supermoular function. Then (a) X C sm X M, (b) ψ(x C ) sl ψ(x M ). Proof As note in Müller (1997), part (a) follows from Theorem 5 in Tchen (1980). Since ψ(x C ) sl ψ(x M ) hols if an only if (18) hols for all non-ecreasing convex functions g : R R for which expectations exists, to prove part (b) it is sufficient to show that for such a function g the function g ψ is supermoular. This follows from Lemma 2.2(b) in Bäuerle (1997). Remark 17 Note that Theorem 16 (b) applies to a large class of interesting functionals, incluing ψ(x) = n i=1 h i (x i ), where the h i s are non-ecreasing (see also Müller (1997)) an ψ(x) = h( n i=1 x i ) for h non-ecreasing an convex. For more examples of supermoular functions or some interesting methos of constructing them, Marshall an Olkin (1979), pp is the stanar reference. Here we want to point out that Theorem 16 (b) oes not apply to (6) because the inicator function of the set {ψ(x) s} is not supermoular. Consier again a portfolio of risks X C = (X 1,...,X n ). In insurance mathematics if ψ(x C ) is to be insure with a retention level, the net premium E[ψ(X C ) ] + is calle the stop-loss premium. A stop-loss premium is etermine once the retention an the multivariate f of X C are given. Hence we set π C,ψ (F 1,...,F n )() := E[ψ(X C ) ] +, P ψ () := sup C C n {π C,ψ (F 1,...,F n )()}. (19) Next proposition states that the univariate stop-loss orer base on the comparison of integrals of non-ecreasing convex functions is equivalent to the stochastic orer base on the comparison of stop-loss premiums. It also explains the name stop-loss for the stochastic orer sl. Proposition 18 Let X an Y be two real ranom variables. Then the following statements are equivalent: (a) X sl Y, (b) For all real E[X ] + E[Y ] +. (20) A simple proof can be foun in Müller an Stoyan (2002), Theorem

16 Accoring to Proposition 18, part (b) of Theorem 16 is equivalent to P ψ () = π M,ψ (F 1,...,F n )() (21) for all non-ecreasing supermoular functions ψ, real retention an arbitrary imension n. Hence π C,ψ (F 1,...,F n )() is maximize over C n when the fixe marginals of the portfolio have a comonotonic joint istribution, provie ψ is a non-ecreasing supermoular function. Observe that this is a much stronger result if confronte with Theorem 6. It is remarkable to note that this solution is not unique pointwise. In Figure 6 (right) we plot the ensity of a f on R 2 which, though iffering from comonotonicity (left), maximizes π C,+ (Φ, Φ)(0) over C n. However, M is the only epenence structures that attains (19) for all real retentions. Fig. 6. Densities of two-imensional istributions obtaine by comonotonic epenence (left) an by maximizing π C,+ (Φ,Φ)(0) over C 2 (right). 4.2 Changing the optimization approach in the VaR problem Theorem 5 provies a lower boun for the probability σ C,ψ (F 1,...,F n )(s), respectively, an upper boun for the VaR α (ψ(x 1,...,X n )). Hence an insurance company holing the risky position ψ(x C ) knows that VaR α (ψ(x C )) τ W,ψ (F 1,...,F n ) 1 (α), C C n. This equation represents the worst-case VaR scenario an can be very expensive for the insurer to hanle. Recalling the efinition of the family of copulas reaching the above level α, i.e. Ĉ n α = {C C n σ C,ψ (F 1,...,F n )(s) = α}, two quite natural conitions the insurer may ask to be satisfie are (a) Ĉn α, (b) Ĉn α oes not epen on s. 16

17 1 0.9 r = r = r = r = r = lower boun comonotonicity C r scenarios Fig. 7. Range for È[X 1 + X 2 < s] uner ifferent epenence scenarios for a stanar uniform portfolio. Inee (a) requires the boun to be sharp, i.e. that there exists a structure of epenence which etermines the worst-case VaR scenario. Conition (b) requires the worst scenario for the insurance company to be inepenent of the threshol s or, equivalently, inepenent of the parameter α that is chosen by the company or by the regulator to evaluate the aggregate risk. In the previous sections we showe that (a) is violate when n > 2 but it hols in the two-imensional case, while (b) is violate even when n = 2. Our aim here is to change the optimization approach so that the solutions satisfy conitions (a) an (b). In orer to o this, we efine the worst-case VaR scenario over a suitable range for the threshol s, rather than on a single value. Figure 7 explains how that can be one. In this graph we plot σ Cr,+(F 1, F 2 )(s) for ifferent values of r (0, 1) in case of two uniform marginals together with the best-possible lower boun m + (s) an the comonotonic curve σ M,+ (F 1, F 2 )(s). As a consequence of Theorem 7, every copula C r gives a lower boun that meets the curve m + (s) at the corresponing threshol an then becomes one. The intuition behin this plot is that the comonotonic copula, though never meeting the boun m + (s), is closer to it than any other copula on average. This iea can now be formalize by introucing a loss function Λ to measure the error committe by evaluating the risky position using a fixe copula C C n rather than the appropriate worst-possible structure of epenence. We then integrate the loss function over a suitable set B an we search for the infimum over the class of all n-copulas. Defining { } r ψ := inf Λ[σ C,ψ (F 1,...,F n )(s) m ψ (s)]s C C n B (22) for some non-ecreasing functional Λ : [0, 1] R + 0, we therefore introuce a measure for the istance between the worst an the comonotonic scenario in Figure 7, where Λ can be viewe as a weighting function. Let Ĉn (B) enote the set of 17

18 copulas leaing to (22). To focus our attention, we choose B = [, ). Theorem 19 Let Λ = I an B = [, ). Then, for every real threshol an non-ecreasing supermoular function ψ satisfying E[ψ(X M )] <, we have that M Ĉn ([, )). Proof Let m ψ(s)s <. Note that m ψ (s) epens only on the fixe marginals, so we obtain { } r ψ = inf [σ C,ψ (F 1,...,F n )(s) m ψ (s)] s C C n { } = sup [È[ψ(X C ) s] m ψ (s)] s C C n { } = m ψ (s) s sup È[ψ(X C ) s] s = m ψ (s) s sup C C n C C n { } È[ψ(X C ) > s] s, where the last step is obtaine since È[ψ(X C ) = s] can be positive at most for countably many values of s, so that the last two integrals containe in the brackets are the same. Finally, recalling that È[ψ(X C ) > s] s = E[ψ(X C ) ] +, it follows that r ψ = = = { } m ψ (s)s sup È[ψ(X C ) > s] s C C n m ψ (s)s sup C C n {E[ψ(X C ) ] + } m ψ (s)s P ψ (). Since E[ψ(X M )] is finite, (21) finally implies that M Ĉn ([, )). If the integral m ψ(s)s =, trivially Ĉn ([, )) = C n. Remark 20 Theorem 19 is vali for right-open intervals but it oes not hol in full generality. For instance, if we fix a trivial interval consisting of a single point, we go back to the original VaR problem. In such case, M oes not lea to the worst-possible scenario. Remark 21 For the above theorem, the only relevant portfolios (X 1,...,X n ) are those for which m ψ(s)s is finite. This technical conition is satisfie for all marginal istributions of interest. As escribe in the previous sections an illustrate in Figure 8 for the ensity of the N(0, 1)-N(0, 1)-portfolio of Example 9, 18

19 uner the optimizing copula, the mass m + (s) is concentrate on the upper triangle. On the remaining strips, where only one marginal takes high values, the mass is zero. Observe that, by the efinition of copula, m + (s) can not excee the sum of the marginal tails relative to these strips, i.e. 2Φ(s/2). In arbitrary imensions, with the same argument, we have that n m ψ (s) s F i (s/n) s < (23) i=1 for marginal fs F 1,...,F n satisfying /n F i(s)s <. The latter trivially hols if the marginal rvs are integrable. In particular, if X 1,..., X n enote losses an assume only positive values, conition (23) is satisfie provie that their expectations are finite. The finiteness of m ψ(s)s is equivalent to the finiteness of the expectation of the rv K in Proposition 23. We prove this result in full generality for non-ecreasing, continuous ψ an increasing marginals in Proposition 24 in the Appenix Fig. 8. Density an contour plot of the joint istribution of (X 1,X 2 ) maximizing È[X 1 + X 2 s] in Example 9. The main issue unerlying Theorem 19 is that, even if the comonotonic epenence structure oes not lea to the worst-case scenario for the original VaR problem, if an insurance company wants to boun well σ C,ψ (F 1,...,F n )(s) for all threshols in [, ) in the sense efine by (22), comonotonicity provies a pruent evaluation for the aggregate risk. The next theorem elivers an insight into an extension of Theorem 19 for general loss functions Λ, i.e. every increasing, convex function Λ : [0, 1] R + 0 satisfying Λ(0) = 0. For a copula C, let e C,ψ (s) := σ C,ψ (F 1,...,F n )(s) m ψ (s) enote the error function in (22). 19

20 Theorem 22 Consier X C = (X 1,..., X n ) with marginal istributionsf 1,...,F n an let ψ be as in Theorem 19 with 0 m ψ (s)s <. Then there exists C [ C, C ] such that for all loss functions Λ: Λ[σ C,ψ (F 1,...,F n )(s) m ψ (s)]s for every C, where Λ[σ M,ψ (F 1,...,F n )(s) m ψ (s)]s C := sup{ R e C,ψ (s) e M,ψ () s an e M,ψ (s) > e C,ψ (s) on some interval ( C, )}, C := inf{ R e C,ψ (s) e M,ψ (s) s }. (24) Proof Theorem 19 yiels e C,ψ(s)s e M,ψ(s)s for all R. The latter integrals are finite since 0 m ψ (s)s <. Denote with B(, ) an m the Borel sets on (, ) an the Lebesgue measure, respectively. Applying Chong (1974), Theorem 1.6, Theorem 2.1 an Corollary 1.2 to e C,ψ an e M,ψ, with Φ = Λ an (X, Λ, µ) = (X, Λ, µ ) = ((, ), B(, ), m), (24) is equivalent to [e C,ψ (s) u] + s [e M,ψ (s) u] + s (25) for all u R. By efinition, C C an e C,ψ e M,ψ on [ C, ) implying C C. Assume now < C < C an let ( C, C]. Choosing u = e M,ψ ( C ) in (25) we have that [e C,ψ (s) u] + C s = [e C,ψ (s) e M,ψ ( C )] + s + = C C [e C,ψ (s) e M,ψ ( C )] + s< [e M,ψ (s) u] + s, C [e C,ψ (s) e M,ψ ( C )] + s [e M,ψ (s) e M,ψ ( C )] + s which conclues the proof. With respect to a copula C an any loss function Λ, comonotonicity is hence a suitable epenence scenario on [ C, ) an it provies a pruent VaR on the corresponing α-set. Note that, for a copula C, the set in the efinition of C may be empty an hence C arbitrarily small. Contrarily, since ψ(x C ) sl ψ(x M ), we have that C < if the two fs cross finitely many times. Unfortunately, in general, C an C may become arbitrary large. For instance, if the function ψ is unboune, Rüschenorf (1981), Theorem 5 yiels the existence of a copula Ĉ epening on ψ, s an the marginals such that m ψ (s) = σĉ,ψ (F 1,...,F n )(s) σ M,ψ (F 1,...,F n )(s), 20

21 where, in general, Ĉ M implying sup C C n C =. Analogously there exist epenence structures leaing to sup C C n C =. We therefore conclue that an extension of Theorem 19 to loss functionals can only exist for suitable subclasses of C n. 5 The Presence of Information Throughout this paper we mainly assume a no-information scenario an put the lower Fréchet boun W as lower boun for the unknown copula C in Theorem 5. From a mathematical point of view this is an unsatisfactory initial assumption, since, for n > 2, W is not a copula anymore an the boun provie in Theorem 6 fails to be sharp. However, we want to warn the reaer from choosing some a priori assumption such as C Π an we emphasize that similar restrictions may lea to a critical unervaluation of the portfolio risk. The assumption C Π, for instance, is justifie by the iea that the worst-var scenarios is implie by the so-calle positive lower orthant epenent risks with È[X 1 x 1,...,X n x n ] È[X 1 x 1 ]... È[X n x n ]. Unfortunately, this is not true an restricting the optimization to the class {C Π} substantially changes the initial problem, since it oes not allow to focus on riskier portfolios, as long as m + > m Π,+. This is a consequence of the fact that the componentwise orering in the class C n is not complete an, putting a lower boun on a copula, exclues all copulas not comparable to such boun. To point out this aspect, we observe that countermonotonicity plays a central role in the efinition of the family Ĉn α arising from (6), whereas every copula shuffle with it is not comparable with the inepenence scenario. 6 Conclusions In this paper we focus on the copulas leaing to the worst-possible VaR for a function of epenent risks an we emphasize that comonotonicity oes not lie in this family. Such worst-case scenarios epen upon the level α where the VaR is evaluate an therefore are not reasonable from a practical point of view. Moreover, these solutions are known only for two-imensional portfolios or in presence of partial information. The investigation of optimal bouns in arbitrary imensions with no prior information remains open. Therefore, we provie an alternative approach supporting the assumption of comonotonicity in a pruent evaluation of the quantiles of the aggregate position. 21

22 Appenix The operator τ C,ψ In this section we exten a claim of Denuit et al. (1999), p. 37 by showing that the operator τ C,ψ (F 1,...,F n )(s) = sup C(F 1 (x 1 ),..., F n 1 (x n 1 ), Fn (ψ x n (s))), x 1,...,x n 1 R with ψ x n (s) = sup{x n R ψ(x n, x n ) < s}, for fixe x n R n 1, is actually the left-continuous version of a f. Proposition 23 For ψ : R n R non-ecreasing, there exists a ranom variable K such that τ C,ψ (F 1,...,F n )(s) = È[K < s]. Proof Since ψ is non-ecreasing in each component, τ C,ψ (s) := τ C,ψ (F 1,..., F n )(s) is a non-ecreasing function. Hence, we have to show that it is also left-continuous an that its right an left limits converge to one an zero, respectively. To prove that lim s τ C,ψ (s) = 1, we fix ε > 0 an efine u ε = (u ε 1,...,uε n ) as a vector satisfying F i (u ε i ) 1 ε n, i = 1,...,n. The existence of such a vector is straightforwar, since F 1,...,F n are non-efective fs. By efinition, the function ψ û is non-ecreasing an its right limit is either ε n finite or infinite. Suppose it is finite. For every real s, it follows that which implies sup{x n R ψ(u ε n, x n ) < s} lim ψ s û (s) =: R <, ε n ψ(u ε n, x n) s for all x n R. Therefore ψ(u ε n, R) =, which contraicts ψ having R as its range. Hence R = an it is always possible to select a real s ε, epening only on ε, such that ψ û ε n (s ε) > u ε n implying For τ C,ψ (s ε ) we also obtain that τ C,ψ (s ε ) = F n (ψ û ε n (s ε)) F n (u ε n ) 1 ε n. sup C(F 1 (x 1 ),...,F n 1 (x n 1 ), Fn (ψ x n (s ε ))) x 1,...,x n 1 R C(F 1 (u ε 1),...,F n 1 (u ε n 1), Fn (ψ (s û ε))) ε n W(F 1 (u ε 1),..., F n 1 (u ε n 1), F n (ψ û ε n (s ε))) (1 ε )n n + 1 = 1 ε, n 22

23 an, since τ C,ψ is non-ecreasing, τ C,ψ (s) τ C,ψ (s ε ) 1 ε for every s s ε. Hence the right limit converges to one. Similarly, for the left limit, we fix ε > 0 an choose u ε satisfying F i (u ε i) < ε, i = 1,..., n for which L := lim s ψ û ε n (s) =. It is always possible to select a real s ε, epening only on ε, such that ψ û ε n (s ε) < u ε n an F n (ψûε n (s ε)) < F n (u ε n ) < ε. Let now A u ε := {x n R n 1 x i u ε i for some i {1,...,n 1}} an A u ε := {x n R n 1 x i > u ε i for all i = 1,..., n 1}. Then sup C(F 1 (x 1 ),...,F n 1 (x n 1 ), Fn (ψ x n (s ε ))) x n A u ε sup M(F 1 (x 1 ),..., F n 1 (x n 1 ), Fn (ψ x n (s ε ))) x n A u ε F i (u ε i ) < ε. (26) If x A u ε, then ψ x n (s ε ) ψ û ε n (s ε) an hence sup C(F 1 (x 1 ),...,F n 1 (x n 1 ), Fn (ψ x n (s ε ))) x n A u ε sup C(F 1 (x 1 ),...,F n 1 (x n 1 ), Fn (s ε))) (ψûε n x n A u ε Fn (ψ (s û ε)) < ε. ε n (27) From (26) an (27) we have that τ C,ψ (s ε ) < ε. Since τ C,ψ is non-ecreasing, τ C,ψ (s) < τ C,ψ (s ε ) < ε for all s < s ε an the left limit goes to zero. It remains to show that τ C,ψ is left-continuous. For non-ecreasing functions f : R R left-continuity is equivalent to lower-semicontinuity. By Ruin (1974), p. 39, the supremum of any collection of lower-semicontinuous function is lowersemicontinuous. It is sufficient then to show that C(F 1 (x 1 ),...,F n 1 (x n 1 ), F n (ψ x n (s))) is left-continuous in s for every x n R n 1. By uniform continuity of C, leftcontinuity an non-ecreasingness of F n, an non-ecreasingness of ψ x n (s), the problem is reuce to show that ψ x n (s) is left-continuous. By efinition, it is nonecreasing an hence, for every real s, there exists l(s) := lim x s ψ x n (x). Assume now that ψ x n (s) is not left-continuous. Then there exists s 1 with l(s 1 ) < ψ x n (s 1 ). Let l(s 1 ) be finite (otherwise there is nothing to prove). Then, for arbitrary positive ε, we have sup{x n R ψ(x n, x n ) < s 1 ε} l(s 1 ) < 23

24 an hence, whenever x n l(s 1 ), it follows that ψ(x n, x n ) s 1 ε. Since ε is arbitrary, it follows that ψ(x n, l(s 1 )) s 1, contraicting the fact that ψ(x n, x n ) < s for every x n < ψ x n (s 1 ), which conclues the proof. The following proposition states that, uner suitable assumptions, the rv K has finite expectation. Proposition 24 Let X C have increasing, continuous marginals F 1,...,F n, ψ : R n R be non-ecreasing continuous an increasing in the last argument. Then, for K as in Proposition 23, E[K] const + m ψ(s)s < for all R if an only if E[ψ(X M )] <. Proof A ranom variable Y with f F has finite expectation if an only if 0 F(x)x < an 0 F(x)x <, which is equivalent to E[Y ] + < for all R. We therefore have to show that m ψ (s) s < E[ψ(X M ) ] + < (28) for all R. By the efinition of m ψ (s), we have that for R, E[ψ(X M ) ] + = an hence immeiately follows. P[ψ(X M ) s] s m ψ (s) s Assume now that the rhs of (28) hols an let U be uniformly istribute on [0, 1]. By Dhaene et al. (2001), Theorem 2 we have that E[φ(U) ] + = E[ψ(X M ) ] + <, where φ : [0, 1] R, φ(u) := ψ(f 1 1 (u),...,f 1 n (u)). Uner the assumptions of the theorem, φ is continuous an φ(φ 1 (s)) = s. We can write m ψ (s) 1 τ W,ψ (F 1,...,F n )(s) = 1 sup x n R n 1 {(F 1 (x 1 ) + + F n 1 (x n 1 ) + F n (ψ x n (s)) n + 1) + } inf x n R n 1{F 1(x 1 ) + + F n 1 (x n 1 ) + È[X n ψ x n (s)]}. (29) 24

25 Choosing x n = (F1 1 (φ 1 (s)),...,fn 1(φ 1 (s))) in (29) an since ψ is increasing in the last argument, (s) = F 1 ψ x n n (φ 1 (s)). Integrating, we finally obtain: m ψ (s)s = = = n n i=1 n i=1 n i=1 È[X i Fi 1 (φ 1 (s))]s È[F 1 i (U) F 1 (φ 1 (s))]s i n È[U φ 1 (s)]s È[φ(U) φ(φ 1 (s))]s i=1 È[φ(U) s]s = ne[φ(x M ) ] + <. References Artzner, P., Delbaen, F., Eber, J. M. an Heath, D., Coherent measures of risk. Mathematical Finance 9, Bäuerle, N., Monotonicity results for MR/GI/1 queues. Journal of Applie Probability 34, Billingsley, P., Probability an Measure. Thir E. John Wiley & Sons, Chichester. Chong, K.-M., Some extensions of a theorem of Hary, Littlewoo an Pólya an their applications. Canaian Journal of Mathematics 26 (6), Denuit, M., Genest, C. an Marceau, É., Stochastic bouns on sums of epenent risks. Insurance: Mathematics an Economics 25, Dhaene, J., Denuit, M., Goovaerts, M. J., Kaas, R. an Vyncke, D., The concept of comonotonicity in actuarial science an finance: theory. Insurance: Mathematics an Economics 31, Embrechts, P., Höing, A. an Juri, A., Using copulae to boun the Valueat-Risk for functions of epenent risks. Finance & Stochastics 7, Embrechts, P. an Puccetti, G., Bouns on Value-at-Risk. Preprint. Frank, M.J., Nelsen, R.B. an Schweizer, B., Best-possible bouns on the istribution of a sum a problem of Kolmogorov. Probability Theory an Relate Fiels 74, Marshall, A. W. an Olkin, I., Inequalities: Theory of Majorization an its Applications. Acaemic press, New York. Mikusiński, P., Sherwoo, F. an Taylor, M. D., Probabilistic interpretations of copulas an their convex sums. Avances in Probability Distributions with Given Marginals (Giorgio Dall Aglio, S. Kotz an G. Salinetti Es.), Kluwer Acaemic Publisher, Amsteram, pp Müller, A., Stop-loss orer for portfolios of epenent risks. Insurance: 25

26 Mathematics an Economics 21, Müller, A. an Stoyan, D., Comparison Methos for Stochastic Moels an Risks. John Wiley an Sons, Chichester. Nelsen, R.B., An Introuction to Copulas. Springer, Berlin. Ruin, W., Real an Complex Analysis. McGraw-Hill, New York. Rüschenorf, L., Sharpness of Fréchet Bouns. Zeitschrift für Wahrscheinlichkeitstheorie un verwante Gebiete 57, Rüschenorf, L., Ranom variables with maximum sums. Avances in Applie Probability 14, Schweizer, B. an Sklar, A., Operations on istribution functions not erivable from operations on ranom variables. Stuia Mathematica 52, Sklar, A., Ranom variables, joint istribution functions, an copulas. Kybernetica 9, Tchen, A. H., Inequalities for istributions with given marginals. The Annals of Probability 8,

A simple tranformation of copulas

A simple tranformation of copulas A simple tranformation of copulas V. Durrleman, A. Nikeghbali & T. Roncalli Groupe e Recherche Opérationnelle Créit Lyonnais France July 31, 2000 Abstract We stuy how copulas properties are moifie after

More information

Modeling of Dependence Structures in Risk Management and Solvency

Modeling of Dependence Structures in Risk Management and Solvency Moeling of Depenence Structures in Risk Management an Solvency University of California, Santa Barbara 0. August 007 Doreen Straßburger Structure. Risk Measurement uner Solvency II. Copulas 3. Depenent

More information

6 General properties of an autonomous system of two first order ODE

6 General properties of an autonomous system of two first order ODE 6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x

More information

Explicit 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 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 information

1. Aufgabenblatt zur Vorlesung Probability Theory

1. Aufgabenblatt zur Vorlesung Probability Theory 24.10.17 1. Aufgabenblatt zur Vorlesung By (Ω, A, P ) we always enote the unerlying probability space, unless state otherwise. 1. Let r > 0, an efine f(x) = 1 [0, [ (x) exp( r x), x R. a) Show that p f

More information

Topic 7: Convergence of Random Variables

Topic 7: Convergence of Random Variables Topic 7: Convergence of Ranom Variables Course 003, 2016 Page 0 The Inference Problem So far, our starting point has been a given probability space (S, F, P). We now look at how to generate information

More information

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions

Computing Exact Confidence Coefficients of Simultaneous Confidence Intervals for Multinomial Proportions and their Functions Working Paper 2013:5 Department of Statistics Computing Exact Confience Coefficients of Simultaneous Confience Intervals for Multinomial Proportions an their Functions Shaobo Jin Working Paper 2013:5

More information

Logarithmic spurious regressions

Logarithmic spurious regressions Logarithmic spurious regressions Robert M. e Jong Michigan State University February 5, 22 Abstract Spurious regressions, i.e. regressions in which an integrate process is regresse on another integrate

More information

Chapter 6: Energy-Momentum Tensors

Chapter 6: Energy-Momentum Tensors 49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.

More information

The total derivative. Chapter Lagrangian and Eulerian approaches

The total derivative. Chapter Lagrangian and Eulerian approaches Chapter 5 The total erivative 51 Lagrangian an Eulerian approaches The representation of a flui through scalar or vector fiels means that each physical quantity uner consieration is escribe as a function

More information

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012

Lecture Introduction. 2 Examples of Measure Concentration. 3 The Johnson-Lindenstrauss Lemma. CS-621 Theory Gems November 28, 2012 CS-6 Theory Gems November 8, 0 Lecture Lecturer: Alesaner Mąry Scribes: Alhussein Fawzi, Dorina Thanou Introuction Toay, we will briefly iscuss an important technique in probability theory measure concentration

More information

Lower Bounds for the Smoothed Number of Pareto optimal Solutions

Lower Bounds for the Smoothed Number of Pareto optimal Solutions Lower Bouns for the Smoothe Number of Pareto optimal Solutions Tobias Brunsch an Heiko Röglin Department of Computer Science, University of Bonn, Germany brunsch@cs.uni-bonn.e, heiko@roeglin.org Abstract.

More information

Equilibrium in Queues Under Unknown Service Times and Service Value

Equilibrium in Queues Under Unknown Service Times and Service Value University of Pennsylvania ScholarlyCommons Finance Papers Wharton Faculty Research 1-2014 Equilibrium in Queues Uner Unknown Service Times an Service Value Laurens Debo Senthil K. Veeraraghavan University

More information

Least-Squares Regression on Sparse Spaces

Least-Squares Regression on Sparse Spaces Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction

More information

Agmon Kolmogorov Inequalities on l 2 (Z d )

Agmon Kolmogorov Inequalities on l 2 (Z d ) Journal of Mathematics Research; Vol. 6, No. ; 04 ISSN 96-9795 E-ISSN 96-9809 Publishe by Canaian Center of Science an Eucation Agmon Kolmogorov Inequalities on l (Z ) Arman Sahovic Mathematics Department,

More information

Calculus of Variations

Calculus of Variations Calculus of Variations Lagrangian formalism is the main tool of theoretical classical mechanics. Calculus of Variations is a part of Mathematics which Lagrangian formalism is base on. In this section,

More information

arxiv:math/ v1 [math.pr] 19 Apr 2001

arxiv:math/ v1 [math.pr] 19 Apr 2001 Conitional Expectation as Quantile Derivative arxiv:math/00490v math.pr 9 Apr 200 Dirk Tasche November 3, 2000 Abstract For a linear combination u j X j of ranom variables, we are intereste in the partial

More information

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson

JUST THE MATHS UNIT NUMBER DIFFERENTIATION 2 (Rates of change) A.J.Hobson JUST THE MATHS UNIT NUMBER 10.2 DIFFERENTIATION 2 (Rates of change) by A.J.Hobson 10.2.1 Introuction 10.2.2 Average rates of change 10.2.3 Instantaneous rates of change 10.2.4 Derivatives 10.2.5 Exercises

More information

Math 342 Partial Differential Equations «Viktor Grigoryan

Math 342 Partial Differential Equations «Viktor Grigoryan Math 342 Partial Differential Equations «Viktor Grigoryan 6 Wave equation: solution In this lecture we will solve the wave equation on the entire real line x R. This correspons to a string of infinite

More information

MODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS ABSTRACT KEYWORDS

MODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS ABSTRACT KEYWORDS MODELLING DEPENDENCE IN INSURANCE CLAIMS PROCESSES WITH LÉVY COPULAS BY BENJAMIN AVANZI, LUKE C. CASSAR AND BERNARD WONG ABSTRACT In this paper we investigate the potential of Lévy copulas as a tool for

More information

Tractability results for weighted Banach spaces of smooth functions

Tractability results for weighted Banach spaces of smooth functions Tractability results for weighte Banach spaces of smooth functions Markus Weimar Mathematisches Institut, Universität Jena Ernst-Abbe-Platz 2, 07740 Jena, Germany email: markus.weimar@uni-jena.e March

More information

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d

A new proof of the sharpness of the phase transition for Bernoulli percolation on Z d A new proof of the sharpness of the phase transition for Bernoulli percolation on Z Hugo Duminil-Copin an Vincent Tassion October 8, 205 Abstract We provie a new proof of the sharpness of the phase transition

More information

Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS"

Mark J. Machina CARDINAL PROPERTIES OF LOCAL UTILITY FUNCTIONS Mark J. Machina CARDINAL PROPERTIES OF "LOCAL UTILITY FUNCTIONS" This paper outlines the carinal properties of "local utility functions" of the type use by Allen [1985], Chew [1983], Chew an MacCrimmon

More information

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments

Time-of-Arrival Estimation in Non-Line-Of-Sight Environments 2 Conference on Information Sciences an Systems, The Johns Hopkins University, March 2, 2 Time-of-Arrival Estimation in Non-Line-Of-Sight Environments Sinan Gezici, Hisashi Kobayashi an H. Vincent Poor

More information

Modelling and simulation of dependence structures in nonlife insurance with Bernstein copulas

Modelling and simulation of dependence structures in nonlife insurance with Bernstein copulas Moelling an simulation of epenence structures in nonlife insurance with Bernstein copulas Prof. Dr. Dietmar Pfeifer Dept. of Mathematics, University of Olenburg an AON Benfiel, Hamburg Dr. Doreen Straßburger

More information

PDE Notes, Lecture #11

PDE Notes, Lecture #11 PDE Notes, Lecture # from Professor Jalal Shatah s Lectures Febuary 9th, 2009 Sobolev Spaces Recall that for u L loc we can efine the weak erivative Du by Du, φ := udφ φ C0 If v L loc such that Du, φ =

More information

arxiv: v4 [math.pr] 27 Jul 2016

arxiv: v4 [math.pr] 27 Jul 2016 The Asymptotic Distribution of the Determinant of a Ranom Correlation Matrix arxiv:309768v4 mathpr] 7 Jul 06 AM Hanea a, & GF Nane b a Centre of xcellence for Biosecurity Risk Analysis, University of Melbourne,

More information

Risk Aggregation with Dependence Uncertainty

Risk Aggregation with Dependence Uncertainty Introduction Extreme Scenarios Asymptotic Behavior Challenges Risk Aggregation with Dependence Uncertainty Department of Statistics and Actuarial Science University of Waterloo, Canada Seminar at ETH Zurich

More information

The Exact Form and General Integrating Factors

The Exact Form and General Integrating Factors 7 The Exact Form an General Integrating Factors In the previous chapters, we ve seen how separable an linear ifferential equations can be solve using methos for converting them to forms that can be easily

More information

Introduction to the Vlasov-Poisson system

Introduction to the Vlasov-Poisson system Introuction to the Vlasov-Poisson system Simone Calogero 1 The Vlasov equation Consier a particle with mass m > 0. Let x(t) R 3 enote the position of the particle at time t R an v(t) = ẋ(t) = x(t)/t its

More information

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x)

The derivative of a function f(x) is another function, defined in terms of a limiting expression: f(x + δx) f(x) Y. D. Chong (2016) MH2801: Complex Methos for the Sciences 1. Derivatives The erivative of a function f(x) is another function, efine in terms of a limiting expression: f (x) f (x) lim x δx 0 f(x + δx)

More information

Acute sets in Euclidean spaces

Acute sets in Euclidean spaces Acute sets in Eucliean spaces Viktor Harangi April, 011 Abstract A finite set H in R is calle an acute set if any angle etermine by three points of H is acute. We examine the maximal carinality α() of

More information

Lower bounds on Locality Sensitive Hashing

Lower bounds on Locality Sensitive Hashing Lower bouns on Locality Sensitive Hashing Rajeev Motwani Assaf Naor Rina Panigrahy Abstract Given a metric space (X, X ), c 1, r > 0, an p, q [0, 1], a istribution over mappings H : X N is calle a (r,

More information

On the number of isolated eigenvalues of a pair of particles in a quantum wire

On the number of isolated eigenvalues of a pair of particles in a quantum wire On the number of isolate eigenvalues of a pair of particles in a quantum wire arxiv:1812.11804v1 [math-ph] 31 Dec 2018 Joachim Kerner 1 Department of Mathematics an Computer Science FernUniversität in

More information

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs

Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Perfect Matchings in Õ(n1.5 ) Time in Regular Bipartite Graphs Ashish Goel Michael Kapralov Sanjeev Khanna Abstract We consier the well-stuie problem of fining a perfect matching in -regular bipartite

More information

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013

Survey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013 Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing

More information

Quantile function expansion using regularly varying functions

Quantile function expansion using regularly varying functions Quantile function expansion using regularly varying functions arxiv:705.09494v [math.st] 9 Aug 07 Thomas Fung a, an Eugene Seneta b a Department of Statistics, Macquarie University, NSW 09, Australia b

More information

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy,

NOTES ON EULER-BOOLE SUMMATION (1) f (l 1) (n) f (l 1) (m) + ( 1)k 1 k! B k (y) f (k) (y) dy, NOTES ON EULER-BOOLE SUMMATION JONATHAN M BORWEIN, NEIL J CALKIN, AND DANTE MANNA Abstract We stuy a connection between Euler-MacLaurin Summation an Boole Summation suggeste in an AMM note from 196, which

More information

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction

FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS. 1. Introduction FLUCTUATIONS IN THE NUMBER OF POINTS ON SMOOTH PLANE CURVES OVER FINITE FIELDS ALINA BUCUR, CHANTAL DAVID, BROOKE FEIGON, MATILDE LALÍN 1 Introuction In this note, we stuy the fluctuations in the number

More information

Qubit channels that achieve capacity with two states

Qubit channels that achieve capacity with two states Qubit channels that achieve capacity with two states Dominic W. Berry Department of Physics, The University of Queenslan, Brisbane, Queenslan 4072, Australia Receive 22 December 2004; publishe 22 March

More information

Lecture 5. Symmetric Shearer s Lemma

Lecture 5. Symmetric Shearer s Lemma Stanfor University Spring 208 Math 233: Non-constructive methos in combinatorics Instructor: Jan Vonrák Lecture ate: January 23, 208 Original scribe: Erik Bates Lecture 5 Symmetric Shearer s Lemma Here

More information

Generalized Tractability for Multivariate Problems

Generalized Tractability for Multivariate Problems Generalize Tractability for Multivariate Problems Part II: Linear Tensor Prouct Problems, Linear Information, an Unrestricte Tractability Michael Gnewuch Department of Computer Science, University of Kiel,

More information

Convergence of Random Walks

Convergence of Random Walks Chapter 16 Convergence of Ranom Walks This lecture examines the convergence of ranom walks to the Wiener process. This is very important both physically an statistically, an illustrates the utility of

More information

Transformations of Copulas and Measures of Concordance

Transformations of Copulas and Measures of Concordance Transformations of Copulas an Measures of Concorance D I S S E R T A T I O N zur Erlangung es akaemischen Graes Doctor rerum naturalium (Dr. rer. nat.) vorgelegt er Fakultät Mathematik un Naturwissenschaften

More information

REAL ANALYSIS I HOMEWORK 5

REAL ANALYSIS I HOMEWORK 5 REAL ANALYSIS I HOMEWORK 5 CİHAN BAHRAN The questions are from Stein an Shakarchi s text, Chapter 3. 1. Suppose ϕ is an integrable function on R with R ϕ(x)x = 1. Let K δ(x) = δ ϕ(x/δ), δ > 0. (a) Prove

More information

Table of Common Derivatives By David Abraham

Table of Common Derivatives By David Abraham Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec

More information

Schrödinger s equation.

Schrödinger s equation. Physics 342 Lecture 5 Schröinger s Equation Lecture 5 Physics 342 Quantum Mechanics I Wenesay, February 3r, 2010 Toay we iscuss Schröinger s equation an show that it supports the basic interpretation of

More information

Extreme Values by Resnick

Extreme Values by Resnick 1 Extreme Values by Resnick 1 Preliminaries 1.1 Uniform Convergence We will evelop the iea of something calle continuous convergence which will be useful to us later on. Denition 1. Let X an Y be metric

More information

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs

Lectures - Week 10 Introduction to Ordinary Differential Equations (ODES) First Order Linear ODEs Lectures - Week 10 Introuction to Orinary Differential Equations (ODES) First Orer Linear ODEs When stuying ODEs we are consiering functions of one inepenent variable, e.g., f(x), where x is the inepenent

More information

Asymptotic Bounds for the Distribution of the Sum of Dependent Random Variables

Asymptotic Bounds for the Distribution of the Sum of Dependent Random Variables Asymptotic Bounds for the Distribution of the Sum of Dependent Random Variables Ruodu Wang November 26, 2013 Abstract Suppose X 1,, X n are random variables with the same known marginal distribution F

More information

LECTURE NOTES ON DVORETZKY S THEOREM

LECTURE NOTES ON DVORETZKY S THEOREM LECTURE NOTES ON DVORETZKY S THEOREM STEVEN HEILMAN Abstract. We present the first half of the paper [S]. In particular, the results below, unless otherwise state, shoul be attribute to G. Schechtman.

More information

How to Minimize Maximum Regret in Repeated Decision-Making

How to Minimize Maximum Regret in Repeated Decision-Making How to Minimize Maximum Regret in Repeate Decision-Making Karl H. Schlag July 3 2003 Economics Department, European University Institute, Via ella Piazzuola 43, 033 Florence, Italy, Tel: 0039-0-4689, email:

More information

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS

ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS ALGEBRAIC AND ANALYTIC PROPERTIES OF ARITHMETIC FUNCTIONS MARK SCHACHNER Abstract. When consiere as an algebraic space, the set of arithmetic functions equippe with the operations of pointwise aition an

More information

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling

Balancing Expected and Worst-Case Utility in Contracting Models with Asymmetric Information and Pooling Balancing Expecte an Worst-Case Utility in Contracting Moels with Asymmetric Information an Pooling R.B.O. erkkamp & W. van en Heuvel & A.P.M. Wagelmans Econometric Institute Report EI2018-01 9th January

More information

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device

Closed and Open Loop Optimal Control of Buffer and Energy of a Wireless Device Close an Open Loop Optimal Control of Buffer an Energy of a Wireless Device V. S. Borkar School of Technology an Computer Science TIFR, umbai, Inia. borkar@tifr.res.in A. A. Kherani B. J. Prabhu INRIA

More information

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1

d dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1 Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of

More information

Lecture 6: Calculus. In Song Kim. September 7, 2011

Lecture 6: Calculus. In Song Kim. September 7, 2011 Lecture 6: Calculus In Song Kim September 7, 20 Introuction to Differential Calculus In our previous lecture we came up with several ways to analyze functions. We saw previously that the slope of a linear

More information

Robustness and Perturbations of Minimal Bases

Robustness and Perturbations of Minimal Bases Robustness an Perturbations of Minimal Bases Paul Van Dooren an Froilán M Dopico December 9, 2016 Abstract Polynomial minimal bases of rational vector subspaces are a classical concept that plays an important

More information

arxiv: v1 [math.ap] 6 Jul 2017

arxiv: v1 [math.ap] 6 Jul 2017 Local an global time ecay for parabolic equations with super linear first orer terms arxiv:177.1761v1 [math.ap] 6 Jul 17 Martina Magliocca an Alessio Porretta ABSTRACT. We stuy a class of parabolic equations

More information

Characterization of Upper Comonotonicity via Tail Convex Order

Characterization of Upper Comonotonicity via Tail Convex Order Characterization of Upper Comonotonicity via Tail Convex Order Hee Seok Nam a,, Qihe Tang a, Fan Yang b a Department of Statistics and Actuarial Science, University of Iowa, 241 Schaeffer Hall, Iowa City,

More information

Lecture 1: Review of Basic Asymptotic Theory

Lecture 1: Review of Basic Asymptotic Theory Lecture 1: Instructor: Department of Economics Stanfor University Prepare by Wenbo Zhou, Renmin University Basic Probability Theory Takeshi Amemiya, Avance Econometrics, 1985, Harvar University Press.

More information

Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane

Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Arch. Rational Mech. Anal. Digital Object Ientifier (DOI) https://oi.org/1.17/s25-18-1336-z Suppressing Chemotactic Blow-Up Through a Fast Splitting Scenario on the Plane Siming He & Eitan Tamor Communicate

More information

LOCAL AND GLOBAL MINIMALITY RESULTS FOR A NONLOCAL ISOPERIMETRIC PROBLEM ON R N

LOCAL AND GLOBAL MINIMALITY RESULTS FOR A NONLOCAL ISOPERIMETRIC PROBLEM ON R N LOCAL AND GLOBAL MINIMALITY RSULTS FOR A NONLOCAL ISOPRIMTRIC PROBLM ON R N M. BONACINI AND R. CRISTOFRI Abstract. We consier a nonlocal isoperimetric problem efine in the whole space R N, whose nonlocal

More information

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection

Online Appendix for Trade Policy under Monopolistic Competition with Firm Selection Online Appenix for Trae Policy uner Monopolistic Competition with Firm Selection Kyle Bagwell Stanfor University an NBER Seung Hoon Lee Georgia Institute of Technology September 6, 2018 In this Online

More information

Interconnected Systems of Fliess Operators

Interconnected Systems of Fliess Operators Interconnecte Systems of Fliess Operators W. Steven Gray Yaqin Li Department of Electrical an Computer Engineering Ol Dominion University Norfolk, Virginia 23529 USA Abstract Given two analytic nonlinear

More information

A Measure of Monotonicity of Two Random Variables

A Measure of Monotonicity of Two Random Variables Journal of Mathematics and Statistics 8 (): -8, 0 ISSN 549-3644 0 Science Publications A Measure of Monotonicity of Two Random Variables Farida Kachapova and Ilias Kachapov School of Computing and Mathematical

More information

Bivariate distributions characterized by one family of conditionals and conditional percentile or mode functions

Bivariate distributions characterized by one family of conditionals and conditional percentile or mode functions Journal of Multivariate Analysis 99 2008) 1383 1392 www.elsevier.com/locate/jmva Bivariate istributions characterize by one family of conitionals an conitional ercentile or moe functions Barry C. Arnol

More information

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations

Diophantine Approximations: Examining the Farey Process and its Method on Producing Best Approximations Diophantine Approximations: Examining the Farey Process an its Metho on Proucing Best Approximations Kelly Bowen Introuction When a person hears the phrase irrational number, one oes not think of anything

More information

Discrete Mathematics

Discrete Mathematics Discrete Mathematics 309 (009) 86 869 Contents lists available at ScienceDirect Discrete Mathematics journal homepage: wwwelseviercom/locate/isc Profile vectors in the lattice of subspaces Dániel Gerbner

More information

On Kusuoka Representation of Law Invariant Risk Measures

On Kusuoka Representation of Law Invariant Risk Measures MATHEMATICS OF OPERATIONS RESEARCH Vol. 38, No. 1, February 213, pp. 142 152 ISSN 364-765X (print) ISSN 1526-5471 (online) http://dx.doi.org/1.1287/moor.112.563 213 INFORMS On Kusuoka Representation of

More information

A simple model for the small-strain behaviour of soils

A simple model for the small-strain behaviour of soils A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:

More information

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations

Characterizing Real-Valued Multivariate Complex Polynomials and Their Symmetric Tensor Representations Characterizing Real-Value Multivariate Complex Polynomials an Their Symmetric Tensor Representations Bo JIANG Zhening LI Shuzhong ZHANG December 31, 2014 Abstract In this paper we stuy multivariate polynomial

More information

Calculus in the AP Physics C Course The Derivative

Calculus in the AP Physics C Course The Derivative Limits an Derivatives Calculus in the AP Physics C Course The Derivative In physics, the ieas of the rate change of a quantity (along with the slope of a tangent line) an the area uner a curve are essential.

More information

On the Aloha throughput-fairness tradeoff

On the Aloha throughput-fairness tradeoff On the Aloha throughput-fairness traeoff 1 Nan Xie, Member, IEEE, an Steven Weber, Senior Member, IEEE Abstract arxiv:1605.01557v1 [cs.it] 5 May 2016 A well-known inner boun of the stability region of

More information

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that

there is no special reason why the value of y should be fixed at y = 0.3. Any y such that 25. More on bivariate functions: partial erivatives integrals Although we sai that the graph of photosynthesis versus temperature in Lecture 16 is like a hill, in the real worl hills are three-imensional

More information

Some Examples. Uniform motion. Poisson processes on the real line

Some Examples. Uniform motion. Poisson processes on the real line Some Examples Our immeiate goal is to see some examples of Lévy processes, an/or infinitely-ivisible laws on. Uniform motion Choose an fix a nonranom an efine X := for all (1) Then, {X } is a [nonranom]

More information

Implicit Differentiation

Implicit Differentiation Implicit Differentiation Thus far, the functions we have been concerne with have been efine explicitly. A function is efine explicitly if the output is given irectly in terms of the input. For instance,

More information

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms

On Characterizing the Delay-Performance of Wireless Scheduling Algorithms On Characterizing the Delay-Performance of Wireless Scheuling Algorithms Xiaojun Lin Center for Wireless Systems an Applications School of Electrical an Computer Engineering, Purue University West Lafayette,

More information

arxiv: v1 [math.co] 15 Sep 2015

arxiv: v1 [math.co] 15 Sep 2015 Circular coloring of signe graphs Yingli Kang, Eckhar Steffen arxiv:1509.04488v1 [math.co] 15 Sep 015 Abstract Let k, ( k) be two positive integers. We generalize the well stuie notions of (k, )-colorings

More information

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k

Robust Forward Algorithms via PAC-Bayes and Laplace Distributions. ω Q. Pr (y(ω x) < 0) = Pr A k A Proof of Lemma 2 B Proof of Lemma 3 Proof: Since the support of LL istributions is R, two such istributions are equivalent absolutely continuous with respect to each other an the ivergence is well-efine

More information

WUCHEN LI AND STANLEY OSHER

WUCHEN LI AND STANLEY OSHER CONSTRAINED DYNAMICAL OPTIMAL TRANSPORT AND ITS LAGRANGIAN FORMULATION WUCHEN LI AND STANLEY OSHER Abstract. We propose ynamical optimal transport (OT) problems constraine in a parameterize probability

More information

Separation of Variables

Separation of Variables Physics 342 Lecture 1 Separation of Variables Lecture 1 Physics 342 Quantum Mechanics I Monay, January 25th, 2010 There are three basic mathematical tools we nee, an then we can begin working on the physical

More information

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp

. ISSN (print), (online) International Journal of Nonlinear Science Vol.6(2008) No.3,pp . ISSN 1749-3889 (print), 1749-3897 (online) International Journal of Nonlinear Science Vol.6(8) No.3,pp.195-1 A Bouneness Criterion for Fourth Orer Nonlinear Orinary Differential Equations with Delay

More information

II. First variation of functionals

II. First variation of functionals II. First variation of functionals The erivative of a function being zero is a necessary conition for the etremum of that function in orinary calculus. Let us now tackle the question of the equivalent

More information

Iterated Point-Line Configurations Grow Doubly-Exponentially

Iterated Point-Line Configurations Grow Doubly-Exponentially Iterate Point-Line Configurations Grow Doubly-Exponentially Joshua Cooper an Mark Walters July 9, 008 Abstract Begin with a set of four points in the real plane in general position. A to this collection

More information

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency

Transmission Line Matrix (TLM) network analogues of reversible trapping processes Part B: scaling and consistency Transmission Line Matrix (TLM network analogues of reversible trapping processes Part B: scaling an consistency Donar e Cogan * ANC Eucation, 308-310.A. De Mel Mawatha, Colombo 3, Sri Lanka * onarecogan@gmail.com

More information

Monotonicity of facet numbers of random convex hulls

Monotonicity of facet numbers of random convex hulls Monotonicity of facet numbers of ranom convex hulls Gilles Bonnet, Julian Grote, Daniel Temesvari, Christoph Thäle, Nicola Turchi an Florian Wespi arxiv:173.31v1 [math.mg] 7 Mar 17 Abstract Let X 1,...,

More information

θ x = f ( x,t) could be written as

θ x = f ( x,t) could be written as 9. Higher orer PDEs as systems of first-orer PDEs. Hyperbolic systems. For PDEs, as for ODEs, we may reuce the orer by efining new epenent variables. For example, in the case of the wave equation, (1)

More information

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors

Math Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+

More information

Extreme Negative Dependence and Risk Aggregation

Extreme Negative Dependence and Risk Aggregation Extreme Negative Dependence and Risk Aggregation Bin Wang and Ruodu Wang January, 015 Abstract We introduce the concept of an extremely negatively dependent (END) sequence of random variables with a given

More information

The Sokhotski-Plemelj Formula

The Sokhotski-Plemelj Formula hysics 25 Winter 208 The Sokhotski-lemelj Formula. The Sokhotski-lemelj formula The Sokhotski-lemelj formula is a relation between the following generalize functions (also calle istributions), ±iǫ = iπ(),

More information

Calculus of Variations

Calculus of Variations 16.323 Lecture 5 Calculus of Variations Calculus of Variations Most books cover this material well, but Kirk Chapter 4 oes a particularly nice job. x(t) x* x*+ αδx (1) x*- αδx (1) αδx (1) αδx (1) t f t

More information

A Sketch of Menshikov s Theorem

A Sketch of Menshikov s Theorem A Sketch of Menshikov s Theorem Thomas Bao March 14, 2010 Abstract Let Λ be an infinite, locally finite oriente multi-graph with C Λ finite an strongly connecte, an let p

More information

Bounds on the value-at-risk for the sum of possibly dependent risks

Bounds on the value-at-risk for the sum of possibly dependent risks Insurance: Mathematics and Economics 37 (2005) 135 151 Bounds on the value-at-risk for the sum of possibly dependent risks Mhamed Mesfioui, Jean-François Quessy Département de Mathématiques et d informatique,

More information

Introduction to variational calculus: Lecture notes 1

Introduction to variational calculus: Lecture notes 1 October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational

More information

Tail Mutual Exclusivity and Tail- Var Lower Bounds

Tail Mutual Exclusivity and Tail- Var Lower Bounds Tail Mutual Exclusivity and Tail- Var Lower Bounds Ka Chun Cheung, Michel Denuit, Jan Dhaene AFI_15100 TAIL MUTUAL EXCLUSIVITY AND TAIL-VAR LOWER BOUNDS KA CHUN CHEUNG Department of Statistics and Actuarial

More information

How to get bounds for distribution convolutions? A simulation study and an application to risk management

How to get bounds for distribution convolutions? A simulation study and an application to risk management How to get bounds for distribution convolutions? A simulation study and an application to risk management V. Durrleman, A. Nikeghbali & T. Roncalli Groupe de Recherche Opérationnelle Crédit Lyonnais France

More information

SHARP BOUNDS FOR EXPECTATIONS OF SPACINGS FROM DECREASING DENSITY AND FAILURE RATE FAMILIES

SHARP BOUNDS FOR EXPECTATIONS OF SPACINGS FROM DECREASING DENSITY AND FAILURE RATE FAMILIES APPLICATIONES MATHEMATICAE 3,4 (24), pp. 369 395 Katarzyna Danielak (Warszawa) Tomasz Rychlik (Toruń) SHARP BOUNDS FOR EXPECTATIONS OF SPACINGS FROM DECREASING DENSITY AND FAILURE RATE FAMILIES Abstract.

More information

Linear First-Order Equations

Linear First-Order Equations 5 Linear First-Orer Equations Linear first-orer ifferential equations make up another important class of ifferential equations that commonly arise in applications an are relatively easy to solve (in theory)

More information

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control

19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control 19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior

More information