Random convolution of O-exponential distributions

Similar documents
On the random max-closure for heavy-tailed random variables

Randomly Weighted Sums of Conditionnally Dependent Random Variables

Discrete uniform limit law for additive functions on shifted primes

The Subexponential Product Convolution of Two Weibull-type Distributions

Randomly weighted sums under a wide type of dependence structure with application to conditional tail expectation

Asymptotic behavior for sums of non-identically distributed random variables

Weak max-sum equivalence for dependent heavy-tailed random variables

Ruin Probabilities of a Discrete-time Multi-risk Model

Journal of Mathematical Analysis and Applications

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

Asymptotic Tail Probabilities of Sums of Dependent Subexponential Random Variables

Finite-time Ruin Probability of Renewal Model with Risky Investment and Subexponential Claims

On an uniqueness theorem for characteristic functions

On the modification of the universality of the Hurwitz zeta-function

5. Some theorems on continuous functions

First In-Class Exam Solutions Math 410, Professor David Levermore Monday, 1 October 2018

ON THE LIMIT POINTS OF THE FRACTIONAL PARTS OF POWERS OF PISOT NUMBERS

Precise Large Deviations for Sums of Negatively Dependent Random Variables with Common Long-Tailed Distributions

Proof. We indicate by α, β (finite or not) the end-points of I and call

LIMIT THEOREMS FOR SHIFT SELFSIMILAR ADDITIVE RANDOM SEQUENCES

Large deviations for weighted random sums

On lower limits and equivalences for distribution tails of randomly stopped sums 1

arxiv: v3 [math.pr] 27 Sep 2016

MATH FINAL EXAM REVIEW HINTS

Rare Events in Random Walks and Queueing Networks in the Presence of Heavy-Tailed Distributions

Research Article Extended Precise Large Deviations of Random Sums in the Presence of END Structure and Consistent Variation

arxiv: v1 [math.pr] 28 Nov 2007

Heavy tails of a Lévy process and its maximum over a random time interval

Scandinavian Actuarial Journal. Asymptotics for ruin probabilities in a discrete-time risk model with dependent financial and insurance risks

C.7. Numerical series. Pag. 147 Proof of the converging criteria for series. Theorem 5.29 (Comparison test) Let a k and b k be positive-term series

ON INTEGERS EXPRESSIBLE BY SOME SPECIAL LINEAR FORM. 1. Introduction

ON THE COMPOUND POISSON DISTRIBUTION

MATH 117 LECTURE NOTES

Homework #2 Solutions Due: September 5, for all n N n 3 = n2 (n + 1) 2 4

AN INEQUALITY FOR TAIL PROBABILITIES OF MARTINGALES WITH BOUNDED DIFFERENCES

Relative growth of the partial sums of certain random Fibonacci-like sequences

Real Analysis Math 131AH Rudin, Chapter #1. Dominique Abdi

An estimate for the probability of dependent events

Approximation of Heavy-tailed distributions via infinite dimensional phase type distributions

Octavio Arizmendi, Ole E. Barndorff-Nielsen and Víctor Pérez-Abreu

Climbing an Infinite Ladder

WEIGHTED SUMS OF SUBEXPONENTIAL RANDOM VARIABLES AND THEIR MAXIMA

Characterization through Hazard Rate of heavy tailed distributions and some Convolution Closure Properties

On the Borel-Cantelli Lemma

Lectures 16-17: Poisson Approximation. Using Lemma (2.4.3) with θ = 1 and then Lemma (2.4.4), which is valid when max m p n,m 1/2, we have

Study of some equivalence classes of primes

Wald for non-stopping times: The rewards of impatient prophets

On Some Combinations of Non-Consecutive Terms of a Recurrence Sequence

EQUICONTINUITY AND SINGULARITIES OF FAMILIES OF MONOMIAL MAPPINGS

On the Bennett-Hoeffding inequality

Global Attractivity of a Higher-Order Nonlinear Difference Equation

Weighted Sums of Orthogonal Polynomials Related to Birth-Death Processes with Killing

MATH 409 Advanced Calculus I Lecture 7: Monotone sequences. The Bolzano-Weierstrass theorem.

ON THE LINEAR INDEPENDENCE OF THE SET OF DIRICHLET EXPONENTS

One-sided power sum and cosine inequalities

arxiv: v2 [math.pr] 4 Sep 2017

GAUSSIAN MEASURE OF SECTIONS OF DILATES AND TRANSLATIONS OF CONVEX BODIES. 2π) n

THE SET OF RECURRENT POINTS OF A CONTINUOUS SELF-MAP ON AN INTERVAL AND STRONG CHAOS

Oscillation Criteria for Certain nth Order Differential Equations with Deviating Arguments

Boundary value problem with integral condition for a Blasius type equation

Semi-self-decomposable distributions on Z +

Ruin Probabilities in a Dependent Discrete-Time Risk Model With Gamma-Like Tailed Insurance Risks

1/12/05: sec 3.1 and my article: How good is the Lebesgue measure?, Math. Intelligencer 11(2) (1989),

The main results about probability measures are the following two facts:

On the Local Profiles of Trees

3 Integration and Expectation

Precise Estimates for the Ruin Probability in Finite Horizon in a Discrete-time Model with Heavy-tailed Insurance and Financial Risks

SPACES ENDOWED WITH A GRAPH AND APPLICATIONS. Mina Dinarvand. 1. Introduction

Brownian Motion. 1 Definition Brownian Motion Wiener measure... 3

Practical approaches to the estimation of the ruin probability in a risk model with additional funds

A Note on the 2 F 1 Hypergeometric Function

Drazin Invertibility of Product and Difference of Idempotents in a Ring

µ (X) := inf l(i k ) where X k=1 I k, I k an open interval Notice that is a map from subsets of R to non-negative number together with infinity

,... We would like to compare this with the sequence y n = 1 n

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

Structure of R. Chapter Algebraic and Order Properties of R

MAJORIZING MEASURES WITHOUT MEASURES. By Michel Talagrand URA 754 AU CNRS

where A is a symmetric nonnegative-definite d d matrix, γ R d and ν is a measure on R d (called the Lévy measure) satisfying

MATH 102 INTRODUCTION TO MATHEMATICAL ANALYSIS. 1. Some Fundamentals

Nonhomogeneous linear differential polynomials generated by solutions of complex differential equations in the unit disc

A derivative-free nonmonotone line search and its application to the spectral residual method

On Approximation of Analytic Functions by Periodic Hurwitz Zeta-Functions

THEOREMS, ETC., FOR MATH 515

Nonnegative k-sums, fractional covers, and probability of small deviations

Series of Error Terms for Rational Approximations of Irrational Numbers

SOME FIXED POINT RESULTS FOR ADMISSIBLE GERAGHTY CONTRACTION TYPE MAPPINGS IN FUZZY METRIC SPACES

2 Lebesgue integration

THIELE CENTRE. Markov Dependence in Renewal Equations and Random Sums with Heavy Tails. Søren Asmussen and Julie Thøgersen

Cullen Numbers in Binary Recurrent Sequences

Mi Ryeong Lee and Hye Yeong Yun

Characterizations on Heavy-tailed Distributions by Means of Hazard Rate

Course 212: Academic Year Section 1: Metric Spaces

On the Bennett-Hoeffding inequality

Mathematische Annalen

Combinatorial Proof of the Hot Spot Theorem

TAIL ASYMPTOTICS FOR A RANDOM SIGN LINDLEY RECURSION

Copyright 2010 Pearson Education, Inc. Publishing as Prentice Hall.

Essential Background for Real Analysis I (MATH 5210)

CHARACTERIZATION OF NONCORRELATED PATTERN SEQUENCES AND CORRELATION DIMENSIONS. Yu Zheng. Li Peng. Teturo Kamae. (Communicated by Xiangdong Ye)

METRIC HEIGHTS ON AN ABELIAN GROUP

Transcription:

Nonlinear Analysis: Modelling and Control, Vol. 20, No. 3, 447 454 ISSN 1392-5113 http://dx.doi.org/10.15388/na.2015.3.9 Random convolution of O-exponential distributions Svetlana Danilenko a, Jonas Šiaulys b a Faculty of Fundamental Sciences, Vilnius Gediminas Technical University Saulėtekio ave. 11, LT-10223 Vilnius, Lithuania svetlana.danilenko@vgtu.lt b Faculty of Mathematics and Informatics, Vilnius University Naugarduko str. 24, LT-03225 Vilnius, Lithuania jonas.siaulys@mif.vu.lt Received: October 15, 2014 / Revised: February 12, 2015 / Published online: March 31, 2015 Abstract. Assume that ξ 1, ξ 2,... are independent and identically distributed non-negative random variables having the O-exponential distribution. Suppose that η is a nonnegative non-degenerate at zero integer-valued random variable independent of ξ 1, ξ 2,.... In this paper, we consider the conditions for η under which the distribution of random sum ξ 1 + ξ 2 + + ξ η remains in the class of O-exponential distributions. Keywords: long tail, random sum, closure property, O-exponential distribution. 1 Introduction Let ξ 1, ξ 2,... be independent copies of a random variable (r.v.) ξ with distribution function (d.f.) F ξ. Let η be a nonnegative non-degenerate at zero integer-valued r.v. independent of {ξ 1, ξ 2,... }. We pose that F ξ is O-exponential and we find minimal conditions under which the d.f. F Sη (x) := P(ξ 1 + ξ 2 + + ξ η x) = P(η = n)p(ξ 1 + ξ 2 + + ξ n x) = n=0 P(η = n)fξ n (x) n=0 belongs to the class of O-exponential distributions as well. Here and elsewhere in this paper, F n denotes the n-fold convolution of d.f. F. Theorem 1 below is the main result of this paper. Before the exact formulation of this theorem, we recall the definition of O-exponential and some related d.f. s classes. In all definitions below, we assume that F (x) = 1 F (x) > 0 for all x R. The authors are ported by a grant (No. MIP-13079) from the Research Council of Lithuania. c Vilnius University, 2015

448 S. Danilenko, J. Šiaulys Definition 1. For γ > 0, by L(γ) we denote the class of exponential d.f.s, i.e. F L(γ) if for any fixed real y, F (x + y) lim = e γy. F (x) In the case γ = 0, class L(0) is called the long-tailed distribution class and is denoted by L. Definition 2. A d.f. F belongs to the dominated varying-tailed class (F D) if for any fixed y (0, 1), F (xy) F (x) <. Definition 3. A d.f. F is O-exponential (F OL) if for any fixed y R, 0 < lim inf F (x + y) F (x) F (x + y) F (x) <. It is easy to see that the following inclusions hold: D OL, L OL, L(γ) OL. In [2, 3], Cline claimed that d.f. F Sη remains in the class L(γ) if F ξ L(γ) and η is any nonnegative non-degenerate at zero integer-valued r.v. Albin [1] observed that Cline s result is false in general. He obtained that d.f. F Sη remains in the class L(γ) if F ξ belongs to the class L(γ) and Ee δη < for each δ > 0. In order to prove this claim, author used the upper estimate F n (x t) (1 + ε)e γt, (1) F n (x) provided that ε > 0, t R, F L(γ), x n(c 1 t) + t and c 1 = c 1 (ε, t) is sufficiently large such that F (x t) (1 + ε)e γt F (x) for x c 1 (see [1, Lemma 1]). Unfortunately, the obtained estimate holds for positive t only. If t is negative, then the above estimate is incorrect in general. This fact was shown by Watanabe and Yamamuro (see [8, Remark 6.1]). Thus, the Cline proposition that P(ξ 1 + ξ 2 + + ξ η x) belongs to the class L(γ) remains not proved. In this paper, we investigate a wider class, OL, instead of the class L(γ). We show that the d.f. of the sum ξ 1 + ξ 2 + + ξ η remains in the class OL, if r.v. η satisfies the conditions similar to that in [1]. The following theorem is the main statement in this paper. Theorem 1. Let ξ 1, ξ 2,... be independent copies of a nonnegative r.v. ξ with d.f. F ξ. Let η be a nonnegative, non-degenerate at zero, integer-valued and independent of {ξ 1, ξ 2,...} r.v. with d.f. F η. If F ξ belongs to the class OL and F η (δx) = O( xf ξ (x)) for each δ (0, 1), then F Sη OL. γ 0 http://www.mii.lt/na

Random convolution of O-exponential distributions 449 A detailed proof of Theorem 1 is presented in Section 3. Note that the proof is similar to that of Theorem 6 in [5]. The following assertion actually shows that Albin s conditions for the counting r.v. η are sufficient for d.f. F Sη to remain in the class OL. The proof of the following corollary is also presented in Section 3. Corollary 1. Let ξ 1, ξ 2,... be a sequence of independent nonnegative r.v.s with common d.f. F ξ OL. (i) D.f. P(ξ 1 + + ξ n x) belongs to the class OL for each fixed n N. (ii) Let η be a r.v. which is nonnegative, non-degenerate at zero, integer-valued and independent of {ξ 1, ξ 2,...}. If Ee εη < for each ε > 0, then F Sη OL. 2 Auxiliary lemmas Before proving our main results, we give three auxiliary lemmas. The first lemma is well known classical estimate for the concentration function of a sum of independent and identically distributed r.v.s. The proof of Lemma 1 can be found in [6] (see Theorem 2.22), for instance. Lemma 1. Let X 1, X 2,..., be a sequence of independent r.v.s with a common nondegenerate d.f. Then there exists a constant c 2, independent of λ and n, such that P(x X 1 + X 2 + + X n x + λ) c 2 (λ + 1)n 1/2 x R for all λ 0 and all n N. The second auxiliary lemma is due to Shimura and Watanabe (see [7, Prop. 2.2]). The lemma describes an important property of a d.f. from the class OL. Lemma 2. Let F be a d.f. from the class OL. Then there exists positive such that lim e x F (x) =. The last auxiliary lemma is crucial in the proof of Theorem 1. The elements of the statement below can be found in [4] (see the proof of Theorem 3(b)). Inequality (1), which is a particular case of the statement below, is proved in [1] (see Lemma 2.1). Leipus and Šiaulys [5] generalized Albin s inequality (1) for an arbitrary d.f. with unbounded port. The analytical proof of Lemma 3 is given in [5] (see proof of Lemma 4). In this paper, we present another, completely probabilistic proof of the lemma below having in mind the importance of the statement. Lemma 3. Let d.f. F be such that F (x) > 0 for all x R. Suppose that F (x t) d 1 x d 2 F (x) Nonlinear Anal. Model. Control, 20(3):447 454

450 S. Danilenko, J. Šiaulys for some positive constants t, d 1 and d 2 > t. Then, for all n = 1, 2,..., we have: F n (x t) d x n(d 2 t)+t F n 1. (x) Proof of Lemma 3. Let X be a r.v. with d.f. F. Then the condition of Lemma 3 says that for some positive t, d 1, d 2 > t, and we need to prove that P(X > x t) d 1 (2) x d 2 P(X > x) P(Sn X > x t) x (nd 2 t)+t P(Sn X d 1 (3) > x) for all n N, where S X n = X 1 + + X n, and X 1, X 2,... are independent copies of X. The proof is proceeded by induction on n. According to condition (2), inequality (3) holds for n = 1. Suppose now that N 1. For arbitrary real x, z and t > 0, we obtain P ( S X N+1 > x ) = P ( S X N + X N+1 > x, X N+1 x z ) + P ( S X N + X N+1 > x, S X N z ) + P(X N+1 > x z) P ( S X N > z ) P ( S X N > x X N+1, x X N+1 z ) + P ( X N+1 > x S X N, x S X N x z + t ) If we replace x by x t and z by z t then we get + P(X N+1 > x z)p ( S X N > z ). (4) P ( S X N+1 > x t ) = P ( S X N + X N+1 > x t, X N+1 x z ) + P ( S X N + X N+1 > x t, S X N z t ) + P(X N+1 > x z)p ( S X N > z t ) = P ( S X N > x X N+1 t, x X N+1 z ) + P ( X N+1 > x S X N t, x S X N x z + t ) + P(X N+1 > x z)p ( S X N > z t ). (5) R.v.s X 1, X 2,... are independent. Therefore, P ( SN X > x X N+1 t, x X N+1 z ) = E ( E ( 1 {S X N >x X N+1 t}1 {x XN+1 z} x X N+1 = y )) = E ( 1 {y z} E ( 1 {S X N >y t} x XN+1 = y )) = E ( 1 {y z} P ( SN X > y t )) http://www.mii.lt/na

Random convolution of O-exponential distributions 451 P(SN X > y t) y z P(SN X > y) E ( 1 {y z} P ( SN X > y )) P(SN X = > y t) y z P(SN X > y) P ( SN X > x X N+1, x X N+1 z ), (6) where 1 A denotes the indicator function of an event A. Similarly, P ( X N+1 > x S X N t, x S X N x z + t ) y x z+t P(X N+1 > y t) P ( X N+1 > x SN X, x SN X x z + t ). (7) P(X N+1 > y) Using estimates (4) (7), we obtain P(S X N+1 > x t) P(S X N+1 > x) { P(SN X max > y t) y z P(SN X > y), y x z+t } P(X > y t) P(X > y) (8) if x, z R, t > 0 and N 1. Suppose now that (3) is satisfied for n = N. We will show that (3) holds for n = N + 1. Condition (2) and estimate (8) imply, taking z = z N = Nx/(N + 1) + t/(n + 1) and w N = x z N + t = x/(n + 1) + Nt/(N + 1), that P(S X N+1 > x t) P(S X N+1 > x) { max P(SN X > y t) y z N P(SN X > y) if x (N + 1)(d 2 t) + t, because, in this case, z N N(d 2 t) + t and w N d 2. } P(X > y t), d 1 y w N P(X > y) So, estimate (3) holds for n = N + 1 and the validity of (3) for all n follows by induction. 3 Proofs of main results In this section, we present detailed proofs of our main results. Proof of Theorem 1. First, we show that F Sη (x a) F Sη (x) = < (9) for each a R. If a 0, then for all x R, and estimate (9) is obvious. Nonlinear Anal. Model. Control, 20(3):447 454

452 S. Danilenko, J. Šiaulys Suppose now that a > 0. Since F ξ OL, we derive that F ξ (x a) F ξ (x) = c 3 (10) for some finite positive quantity c 3 maybe depending on a. So, there exists some K = K a > a + 1 such that Applying Lemma 3, we obtain that F ξ (x a) 2c 3. (11) x K F ξ (x) P(S n > x a) F n = x n(k a)+a P(S n > x) x n (K a)+a where and below S n = ξ 1 + ξ 2 + + ξ n if n N. For an arbitrarily chosen positive x, we have n=1 ξ (x a) 2c 3, (12) (x) F n ξ = P(S n > x)p(η = n) P(ξ 1 > x)p(η = n) If x K, then, using (12), we get: n=1 = F ξ (x)p(η 1). (13) ( = P S η > x a, η x a ) ( + P S η > x a, η > x a ) = P(S n > x a) P(η = n) n (x a)/(k a) + 2c 3 + + n (x a)/(k a) P(S n > x a)p(η = n) P(S n > x)p(η = n) P(S n > x a)p(η = n) P(S n > x) P(η = n) P(S n > x)p(η = n) http://www.mii.lt/na

Random convolution of O-exponential distributions 453 c 4 + n=1 P(S n > x)p(η = n) with c 4 = max{2c 3, 1}. According to Lemma 1, we obtain P(x a < S n x) P(η = n) (14) P(x a < S n x) c 5 (a + 1) 1, x R n where the constant c 5 is independent of a and n. Thus, inequality (14) implies c 4 + c 5 (a + 1) c 4 + c 5 (a + 1)P x a ( P(η = n) n η > x a ) (15) provided that x K. Inequalities (13) and (15) imply that, for x K, it holds c 4 + c ( ) 5 (a + 1) x a x a P(η 1)Fξ (x) F η. Consequently, c 4 + c 5 (a + 1) P(η 1) = c 4 + c 3 c 5 (a + 1) P(η 1) F η ((x a)/()) x a Fξ (x a) F η (x/()) x Fξ (x) < F ξ (x a) F ξ (x) due to equality (10) and requirement F η (δx) = O( x F ξ (x)) which holds for arbitrary δ (0, 1). Therefore, relation (9) is satisfied for for all a R. It remains to prove that lim inf F Sη (x a) F Sη (x) = lim inf > 0, where a is an arbitrarily chosen real number. But this relation follows from the proved estimate (9), because F ξ (x)p(η 1) > 0 Nonlinear Anal. Model. Control, 20(3):447 454

454 S. Danilenko, J. Šiaulys for each positive number x, and so lim inf = ( ) 1 P(S η > x + a) > 0. The last inequality, together with estimate (9), implies that d.f. F Sη belongs to the class OL. Theorem 1 is proved. Proof of Corollary 1. Part (i) of Corollary 1 is evident. So we only prove part (ii). Let δ (0, 1). According to the Markov inequality, we have F η (δx) = P(η > δx) = P ( e yη > e yδx) e δyx Ee yη (16) for each y > 0. The d.f. F ξ belongs to the class OL. Therefore, Lemma 2 implies that e x F ξ (x) as x. for some positive. Choosing y = /δ > 0 in (16), we obtain: F η (δx) x Fξ (x) Ee yη e δyx x F ξ (x) = 1 1 x e x F ξ (x) Ee( /δ)η 0 because Ee εη is finite for an arbitrarily positive ε according to the main condition of Corollary 1. The statement of Corollary 1 follows now from Theorem 1. Acknowledgment. helpful comments. We would like to thank anonymous referees for the detailed and References 1. J.M.P. Albin, A note on the closure of convolution power mixtures (random sums) of exponential distributions, J. Aust. Math. Soc., 84:1 7, 2008. 2. D.B.H. Cline, Convolutions of the distributions with exponential tails, J. Aust. Math. Soc., 43:347 365, 1987. 3. D.B.H. Cline, Convolutions of the distributions with exponential tails: Corrigendum, J. Aust. Math. Soc., 48:152 153, 1990. 4. P. Embrechts, C.M. Goldie, On closure and factorization properties of subexponential and related distributions, J. Aust. Math. Soc., Ser. A, 29:243 256, 1980. 5. R. Leipus, J. Šiaulys, Closure of some heavy-tailed distribution classes under random convolution, Lith. Math. J., 52:249 258, 2012. 6. V.V. Petrov, Limit Theorems of Probability Theory: Sequences of Independent Random Variables, Oxford Studies in Probability, Vol. 4, Clarendon Press, Oxford, 1995. 7. T. Shimura, T. Watanabe, Infinite divisibility and generalized subexponentiality, Bernoulli, 11:445 469, 2005. 8. T. Watanabe, K. Yamamuro, Ratio of the tail of an infinitely divisible distribution on the line to that of its Lévy measure, Electron. J. Probab., 15:44 74, 2010. http://www.mii.lt/na