ON BOUNDS FOR THE DIFFERENCE BETWEEN THE STOP-LOSS TRANSFORMS OF TWO COMPOUND DISTRIBUTIONS

Similar documents
DEP ARTEMENT TOEGEP ASTE ECONOMISCHE WETENSCHAPPEN

Comonotonicity and Maximal Stop-Loss Premiums

Two binomial methods for evaluating the aggregate claims distribution in De Pril s individual risk model

EVALUATING THE BIVARIATE COMPOUND GENERALIZED POISSON DISTRIBUTION

Mathematical models in regression credibility theory

The compound Poisson random variable s approximation to the individual risk model

Recursive Calculation of Finite Time Ruin Probabilities Under Interest Force

Probability Transforms with Elliptical Generators

Remarks on quantiles and distortion risk measures

Characterization of Upper Comonotonicity via Tail Convex Order

Asymptotics of random sums of heavy-tailed negatively dependent random variables with applications

P i [B k ] = lim. n=1 p(n) ii <. n=1. V i :=

Ruin Probabilities of a Discrete-time Multi-risk Model

LECTURE 15, WEDNESDAY

Estimating VaR in credit risk: Aggregate vs single loss distribution

Discounted probabilities and ruin theory in the compound binomial model

Tail Mutual Exclusivity and Tail- Var Lower Bounds

Discrete Distributions Chapter 6

The Ruin Probability of a Discrete Time Risk Model under Constant Interest Rate with Heavy Tails

MONOTONICITY OF RATIOS INVOLVING INCOMPLETE GAMMA FUNCTIONS WITH ACTUARIAL APPLICATIONS

Applications of axiomatic capital allocation and generalized weighted allocation

a 11 x 1 + a 12 x a 1n x n = b 1 a 21 x 1 + a 22 x a 2n x n = b 2.

On a compound Markov binomial risk model with time-correlated claims

On a Class of Approximative Computation Methods in the Individual Risk Model

Math 110 (S & E) Textbook: Calculus Early Transcendentals by James Stewart, 7 th Edition

ELEMENTARY LINEAR ALGEBRA

THREE METHODS TO CALCULATE THE PROBABILITY OF RUIN. University of Lausanne, Switzerland

A new approach for stochastic ordering of risks

Some Approximations on the Probability of Ruin and the Inverse Ruin Function

417 P. ERDÖS many positive integers n for which fln) is l-th power free, i.e. fl n) is not divisible by any integral l-th power greater than 1. In fac

Monotonicity and Aging Properties of Random Sums

A multivariate dependence measure for aggregating risks

THE DENOMINATORS OF POWER SUMS OF ARITHMETIC PROGRESSIONS. Bernd C. Kellner Göppert Weg 5, Göttingen, Germany

Additional Practice Lessons 2.02 and 2.03

Finding Divisors, Number of Divisors and Summation of Divisors. Jane Alam Jan

A Gentle Introduction to Gradient Boosting. Cheng Li College of Computer and Information Science Northeastern University

A UNIFIED APPROACH TO GENERATE RISK MEASURES ABSTRACT KEYWORDS. Insurance premium principle, Risk measure, Markov inequality 1.

Upper stop-loss bounds for sums of possibly dependent risks with given means and variances

On Finite-Time Ruin Probabilities for Classical Risk Models

Point sets and certain classes of sets

CHAPTER 8: EXPLORING R

CSC 344 Algorithms and Complexity. Proof by Mathematical Induction

Links Between Sums Over Paths in Bernoulli s Triangles and the Fibonacci Numbers

Probability Generating Functions

A Comparison: Some Approximations for the Aggregate Claims Distribution

Tail Conditional Expectations for Extended Exponential Dispersion Models

CHAPTER 4 THE RIEMANN HYPOTHESIS FOR CURVES OVER FINITE FIELDS

An estimate for the probability of dependent events

Nonlife Actuarial Models. Chapter 5 Ruin Theory

Applying the proportional hazard premium calculation principle

SCTT The pqr-method august 2016

ELEMENTARY LINEAR ALGEBRA

A New Shuffle Convolution for Multiple Zeta Values

Poisson Processes. Stochastic Processes. Feb UC3M

Katholieke Universiteit Leuven, Belgium. University of Bergen, Norway

FIBONACCI NUMBERS AND DECIMATION OF BINARY SEQUENCES

A COMPOUND POISSON APPROXIMATION INEQUALITY

A matrix over a field F is a rectangular array of elements from F. The symbol

P. ERDÖS 2. We need several lemmas or the proo o the upper bound in (1). LEMMA 1. { d ( (k)) r 2 < x (log x)c5. This result is due to van der Corput*.

A well-quasi-order for tournaments

On Finite-Time Ruin Probabilities in a Risk Model Under Quota Share Reinsurance

Remarks on quantiles and distortion risk measures

Binary Convolutional Codes of High Rate Øyvind Ytrehus

Some Congruences for the Partial Bell Polynomials

CSI Mathematical Induction. Many statements assert that a property of the form P(n) is true for all integers n.

Bounds for Sums of Non-Independent Log-Elliptical Random Variables

Definitions, Theorems and Exercises. Abstract Algebra Math 332. Ethan D. Bloch

arxiv: v2 [math.pr] 22 Apr 2014 Dept. of Mathematics & Computer Science Mathematical Institute

4 Sums of Independent Random Variables

HIGHER-ORDER DIFFERENCES AND HIGHER-ORDER PARTIAL SUMS OF EULER S PARTITION FUNCTION

Non Uniform Bounds on Geometric Approximation Via Stein s Method and w-functions

Quadratic Forms of Skew Schur Functions

MATH 56A: STOCHASTIC PROCESSES CHAPTER 2

Technical Report No. 13/04, December 2004 INTRODUCING A DEPENDENCE STRUCTURE TO THE OCCURRENCES IN STUDYING PRECISE LARGE DEVIATIONS FOR THE TOTAL

Sequences and Series. College Algebra

ELEMENTARY LINEAR ALGEBRA

GENERALIZED ANNUITIES AND ASSURANCES, AND INTER-RELATIONSHIPS. BY LEIGH ROBERTS, M.Sc., ABSTRACT

EXTENSIONS OF KATZ PANJER FAMILIES OF DISCRETE DISTRIBUTIONS

Ruin Theory. A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science at George Mason University

Non-Life Insurance: Mathematics and Statistics

Physics 342 Lecture 23. Radial Separation. Lecture 23. Physics 342 Quantum Mechanics I

arxiv:math/ v1 [math.st] 9 Nov 2005

The Littlewood-Richardson Rule

EE376A: Homework #3 Due by 11:59pm Saturday, February 10th, 2018

n! (k 1)!(n k)! = F (X) U(0, 1). (x, y) = n(n 1) ( F (y) F (x) ) n 2

Random Variables. P(x) = P[X(e)] = P(e). (1)

SOLUTIONS TO PROBLEMS

A New Characterization of Boolean Rings with Identity

ON PARTITION FUNCTIONS OF ANDREWS AND STANLEY

A random variable is a quantity whose value is determined by the outcome of an experiment.

Upper Bounds for the Ruin Probability in a Risk Model With Interest WhosePremiumsDependonthe Backward Recurrence Time Process

Chapter 5: Integer Compositions and Partitions and Set Partitions

Math Treibergs. Solutions to Fourth Homework February 13, 2009

CS 6820 Fall 2014 Lectures, October 3-20, 2014

Characterizing Cycle Partition in 2-Row Bulgarian Solitaire

1: PROBABILITY REVIEW

Alternative Characterization of Ergodicity for Doubly Stochastic Chains

Estimators for the binomial distribution that dominate the MLE in terms of Kullback Leibler risk

Seunghee Ye Ma 8: Week 2 Oct 6

Factorization in Integral Domains II

Transcription:

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 the stoploss transforms of two compound distributions with the same severity distribution. The class contains bounds of any degree of accuracy in the sense that the bounds can be chosen as close to the exact value as desired; the time required to compute the bounds increases with the accuracy. 1. INTRODUCTION During the last twenty years, there has grown up a large literature on approximations and inequalities for stop-loss premiums under various assumptions. One way of approximation is to approximate the original distribution with another distribution that makes the evaluation simpler. In such cases it is useful to have bounds for the difference between the exact stop-loss premium and the approximation, that is, it is of interest to have bounds for the difference between the stop-loss transforms of two distributions. When approximating the stop-loss transform of a compound distribution, it is sometimes convenient to replace the counting distribution with another distribution, e.g. a Bernoulli distribution or a Poisson distribution, and keep the severity distribution unchanged. Such approximations are discussed by i.a. Dhaene & Sundt (1996). In the present note we deduce classes of bounds for the difference between the stoploss transforms of two compound distributions with the same severity distribution. The classes contain bounds of any degree of accuracy in the sense that the bounds can be chosen as close to the exact value as desired; the time required to compute the bounds increases with the accuracy. 2. NOTATION AND CONVENTIONS Let ~ and Z denote the sets of respectively the non-negative real numbers and the non-negative integers, and Pn and it::'z+ the classes of probability distributions with t University of Bergen, Norway 2 Katholieke Universiteit Leuven, Belgium ASTIN BULLETIN, Vol 26, No. 2. 1996, pp. 225-231

226 BJORN SUNDT - JAN DHAENE finite mean on respectively ~+ and &. For distributions in Re+, we shall denote the cumulative distribution function with a capital letter and, for distributions in Pz, the discrete density function with the corresponding lower case letter. The stop-loss transform of a distribution will be denoted by a horizontal bar on the top of the symbol of the distribution, that is, for a distribution F E Pe+, we have the stop-loss transform F given by F(x)=f?(y-x)dF(y)=f?(1-F(y))dy. (x_>o) The mean of F is denoted by #~-, that is, PF = F(O)= ~ ydf(y)=j'~ (1- F(y))dy. We shall denote a compound distribution with counting distribution P E Pz and severity distribution H E Pe by P v H, that is, p v H = n=0 p(.)h"', where/-/"*denotes the n-fold convolution of H. For FE PR+ and re K'+, we define the approximation Fr by Fr(x)=(lF.(x) (x(o<-x<r)_>r) This approximation can be interpreted as the distribution obtained by setting all observations greater than r equal to r. The limiting cases r = 0 and r = ~ correspond to respectively the distribution concentrated in zero and the original distribution F. We shall interpret Z i=a b I)i = 0 when b < a. 3. RESULTS 3,1. The following lemma is proved as formula (38) in De Pril & Dhaene (1992) for the special case r = I ; the proof is easily extended to the general case. Lemma 1 For H E P~+, r, m E Z such that r _<m, and x E 7~+, we have (m - r)h (x) _< H'*(x)- Hr*(x) _<(m- r)12 H. Lemma 2 For P E Pz, HE P~+, re Z+, and x E 7~, we have H (x)p (r) _< P v H (x)- Prv H(x) _< #H P(r). (1)

DIFFERENCE STOP-LOSS ORDER AND TWO COMPOUND DISTRIBUTIONS 227 Proof. We have Application of Lemma 1 gives oo e v n(x)- ~ v 14(x) = y~ p(n) (n"'(x)-/4"(x)). n=r n=r p(~) (n-r)n(x)_< ev H(x)-er v n(x)~_ ~p(n) (n -,')~., liar from which we obtain (1). Q.E.D The second inequality in (1) was proved under more general assumptions by Sundt (1991), who also showed that 0 _< Pv H(x)- Prv H(x), which is weaker than the first inequality in (1). If Pr = P, that is, P(r) = 1, then the bounds in (I) become equal to zero. Lemma 1 appears as a special case of Lemma 2 by letting P be the distribution concentrated in m. 3.2. For P, Q ~z+, H ~P~+, r Z+, and x "R+, we introduce B r (x;p,q,h)= Prv H(x)-Qr v H(x)+PnP(r )- H(x)Q(r), (2) which can also be written as r-i O r (X; P, Q, H)= Z(p(n)-q(nl)H"*(x)- n=l ( P(r - 1) - Q(r - 1))H r (x) +/'/n P(r ) - H (x ) Q (r ). (3) Theorem 1 For P, Q Pz+, H PR+, r Z, and x 7~, we have -Br(x;Q,P,H) _< P v H(x)-a v H(x) _< Br(X;P,Q,H ). (4) Proof. Application of Lemma 2 gives P v H(x) - O v H(x) _< Pr V H(x) + pt.ip(r) - Or v H(x) - H(x)O(r) = Br(X;P,Q,H), which proves the second inequality in (4). The first inequality follows by symmetry. This completes the proof of Theorem 1. Q.E.D. We shall look at some special cases of Theorem 1: 1. As Br(x;P, Pr,H)=I-IHP(r) Br(X;Pr,P,H)=-H(x)P(r), we see that Lemma 2 (and thus also Lemma I) is a special case of Theorem 1.

228 BJORN SUNDT-JAN DHAENE 2. From (3) we obtain B, (x;p,a, H) = -(p(0)- q(0)) H(x) +/1H P(1)- H(x)O(1) = For H ~ Pz+ this case is discussed in Dhaene & Sundt (1996). 3. From (2) we obtain B 0 (x; P,Q, H) = ~HI.te - H(x)lJa. (5) 4. If P(r) = Q(r) = 1, then Pr = P and Qr = Q, and from (2) we obtain B~ (x; P,a, H) = -B~ (x;a, P, H) = P v H(x)- a v H(x), that is, in this case Theorem 1 becomes trivial. 5. From (2) we obtain Br(x;P,P,H)=(II H - H(x))P(r), (6) that is, unfortunately the bounds in Theorem 1 do not in general become equal to zero when comparing two identical compound distributions. 3,3. Let D/x,'P,Q,H) denote the difference between the upper and lower bound in Theorem 1, that is, Dr(x;P,Q,H)= Br(x;P,Q,H)+ B~(x;Q,P,H). (7) Then Or(x ~ P, Q, H): (]'/H -- "n(x))('p(~')'~ ~)(F)). (8) We see that Dr(x;P,Q,H) decreases to zero when rincreases to infinity, that is, we can make the difference between the upper and lower bound in Theorem 1 as small as desired by making r sufficiently large. We see that Or(X;P,Q,I-I ) increases from zero to I.tH(-P(r)+O.(r)) when x increases from zero to infinity. Thus our bounds are most accurate for low values of x. Furthermore, if for some e > 0 we choose r such that then Dr(x:P,Q,H) <e for all x E ;P+. 3.4. Let _ m E P(r)+Q(r) < -- lzn br (x; P,Q, H) = B r (x; P, Q, H) - Br+ 1 (x; P, Q, H).

DIFFERENCE STOP-LOSS ORDER AND TWO COMPOUND DISTRIBUTIONS 229 From (3) and trivial calculus we obtain br(x;p,q,n)=(p(r)-q(r))(h(r+l)*(x) - Hr*(x)l+ /2H(1- P(r))- H(x)(1-O(r)). (9) By rewriting (9) as (i-q(r))ih(r+i'*(x) - Hr*(x)- n(x)l and application of Lemma 1, we see that br(x;p,q,h) _> O. Thus Br(x;P,Q,H ) is nonincreasing in r. This implies that in (4), the upper bound is non-increasing and the lower bound is non-decreasing in r, and as Dr(x;P,Q,H) goes to zero when rgoes to infinity, both bounds converge to P v H(x) - Q v H(x). Formula (9) can be applied for recursive evaluation of &(x;p,q.h). Furthermore, when we have found Br(x;P,Q,H), we easily obtainbr(x;q,p,h) from (7) and (8). 3.5. The main purpose of the present subsection is to deduce an improvement of the bounds in Theorem 1. For doing that, we shall need the following lemma. Lemma 3 For P, Q ~ Pz+, H ~ ~+, r ~ Z and x ~ R+, we have br(x;p,q,h)_>(~ H -H(x))(l-max(P(r),Q(r))). (10) Proof. We apply Lemma 1 in (9). If P(r) _> Q(r), then b r (x; P,Q, H) _>(P(r) - Q(r))H(x) + IZn (1 - P(r)) - H(x)(1 - Q(r)), that is, br(x;p,q,h )_>(u H - H(x))(I- P(r)). (I I) Analogously, if P(r) < Q(r), then b, (x; P, Q, H) _>( P(r ) - Qfr )) la H + #H (l - P(r )) - H(x)(1 - e(r )), that is, b,(x; P,Q,H)_> (#h - H(x))(l-Q(r)). (12) From (I I) and (12) we obtain (I0). Q.E.D.

230 BJORNSUNDT-JANDHAENE Theorem 2 For P, Q E Pz, H ~ PR+, r ~ Z, and x E 7;?+, we have Br(X;Q,P,H)+(p H -H(x))~(1-max(P(k),Q(k)))_< Pv H(x)-Ov H(x)< k=r nr(x;e,a,")-(]l H - H(x))~(1-max(P(k ),Q(k ))). k=r (13) Proof. For s ~ Z such that s _> r, we obtain by applying successively Theorem 1 and Lemma 3 P v H(x) - Q v H(x) _< Bs(x; P, Q, H) = s-i Br(x; P, Q, H) - Z bk (x; P, Q, H) _< k=r s-i Br(X;P,Q,H)-(I.tH - H(x))Z(I-max(P(k),Q(k)) ). By letting s go to infinity, we obtain the second inequality in (13); the first inequality follows by symmetry. This completes the proof of Theorem 2. Q.E.D. We see that in (13), like in (4), the lower bound is non-decreasing in rand the upper bound is non-increasing in r. The infinite summation in (13) may seem complicated. However, as all the terms are non-negative, we obtain weaker bounds by including only a finite number of terms. Furthermore, if P or Q has a finite support, then only a finite number of terms are non-zero. In the following corollary we consider another case where the summation obtains a particularly simple form. Corollary 1 Let P, Q E Pz+, H ~ Pe+, and x E 7~+. If there exists a non-negative integer s (possibly equal to infinity) such that then Q(y)<P(y) k=r (y= 0,1... s- 1)'~ t- Q(y)_> P(y), (y= s,s + l... ) J (14) - 8/x;Q,P,H)+(u. _< Pv H(x)-Qv H(x) < Be(x; P, Q, H) - (,tt H - H(x))('P(r) - "P(s)+ Q(s)) (r=0,1... s-i) (15)

DIFFERENCE STOP-LOSS ORDER AND TWO COMPOUND DISTRIBUTIONS 231 -B,. (x;q, P, H) + (11 H - H(x ))Q(r) _< P v H(x ) - Q v H(x ) <_ Br(x; P,Q,H)-(I11 t ---H(x))Q)(r). (r = s,s + 1... ) (16) Proof. For r = 0, 1... s-1, we have '~ (1 - max(p(k),q(k)))= ~_~(1- P(k))+ ~ (1- Q(k ))= k=r k=r k=s 7(r)-- 7(s) + 0(s), and insertion in (13) gives (15). The inequalities (16) are proved analogously. This completes the proof of Corollary 1. Q.E.D. If we in addition to (14) assume that/.to <,up, then we have the stop-loss orderings Q < P and Q v H < P v H; for proofs cf. e.g. Goovaerts, Kaas, van Heerwaarden, & Bauwelinckx (1990). At the end of subsection 3.2. we pointed out that unfortunately the bounds in Theorem 1 do not become equal to zero when Q = P. From (6) and (15) we see that the improved bounds of Theorem 2 do not have this deficiency. ACKNOWLEDGEMENT The authors are grateful to an anonymous referee for useful suggestions. REFERENCES DE PRIL, N. & DHAENE, J. (1992). Error bounds for compound Poisson approximations to the individual risk model. ASTIN Bulletin 22, 135-148. DHAENE, J. ~.. SUNDT, B. (1996). On error bounds for approximations to aggregate claims distributions. Submitted for publication in ASTIN Bulletin. GOOVAERTS M.J., KAAS, R., VAN HEERWAARDEN, A.E., & BAUWEL[NCKX, T. (1990). Effective Actuarial Methods. Norlh-Holland, Amsterdam. StmoT, B. (1991). On approximating aggregate claims distributions and stop-loss premiums by truncation. Insurance: Mathematics and Economics I0, 133-136. BJORN SUNDT Department of Mathematics University of Bergen Alldgaten 55 N-5007 Bergen Norway JAN DHAENE Departement Toegepaste Economische Wetenschappen Katholieke Universiteit Leuven Huis Eighen Heerd Minderbroederstraat 5 B-3000 Leuven Belgium