Next Article in Journal
Numerical Investigation of the Magnetized Reactive Viscous Couple Stress Fluid Flow Down an Inclined Riga Plate with Variable Viscosity
Next Article in Special Issue
Delicate Comparison of the Central and Non-Central Lyapunov Ratios with Applications to the Berry–Esseen Inequality for Compound Poisson Distributions
Previous Article in Journal
Linear and Energy-Stable Method with Enhanced Consistency for the Incompressible Cahn–Hilliard–Navier–Stokes Two-Phase Flow Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparing Compound Poisson Distributions by Deficiency: Continuous-Time Case

by
Vladimir Bening
1,2 and
Victor Korolev
1,2,3,4,*
1
Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, 119991 Moscow, Russia
2
Moscow Center for Fundamental and Applied Mathematics, 119991 Moscow, Russia
3
Federal Research Center “Computer Science and Control”, Russian Academy of Sciences, 119333 Moscow, Russia
4
Department of Mathematics, School of Science, Hangzhou Dianzi University, Hangzhou 310018, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(24), 4712; https://doi.org/10.3390/math10244712
Submission received: 17 November 2022 / Revised: 6 December 2022 / Accepted: 8 December 2022 / Published: 12 December 2022

Abstract

:
In the paper, we apply a new approach to the comparison of the distributions of sums of random variables to the case of Poisson random sums. This approach was proposed in our previous work (Bening, Korolev, 2022) and is based on the concept of statistical deficiency. Here, we introduce a continuous analog of deficiency. In the case under consideration, by continuous deficiency, we will mean the difference between the parameter of the Poisson distribution of the number of summands in a Poisson random sum and that of the compound Poisson distribution providing the desired accuracy of the normal approximation. This approach is used for the solution of the problem of determination of the distribution of a separate term in the Poisson sum that provides the least possible value of the parameter of the Poisson distribution of the number of summands guaranteeing the prescribed value of the ( 1 α ) -quantile of the normalized Poisson sum for a given α ( 0 , 1 ) . This problem is solved under the condition that possible distributions of random summands possess coinciding first three moments. The approach under consideration is applied to the collective risk model in order to determine the distribution of insurance payments providing the least possible time that provides the prescribed Value-at-Risk. This approach is also used for the problem of comparison of the accuracy of approximation of the asymptotic ( 1 α ) -quantile of the sum of independent, identically distributed random variables with that of the accompanying infinitely divisible distribution.

1. Introduction

This paper is a complement to our previous work [1], where we considered a version of the problem of stochastic ordering and proposed an approach based on the concept of deficiency that is well-known in asymptotic statistics; see, e.g., [2] and later publications [3,4,5,6]. In the paper [1], we used the approach mentioned above in order to establish a kind of stochastic order for the distributions of sums of independent random variables (r.v.s) based on the comparison of the number of summands required for the distribution of the sum to have the desired asymptotic properties (for the problems and methods related to stochastic ordering, see, e.g., [7]). Here, we apply this approach to the comparison of the distributions of sums of random variables to the case of Poisson random sums.
In statistics, as well as in [1], the deficiency is measured in integer units and correspondingly has the meaning of either the number of additional observations required for a statistical procedure to attain the same quality as the ‘optimal’ procedure in statistics or the number of additional summands in the sum required to attain the desired accuracy of the normal approximation in [1]. Unlike these cases, in the present paper, we deal with the compound Poisson distributions and introduce a continuous analog of deficiency. The extension of the approach proposed in [1] for non-random sums of independent r.v.s to Poisson random sums is possible due to the asymptotic normality of the latter as the parameter of the Poisson distribution of the number of summands infinitely grows. In the case under consideration, by continuous deficiency, we mean the difference between the parameter of the Poisson distribution of the number of summands in a Poisson random sum and that of the compound Poisson distribution providing the desired accuracy of the normal approximation. This approach is used for the solution of the problem of determination of the distribution of a separate term in the Poisson sum that provides the least possible value of the parameter of the Poisson distribution of the number of summands guaranteeing the prescribed value of the ( 1 α ) -quantile of the normalized Poisson sum for a given α ( 0 , 1 ) .
This problem is solved under the condition that possible distributions of random summands possess coinciding first three moments. Therefore, we can say that, in this problem, we deal with ‘fine tuning’ of the distribution of a separate summand since we assume that different possible distributions of random summands may differ only by their kurtosis. In the setting under consideration, the best distribution delivers the smallest value of the parameter of the compounding Poisson distribution. This problem is actually a particular case of the problem of quantification of the accuracy of approximations of the compound Poisson distributions provided by limit theorems of probability theory. The main mathematical tools used in the paper are asymptotic expansions for the compound Poisson distributions and their quantiles.
The formal setting mentioned above can be applied to solving some practical problems dealing with the collective risk insurance models where it is traditional to describe the cumulative insurance payments by the compound Poisson process. The approach under consideration makes it possible to determine the distribution of insurance payments providing the least possible time that provides the prescribed Value-at-Risk.
To make the above-mentioned more clear, consider an insurance company that starts its activity at time t 0 = 0 . Within the classical collective risk model [8], the total insurance payments at some time t have the form of a sum of a random number (number of payments by the time t) of independent identically distributed r.v.s (insurance payments), that is, of a Poisson random sum. In this model, the number of insurance payments by time t follows the Poisson process N λ ( t ) with some intensity λ > 0 . We assume that the parameter λ is uncontrollable and fixed. Since N λ ( t ) has the same distribution as N 1 ( λ t ) and the parameter λ is assumed fixed, the setting under consideration concerns the problem of determination of the distribution of an individual insurance payment providing the least possible t guaranteeing the prescribed Value-at-Risk for the average losses of the insurance company within the time interval [ 0 , t ] .
The approach considered in the paper can be used when the distributions of the summands (possible losses) are known only up to their first three moments, and the exact Value-at-Risk is not known for sure.
Within the framework of the collective risk model in the setting under consideration, the problem consists in the description of the best strategy of the insurance company. Here, the choice of the terms of a contract (e.g., the amount of insurance payment related to each possible insurance event) is meant as a strategy. That is, a strategy consists in the determination of the distribution of an insurance payment. Briefly, the problem is to choose an optimal distribution of a separate insurance payment among the distributions that have the same first three moments so that the desired goal is achieved within the least possible time interval.
We also consider the application of the proposed approach to the study of the asymptotic properties of non-random sums of independent identically distributed r.v.s as compared to those of the compound Poisson distributions with the same expectation. It is well-known that, in many respects, these properties coincide. This phenomenon manifests itself, for example, in the form of the method of accompanying infinitely divisible distributions (see, e.g., [9], Chapter 4, Section 24). Therefore, it is of certain interest to investigate the accuracy of the approximation of the characteristics of sums of independent r.v.s as compared to that of the accompanying infinitely divisible (that is, corresponding compound Poisson) laws. This problem was studied by many specialists; see, e.g., [10,11,12,13] and the references therein. Unlike most preceding works where the approximation of distribution functions was discussed, here we consider the application of accompanying laws to a somewhat inverse problem of approximation of quantiles.
The paper is organized as follows. Section 2 contains a short overview of the results concerning the asymptotic expansions for compound Poisson distributions. Here we also formulate basic lemmas to be used in the next sections. The main results are presented in Section 3. In Section 3.1, we introduce the notion of the α -reserve in the collective risk model and present some asymptotic expansions for this quantity. In Section 3.2, a continuous-time analog of the notion of deficiency is introduced. Here we also prove some general results concerning the continuous-time deficiency. In Section 3.3, we consider the problem of comparison of compound Poisson distributions by deficiency and present the asymptotic formula for the deficiency of one compound Poisson distribution with respect to the other. In Section 3.4, we deal with the problem of comparison of the distributions of Poisson random sums with those of non-random sums. Actually, this problem consists in the comparison of the accuracy of approximation of the asymptotic ( 1 α ) -quantile of the sum of independent identically distributed random variables with that of the accompanying infinitely divisible distribution.

2. Notation and Auxiliary Results

Throughout what follows, we will assume that all the random variables and processes are defined on the same probability space ( Ω , F , P ) . The expectation and variance with respect to the measure P will be, respectively, denoted E and D . The set of real numbers and natural numbers will be, respectively, denoted R and N . The distribution function of the standard normal law will be denoted Φ ( x ) ,
Φ ( x ) = 1 2 π x φ ( y ) d y , φ ( x ) = 1 2 π exp x 2 2 , x R .
The distribution of a random variable X will be denoted L ( X ) .
Let X 1 , X 2 , be independent identically distributed random variables. Let N λ be the random variable with the Poisson distribution with parameter λ . Assume that for each λ > 0 , the random variables N λ , X 1 , X 2 , are independent. Let S λ be the Poisson random sum, S λ = X 1 + + X N λ . If N λ = 0 , then S λ is assumed to equal to zero. Assume that E X 1 = a and D X 1 = σ 2 > 0 exist. For integer k 0 , denote E X 1 k = α k . Of course, α 0 = 1 , α 1 = a and α 2 = σ 2 + a 2 .
Recall some facts concerning the asymptotic expansions for the compound Poisson distributions (sf. [8,14,15]).
Denote the characteristic functions of the random variables X 1 and S λ as f ( t ) and h λ ( t ) , respectively. It is well-known that if f ( t ) has r continuous derivatives, then, as t 0 , we have
f ( t ) = 1 + i a t α 2 t 2 2 + ( i t ) 2 k = 1 r 2 ( i t ) k α k + 2 ( k + 2 ) ! + o ( t r ) .
A random variable X 1 is said to satisfy the Cramér condition (C), if
lim sup | t | | f ( t ) | < 1 .
For k = 0 , 1 , 2 , define the function H k ( x ) : R R as
H k ( x ) ( 1 ) k ϕ ( k ) ( x ) ϕ ( x ) .
The function H k ( x ) , x R , so defined, is a polynomial of degree k and is called the Hermite polynomial of degree k.
It is easy to calculate that
H 0 ( x ) = 1 , H 1 ( x ) = x , H 2 ( x ) = x 2 1 , H 3 ( x ) = x 3 3 x , H 4 ( x ) = x 4 6 x 2 + 3 ,
H 5 ( x ) = x 5 10 x 3 + 15 x , H 6 ( x ) = x 6 15 x 4 + 45 x 2 15 .
Let m be a nonnegative integer and q k R , k = 0 , , m . Consider the polynomial
q ( x ) = k = 0 m q k x k .
Let H 0 ( x ) , , H m ( x ) be Hermite polynomials. Let
Q ( x ) = k = 0 m q k H k ( x ) .
Then it is easy to make sure that the function v ( t ) = q ( i t ) exp { t 2 / 2 } is the Fourier transform of the function V ( x ) = Q ( x ) ϕ ( x ) . Throughout what follows, we will assume that r 3 is a fixed integer number.
For a complex z, let
f ˜ ( z ) = k = 1 r 2 α k + 2 z k ( k + 2 ) ! .
Obviously, f ˜ ( z ) is a polynomial of degree r 2 with real coefficients; moreover, f ˜ ( 0 ) = 0 . From (1), it follows that
f ( t ) 1 i a t + α 2 t 2 2 = ( i t ) 2 f ˜ ( i t ) + o ( t r )
as t 0 . For λ > 0 and a complex z let
p λ ( z ) = k = 1 r 2 1 k ! z 2 α 2 f ˜ z λ α 2 k .
It can be easily made sure that there exist integer m 3 and polynomials q k ( z ) with real coefficients, k = 3 , , m , not depending on λ such that
p λ ( z ) = k = 3 m λ k / 2 + 1 q k ( z )
for all λ > 0 and complex z. Moreover, these polynomials q k ( z ) are uniquely determined by (3) and (4). Let
q k ( z ) = j = 3 L k q k , j z j
be the corresponding representation of q k ( z ) with q k , j R ( j = 3 , , L k ), L k 3 ( k = 3 , , m ). Let H j ( x ) be the Hermite polynomials. For x R and k = 3 , , m let
R k ( x ) = j = 3 L k q k , j H j 1 ( x ) .
The function R k ( x ) is called the Edgeworth polynomial of degree k.
For λ > 0 and complex z from (3) and (4), we easily obtain
p λ ( z ) = k = 3 ( r 2 ) 2 + 2 λ k / 2 + 1 k 2 r 2 j k 2 α k , j z k + 2 ( j 1 ) ,
where
j ! α k , j = 3 n 1 n j r n 1 + + n j = k + 2 ( j 1 ) α n 1 · · α n j n 1 ! · · n j ! α 2 k / 2 j + 1 .
Therefore, in (4) and (5), we should set m = ( r 2 ) 2 + 2 and L k = 3 ( k 2 ) ( k = 3 , , m ).
For x R , λ > 0 and r N define the functions G λ , r ( x ) as
G λ , r ( x ) = Φ ( x ) + ϕ ( x ) k = 3 r λ k / 2 + 1 R k ( x ) .
In particular, for r = 3 , we have
R 3 ( x ) = α 3 6 α 2 3 / 2 H 2 ( x )
and
G λ , 3 ( x ) = Φ ( x ) α 3 6 α 2 3 / 2 λ ( x 2 1 ) ϕ ( x ) .
For r = 4 , we have
R 4 ( x ) = α 4 24 α 2 2 H 3 ( x ) α 3 2 72 α 2 3 H 5 ( x )
and
G λ , 4 ( x ) = Φ ( x ) α 3 6 α 2 3 / 2 λ ( x 2 1 ) ϕ ( x ) ϕ ( x ) λ α 4 24 α 2 2 ( x 3 3 x ) α 3 2 72 α 2 3 ( x 5 10 x 3 + 15 x ) .
Moreover, if ϰ 3 ( S λ ) and ϰ 4 ( S λ ) are the skewness and kurtosis of the random variable S λ ,
ϰ 3 ( S λ ) E S λ E S λ D S λ 3 = E S λ α 1 λ λ α 2 3 = α 3 λ α 2 3 / 2 ,
ϰ 4 ( S λ ) E S λ E S λ D S λ 4 3 = E S λ α 1 λ λ α 2 4 3 = α 4 λ α 2 2 ,
then (7) and (8) can be rewritten as
G λ , 3 ( x ) = Φ ( x ) ϰ 3 ( S λ ) 6 Φ ( 3 ) ( x )
and
G λ , 4 ( x ) = Φ ( x ) ϰ 3 ( S λ ) 6 Φ ( 3 ) ( x ) + ϰ 4 ( S λ ) 24 Φ ( 4 ) ( x ) + ϰ 3 2 ( S λ ) 72 Φ ( 6 ) ( x ) .
Lemma 1. 
Let r > 3 . Assume that the distribution of the random variable X 1 satisfies the Cramér condition ( C ) (see (2)). Then
sup x P S λ a λ λ ( a 2 + σ 2 ) < x G λ , r ( x ) = o λ r / 2 + 2 ,
that is,
lim λ λ r / 2 1 sup x P S λ a λ λ ( a 2 + σ 2 ) < x G λ , r ( x ) = 0 .
This statement is a particular case of Theorem 4.4.1 in [15].
Our further reasoning is based on the following general statement dealing with the asymptotic behavior of the quantiles of univariate distributions of a random process.
Let Z ( t ) , t 0 , be a random process. Assume that for each t 0 the distribution of the random variable Z ( t ) is continuous. For β ( 0 , 1 ) and t 0 , the left β -quantile of the random variable Z ( t ) will be denoted u β ( t ) :
u β ( t ) = inf { u : P ( Z ( t ) < u ) β } .
Lemma 2. 
Assume that, as t , the distribution function of the random process Z ( t ) admits the asymptotic expansion of the form
P ( Z ( t ) < x ) = Ψ 0 ( x ) + t 1 / 2 Ψ 1 ( x ) + t 1 Ψ 2 ( x ) + o ( t 1 ) .
Moreover, let the functions Ψ 0 ( x ) , Ψ 1 ( x ) and Ψ 2 ( x ) be continuous and Ψ 0 ( x ) > 0 . Then for any β ( 0 , 1 ) , we have
u β ( t ) = u β Ψ 1 ( u β ) Ψ 0 ( u β ) t + Ψ 0 ( u β ) Ψ 1 ( u β ) Ψ 1 ( u β ) Ψ 0 ( u β ) 2 Ψ 2 ( u β ) 1 2 Ψ 1 2 ( u β ) Ψ 0 ( u β ) Ψ 0 ( u β ) 3 t + o ( t 1 ) ,
where u β is the left β-quantile of the distribution function Ψ 0 ( x ) : Ψ 0 ( u β ) = β .
For the proof of this statement, see [15], Section 4.5.
Remark 1. 
If we set
u ¯ β ( t ) = u β Ψ 1 ( u β ) Ψ 0 ( u β ) t + Ψ 0 ( u β ) Ψ 1 ( u β ) Ψ 1 ( u β ) Ψ 0 ( u β ) 2 Ψ 2 ( u β ) 1 2 Ψ 1 2 ( u β ) Ψ 0 ( u β ) Ψ 0 ( u β ) 3 t ,
then it is not difficult to make sure that under the conditions of Lemma 2, we have
P ( Z ( t ) < u ¯ β ( t ) ) = β + o ( t 1 ) .
From Lemmas 1 and 2, it follows that if α 4 = E X 1 4 < and the random variable X 1 satisfies the Cramér (C) condition (2), then
P S λ a λ λ ( a 2 + σ 2 ) < x = Φ ( x ) + Ψ 1 ( x ) λ + Ψ 2 ( x ) λ + o ( λ 1 )
where
Ψ 1 ( x ) = α 3 6 α 2 3 / 2 ϕ ( x ) H 2 ( x ) , Ψ 2 ( x ) = ϕ ( x ) α 4 24 α 2 2 H 3 ( x ) + α 3 2 72 α 2 3 H 5 ( x ) .
Therefore, setting t = λ , Z ( t ) = S λ , Ψ 0 ( x ) = Φ ( x ) , from Lemma 2, we obtain the following result. For β ( 0 , 1 ) , let w β ( λ ) and u β be the β -quantiles of the random variable S λ and of the standard normal distribution, respectively.
Lemma 3. 
Let E X 1 4 < , and let the random variable X 1 satisfy the Cramér ( C ) condition (2). Then, as λ , we have
w β ( λ ) = a λ + u β λ α 2 + α 3 H 2 ( u β ) 6 α 2 +
+ 1 λ α 2 5 / 2 α 3 2 72 H 5 ( u β ) 2 H 2 ( u β ) H 3 ( u β ) + 4 u β H 2 2 ( u β ) + α 4 α 2 24 H 3 ( u β ) + o ( λ 1 / 2 )
where H k ( x ) are the Hermite polynomials.

3. Main Results

3.1. The Asymptotic Expansions for the α -Reserve in the Collective Risk Model

Let X 1 , X 2 , be independent identically distributed r.v.s such that
X 1 2 > 0 , | X 1 | 4 + δ < , δ > 0 .
Assume that the r.v. X 1 satisfies the Cramér ( C ) condition (2). For t > 0 , let the r.v. N t have the Poisson distribution with parameter λ t , where λ > 0 is a fixed parameter. Assume that for each t > 0 the r.v.s N t , X 1 , X 2 , are independent. Consider the Poisson random sum
S t = X 1 + X N t .
In terms of the collective risk model, the r.v.s X j can be interpreted as individual insurance claims, and the r.v. S t can be interpreted as the total insurance payment of an insurance company by the time t.
Let α ( 0 , 1 ) . Define the standardized α-reserve C α * ( t ) by the formula
P S t λ t E X 1 λ t E X 1 2 C α * ( t ) = α + o ( t 1 ) , t .
Along with the set X 1 , X 2 , consider another set Y 1 , Y 2 , of independent identically distributed r.v.s such that
Y 1 2 > 0 , | Y 1 | 4 + δ < , δ > 0 .
Assume that the r.v. Y 1 satisfies the Cramér ( C ) condition (2). Also assume that for each t > 0 , the r.v. N t having the Poisson distribution with parameter λ t is independent of the set Y 1 , Y 2 , Denote
T t = Y 1 + + Y N t .
In the same way as (11), define the standardized α-reserve C α * * ( t ) for the sequence Y 1 , Y 2 , as
P T t λ t E Y 1 λ t E Y 1 2 C α * * ( t ) = α + o ( t 1 ) , t .
Lemmas 2 and 3 directly imply the following statement. For α ( 0 , 1 ) let u α be the 1 α -quantile of the standard normal distribution, that is, Φ ( u α ) = 1 α .
Theorem 1. 
Let α ( 0 , 1 ) and the r.v.s X 1 , X 2 , and Y 1 , Y 2 , satisfy conditions (10), (12) and (2). Then, as t ,
C α * ( t ) = u α + E X 1 3 ( u α 2 1 ) 6 λ t ( E X 1 2 ) 3 / 2 + 1 12 λ t E X 1 2 ( E X 1 3 ) 2 E X 1 2 ( 5 u α 2 u α 3 ) + E X 1 4 2 ( E X 1 2 ) 2 ( u α 3 3 u α ) + o ( t 1 ) ,
C α * * ( t ) = u α + E Y 1 3 ( u α 2 1 ) 6 λ t ( E Y 1 2 ) 3 / 2 + 1 12 λ t E Y 1 2 ( E Y 1 3 ) 2 E Y 1 2 ( 5 u α 2 u α 3 ) + E Y 1 4 2 ( E Y 1 2 ) 2 ( u α 3 3 u α ) + o ( t 1 ) .
We see that if the first three moments of X 1 and Y 1 coincide, then C α * ( t ) and C α * * ( t ) differ only by the terms of order O ( t 1 ) .
Now if we define the α-reserves C ˜ α * ( t ) and C ˜ α * * ( t ) as
P S t C ˜ α * ( t ) = α + o ( t 1 ) , and P T t C ˜ α * * ( t ) = α + o ( t 1 ) , t ,
then
C ˜ α * ( t ) = λ t E X 1 2 · C α * ( t ) + λ t E X 1 and C ˜ α * * ( t ) = λ t E Y 1 2 · C α * ( t ) + λ t E Y 1 .

3.2. A Continuous-Time Analog of Deficiency

In this section, we will propose an approach to the comparison of the two compound Poisson distributions in terms of the ‘continuous’ analog of deficiency. For the traditional definition of deficiency as the number of additional observations required for a statistical procedure to attain the desired quality, we refer the reader to the papers [1,2,3,5,6]. Here, we will introduce its continuous-time analog.
Consider two stochastic processes X ( t ) and Y ( t ) , t 0 . We will be interested in the asymptotic behavior of the probabilities of X ( t ) and Y ( t ) to exceed a given threshold.
For α ( 0 , 1 ) let c α ( t ) be the asymptotic ( 1 α ) -quantile of X ( t ) :
P X ( t ) c α ( t ) = α + o ( t 1 ) , t .
Lemma 2 directly implies the following statement.
Proposition 1. 
Assume that there exist distribution function G ( x ) and the functions g 1 ( x ) and g 2 ( x ) such that
sup x R | P X ( t ) < x G ( x ) 1 t g 1 ( x ) 1 t g 2 ( x ) | = o ( t 1 ) ,
where the functions G ( x ) , g 1 ( x ) and g 2 ( x ) are smooth enough. Then the the asymptotic ( 1 α ) -quantile of X ( t ) admits the asymptotic expansion
c α ( t ) = c α g 1 ( c α ) G ( c α ) t 1 t G ( c α ) g 1 2 ( c α ) 2 ( G ( c α ) ) 3 + G ( c α ) g 2 ( c α ) g 1 ( c α ) g 1 ( c α ) ( G ( c α ) ) 2 + o ( t 1 ) ,
where c α is the ( 1 α ) -quantile of the distribution function G ( x ) , that is, G ( c α ) = 1 α .
Assume that the asymptotic expansion for the distribution function of Y ( t ) has the form
P Y ( t ) < x = G ( x ) + 1 t g 1 ( x ) + 1 t g ¯ 2 ( x ) + o ( t 1 ) ,
where the functions G ( x ) , g 1 ( x ) and g ¯ 2 ( x ) are smooth enough. The asymptotic expansion (14) differs from that for the distribution function of X ( t ) in Proposition 1 only by the term of order t 1 , that is, the two distributions are close enough.
Define the positive function m ( t ) , t > 0 , by the equality
P t Y ( m ( t ) ) c α ( m ( t ) ) = α + o ( t 1 ) , t .
If m ( t ) t = d + o ( 1 ) , d R , t , then the number d is called the asymptotic deficiency of the distribution L ( Y ( t ) ) with respect to the distribution L ( X ( t ) ) . In other words, d is the asymptotic ‘additional’ time required for the process Y ( t ) to attain the quantile of the same order as that of X ( t ) .
Theorem 2. 
Assume that conditions (13) and (14) hold. Then the asymptotic deficiency d of the distribution L ( Y ( t ) ) with respect to the distribution L ( X ( t ) ) has the form
d = 2 g 2 ( c α ) g ¯ 2 ( c α ) G ( c α ) c α + o ( 1 ) .
The proof of this statement repeats that of Theorem 3.1 in [1] up to notation (furthermore, unfortunately, in formula (16) of [1], the coefficient n analogous to t in (15) of the present paper was erroneously omitted).

3.3. The Comparison of Compound Poisson Distributions by Deficiency

In this section, we will discuss the asymptotic deficiency of the compound Poisson distributions providing a given ( 1 α ) -quantile of the normalized Poisson random sums. For this purpose, we will use Theorem 2.
Define the average Poisson random sums  S ¯ t and T ¯ t by the formulas
S ¯ t = S t λ t E X 1 t λ E X 1 2 , T ¯ t = T t λ t E Y 1 t λ E Y 1 2 .
Define the asymptotic deficiency  d R of T ¯ t with respect to S ¯ t by the formula
P t · T ¯ t ¯ C α * ( t ¯ ) = α + o ( t ) , t ,
where t ¯ = t + d + o ( 1 ) , that is, d is the ‘additional time’ required for the normalized average Poisson random sum t · T ¯ t to exceed the asymptotic α -reserve C α * ( t ) of the normalized average Poisson random sum t · S ¯ t .
To apply Theorem 2, assume that
E X 1 3 ( E X 1 2 ) 3 / 2 = E Y 1 3 ( E Y 1 2 ) 3 / 2 .
Condition (16) holds, e.g., if the first three moments of X 1 and Y 1 coincide.
Theorem 2 directly implies the following statement.
Theorem 3. 
Assume that the r.v.s N t , X 1 , X 2 , ; Y 1 , Y 2 , satisfy conditions (2), (10) and (16). Then, as t , the ‘additional time’ d has the form
d = ( 3 u α 2 ) 12 E X 1 4 ( E X 1 2 ) 2 E Y 1 4 ( E Y 1 2 ) 2 + o ( 1 ) .
Remark 2. 
If E X 1 = E Y 1 = 0 , then (17) can be rewritten as
d = 1 12 ( 3 u α 2 ) ϰ 4 ( X 1 ) ϰ 4 ( Y 1 ) + o ( 1 ) ,
That is, in this case, the continuous-time analog of asymptotic deficiency is determined by the difference of kurtoses.

3.4. Comparing the Distributions of Poisson Random Sums with Those of Non-Random Sums

It is well-known that the asymptotic properties of non-random sums of independent identically distributed r.v.s coincide with those of the compound Poisson distributions with the same expectation. This phenomenon manifests itself, for example, in the form of the method of accompanying infinitely divisible distributions (see, e.g., [9], Chapter 4, Section 24). Therefore, it is of certain interest to investigate the accuracy of the approximation of the characteristics of sums of independent r.v.s as compared to that of the accompanying infinitely divisible (that is, corresponding compound Poisson) laws. This problem was studied by many specialists, see, e.g., [10,11,12,13]. Unlike most preceding works where the approximation of distribution functions was discussed, here we consider the application of accompanying laws to a somewhat inverse problem of approximation of quantiles.
Here, we will not assume the possibility of the interpretation of the presented results in terms of a collective risk model where at least the expectations of X j should be positive. Assume that the independent identically distributed r.v.s X 1 , X 2 , are standardized:
E X 1 = 0 , E X 1 2 = 1 .
Again, let N t be an r.v. with the Poisson distribution with parameter λ t , where λ > 0 is fixed. Assume that for each t > 0 the random variables N t , X 1 , X 2 , are independent. Consider the problem of comparison of the distribution of a normalized Poisson random sum
S t * = X 1 + + X N t λ t
with the distribution of the corresponding non-random sum
U t * = X 1 + + X [ λ t ] [ λ t ]
as t , where the symbol [ a ] denotes the integer part of a real number a. For definiteness, if N t = 0 , then S t * is assumed to be equal to zero.
If conditions (18), (10) and (2), then Lemmas 1 and 2 imply (see (9)) that, as t ,
P ( S t * < x ) = Φ ( x ) E X 1 3 6 λ t φ ( x ) ( x 2 1 )
φ ( x ) 24 λ t E X 1 4 ( x 3 3 x ) + ( E X 1 3 ) 2 3 ( x 5 10 x + 15 x ) + o ( t 1 ) ,
whereas the classical theory of asymptotic expansions in the central limit theorem (e.g., see [16]) yields that
P ( U t * < x ) = Φ ( x ) E X 1 3 6 [ λ t ] φ ( x ) ( x 2 1 )
φ ( x ) 24 [ λ t ] ( E X 1 4 3 ) ( x 3 3 x ) + ( E X 1 3 ) 2 3 ( x 5 10 x + 15 x ) + o ( t 1 ) .
Note that (19) and (20) differ in that, in (19), the kurtosis of X 1 is present in the non-normalized form ϰ 4 * ( X 1 ) = E X 1 4 , whereas in (20), there stands the normalized kurtosis ϰ 4 ( X 1 ) = E X 1 4 3 .
From the obvious inequalities
λ t 1 [ λ t ] λ t
it follows that, as t ,
1 λ t 1 [ λ t ] 1 λ t 1 = 1 λ t 1 + 1 λ t + O ( t 2 )
and
1 [ λ t ] = 1 λ t + O ( t 3 / 2 ) .
Therefore, relation (20) can be rewritten as
P ( U t * < x ) = Φ ( x ) E X 1 3 6 λ t φ ( x ) ( x 2 1 )
φ ( x ) 24 λ t E ( X 1 4 3 ) ( x 3 3 x ) + ( E X 1 3 ) 2 3 ( x 5 10 x + 15 x ) + o ( t 1 ) .
Denote U ¯ t * = U t * / t . Let α ( 0 , 1 ) . Define the asymptotic ( 1 α ) -quantile C α ( t ) of S t * by the relation
P S t * C α ( t ) = α + o ( t 1 ) , t .
Define the number d R by the formula
P t U ¯ t ¯ * C α ( t ¯ ) = α + o ( t 1 ) , t ,
where t ¯ = t + d + o ( 1 ) . Now relations (19), (21) and Theorem 2 directly imply the following statement.
Theorem 4. 
Let α ( 0 , 1 ) . Assume that the r.v.s N t , X 1 , X 2 , satisfy conditions (18), (10) and (2). Then
d = 3 u α 2 4 + o ( 1 )
as t , where Φ ( u α ) = 1 α .
Remark 3. 
The quantity d can be interpreted as the asymptotic deficiency of the distribution of a non-random sum with respect to the corresponding accompanying compound Poisson distribution. Note that under the conditions of Theorem 4, d does not depend on the distribution of X 1 . If α > 0.0417 . . . , then d is asymptotically positive, that is, the (accompanying) compound Poisson distribution of the r.v. S t * provides better accuracy for the approximation of the asymptotic ( 1 α ) -quantile of U t * .

4. Conclusions

This paper is a continuation of our previous work [1] and deals with a version of the problem of stochastic ordering. We follow an approach based on the concept of deficiency, which is well-known in asymptotic statistics. In the present paper, we considered compound Poisson distributions and introduced a continuous analog of deficiency. It was suggested to understand the continuous deficiency as the difference between the parameter of the compounding distribution of a Poisson random sum and that of the compound Poisson distribution providing the desired accuracy of the normal approximation. The asymptotic representations for the continuous deficiency were obtained under the condition that possible distributions of random summands possess coinciding first three moments. Therefore, we can say that, in this problem, we deal with ‘fine tuning’ of the distribution of a separate summand since we assume that different possible distributions of random summands can differ only by their kurtosis. In the setting under consideration, the best distribution delivers the smallest value of the parameter of the compounding Poisson distribution. The main mathematical tools used in the paper are asymptotic expansions for the compound Poisson distributions and their quantiles. The formal setting mentioned above was applied to solving some practical problems dealing with the collective risk insurance models where it is traditional to describe the cumulative insurance payments by the compound Poisson process. The approach under consideration makes it possible to determine the distribution of insurance payments providing the least possible time that provides the prescribed Value-at-Risk. We also considered the application of the proposed approach to the study of the asymptotic properties of non-random sums of independent identically distributed r.v.s as compared to those of the compound Poisson distributions with the same expectation. We investigate the accuracy of the approximation of the characteristics of sums of independent r.v.s as compared to that of the accompanying infinitely divisible (that is, corresponding compound Poisson) laws. Unlike most preceding works where the approximation of distribution functions was discussed, here we considered the application of accompanying laws to a somewhat inverse problem of approximation of quantiles.

Author Contributions

Conceptualization, V.B. and V.K.; methodology, V.B. and V.K.; validation, V.B. and V.K.; formal analysis, V.B. and V.K.; investigation, V.B. and V.K.; writing—original draft preparation, V.K.; writing—review and editing, V.K.; supervision, V.K.; project administration, V.K.; funding acquisition, V.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Russian Science Foundation, grant 22-11-00212.

Acknowledgments

The authors thank Alexander Zeifman for his help in the final preparation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bening, V.E.; Korolev, V.Y. Comparing distributions of sums of random variables by deficiency: Discrete case. Mathematics 2022, 10, 454. [Google Scholar] [CrossRef]
  2. Hodges, J.L.; Lehmann, E.L. Deficiency. Ann. Math. Stat. 1970, 41, 783–801. [Google Scholar] [CrossRef]
  3. Xiang, X. Deficiency of the sample quantile estimator with respect to kernel quantile estimators for censored data. Ann. Stat. 1995, 23, 836–854. [Google Scholar] [CrossRef]
  4. Bening, V.E. Asymptotic Theory of Testing Statistical Hypotheses: Efficient Statistics, Optimality, Power Loss, and Deficiency; Walter de Gruyter: Berlin, Germany, 2011. [Google Scholar]
  5. Torgersen, E. Comparison of Statistical Experiments; Cambridge University Press: Cambridge, UK, 1991. [Google Scholar]
  6. Bening, V.E.; Korolev, V.Y.; Zeifman, A.I. Calculation of the deficiency of some statistical estimators constructed from samples with random sizes. Colloq. Math. 2019, 157, 157–171. [Google Scholar] [CrossRef]
  7. Müller, A.; Stoyan, D. Comparison Methods for Stochastic Models and Risks; John Wiley & Sons: Chichester, UK, 2002. [Google Scholar]
  8. Cramér, H. Collective Risk Theory; Skandia Jubilee Volume: Stockholm, Sweden, 1955. [Google Scholar]
  9. Gnedenko, B.V.; Kolmogorov, A.N. Limit Distributions for Sums of Independent Random Variables; Addison-Wesley: Reading, MA, USA, 1954. [Google Scholar]
  10. Prokhorov, Y.V. Strong stability of sums and infinitely divisible distributions. Theory Probab. Its Appl. 1958, 4, 141–153. [Google Scholar] [CrossRef]
  11. Arak, T.V. On the approximation of n-fold convolutions of distributions having non-negative characteristic functions with accompanying laws. Theory Probab. Its Appl. 1981, 25, 221–245. [Google Scholar] [CrossRef]
  12. Čekanavičius, V. Approximation by accompanying distributions and asymptotic expansions. I; II. Lith. Math. J. (Liet. Mat. Rink.) 1989, 29, 171–178, 402–415. [Google Scholar]
  13. Zaitsev, A.Y. Approximation of convolutions of multi-dimensional symmetric distributions by accompanying laws. J. Sov. Math. 1992, 61, 1859–1872. [Google Scholar] [CrossRef]
  14. Von Chossy, R.; Rappl, G. Some approximation methods for the distribution of random sums. Insur. Math. Econ. 1983, 2, 251–270. [Google Scholar] [CrossRef]
  15. Bening, V.; Korolev, V. Generalized Poisson Models and Their Applications in Insurance and Finance; Walter de Gruyter: Berlin, Germany, 2012. [Google Scholar]
  16. Bhattacharya, R.N.; Ranga Rao, R. Normal Approximation and Asymptotic Expansions; John Wiley & Sons: New York, NY, USA; London, UK; Sydney, Australia, 1976. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bening, V.; Korolev, V. Comparing Compound Poisson Distributions by Deficiency: Continuous-Time Case. Mathematics 2022, 10, 4712. https://doi.org/10.3390/math10244712

AMA Style

Bening V, Korolev V. Comparing Compound Poisson Distributions by Deficiency: Continuous-Time Case. Mathematics. 2022; 10(24):4712. https://doi.org/10.3390/math10244712

Chicago/Turabian Style

Bening, Vladimir, and Victor Korolev. 2022. "Comparing Compound Poisson Distributions by Deficiency: Continuous-Time Case" Mathematics 10, no. 24: 4712. https://doi.org/10.3390/math10244712

APA Style

Bening, V., & Korolev, V. (2022). Comparing Compound Poisson Distributions by Deficiency: Continuous-Time Case. Mathematics, 10(24), 4712. https://doi.org/10.3390/math10244712

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop