Approximation of the conditional number of exceedances

Transcription

1 Approximation of the conditional number of exceedances Han Liang Gan University of Melbourne July 9, 2013 Joint work with Aihua Xia Han Liang Gan Approximation of the conditional number of exceedances 1 / 16

2 Conditional exceedances Han Liang Gan Approximation of the conditional number of exceedances 2 / 16

3 Conditional exceedances Given a sequence of identical but not independent random variables X 1,..., X n, define n N s,n = 1 {Xi >s}. i=1 Han Liang Gan Approximation of the conditional number of exceedances 2 / 16

4 Conditional exceedances Given a sequence of identical but not independent random variables X 1,..., X n, define n N s,n = 1 {Xi >s}. i=1 The fragility index of order m, denoted FI n (m) is defined as FI n (m) = lim s E(N s,n N s,n m). Han Liang Gan Approximation of the conditional number of exceedances 2 / 16

5 Conditional exceedances Given a sequence of identical but not independent random variables X 1,..., X n, define n N s,n = 1 {Xi >s}. i=1 The fragility index of order m, denoted FI n (m) is defined as FI n (m) = lim s E(N s,n N s,n m). Applications: Insurance. Han Liang Gan Approximation of the conditional number of exceedances 2 / 16

6 Compound Poisson limit Under certain conditions, Tsing et. al. (1988) showed that the exceedance process converges to a Compound Poisson limit. Han Liang Gan Approximation of the conditional number of exceedances 3 / 16

7 Assumptions Han Liang Gan Approximation of the conditional number of exceedances 4 / 16

8 Assumptions Existence of the limits. Han Liang Gan Approximation of the conditional number of exceedances 4 / 16

9 Assumptions Existence of the limits. Tail behaviour is nice enough. Han Liang Gan Approximation of the conditional number of exceedances 4 / 16

10 Assumptions Existence of the limits. Tail behaviour is nice enough. Mixing condition. Han Liang Gan Approximation of the conditional number of exceedances 4 / 16

11 Assumptions Existence of the limits. Tail behaviour is nice enough. Mixing condition. Uniform convergence. Han Liang Gan Approximation of the conditional number of exceedances 4 / 16

12 Assumptions Existence of the limits. Tail behaviour is nice enough. Mixing condition. Uniform convergence. Uniform integrability. Han Liang Gan Approximation of the conditional number of exceedances 4 / 16

13 Alternatives In contrast to limit theorems, we can use approximation theory. Our tool of choice will be Stein s method. Han Liang Gan Approximation of the conditional number of exceedances 5 / 16

14 Alternatives In contrast to limit theorems, we can use approximation theory. Our tool of choice will be Stein s method. For the next section of the talk, we shall stick to conditioning on just one exceedance, Poisson approximation, and total variation distance. Han Liang Gan Approximation of the conditional number of exceedances 5 / 16

15 A simple example Let X 1,..., X n be a sequence of i.i.d. exponential random variables, p = p s = P(X 1 > s) and Z λ Po(λ). Han Liang Gan Approximation of the conditional number of exceedances 6 / 16

16 A simple example Let X 1,..., X n be a sequence of i.i.d. exponential random variables, p = p s = P(X 1 > s) and Z λ Po(λ). d TV (L(N 1 ), Po 1 (λ)) = 1 2 P(N 1 = j) P(Zλ 1 = j). j=1 Han Liang Gan Approximation of the conditional number of exceedances 6 / 16

17 A simple example (2) Set λ such that e λ = (1 p) n. Han Liang Gan Approximation of the conditional number of exceedances 7 / 16

18 A simple example (2) Set λ such that e λ = (1 p) n. So now we have d TV (L(N 1 ), Po 1 (λ)) = 1 2(1 P(N = 0)) P(N = j) P(Z λ = j). j=1 Han Liang Gan Approximation of the conditional number of exceedances 7 / 16

19 A simple example (3) Set λ = np, and we can reduce the problem to known results. d TV (L(N 1 ), Po 1 (λ)) = 1 2(1 P(N = 0)) [ P(N = j) P(Z λ = j) j=1 + P(Z λ = j) P(Z λ = j) ]. Han Liang Gan Approximation of the conditional number of exceedances 8 / 16

20 A simple example (3) Set λ = np, and we can reduce the problem to known results. d TV (L(N 1 ), Po 1 (λ)) = 1 2(1 P(N = 0)) [ P(N = j) P(Z λ = j) j=1 + P(Z λ = j) P(Z λ = j) ]. Using results from Barbour, Holst and Janson (1992), it can be shown that: d TV (L(N 1 ), Po 1 (λ)) = p + o(p). Han Liang Gan Approximation of the conditional number of exceedances 8 / 16

21 The conditional Poisson Stein Identity Lemma W Po 1 (λ) if and only if for all bounded functions g : Z + R, E [ λg(w + 1) Wg(W ) 1 {W 2} ] = 0. Han Liang Gan Approximation of the conditional number of exceedances 9 / 16

22 Stein s method in one slide We construct a function g for any set A, that satisfies: 1 {j A} Po 1 (λ){a} = λg(j + 1) jg(j) 1 {j 2}. Han Liang Gan Approximation of the conditional number of exceedances 10 / 16

23 Stein s method in one slide We construct a function g for any set A, that satisfies: 1 {j A} Po 1 (λ){a} = λg(j + 1) jg(j) 1 {j 2}. It now follows that, P(W A) Po(λ){A} = E [ ] λg(w + 1) Wg(W ) 1 {W 2}. Han Liang Gan Approximation of the conditional number of exceedances 10 / 16

24 Stein s method in one slide We construct a function g for any set A, that satisfies: 1 {j A} Po 1 (λ){a} = λg(j + 1) jg(j) 1 {j 2}. It now follows that, P(W A) Po(λ){A} = E [ ] λg(w + 1) Wg(W ) 1 {W 2}. We then only need to bound the right hand side. Han Liang Gan Approximation of the conditional number of exceedances 10 / 16

25 Generator interpretation If we set g(w) = h(w) h(w 1), the equation becomes: Han Liang Gan Approximation of the conditional number of exceedances 11 / 16

26 Generator interpretation If we set g(w) = h(w) h(w 1), the equation becomes: E [ λ (h(w + 1) h(w )) W (h(w ) h(w 1))1 {W 2} ]. Han Liang Gan Approximation of the conditional number of exceedances 11 / 16

27 Generator interpretation If we set g(w) = h(w) h(w 1), the equation becomes: E [ λ (h(w + 1) h(w )) W (h(w ) h(w 1))1 {W 2} ]. This looks like the generator for a immigration death process. Han Liang Gan Approximation of the conditional number of exceedances 11 / 16

28 Generator interpretation If we set g(w) = h(w) h(w 1), the equation becomes: E [ ] λ (h(w + 1) h(w )) W (h(w ) h(w 1))1 {W 2}. This looks like the generator for a immigration death process. Moreover, the stationary distribution for this generator is Po 1 (λ). Han Liang Gan Approximation of the conditional number of exceedances 11 / 16

29 Stein Factors We can use theorem 2.10 from Brown & Xia (2001) to get: Han Liang Gan Approximation of the conditional number of exceedances 12 / 16

30 Stein Factors We can use theorem 2.10 from Brown & Xia (2001) to get: Theorem The solution g that satisfies the Stein Equation for the total variation distance satisfies: g 1 e λ λe λ λ(1 e λ. ) Han Liang Gan Approximation of the conditional number of exceedances 12 / 16

31 The example revisited Han Liang Gan Approximation of the conditional number of exceedances 13 / 16

32 The example revisited Recall that from before: d TV (L(N 1 ), Po 1 (λ)) = p + o(p). Han Liang Gan Approximation of the conditional number of exceedances 13 / 16

33 The example revisited Recall that from before: d TV (L(N 1 ), Po 1 (λ)) = p + o(p). Using Stein s method directly for conditional Poisson approximation, it can be shown that: d TV (L(N 1 ), Po 1 (λ )) = p 2 + o(p). Han Liang Gan Approximation of the conditional number of exceedances 13 / 16

34 Generalisations The conditional Poisson Stein Identity generalises for conditioning on multiple exceedances: Han Liang Gan Approximation of the conditional number of exceedances 14 / 16

35 Generalisations The conditional Poisson Stein Identity generalises for conditioning on multiple exceedances: Lemma W Po a (λ) if and only if for all bounded functions g : {a, a + 1,...} R, E [ λg(w + 1) Wg(W ) 1 {W >a} ] = 0. Han Liang Gan Approximation of the conditional number of exceedances 14 / 16

36 Negative Binomial Approximation We can do exactly the same for negative binomial random variables. Lemma W Nb a (r, p) if and only if for all bounded functions g : {a, a + 1,...} R, E [ (1 p)(r + W )g(w ) Wg(W )) 1 {W >a} ] = 0. Han Liang Gan Approximation of the conditional number of exceedances 15 / 16

37 Different metrics Have we solved our original problem? Can we calculate errors when estimating fragility indicies? Han Liang Gan Approximation of the conditional number of exceedances 16 / 16

38 Different metrics Have we solved our original problem? Can we calculate errors when estimating fragility indicies? No, we still have the uniform integrability problem. Han Liang Gan Approximation of the conditional number of exceedances 16 / 16

39 Different metrics Have we solved our original problem? Can we calculate errors when estimating fragility indicies? No, we still have the uniform integrability problem. We need to use a stronger metric, such as the Wasserstein distance. Han Liang Gan Approximation of the conditional number of exceedances 16 / 16