The compound Poisson random variable s approximation to the individual risk model
|
|
- Adrian Mitchell
- 5 years ago
- Views:
Transcription
1 Insurance: Mathematics and Economics 36 (25) The compound Poisson random variable s approximation to the individual risk model Jingping Yang, Shulin Zhou, Zhenyong Zhang LMAM, School of Mathematical Sciences, Peking University, Beijing 1871, China Received December 23; received in revised form October 24; accepted 14 October 24 Abstract In this paper we study approximating the total loss associated with the individual insurance risk model by a compound Poisson random variable. By minimizing the expectation of the absolute deviation of the compound Poisson random variable from the true total loss, we investigate not only the optimal compound Poisson random variable but also the numerical calculation of the approximation error. We also discuss the influence of the Poisson parameter on the approximation error. 24 Elsevier B.V. All rights reserved. JEL classification: C63; O29 2 Subject and Insurance Branch Codes: IM2 Keywords: Compound Poisson approximation; Individual risk model; Wasserstein distance; Comonotonic 1. Introduction In actuarial science, the approximation of the individual risk model by a compound Poisson model plays an important role. This is mainly due to the compound Poisson model s advantages in recursive calculation (Panjer, 1981), combination and decomposition (Panjer and Willmot, 1992, Chapter 6 or Kaas et al., 21, Chapter 3). In the development on this topic, there are a lot of papers concerned with the approximation. Bühlmann et al. (1977) illustrated why a cautious insurer should prefer the compound Poisson model to the individual risk model in the sense of stop-loss order. Gerber (1979, Chapter 4) gave a description of the choice of the Poisson parameter and introduced two cases which are often used in the later discussions. Gerber (1984), Hipp (1985, 1986), Michel (1987), De Pril and Dhaene (1992), Sundt (1993) and Dhaene and Sundt (1997) investigated the error bounds for approximation in Corresponding author. addresses: yangjp@math.pku.edu.cn (J. Yang), szhou@math.pku.edu.cn (S. Zhou) /$ see front matter 24 Elsevier B.V. All rights reserved. doi:1.116/j.insmatheco
2 58 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) terms of distribution or stop-loss premium. Kaas et al. (1988b) discussed the approximation of the aggregate claims and the stop-loss premiums by approximating the aggregate claims by the sum of a compound Poisson random variable (r.v.) and another r.v. determined by stop-loss order. Kuon et al. (1993) studied the approximation quality when the portfolio keeps growing. The former papers are mainly focused on approximation in the aggregate claims distribution and related functions, such as stop-loss premiums. However, in practice one initially concerns about the individual risk model. Thus one natural question arises: given the observation data from the individual risk model, how should one determine the r.v. s in the corresponding compound Poisson model? From the practical viewpoint, the r.v. s in the approximation model have to be determined from the observation data. The aim of this paper is to develop a method for carrying out such an approximation. By minimizing the expectation of the absolute deviation of compound Poisson r.v. s from the total loss associated with the individual risk model, we present an optimal approximation model. We also give a numerical method to evaluate the approximation error. Finally we discuss the influence of the Poisson parameter on the approximation error. Throughout this paper, it is always assumed that the individual risks are independent. We first consider the case that the individual risks are homogenous, then apply the homogenous results to approximate the heterogenous risk model. Consider a portfolio containing n homogenous insurance risks. Let X i denote the loss associated with the ith risk, i = 1, 2,...,n. Assume that X 1,X 2,...,X n are independent identically distributed with common distribution F, where the variance Var(X 1 ) is finite and <F() < 1. The total loss for the portfolio {X 1,X 2,...,X n } can be written as S n = X i. (1.1) The number of claims for the portfolio is denoted as N n, i.e., N n = #{i : X i >,i n}. For the technical reason, an auxiliary sequence of risks {X n+1,x n+2,...,} is introduced, where X 1,X 2,... are independent identically distributed. Denote M 1 = inf{i 1:X i > }, M n = inf{i >M n 1 : X i > }, n 2. It is known that P(M i < ) = 1. M i is the index of the risk that a claim occurs and M i is strictly increasing with respect to i. By using the auxiliary risks {X n+1,x n+2,...}, we can define a claim sequence Y i = X Mi, i 1, which is a subsequence of {X i,i 1} and is well-defined. It is proved that Y i,i 1 are independent with common distribution function F Y (x) = F(x) F(), x. 1 F() It can also be shown that N n is Binomial(n, q) with parameter q = 1 F(), P(N n = m) = C m n qm (1 q) n m, m =, 1,...,n,
3 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) and N n is independent of the sequence Y i,i 1. Here Cn m denotes (n m ). Based on the above method, the link between the homogenous individual risk model and the compound binomial model has been established. We refer to Li and Yang (21) for the above theoretical results. By the definition of {Y i }, we know that S n in (1.1) can be expressed as N n S n = Y i. (1.2) Each term Y i of S n in (1.2) corresponds to an actual claim while there are many terms X i in (1.1) which may equal zero. Among these terms in (1.2), the individual claims Y i, i 1 are independent and identically distributed, N n is the number of claims in the portfolio {X 1,X 2,...,X n }, and the independence between N n and Y i is fulfilled. Summing over the first N n claims of the sequence Y i,i 1, we thus obtain the total loss S n associated with the portfolio {X 1,X 2,...,X n }. Now we consider the approximation model. It is natural to assume that the approximation model has the same claim sequence as the individual risk model. The number of claims N(θ), as a Poisson r.v., should be estimated from the observation data. Furthermore, the independence between N(θ) and Y i,i 1 is assumed. An approximation to S n is defined by N(θ) S = Y i. (1.3) In the above approximation, the main problem is to determine the Poisson r.v. N(θ). Let F poi(θ) denote Poisson distribution with mean θ, and F bin(n,q) denote binomial distribution with parameters (n, q), where q = 1 F(). The family of all Poisson r.v. s, which have common mean θ and are independent of the claim sequence Y i,i 1, is denoted as R(F poi(θ) ). As mentioned above, the optimal Poisson r.v. N n (θ) R(F poi(θ)) is determined by the following minimizing principle: H n (θ) =: E S n N n (θ) Y i = inf E S n S, (1.4) N(θ) R(F poi(θ) ) where S n and S are given, respectively, in (1.2) and (1.3). Here H n (θ) represents the approximation error. We organize our paper in the following frame. In Sections 2 6 we investigate the homogenous individual risk models. We first concentrate on the case that the Poisson parameter θ is fixed, then study the influence of the Poisson parameter on the approximation error. In Section 2 we show the existence of the optimal Poisson r.v. Nn (θ), and give an explicit expression for Nn (θ). In Section 3 we prove some optimal properties of N n (θ). In Section 4 we present a numerical method to evaluate the approximation error. In Section 5 we apply our approximation argument to the functions of the total loss. In Section 6 we demonstrate that H n (θ) has a unique minimum point, and make a comparison between the minimum Poisson parameter and those parameters normally used in literatures. In Section 7 we apply our main results to the heterogenous individual risk model. In Appendix A we give the detailed proof of the main theorems in Section The existence of the optimal Poisson r.v. Given the risk sequence {X 1,X 2,...}, ar.v.u can be defined as below. Construct r.v. s U m,m=, 1,...,n such that U m,m=, 1,...,nand X i,i 1 are independent. It is assumed that U m is uniformly distributed over
4 6 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) (F bin(n,q) (m 1)),F bin(n,q) (m)] with probability density function 1 F bin(n,q) (m) F bin(n,q) (m 1). Then define U = U m I {Nn =m}. m= Under the above construction, U m,m=, 1, 2,...,nare independent of the number of claims N n and the claim sequence Y i,i= 1, 2,... Thus by the independence between N n and Y i, U is also independent of Y i. Moreover, 1, x F bin(n,q) (m), x F bin(n,q) (m 1) P(U x N n = m) = F bin(n,q) (m) F bin(n,q) (m 1),F bin(n,q)(m 1) x<f bin(n,q) (m), (2.1), x<f bin(n,q) (m 1). The inverse function F 1 bin(n,q) (y)off bin(n,q)(y) is defined as F 1 bin(n,q) (y) = inf{x : F bin(n,q)(x) y}, y [, 1] and the inverse function F 1 poi(θ) (y)off poi(θ)(y) is defined as F 1 poi(θ) (y) = inf{x : F poi(θ)(x) y}. Lemma 1. The r.v. U is uniformly distributed over [, 1], and N n = F 1 bin(n,q) (U). Proof. For fixed s (, 1), denote M = sup{k : F bin(n,q) (k) s, k =, 1,...}. It follows from (2.1) that: and Thus P(U s N n = M + 1) = s F bin(n,q) (M) F bin(n,q) (M + 1) F bin(n,q) (M) P(U s N n = m) =, m>m+ 1; P(U s N n = m) = 1, m M. P(U s) = = P(U s N n = m)p(n n = m) m= m=m+1 P(U s N n = m)p(n n = m) + M P(U s N n = m)p(n n = m) m=
5 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) = P(U s N n = M + 1)P(N n = M + 1) + M P(N n = m) s F bin(n,q) (M) = F bin(n,q) (M + 1) F bin(n,q) (M) {F bin(n,q)(m + 1) F bin(n,q) (M)}+F bin(n,q) (M) = s. That is, the r.v. U is uniformly distributed over [, 1]. Next we prove that N n = F 1 bin(n,q) (U) holds. Fix non-negative integer m, m n. Suppose that N n = m. Thus F bin(n,q) (m 1) <U F bin(n,q) (m) follows. Equivalently, F 1 bin(n,q) (U) = m. Hence N n = F 1 bin(n,q) (U) = m. The lemma is proved. In the following theorem we prove the existence of the optimal Poisson r.v. Nn (θ), and then give a mathematical expression for the corresponding approximation error. Theorem 1. The Poisson r.v. Nn (θ) defined by N 1 1 n (θ) = Fpoi(θ) (U) satisfies (1.4), i.e., Fpoi(θ) (U) is a possible choice for the optimal Poisson r.v. with mean θ. Moreover, the approximation error H n (θ) satisfies H n (θ) = E(Y 1 )E F 1 1 bin(n,q) (U) Fpoi(θ) (U) =E(Y 1) m= F bin(n,q) (x) F poi(θ) (x) dx. (2.2) Proof. Let R(F bin(n,q),f poi(θ) ) denote the family containing all random vectors (W, V ), where W has distribution F bin(n,q) and V has distribution F poi(θ), and (W, V ) is independent of the claim sequence Y i,i 1. Using the independence relation between Y i,i 1 and R(F bin(n,q),f poi(θ) ), we find that inf E W V Y (W,V ) R(F bin(n,q),f poi(θ) ) i Y i = inf E W V (W,V ) R(F bin(n,q),f poi(θ) ) I {V<W} Y i + I {V>W} Y i i=v +1 i=w+1 = E(Y 1 ) inf E W V. (2.3) (W,V ) R(F bin(n,q),f poi(θ) ) From the result of L 1 -Wassertein distance (Barrio et al., 1999, orshorack and Wellner, 1986), it holds that 1 1 inf E W V =E F (W,V ) R(F bin(n,q),f poi(θ) ) bin(n,q) (U) Fpoi(θ) (U) = F bin(n,q) (x) F poi(θ) (x) dx. (2.4) It is obvious that H n (θ) = inf E N n N Y i Y i = inf N R(F poi(θ) ) E W V Y i Y i. (2.5) (W,V ) R(F bin(n,q),f poi(θ) ) Combining Eqs. (2.3) (2.5), we obtain (2.2) and complete the proof. In this paper we choose the optimal Poisson r.v. as Nn 1 (θ) = Fpoi(θ) (U). (2.6) In modern risk theory, comonotonicity is an important concept. The two risks P and Q are comonotonic if there exist two non-decreasing real-valued functions u, v and a risk Z such that P = u(z), Q = v(z) (see Wang et al., 1997). By Lemma 1 and (2.6), Nn (θ) and N n are both non-decreasing functions of the r.v. U. Thus they are comonotonic. It means that with the varying of U, Nn (θ) increases when N n increases, and Nn (θ) decreases when N n decreases.
6 62 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Some results on N n (θ) The joint distribution F Nn,N(θ) of (N n,n(θ)) satisfies max{f bin(n,q) (x) + F poi(θ) (y) 1, } F Nn,N(θ)(x, y) min{f bin(n,q) (x),f poi(θ) (y)}. (3.1) See Joe (1997) for its related discussion. Note that the distribution function of (N n,nn (θ)) equals the upper bound min{f bin(n,q) (x),f poi(θ) (y)} Some optimal properties We will show some optimal properties of Nn (θ) in the following theorem. Theorem 2. The optimal r.v. Nn (θ) in (2.6) is optimal in the following senses: and E(N n N n (θ))2 = Var(N n N n (θ)) = inf E(N n N(θ)) 2 (3.2) N(θ) R(F poi(θ) ) inf Var(N n N(θ)). (3.3) N(θ) R(F poi(θ) ) Proof. Recalling inequality (3.1) we have that for x, y [, ), P(N n >x,n(θ) >y) min{p(n n >x),p(n(θ) >y)}. Applying Hoeffding s identity (Joe, 1997, p. 23), we obtain E{N n N(θ)} = Thus we have = P(N n >x,n(θ) >y)dx dy P(N n >x,n n (θ) >y)dx dy = E{N nn n (θ)}. min{p(n n >x),p(n(θ) >y)} dx dy E(N n N(θ)) 2 = EN 2 n + E(N(θ))2 2E{N n N(θ)} EN 2 n + E(N(θ))2 2E{N n N n (θ)} which implies (3.2). By using (3.2) and = E(N n N n (θ))2, Var(N n N(θ)) = E(N n N(θ)) 2 (EN n EN(θ)) 2 = E(N n N(θ)) 2 (nq θ) 2, we obtain (3.3). Thus we complete the proof of the theorem. The above theorem asserts the optimality of Nn (θ) in the senses of (3.2) and (3.3) The joint probability of (N n,n n (θ)) In this section we give several formulas to evaluate the joint distribution of (N n,n n (θ)).
7 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Given N n = m, m =, 1,...,n, the conditional probability P(Nn (θ) = k N n = m) = P(F 1 poi(θ) (U) = k U (F bin(n,q)(m 1),F bin(n,q) (m)]) = P(U (F poi(θ) (k 1),F poi(θ) (k)) U (F bin(n,q) (m 1),F bin(n,q) (m)]) = min{f poi(θ)(k),f bin(n,q) (m)} max{f poi(θ) (k 1),F bin(n,q) (m 1)}. F bin(n,q) (m) F bin(n,q) (m 1) Then the joint probability P(N n = m, Nn (θ) = k) = P(N n (θ) = k N n = m)p(n n = m) = (min{f poi(θ) (k),f bin(n,q) (m)} max{f poi(θ) (k 1),F bin(n,q) (m 1)}). (3.4) Furthermore, the probability P(N n (θ) = N n) can be evaluated by P(N n (θ) = N n) = = P(Nn (θ) = N n = i) (min{f poi(θ) (i),f bin(n,q) (i)} max{f poi(θ) (i 1),F bin(n,q) (i 1)}). i= Table 1 gives some values of the joint probabilities of (N n,n n (θ)) in the case n = 1,θ = 1,p=.1,EY 1 = 1, calculated by (3.4). Table 2 gives some results in the case n = 1,θ =.1,p=.1,EY 1 = 1. The values in the two tables show the comonotonic trend between N n and N n (θ). Table 1 The probability P(N n = i, N n (θ) = j) when n = 1,θ = 1,p=.1 and EY 1 = 1 j = j = 1 j = 2 j = 3 j = 4 j = 5 j = 6 i =.3677 i = E 8 i = E 5 i = E 5 i = E 5 i = E 6 i = 6.5 i = 7 Table 2 The probability P(N n = i, N n (θ) = j) when n = 1,θ =.1,p=.1 and EY 1 = 1 j = j = 1 j = 2 j = 3 j = 4 j = 5 j = 6 i =.99 i = E E 6 i = E E 8 i = E E 1 i = E E 13 i = E E 16 i = E 16
8 64 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Evaluating the approximation error H n (θ) measures the expectation of the absolute approximation error. In this section we will show how to compute the error. At the beginning we will list a useful lemma. Lemma 2. If θ nlog(1 q), then F bin(n,q) (t) F poi(θ) (t) for all t. Conversely, if F bin(n,q) () F poi(θ) (), then θ nlog(1 q). Lemma 2 can be obtained directly by Theorem 1 of De Pril and Dhaene (1992). Theorem 3. (a) If θ n log(1 q), then H n (θ) = (θ nq)ey 1. (b) If θ< n log(1 q), then i H n (θ) = 2 (i j)p(n n = i, Nn (θ) = j) + θ nq EY 1. (4.2) Proof. i= j= (a) In the case θ n log(1 q), by Lemma 2 we have (4.1) P(N n >t) P(Nn (θ) >t), t [, ]. Then using (2.2), we obtain H n (θ) = EY 1 P(N n t) P(Nn (θ) t) dt = EY 1 = EY 1 (EN n (θ) EN n) = EY 1 (θ nq). (b) In the case θ< n log(1 q), by Theorem 1 we have {P(N n (θ) >t) P(N n >t)} dt H n (θ) = EY 1 E N n N n (θ) =EY 1 {2E(N n N n (θ)) + E(N n N n (θ)} = EY 1 {2E(N n Nn (θ)) + + θ nq} i = 2 (i j)p(n n = i, Nn (θ) = j) + θ nq EY 1, i= j= where the probability function of (N n,nn (θ)) can be calculated by (3.4). The Poisson parameter θ is often chosen as θ = nq or θ = nlog(1 q) (Gerber, 1979, Chapter 4). Tables 3 5 provide some numerical values of H n (nq) and H n ( n log(1 q)) when EY 1 = 1. By comparing H n (nq) with H n ( n log(1 q)), one can analyze the accuracy of the approximation. For instance, in the case nq =.5, H 1 (nq) >H 1 ( n log(1 q)),h 1 (nq) <H 1 ( n log(1 q)) and H 1 (nq) <H 1 ( n log(1 q)). Generally speaking, one cannot assert that the parameter nq is optimal than the parameter n log(1 q).
9 Table 3 The expected errors when n = 1 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) q nq n log(1 q) H n (nq) H n ( n log(1 q)) Table 4 The expected errors when n = 1 q nq n log(1 q) H n (nq) H n ( n log(1 q)) E E Table 5 The expected errors when n = 1 q nq n log(1 q) H n (nq) H n ( n log(1 q)) Approximation to the total loss and related functions Given the values {N n,y 1,Y 2,...,Y Nn } determined from the portfolio {X 1,X 2,...,X n }, the optimal Poisson r.v. in (2.6) can be determined as follows. First simulate the r.v. U by (2.1), then calculate the Poisson r.v. Nn (θ) = F 1 poi(θ) (U). In the case N n Nn (θ), the values of the claim sequence Y 1,Y 2,...,Y N n (θ) can be obtained from the set {Y 1,Y 2,...,Y Nn }. Otherwise, simulation is needed to obtain the r.v. s {Y Nn +1,...,Y N n (θ) }. Now we consider the total amount of the type n g(x i ), where g is a non-negative measurable function and g() =. An approximation for n g(x i )is N n (θ) g(y i ). The corresponding approximation error equals h n (θ, g) =: N n (θ) g(x i ) g(y i ).
10 66 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Note that when g(x) = x, S n = N n (θ) Y i + h n (θ, g). Theorem 4. For the non-negative function g with g() =, Eh n (θ, g) = (nq θ)eg(y 1 ) (5.1) and E h n (θ, g) = Eg(Y 1) H n (θ). EY 1 Proof. It is easy to prove (5.1). Here we only prove (5.2) in detail. Since g() =, we have (5.2) Thus N n g(x i ) = g(y i ). N n N n (θ) E h n (θ, g) =E g(y i ) g(y i ) = Eg(Y 1) H n (θ). EY 1 Therefore, the theorem is proved. Then We would like to mention one interesting fact when θ = n log(1 q). In this case, by Lemma 2 we have Hence F 1 1 poi(θ) (s) Fbin(n,q) (s), s (, 1) N n = F 1 1 bin(n,q) (U) Fpoi(θ) (U) = N n (θ). N n g(x i ) = g(y i ) N n (θ) g(y i ). Thus the approximated total amount is always greater than the actual aggregate claims. For the portfolio {X 1,X 2,...,X n }, we consider the excess-of-loss reinsurance treaty with retention M. The direct insurer pays the amount S D,n (M) = X i M
11 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) and the reinsurer pays the amount S R,n (M) = (X i M) +, where X i M = min{x i,m} and (X i M) + = max{x i M, }. Use the two compound Poisson r.v. s S AP D,n N n (θ) (M) = Y i M and S AP R,n N n (θ) (M) = (Y i M) + to approximate S D,n (M) and S R,n (M), respectively. By Theorem 4, the expected approximation errors are E S R,n (M) SR,n AP (M) =E(Y 1 M) H n (θ), E S D,n (M) SD,n AP EY (M) =E(Y 1 M) + H n (θ), 1 EY 1 respectively. Note that the sum of the expected approximation errors equals H n (θ), i.e., E S R,n (M) S AP R,n (M) +E S D,n(M) S AP D,n (M) =H n(θ). 6. The uniqueness of the Poisson parameter to minimizing H n (θ) In this section we investigate whether there exists a unique θn, such that H n (θ n ) = min θ H n(θ). (6.1) The answer to this problem is affirmative. For k =, 1, 2,...,n 1, denote θ (k) to be the solution of equation F bin(n,q) (k) = F poi(θ (k) ) (k). Then we have the following result. Theorem 5. It holds that θ (n 1) >θ (n 2) > >θ (1) >θ (). Denote 1 iff bin(n,q) (x) >F poi(θ) (x), sign(f bin(n,q) (x) F poi(θ) (x)) = iff bin(n,q) (x) = F poi(θ) (x), 1 iff bin(n,q) (x) <F poi(θ) (x). In the following we give a theorem on the uniqueness of θ n and present an approach to solve θ n. Here for simplicity we assume EY 1 = 1.
12 68 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Theorem 6. There exists a unique θn with θ n nlog(1 q), such that (6.1) holds. Further, for θ (θ (k+1),θ (k) ),k =, 1, 2,...,n 1, the derivative function of H n with respect to θ satisfies n 1 H n (θ) = (sign(f bin(n,q) (k) F poi(θ) (k)) 1) e θ θ k k= k! + 1. (6.2) H n (θ) is strictly increasing for θ θ(k),k =, 1, 2,...,n 1 and H n (θ n ), H n (θ n +), (6.3) where H n (θ n ) and H n (θ n +) denote, respectively, the left and right limits of H n at point θ n. We will prove Theorems 5 and 6 in Appendix A. In the case n = 1, (6.2) can be used to find θ 1.Forθ< log(1 q), F bin(1,q)() <F poi(θ) (). Then by (6.2) we have H 1 (θ) = 2e θ + 1. Solving H 1 (θ) =, we have θ = log 2. Moreover, H 1 (θ) < for θ<log 2. Then we can conclude that in the case q>.5, H n (θ) achieves its minimum log 2 + q 1 at point θ1 = log 2; when q.5, H n(θ) achieves its minimum q log(1 q) at point θ1 = log(1 q). In general, θn can be solved numerically. By Theorem 6, θ n is the solution of H n (θ) = if the solution exits. Otherwise, it satisfies H n (θ ) and H n (θ+). Thus by Theorem 6, we can focus on the strictly increasing function h(θ) defined as n 1 h(θ) = (sign(f bin(n,q) (k) F poi(θ) (k)) 1) e θ θ k k! k= + 1, θ [, ) and use the Newtonian method to locate θn. Table 6 lists some numerical results. It can be seen that when n is large, the difference between θn and nq is very small. Thus nq is a good approximation for θn when nq is large while n log(1 q) is a good approximation when nq is small. 7. Approximation to the heterogenous individual risk model In practice, the individual risks X i of a portfolio are often independent, but not identically distributed. In this section, we will apply the above results to the heterogenous individual risk model. We can divide the heterogenous portfolio into several independent homogenous portfolios, then approximate every homogenous portfolio separately by our method. Assume that the portfolio can be divided into m homogenous portfolios and the aggregate claims from the ith portfolio can be represented by IS(i) for i = 1, 2,...,m. For the ith homogenous portfolio, N (i) IS(i) = Y i,j, j=1
13 Table 6 The optimal θ n and the corresponding error H n(θ n ) J. Yang et al. / Insurance: Mathematics and Economics 36 (25) n q nq n log(1 q) θ n H n (nq) H n ( n log(1 q)) H n (θ n ) E 6 1.1E 6 1.1E where N (i) is the number of claims and the claim sequences Y i,j,j = 1, 2,..., are determined as in Section 1, i = 1, 2,...,m. Then the aggregate claims for the heterogenous portfolio is IS = IS(i) = N (i) Y i,j. j=1 Note that the following properties hold: 1. N (i),i= 1, 2,...,m, are independent. 2. For a given i, Y i,j,j = 1, 2,...are independent identically distributed, here the common distribution is denoted as F i. Further, Y i,j,i= 1, 2,...,m,j = 1, 2,...,are independent. 3. Claim sequences are independent of the numbers of claims. Let CS(i, θ i ) be compound Poisson with Poisson r.v. N(i, θ i ) and claim sequence Y i,j, j = 1, 2,..., N(i, θ i ) and Y i,j, j = 1, 2,...are independent. Here θ i is the Poisson parameter. Then CS(i, θ i ) = N(i,θ i ) j=1 Y i,j. The ith group s optimal Poisson r.v. determined by (2.6) is denoted as N (i, θ i ), the optimal Poisson parameter satisfying (6.1) is denoted as θ (i). The independence between the Poisson r.v. s N(i, θ i ),i= 1, 2,...,mis assumed. Denote CS (i, θ i ) = N (i,θ i ) j=1 CS(θ 1,θ 2,...,θ m ) = Y i,j, CS(i, θ i ) j=1
14 7 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) and CS (θ 1,θ 2,...,θ m ) = CS (i, θ i ). j=1 CS(θ 1,θ 2,...,θ m ) and CS (θ 1,θ 2,...,θ m ) are compound Poisson r.v. s, each with Poisson parameter m θ i and claim distribution j=1 θ j m θ i F j. Among all CS(θ 1,θ 2,...,θ m ) s, we will prove that CS (θ 1,θ 2,...,θ m ) is optimal in the following sense. Theorem 7. It holds that and E(IS CS (θ 1,θ 2,...,θ m )) 2 = E IS CS (θ 1,θ 2,...,θ m ) inf E(IS CS(θ 1,θ 2,...,θ m )) 2 (7.1) N(i,θ i ),i m E IS(i) CS (i, θ i ) Proof. According to the above notations, we have E IS(i) CS(i, θ i ). (7.2) E(IS CS(θ 1,θ 2,...,θ m )) 2 = Var(IS CS(θ 1,θ 2,...,θ m )) +{E(IS CS(θ 1,θ 2,...,θ m ))} 2 { m } 2 = Var(IS(i) CS(i, θ i )) + EY i,1 E(N (i) θ i ) Thus by Theorems 1 and 2, it holds that E(IS CS (θ 1,θ 2,...,θ m )) 2 = = E N (i) N(i, θ i ) Var(Y i,1 ) + Var(N (i) N(i, θ i ))E(Yi,1 2 ) { m 2 + E(Y i,1 )E(N (i) θ i )}. { m E N (i) N (i, θ i ) Var(Y i,1 ) + Var(N (i) N (i, θ i ))E(Yi,1 2 ) + E(Y i,1 )E(N (i) θ i ) { m E N (i) N(i, θ i ) Var(Y i,1 ) + Var(N (i) N(i, θ i ))E(Yi,1 2 ) + E(Y i,1 )E(N (i) θ i ) = E(IS CS(θ 1,θ 2,...,θ m )) 2, } 2 } 2 thus (7.1) holds. It is obvious that (7.2) holds. Now we complete the proof.
15 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) When every homogenous portfolio s size is large enough, it is advised to choose θ i = EN (i),i= 1, 2,...,m.For this choice of the Poisson parameter, the approximation CS (θ 1,θ 2,...,θ m ) produces good approximation effects. It can be seen from the following two aspects: (1) First, from Table 6 we can find that θ (i) EN (i). Thus EN (i) is a good approximation for the optimal Poisson parameter θ (i). (2) Second, (7.2) shows that the approximation error is bounded by the sum of the approximation errors for each homogenous portfolio. Thus the corresponding error is small when every homogenous portfolio s size is large enough. 8. Conclusions In this paper, we presented a new method to approximate the individual risk model by a compound Poisson r.v. We investigated the determination of the r.v. s in the compound Poisson model and the calculation of the approximation error. During our discussion, we first focused on the homogenous individual risk models, then applied the results to the heterogenous individual risk models. Numerical results showed that our approximation model provides a good approximation. Acknowledgements The authors would like to thank the referees for their helpful suggestions. The first author was supported by NSFC grants , and the second author was supported by NSFC grant Appendix A In what follows, we will prove Theorems 5 and 6. First we introduce several lemmas. Lemma A.1. There exists a non-negative integer C n 1 such that P(N n = k + 1) P(N n = k) P(N(θ) = k + 1) >,k<c, P(N(θ) = k) P(N(θ) = k + 1) <,n>k>c, P(N(θ) = k) P(N(θ) = k + 1),k= C. P(N(θ) = k) Proof. It is easy to verify that for k + 1 n, P(N n = k + 1) P(N n = k) = Ck+1 n q k+1 (1 q) n k 1 Cn kqk (1 q) n k = n k q k q (A.1)
16 72 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) and P(N(θ) = k + 1) P(N(θ) = k) = θ k + 1. (A.2) Comparing (A.1) with (A.2), we observe that P(N n = k + 1) P(N n = k) P(N(θ) = k + 1) P(N(θ) = k) if and only if (n k){q/(1 q)} θ, i.e., Denote k n θ(1 q). q { [ C = min n 1, n ]} θ(1 q). q Here [x] denotes the integral part of x. Thus if C<n 1, for the integer n>k>cit holds that P(N n = k + 1) P(N n = k) And if k<c, P(N n = k + 1) P(N n = k) < > P(N(θ) = k + 1). P(N(θ) = k) P(N(θ) = k + 1) P(N(θ) = k) holds. In the case k = C, it holds that P(N n = k + 1) P(N n = k) P(N(θ) = k + 1) P(N(θ) = k) Hence the lemma is proved. It is easy to verify the following lemma. We omit its proof here. Lemma A.2. Suppose a n >,b n > and a n+1 /a n b n+1 /b n,n. Then for n 1, a n + a n 1 + +a a n 1 + +a b n + b n 1 + +b b n 1 + +b. Lemma A.3. Suppose that for non-negative integer M, F bin(n,q) (M) F poi(θ) (M). Then F bin(n,q) (M + 1) >F poi(θ) (M + 1). (A.3) Proof. It is trivial for the case M n. Now we consider the case M<n. Assume θ< n log(1 q). Otherwise, by Lemma 4.1 we see that F bin(n,q) (k) F poi(θ) (k),k =, 1, 2,..., thus (A.3) holds. We will divide our proof into two cases M C and M>C, where C is the integer defined in Lemma A.1.
17 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Case A.1 (M C). Since θ< n log(1 q), it leads to P(N n = ) <P(N(θ) = ). Then according to the definition of C and F bin(n,q) (M) F poi(θ) (M), by Lemmas A.1 and A.2 we conclude that P(N n = M) >P(N(θ) = M). Also by Lemma A.1 we have P(N n = M + 1) P(N n = M) Then it follows that: Thus P(N(θ) = M + 1). P(N(θ) = M) P(N n = M + 1) >P(N(θ) = M + 1). F bin(n,q) (M + 1) = F bin(n,q) (M) + P(N n = M + 1) >F poi(θ) (M) + P(N(θ) = M + 1) = F poi(θ) (M + 1), which implies (A.3). Case A.2 (M >C). If P(N n = M + 1) P(N(θ) = M + 1), then it is obvious that (A.3) holds. Next we consider the case P(N n = M + 1) <P(N(θ) = M + 1). Since P(N n = M + 1) <P(N(θ) = M + 1) and M>C,byLemma A.1 we have Then P(N n = k) <P(N(θ) = k), k M + 1. P(N n M + 1) = 1 > 1 Thus the lemma is proved. k=m+2 k=m+2 P(N n = k) > 1 k=m+2 P(N(θ) = k) P(N(θ) = k) = P(N(θ) M + 1). Proof of Theorem 5. Fix k<n. From the definition F bin(n,q) (k) = F poi(θ (k) )(k) and Lemma A.3, wehave F bin(n,q) (k + 1) >F poi(θ (k) )(k + 1). Since for fixed k the function F poi(θ) (k + 1) is strictly decreasing about θ, wehave θ (k+1) <θ (k). Now we complete the proof of the theorem.
18 74 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Lemma A.4. For θ θ (i),i=, 1,...,n 1, the derivative H n (θ) of H n(θ) exists and n 1 H n (θ) = (sign(f bin(n,q) (k) F poi(θ) (k)) 1) e θ θ k k= k! + 1. (A.4) Proof. Fix θ > and θ θ (i),i=, 1,...,n 1. Denote g(θ, x) = F bin(n,q) (x) F poi(θ) (x). First consider the case θ>θ. From (2.2) we have H n (θ) H n (θ ) θ θ = 1 θ θ = 1 θ θ + 1 θ θ ( g(θ, x) g(θ,x) )dx (g(θ, x)i {g(θ,x)>} g(θ,x)i {g(θ,x)>})dx ( g(θ, x)i {g(θ,x)<} + g(θ,x)i {g(θ,x)<})dx =: A + B. By using the fact that g(θ, x) >g(θ,x) for θ>θ and g(θ,x), we have A = 1 (g(θ, x)i {g(θ,x)>} g(θ,x)i {g(θ,x)>} )dx θ θ + 1 θ θ = 1 θ θ + 1 θ θ (g(θ,x)(i {g(θ,x)>} I {g(θ,x)>})dx (F poi(θ )(x) F poi(θ) (x))i {g(θ,x)>} dx (A.5) (F bin(n,q) (x) F poi(θ )(x))i {g(θ,x)>,g(θ,x)<} dx. (A.6) For the case g(θ, x) > and g(θ,x) <, it holds that F bin(n,q) (x) >F poi(θ) (x) and F bin(n,q) (x) <F poi(θ )(x). Thus in this case it follows that: F bin(n,q) (x) F poi(θ )(x) θ θ I {g(θ,x)>,g(θ,x)<} dx F poi(θ )(x) F poi(θ) (x) I {g(θ,x)>,g(θ,x)<} dx. (A.7) θ θ For x [, ), θ F poi(θ)(x) = e θ θ [x] [x]! (2θ ) [x], θ < 2θ, [x]! where [x] is the integral part of x. Then by the dominated convergence theorem and g(θ, x) g(θ,x)wehave 1 θ θ (F poi(θ )(x) F poi(θ) (x))i {g(θ,x)>} dx e θ θ [x] I {g(θ,x)>} dx, [x]! (A.8)
19 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) and by (A.7) we obtain F bin(n,q) (x) F poi(θ) (x) I {g(θ,x)>,g(θ,x)<} dx θ θ F poi(θ )(x) F poi(θ) (x) I {g(θ,x)>,g(θ,x)<} dx. (A.9) θ θ Combining (A.8) and (A.9), we conclude that lim A = θ θ Similarly, we have lim B = θ θ e θ θ [x] I {g(θ,x)>} dx. [x]! e θ θ [x] I {g(θ,x)<} dx. [x]! Summarizing the above two equations and (A.5), we deduce that H n (θ) H n (θ ) lim = θ θ θ θ Similarly, H n (θ) H n (θ ) lim = θ θ θ θ e θ θ [x] sign(f bin(n,q) (x) F poi(θ )(x)) dx. [x]! e θ θ [x] sign(f bin(n,q) (x) F poi(θ )(x)) dx [x]! holds. Thus combining (A.1) and (A.11), we obtain H n (θ e θ θ [x] ) = sign(f bin(n,q) (x) F poi(θ )(x)) dx [x]! n 1 = sign(f bin(n,q) (k) F poi(θ )(k)) e θ θ k k! k= n 1 + k=n e θ θ k k! = (sign(f bin(n,q) (k) F poi(θ )(k)) 1) e θ θ k + 1. k! k= (A.1) (A.11) Proof of Theorem 6. Fix θ (θ (i+1),θ (i) ). Since we have F bin(n,q) (i) = F poi(θ (i) ) (i), F bin(n,q)(i + 1) = F poi(θ (i+1) )(i + 1), F bin(n,q) (i) <F poi(θ) (i), F bin(n,q) (i + 1) >F poi(θ) (i + 1). Thus by Lemma A.3 we assert that F bin(n,q) (k) <F poi(θ) (k), k i
20 76 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) and F bin(n,q) (k) >F poi(θ) (k), k i + 1. Then by Lemma A.4 we have i H e θ θ k (θ) = 2 k! k= + 1, which is strictly increasing with respect to θ (θ (i+1),θ (i) ). From the above equation it is easy to verify H (θ (i) ) <H (θ (i) +), i =, 1, 2,...,n 1. Thus we conclude that H (θ) is strictly increasing with respect to θ θ (i). Since H n (θ) is continuous, from Theorem 3 we assert that the minimum point θn exists and θ n n log(1 q). By the monotonicity of H n (θ) with respect to θ, we assert that H n (θ n +) and H n (θ n ). Since H n (θ) is strictly increasing, thus θn is unique. Hence the uniqueness is proved. As (6.2) has been proved in Lemma A.4, we now complete the proof of Theorem 6. References Barrio, E.D., Gine, E., Matran, C., Central limit theorems for the Wasserstein distance between the empirical and the true distributions. The Annals of Probability 27 (2), Bühlmann, H., Gagliardi, B., Gerber, H.U., Straub, E., Some inequalities for stop-loss premiums. Astin Bulletin 9, Dhaene, J., Sundt, B., On error bounds for approximations to aggregate claims distributions. Astin Bulletin 27 (2), De Pril, N., Dhaene, J., Error bounds for the compound Poisson approximations of the individual risk model. Astin Bulletin 22 (2), Gerber, H.U., An Introduction to Mathematical Risk Theory. S.S. Hueber Foundation Monograph No. 8. Richard D. Irwin, Homewood, IL. Gerber, H., Error bounds for the compound Poisson approximation. Insurance: Mathematics and Economics 3, Hipp, C., Approximation of aggregate claims distributions by compound Poisson approximation. Insurance: Mathematics and Economics 4, Hipp, C., Improved approximations for the aggregate claims distribution in the individual model. Astin Bulletin 16 (2), Joe, H., Multivariate Models and Dependence Concepts. Chapman & Hall, New York, pp Kaas, R., Van Heerwaarden, Goovaerts, M.J., 1988b. Between individual and collective model for the total claims. Astin Bulletin 18 (2), Kaas, R., Goovaerts, M., Dhaene, J., Denuit, M., 21. Modern Actuarial Risk Theory. Kluwer Academic Publishers, London. Kuon, S., Radtke, M., Reich, A., An approximate way to switch from the individual risk model to the collective one. Astin Bulletin 23 (1), Li, X.D., Yang, J.P., 21. Compound Poisson approximation to individual risk models. Acta Scientiarum Naturalium, Universitatis Pekinensis 37 (5), (in Chinese). Michel, R., An improved error bound for the compound Poisson approximation of a family of a nearly homogeneous portfolio. Astin Bulletin 17 (2), Panjer, H.H., Recursive evaluation of a family of compound distributions. Astin Bulletin 12 (1), Panjer, H.H., Willmot, G.E., Insurance Risk Models. Society of Actuaries, Schaumburg. Shorack, G.R., Wellner, J.A., Empirical Processes with Applications to Statistics. Wiley, New York, p. 64. Sundt, B., Improved error bounds for Bertram s method. Astin Bulletin 23 (2),
21 J. Yang et al. / Insurance: Mathematics and Economics 36 (25) Wang, S., Young, V.R., Panjer, H.H., Axiomatic characterization of insurance prices. Insurance: Mathematics and Economics 21, Further reading Hesselager, O., Order relations for some distributions. Insurance: Mathematics and Economics 16, Kaas, R., Van Heerwaarden, Goovaerts, M.J., 1988a. On stop-loss premiums for the individual risk model. Astin Bulletin 18, Sundt, B., Dhaene, J., On approximating claims distributions and stop-loss premiums by truncation. Insurance: Mathematics and Economics 1, Sundt, B., Dhaene, J., On bounds for the difference between the stop-loss transforms of the two compound distributions. Astin Bulletin 26 (2),
EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION
An. Şt. Univ. Ovidius Constanţa Vol. 92), 2001, 181 192 EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION Raluca Vernic Dedicated to Professor Mirela Ştefănescu on the occasion of her
More informationComonotonicity and Maximal Stop-Loss Premiums
Comonotonicity and Maximal Stop-Loss Premiums Jan Dhaene Shaun Wang Virginia Young Marc J. Goovaerts November 8, 1999 Abstract In this paper, we investigate the relationship between comonotonicity and
More informationApplying the proportional hazard premium calculation principle
Applying the proportional hazard premium calculation principle Maria de Lourdes Centeno and João Andrade e Silva CEMAPRE, ISEG, Technical University of Lisbon, Rua do Quelhas, 2, 12 781 Lisbon, Portugal
More informationUpper stop-loss bounds for sums of possibly dependent risks with given means and variances
Statistics & Probability Letters 57 (00) 33 4 Upper stop-loss bounds for sums of possibly dependent risks with given means and variances Christian Genest a, Etienne Marceau b, Mhamed Mesoui c a Departement
More informationCharacterization 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 informationON BOUNDS FOR THE DIFFERENCE BETWEEN THE STOP-LOSS TRANSFORMS OF TWO COMPOUND DISTRIBUTIONS
ON BOUNDS FOR THE DIFFERENCE BETWEEN THE STOP-LOSS TRANSFORMS OF TWO COMPOUND DISTRIBUTIONS BJORN SUNDT I AND JAN DHAENE 2 ABSTRACT In the present note we deduce a class of bounds for the difference between
More informationA new approach for stochastic ordering of risks
A new approach for stochastic ordering of risks Liang Hong, PhD, FSA Department of Mathematics Robert Morris University Presented at 2014 Actuarial Research Conference UC Santa Barbara July 16, 2014 Liang
More informationCalculation of Bayes Premium for Conditional Elliptical Risks
1 Calculation of Bayes Premium for Conditional Elliptical Risks Alfred Kume 1 and Enkelejd Hashorva University of Kent & University of Lausanne February 1, 13 Abstract: In this paper we discuss the calculation
More informationTwo binomial methods for evaluating the aggregate claims distribution in De Pril s individual risk model
Two binomial methods for evaluating the aggregate claims distribution in De Pril s individual risk model Sundt, Bjørn Vital Forsikring ASA P.O.Box 250, N-1326 Lysaker, Norway Phone: (+47) 67 83 44 71 Fax:
More informationRuin Probabilities of a Discrete-time Multi-risk Model
Ruin Probabilities of a Discrete-time Multi-risk Model Andrius Grigutis, Agneška Korvel, Jonas Šiaulys Faculty of Mathematics and Informatics, Vilnius University, Naugarduko 4, Vilnius LT-035, Lithuania
More informationType II Bivariate Generalized Power Series Poisson Distribution and its Applications in Risk Analysis
Type II Bivariate Generalized Power Series Poisson Distribution and its Applications in Ris Analysis KK Jose a, and Shalitha Jacob a,b a Department of Statistics, St Thomas College, Pala, Arunapuram, Kerala-686574,
More informationDiscrete Distributions Chapter 6
Discrete Distributions Chapter 6 Negative Binomial Distribution section 6.3 Consider k r, r +,... independent Bernoulli trials with probability of success in one trial being p. Let the random variable
More informationCompetitive Equilibria in a Comonotone Market
Competitive Equilibria in a Comonotone Market 1/51 Competitive Equilibria in a Comonotone Market Ruodu Wang http://sas.uwaterloo.ca/ wang Department of Statistics and Actuarial Science University of Waterloo
More informationPolynomial approximation of mutivariate aggregate claim amounts distribution
Polynomial approximation of mutivariate aggregate claim amounts distribution Applications to reinsurance P.O. Goffard Axa France - Mathematics Institute of Marseille I2M Aix-Marseille University 19 th
More informationAsymptotics of random sums of heavy-tailed negatively dependent random variables with applications
Asymptotics of random sums of heavy-tailed negatively dependent random variables with applications Remigijus Leipus (with Yang Yang, Yuebao Wang, Jonas Šiaulys) CIRM, Luminy, April 26-30, 2010 1. Preliminaries
More informationAlmost Sure Convergence of the General Jamison Weighted Sum of B-Valued Random Variables
Acta Mathematica Sinica, English Series Feb., 24, Vol.2, No., pp. 8 92 Almost Sure Convergence of the General Jamison Weighted Sum of B-Valued Random Variables Chun SU Tie Jun TONG Department of Statistics
More informationA multivariate dependence measure for aggregating risks
A multivariate dependence measure for aggregating risks Jan Dhaene 1 Daniël Linders 2 Wim Schoutens 3 David Vyncke 4 December 1, 2013 1 KU Leuven, Leuven, Belgium. Email: jan.dhaene@econ.kuleuven.be 2
More informationRemarks on quantiles and distortion risk measures
Remarks on quantiles and distortion risk measures Jan Dhaene Alexander Kukush y Daniël Linders z Qihe Tang x Version: October 7, 202 Abstract Distorted expectations can be expressed as weighted averages
More informationTail 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 informationOn Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance
Applied Mathematical Sciences, Vol. 11, 217, no. 53, 269-2629 HIKARI Ltd, www.m-hikari.com https://doi.org/1.12988/ams.217.7824 On Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance
More informationThe Credibility Estimators with Dependence Over Risks
Applied Mathematical Sciences, Vol. 8, 2014, no. 161, 8045-8050 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ams.2014.410803 The Credibility Estimators with Dependence Over Risks Qiang Zhang
More informationRuin probabilities of the Parisian type for small claims
Ruin probabilities of the Parisian type for small claims Angelos Dassios, Shanle Wu October 6, 28 Abstract In this paper, we extend the concept of ruin in risk theory to the Parisian type of ruin. For
More informationDistortion Risk Measures: Coherence and Stochastic Dominance
Distortion Risk Measures: Coherence and Stochastic Dominance Dr. Julia L. Wirch Dr. Mary R. Hardy Dept. of Actuarial Maths Dept. of Statistics and and Statistics Actuarial Science Heriot-Watt University
More informationST5215: Advanced Statistical Theory
Department of Statistics & Applied Probability Wednesday, October 19, 2011 Lecture 17: UMVUE and the first method of derivation Estimable parameters Let ϑ be a parameter in the family P. If there exists
More informationReinsurance and ruin problem: asymptotics in the case of heavy-tailed claims
Reinsurance and ruin problem: asymptotics in the case of heavy-tailed claims Serguei Foss Heriot-Watt University, Edinburgh Karlovasi, Samos 3 June, 2010 (Joint work with Tomasz Rolski and Stan Zachary
More informationCOMPARING APPROXIMATIONS FOR RISK MEASURES RELATED TO SUMS OF CORRELATED LOGNORMAL RANDOM VARIABLES MASTERARBEIT
COMPARING APPROXIMATIONS FOR RISK MEASURES RELATED TO SUMS OF CORRELATED LOGNORMAL RANDOM VARIABLES MASTERARBEIT Technische Universität Chemnitz Fakultät für Mathematik Integrierter Internationaler Master-
More informationSome results on Denault s capital allocation rule
Faculty of Economics and Applied Economics Some results on Denault s capital allocation rule Steven Vanduffel and Jan Dhaene DEPARTMENT OF ACCOUNTANCY, FINANCE AND INSURANCE (AFI) AFI 0601 Some Results
More informationA 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 informationChp 4. Expectation and Variance
Chp 4. Expectation and Variance 1 Expectation In this chapter, we will introduce two objectives to directly reflect the properties of a random variable or vector, which are the Expectation and Variance.
More informationA Comparison: Some Approximations for the Aggregate Claims Distribution
A Comparison: ome Approximations for the Aggregate Claims Distribution K Ranee Thiagarajah ABTRACT everal approximations for the distribution of aggregate claims have been proposed in the literature. In
More informationTail Conditional Expectations for. Extended Exponential Dispersion Models
Tail Conditional Expectations for Extended Exponential Dispersion Models A Thesis Presented to the Faculty of the Department of Mathematical Sciences Middle Tennessee State University In Partial Fulfillment
More informationP (A G) dp G P (A G)
First homework assignment. Due at 12:15 on 22 September 2016. Homework 1. We roll two dices. X is the result of one of them and Z the sum of the results. Find E [X Z. Homework 2. Let X be a r.v.. Assume
More informationUniversity Of Calgary Department of Mathematics and Statistics
University Of Calgary Department of Mathematics and Statistics Hawkes Seminar May 23, 2018 1 / 46 Some Problems in Insurance and Reinsurance Mohamed Badaoui Department of Electrical Engineering National
More informationChapter 3: Unbiased Estimation Lecture 22: UMVUE and the method of using a sufficient and complete statistic
Chapter 3: Unbiased Estimation Lecture 22: UMVUE and the method of using a sufficient and complete statistic Unbiased estimation Unbiased or asymptotically unbiased estimation plays an important role in
More informationApplications of axiomatic capital allocation and generalized weighted allocation
6 2010 11 ( ) Journal of East China Normal University (Natural Science) No. 6 Nov. 2010 Article ID: 1000-5641(2010)06-0146-10 Applications of axiomatic capital allocation and generalized weighted allocation
More information40.530: Statistics. Professor Chen Zehua. Singapore University of Design and Technology
Singapore University of Design and Technology Lecture 9: Hypothesis testing, uniformly most powerful tests. The Neyman-Pearson framework Let P be the family of distributions of concern. The Neyman-Pearson
More informationASSOCIATIVE n DIMENSIONAL COPULAS
K Y BERNETIKA VOLUM E 47 ( 2011), NUMBER 1, P AGES 93 99 ASSOCIATIVE n DIMENSIONAL COPULAS Andrea Stupňanová and Anna Kolesárová The associativity of n-dimensional copulas in the sense of Post is studied.
More informationMonotonicity and Aging Properties of Random Sums
Monotonicity and Aging Properties of Random Sums Jun Cai and Gordon E. Willmot Department of Statistics and Actuarial Science University of Waterloo Waterloo, Ontario Canada N2L 3G1 E-mail: jcai@uwaterloo.ca,
More informationDiscounted probabilities and ruin theory in the compound binomial model
Insurance: Mathematics and Economics 26 (2000) 239 250 Discounted probabilities and ruin theory in the compound binomial model Shixue Cheng a, Hans U. Gerber b,, Elias S.W. Shiu c,1 a School of Information,
More informationPANJER CLASS UNITED One formula for the probabilities of the Poisson, Binomial, and Negative Binomial distribution.
PANJER CLASS UNITED One formula for the probabilities of the Poisson, Binomial, and Negative Binomial distribution Michael Facler 1 Abstract. This paper gives a formula representing all discrete loss distributions
More informationACM 116: Lectures 3 4
1 ACM 116: Lectures 3 4 Joint distributions The multivariate normal distribution Conditional distributions Independent random variables Conditional distributions and Monte Carlo: Rejection sampling Variance
More informationOn a compound Markov binomial risk model with time-correlated claims
Mathematica Aeterna, Vol. 5, 2015, no. 3, 431-440 On a compound Markov binomial risk model with time-correlated claims Zhenhua Bao School of Mathematics, Liaoning Normal University, Dalian 116029, China
More informationThe finite-time Gerber-Shiu penalty function for two classes of risk processes
The finite-time Gerber-Shiu penalty function for two classes of risk processes July 10, 2014 49 th Actuarial Research Conference University of California, Santa Barbara, July 13 July 16, 2014 The finite
More informationExperience Rating in General Insurance by Credibility Estimation
Experience Rating in General Insurance by Credibility Estimation Xian Zhou Department of Applied Finance and Actuarial Studies Macquarie University, Sydney, Australia Abstract This work presents a new
More informationNON-MONOTONICITY HEIGHT OF PM FUNCTIONS ON INTERVAL. 1. Introduction
Acta Math. Univ. Comenianae Vol. LXXXVI, 2 (2017), pp. 287 297 287 NON-MONOTONICITY HEIGHT OF PM FUNCTIONS ON INTERVAL PINGPING ZHANG Abstract. Using the piecewise monotone property, we give a full description
More informationA Note On The Erlang(λ, n) Risk Process
A Note On The Erlangλ, n) Risk Process Michael Bamberger and Mario V. Wüthrich Version from February 4, 2009 Abstract We consider the Erlangλ, n) risk process with i.i.d. exponentially distributed claims
More informationSome Approximations on the Probability of Ruin and the Inverse Ruin Function
MATIMYÁS MATEMATIKA Journal of the Mathematical Society of the Philippines ISSN 115-6926 Vol. 38 Nos. 1-2 (215) pp. 43-5 Some Approximations on the Probability of Ruin and the Inverse Ruin Function Lu
More informationMOMENT CONVERGENCE RATES OF LIL FOR NEGATIVELY ASSOCIATED SEQUENCES
J. Korean Math. Soc. 47 1, No., pp. 63 75 DOI 1.4134/JKMS.1.47..63 MOMENT CONVERGENCE RATES OF LIL FOR NEGATIVELY ASSOCIATED SEQUENCES Ke-Ang Fu Li-Hua Hu Abstract. Let X n ; n 1 be a strictly stationary
More informationRisk 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 informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4. (*. Let independent variables X,..., X n have U(0, distribution. Show that for every x (0,, we have P ( X ( < x and P ( X (n > x as n. Ex. 4.2 (**. By using induction or otherwise,
More informationJackknife Euclidean Likelihood-Based Inference for Spearman s Rho
Jackknife Euclidean Likelihood-Based Inference for Spearman s Rho M. de Carvalho and F. J. Marques Abstract We discuss jackknife Euclidean likelihood-based inference methods, with a special focus on the
More informationExam P Review Sheet. for a > 0. ln(a) i=0 ari = a. (1 r) 2. (Note that the A i s form a partition)
Exam P Review Sheet log b (b x ) = x log b (y k ) = k log b (y) log b (y) = ln(y) ln(b) log b (yz) = log b (y) + log b (z) log b (y/z) = log b (y) log b (z) ln(e x ) = x e ln(y) = y for y > 0. d dx ax
More informationMinimization of the root of a quadratic functional under a system of affine equality constraints with application in portfolio management
Minimization of the root of a quadratic functional under a system of affine equality constraints with application in portfolio management Zinoviy Landsman Department of Statistics, University of Haifa.
More informationPM functions, their characteristic intervals and iterative roots
ANNALES POLONICI MATHEMATICI LXV.2(1997) PM functions, their characteristic intervals and iterative roots by Weinian Zhang (Chengdu) Abstract. The concept of characteristic interval for piecewise monotone
More informationLecture 5: Expectation
Lecture 5: Expectation 1. Expectations for random variables 1.1 Expectations for simple random variables 1.2 Expectations for bounded random variables 1.3 Expectations for general random variables 1.4
More informationMultivariate negative binomial models for insurance claim counts
Multivariate negative binomial models for insurance claim counts Peng Shi (Northern Illinois University) and Emiliano A. Valdez (University of Connecticut) 9 November 0, Montréal, Quebec Université de
More informationConvergence in Distribution
Convergence in Distribution Undergraduate version of central limit theorem: if X 1,..., X n are iid from a population with mean µ and standard deviation σ then n 1/2 ( X µ)/σ has approximately a normal
More informationConditional Tail Expectations for Multivariate Phase Type Distributions
Conditional Tail Expectations for Multivariate Phase Type Distributions Jun Cai Department of Statistics and Actuarial Science University of Waterloo Waterloo, ON N2L 3G1, Canada Telphone: 1-519-8884567,
More informationOptimal Reinsurance Strategy with Bivariate Pareto Risks
University of Wisconsin Milwaukee UWM Digital Commons Theses and Dissertations May 2014 Optimal Reinsurance Strategy with Bivariate Pareto Risks Evelyn Susanne Gaus University of Wisconsin-Milwaukee Follow
More informationRecall that in order to prove Theorem 8.8, we argued that under certain regularity conditions, the following facts are true under H 0 : 1 n
Chapter 9 Hypothesis Testing 9.1 Wald, Rao, and Likelihood Ratio Tests Suppose we wish to test H 0 : θ = θ 0 against H 1 : θ θ 0. The likelihood-based results of Chapter 8 give rise to several possible
More informationEstimating VaR in credit risk: Aggregate vs single loss distribution
Estimating VaR in credit risk: Aggregate vs single loss distribution M. Assadsolimani and D. Chetalova arxiv:172.4388v1 [q-fin.cp] 14 Feb 217 Abstract Using Monte Carlo simulation to calculate the Value
More informationMathematical models in regression credibility theory
BULETINUL ACADEMIEI DE ŞTIINŢE A REPUBLICII MOLDOVA. MATEMATICA Number 3(58), 2008, Pages 18 33 ISSN 1024 7696 Mathematical models in regression credibility theory Virginia Atanasiu Abstract. In this paper
More informationSurvival probabilities in bivariate risk models, with application to reinsurance
Survival probabilities in bivariate risk models, with application to reinsurance A. Castañer a, M.M. Claramunt a,, C. Lefèvre b a Dept. Matemàtica Econòmica, Financera i Actuarial, Universitat de Barcelona,
More informationOn the positivity of linear weights in WENO approximations. Abstract
On the positivity of linear weights in WENO approximations Yuanyuan Liu, Chi-Wang Shu and Mengping Zhang 3 Abstract High order accurate weighted essentially non-oscillatory (WENO) schemes have been used
More informationA Lower Bound Theorem. Lin Hu.
American J. of Mathematics and Sciences Vol. 3, No -1,(January 014) Copyright Mind Reader Publications ISSN No: 50-310 A Lower Bound Theorem Department of Applied Mathematics, Beijing University of Technology,
More informationChapter 7. Confidence Sets Lecture 30: Pivotal quantities and confidence sets
Chapter 7. Confidence Sets Lecture 30: Pivotal quantities and confidence sets Confidence sets X: a sample from a population P P. θ = θ(p): a functional from P to Θ R k for a fixed integer k. C(X): a confidence
More informationCase study: stochastic simulation via Rademacher bootstrap
Case study: stochastic simulation via Rademacher bootstrap Maxim Raginsky December 4, 2013 In this lecture, we will look at an application of statistical learning theory to the problem of efficient stochastic
More informationd(x n, x) d(x n, x nk ) + d(x nk, x) where we chose any fixed k > N
Problem 1. Let f : A R R have the property that for every x A, there exists ɛ > 0 such that f(t) > ɛ if t (x ɛ, x + ɛ) A. If the set A is compact, prove there exists c > 0 such that f(x) > c for all x
More informationMultivariate Measures of Positive Dependence
Int. J. Contemp. Math. Sciences, Vol. 4, 2009, no. 4, 191-200 Multivariate Measures of Positive Dependence Marta Cardin Department of Applied Mathematics University of Venice, Italy mcardin@unive.it Abstract
More informationn! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2
Order statistics Ex. 4.1 (*. Let independent variables X 1,..., X n have U(0, 1 distribution. Show that for every x (0, 1, we have P ( X (1 < x 1 and P ( X (n > x 1 as n. Ex. 4.2 (**. By using induction
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 informationLecture 11: Introduction to Markov Chains. Copyright G. Caire (Sample Lectures) 321
Lecture 11: Introduction to Markov Chains Copyright G. Caire (Sample Lectures) 321 Discrete-time random processes A sequence of RVs indexed by a variable n 2 {0, 1, 2,...} forms a discretetime random process
More information4 Invariant Statistical Decision Problems
4 Invariant Statistical Decision Problems 4.1 Invariant decision problems Let G be a group of measurable transformations from the sample space X into itself. The group operation is composition. Note that
More informationA Ruin Model with Compound Poisson Income and Dependence Between Claim Sizes and Claim Intervals
Acta Mathematicae Applicatae Sinica, English Series Vol. 3, No. 2 (25) 445 452 DOI:.7/s255-5-478- http://www.applmath.com.cn & www.springerlink.com Acta Mathema cae Applicatae Sinica, English Series The
More informationSpring 2012 Math 541B Exam 1
Spring 2012 Math 541B Exam 1 1. A sample of size n is drawn without replacement from an urn containing N balls, m of which are red and N m are black; the balls are otherwise indistinguishable. Let X denote
More informationA polynomial expansion to approximate ruin probabilities
A polynomial expansion to approximate ruin probabilities P.O. Goffard 1 X. Guerrault 2 S. Loisel 3 D. Pommerêt 4 1 Axa France - Institut de mathématiques de Luminy Université de Aix-Marseille 2 Axa France
More informationRecursive Calculation of Finite Time Ruin Probabilities Under Interest Force
Recursive Calculation of Finite Time Ruin Probabilities Under Interest Force Rui M R Cardoso and Howard R Waters Abstract In this paper we consider a classical insurance surplus process affected by a constant
More informationλ(x + 1)f g (x) > θ 0
Stat 8111 Final Exam December 16 Eleven students took the exam, the scores were 92, 78, 4 in the 5 s, 1 in the 4 s, 1 in the 3 s and 3 in the 2 s. 1. i) Let X 1, X 2,..., X n be iid each Bernoulli(θ) where
More informationResearch Article Strong Convergence Bound of the Pareto Index Estimator under Right Censoring
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 200, Article ID 20956, 8 pages doi:0.55/200/20956 Research Article Strong Convergence Bound of the Pareto Index Estimator
More informationA filled function method for finding a global minimizer on global integer optimization
Journal of Computational and Applied Mathematics 8 (2005) 200 20 www.elsevier.com/locate/cam A filled function method for finding a global minimizer on global integer optimization You-lin Shang a,b,,,
More informationProbability Transforms with Elliptical Generators
Probability Transforms with Elliptical Generators Emiliano A. Valdez, PhD, FSA, FIAA School of Actuarial Studies Faculty of Commerce and Economics University of New South Wales Sydney, AUSTRALIA Zurich,
More informationUniversity of Mannheim, West Germany
- - XI,..., - - XI,..., ARTICLES CONVERGENCE OF BAYES AND CREDIBILITY PREMIUMS BY KLAUS D. SCHMIDT University of Mannheim, West Germany ABSTRACT For a risk whose annual claim amounts are conditionally
More informationThe Degree of the Splitting Field of a Random Polynomial over a Finite Field
The Degree of the Splitting Field of a Random Polynomial over a Finite Field John D. Dixon and Daniel Panario School of Mathematics and Statistics Carleton University, Ottawa, Canada {jdixon,daniel}@math.carleton.ca
More informationOn the convergence of sequences of random variables: A primer
BCAM May 2012 1 On the convergence of sequences of random variables: A primer Armand M. Makowski ECE & ISR/HyNet University of Maryland at College Park armand@isr.umd.edu BCAM May 2012 2 A sequence a :
More informationOn 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 informationBounds for Sums of Non-Independent Log-Elliptical Random Variables
Bounds for Sums of Non-Independent Log-Elliptical Random Variables Emiliano A. Valdez The University of New South Wales Jan Dhaene Katholieke Universiteit Leuven May 5, 003 ABSTRACT In this paper, we construct
More information1 + lim. n n+1. f(x) = x + 1, x 1. and we check that f is increasing, instead. Using the quotient rule, we easily find that. 1 (x + 1) 1 x (x + 1) 2 =
Chapter 5 Sequences and series 5. Sequences Definition 5. (Sequence). A sequence is a function which is defined on the set N of natural numbers. Since such a function is uniquely determined by its values
More informationX n D X lim n F n (x) = F (x) for all x C F. lim n F n(u) = F (u) for all u C F. (2)
14:17 11/16/2 TOPIC. Convergence in distribution and related notions. This section studies the notion of the so-called convergence in distribution of real random variables. This is the kind of convergence
More informationDEP ARTEMENT TOEGEP ASTE ECONOMISCHE WETENSCHAPPEN
DEP ARTEMENT TOEGEP ASTE ECONOMISCHE WETENSCHAPPEN ONDERZOEKSRAPPORT NR 9513 Some Moment Relations for the Hipp Approximation by JanDHAENE Bj0mSUNDT Nelson DE PRIL Katholieke Universiteit Leuven Naamsestraat
More informationAsymptotic 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 informationOn the distance signless Laplacian spectral radius of graphs and digraphs
Electronic Journal of Linear Algebra Volume 3 Volume 3 (017) Article 3 017 On the distance signless Laplacian spectral radius of graphs and digraphs Dan Li Xinjiang University,Urumqi, ldxjedu@163.com Guoping
More information18Ï È² 7( &: ÄuANOVAp.O`û5 571 Based on this ANOVA model representation, Sobol (1993) proposed global sensitivity index, S i1...i s = D i1...i s /D, w
A^VÇÚO 1 Êò 18Ï 2013c12 Chinese Journal of Applied Probability and Statistics Vol.29 No.6 Dec. 2013 Optimal Properties of Orthogonal Arrays Based on ANOVA High-Dimensional Model Representation Chen Xueping
More informationMAT 257, Handout 13: December 5-7, 2011.
MAT 257, Handout 13: December 5-7, 2011. The Change of Variables Theorem. In these notes, I try to make more explicit some parts of Spivak s proof of the Change of Variable Theorem, and to supply most
More informationSmooth simultaneous confidence bands for cumulative distribution functions
Journal of Nonparametric Statistics, 2013 Vol. 25, No. 2, 395 407, http://dx.doi.org/10.1080/10485252.2012.759219 Smooth simultaneous confidence bands for cumulative distribution functions Jiangyan Wang
More informationDoubly Indexed Infinite Series
The Islamic University of Gaza Deanery of Higher studies Faculty of Science Department of Mathematics Doubly Indexed Infinite Series Presented By Ahed Khaleel Abu ALees Supervisor Professor Eissa D. Habil
More informationLecture Notes on Risk Theory
Lecture Notes on Risk Theory February 2, 21 Contents 1 Introduction and basic definitions 1 2 Accumulated claims in a fixed time interval 3 3 Reinsurance 7 4 Risk processes in discrete time 1 5 The Adjustment
More information3 Measurable Functions
3 Measurable Functions Notation A pair (X, F) where F is a σ-field of subsets of X is a measurable space. If µ is a measure on F then (X, F, µ) is a measure space. If µ(x) < then (X, F, µ) is a probability
More informationCounts using Jitters joint work with Peng Shi, Northern Illinois University
of Claim Longitudinal of Claim joint work with Peng Shi, Northern Illinois University UConn Actuarial Science Seminar 2 December 2011 Department of Mathematics University of Connecticut Storrs, Connecticut,
More informationStatistics Ph.D. Qualifying Exam: Part II November 9, 2002
Statistics Ph.D. Qualifying Exam: Part II November 9, 2002 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. 1 2 3 4 5 6 7 8 9 10 11 12 2. Write your
More informationAsymptotic of Enumerative Invariants in CP 2
Peking Mathematical Journal https://doi.org/.7/s4543-8-4-4 ORIGINAL ARTICLE Asymptotic of Enumerative Invariants in CP Gang Tian Dongyi Wei Received: 8 March 8 / Revised: 3 July 8 / Accepted: 5 July 8
More information