Next Article in Journal
Modeling and Simulation Techniques Used in High Strain Rate Projectile Impact
Next Article in Special Issue
Asymptotically Exact Constants in Natural Convergence Rate Estimates in the Lindeberg Theorem
Previous Article in Journal
Fuzzy Automata as Coalgebras
Previous Article in Special Issue
A Random Walk Model for Spatial Galaxy Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotically Normal Estimators for the Parameters of the Gamma-Exponential Distribution

by
Alexey Kudryavtsev
1,2,* and
Oleg Shestakov
1,2,3,*
1
Faculty of Computational Mathematics and Cybernetics, M. V. Lomonosov Moscow State University, Moscow 119991, Russia
2
Moscow Center for Fundamental and Applied Mathematics, Moscow 119991, Russia
3
Institute of Informatics Problems, Federal Research Center “Computer Science and Control” of the Russian Academy of Sciences, Moscow 119333, Russia
*
Authors to whom correspondence should be addressed.
Submission received: 24 December 2020 / Revised: 22 January 2021 / Accepted: 28 January 2021 / Published: 30 January 2021
(This article belongs to the Special Issue Analytical Methods and Convergence in Probability with Applications)

Abstract

:
Currently, much research attention has focused on generalizations of known mathematical objects in order to obtain adequate models describing real phenomena. An important role in the applied theory of probability and mathematical statistics is the gamma class of distributions, which has proven to be a convenient and effective tool for modeling many real processes. The gamma class is quite wide and includes distributions that have useful properties such as, for example, infinite divisibility and stability, which makes it possible to use distributions from this class as asymptotic approximations in various limit theorems. One of the most important tasks of applied statistics is to obtain estimates of the parameters of the model distribution from the available real data. In this paper, we consider the gamma-exponential distribution, which is a generalization of the distributions from the gamma class. Estimators for some parameters of this distribution are given, and the asymptotic normality of these estimators is proven. When obtaining the estimates, a modified method of moments was used, based on logarithmic moments calculated on the basis of the Mellin transform for the generalized gamma distribution. On the basis of the results obtained, asymptotic confidence intervals for the estimated parameters are constructed. The results of this work can be used in the study of probabilistic models based on continuous distributions with an unbounded non-negative support.

1. Introduction

Estimating unknown parameters is an important problem in applied mathematical statistics. At the same time, in order to improve the consistency of mathematical models and analyzed real processes, researchers consider increasingly complex mathematical abstractions. Many models are traditionally described using continuous distributions with unbounded non-negative supports. For these purposes, special cases of the generalized gamma distribution and beta prime distributions are usually used. The paper considers the problem of estimating the parameters of the distribution proposed in [1], which is closely related to the listed popular distributions.
Definition 1.
We say that the random variable ζ has the gamma-exponential distribution G E ( r , ν , s , t , δ ) with the parameters of bent 0 r < 1 , shape ν 0 , concentration s , t > 0 , and scale δ > 0 if its density at x > 0 is
g E ( x ) = | ν | x t ν 1 δ t ν Γ ( s ) Γ ( t ) Ge r , t r + s ( ( x / δ ) ν ) ,
where E = ( r , ν , s , t , δ ) , and Ge α , β ( x ) is the gamma-exponential function [2]:
Ge α , β ( x ) = k = 0 x k k ! Γ ( α k + β ) , x R , 0 α < 1 , β > 0 .
Function (2) generalizes to the case β 1 , the transformation introduced by Le Roy [3] to study generating functions of a special form. In addition, Function (2) can be considered (under some assumptions) as a special case of the Srivastava–Tomovski function [4], that generalizes the Mittag–Leffler function [5].
In [1] it was shown that the distribution (1) adequately describes Bayesian balance models [6]. This is primarily due to the fact that the distribution with the density (1) can be represented as a scaled mixture of two random variables with generalized gamma distributions.
In turn, the generalized gamma distribution G G ( v , q , θ ) with the density
f ( x ) = | v | x v q 1 e ( x / θ ) v θ v q Γ ( q ) , v 0 , q > 0 , θ > 0 , x > 0 ,
proposed in 1925 by the Italian economist L. Amoroso [7] and often associated with E. W. Stacy [8], who considered in 1962 a particular case of the Amoroso distribution, has proven its validity in many applied problems that use continuous distributions with unbounded non-negative support for modeling. The class of distributions (3) is wide enough and includes exponential distribution; χ 2 -distribution; Erlang distribution; gamma distribution; semi-normal distribution, or distribution of the maximum of the Brownian motion process; Rayleigh distribution; Maxwell–Boltzmann distribution; χ -distribution; Nakagami m-distribution; Wilson–Hilferty distribution; Weibull–Gnedenko distribution and many others, including scaled and inverse analogs of the above.
The problem of estimating parameters of the distribution (3), its special types, and mixtures has a rich history and is still relevant [9,10,11,12].
In [1], it was shown that the gamma-exponential distribution has the following properties.
Lemma 1.
1. Let the independent random variables λ and μ have the distributions G G ( v , q , θ ) and G G ( u , p , α ) , u v > 0 , respectively. Then the distribution of λ coincides with G E ( 0 , v , · , q , θ ) ; the distribution of λ / μ for | u | > | v | coincides with G E ( v / u , v , p , q , θ / α ) ; the distribution of λ / μ for | v | > | u | coincides with G E ( u / v , u , q , p , θ / α ) .
2. For 0 < r < 1 , the density g E ( x ) , E = ( r , ν , s , t , δ ) , coincides with the density of the ratio of independent random variables with generalized gamma distributions G G ( ν , t , δ ) and G G ( ν / r , s , 1 ) .
The possibility of representing the gamma-exponential distribution as a ratio of random variables having the generalized gamma distribution allows it to be used in a wide range of applied problems.
Thus, in demography, the infant mortality rate is defined as the ratio of the number of deaths under the one year age to the number of births over a certain period of time, and the divorce index is defined as the ratio of the total divorce rate to the total marriage rate [13]. In physics, the transformation ratio is the ratio of the output voltage to the input voltage [14], and the universal Kirchhoff function is the ratio of the emissivity to the absorptivity of the body. In queuing theory, the ratio of the intensity of the incoming flow to the intensity of the service determines the system load factor [15]. When simulating emergency situations, the fire hazard of an object is determined by the ratio of the fire threat to the fire protection factor [16]. In reliability theory, the expected uptime is represented as the ratio of the average uptime to the average recovery time [17]. A number of examples can be continued. Each of these characteristics can be considered as the system balance index [6]. The application of a randomized Bayesian approach to the described models makes it possible to study the characteristics of the balance index as a scale mixture of probabilistic laws.
In addition, the five-parameter gamma-exponential distribution can be used to model a wide range of real phenomena, due to the wide variety of its possible densities (see Figure 1).
In practice, the researcher deals with observable quantities that reflect the evolution of the analyzed real process. In relation to these quantities, some model assumptions are made about the form of their distribution. The problem of estimating unknown parameters from real data also arises in the case of modeling a real process using the gamma-exponential distribution. Due to the representation of the density g E ( x ) in terms of a special gamma-exponential function (2), the maximum likelihood method seems to be too complicated. The same can be said about the direct method of moments. For this reason, in [18] it was proposed to estimate the parameters of the gamma-exponential distribution using a modified method based on logarithmic moments.

2. Estimators for the Parameters of the Gamma-Exponential Distribution

Let us introduce the estimators for the parameters of the gamma-exponential distribution. To do this, we define the di-gamma function ψ ( z ) = Γ ( z ) / Γ ( z ) and the functions
R ( x ) = ν 2 x ψ ( t ) ψ ( s ) ;
D r ( x , y ) = exp x ψ ( t ) ν + ψ ( s ) ν 2 y ψ ( t ) ν 2 ψ ( s ) ;
V ( x ) = ψ ( t ) r 3 ψ ( s ) x 3 ;
D ν ( x , y ) = exp x ( ψ ( t ) r ψ ( s ) ) y ψ ( t ) r 3 ψ ( s ) 3 .
Let us introduce a notation for the sample logarithmic moments of the random variable ζ :
L k ( X ) = 1 n i = 1 n ln k X i ,
where X = ( X 1 , , X n ) is a sample from the distribution of ζ .
In [18], the following statements were proved.
Lemma 2.
For fixed parameters ν, t and s of the distribution G E ( r , ν , s , t , δ ) the estimators
r ^ ( X ) = R ( L 2 ( X ) L 1 2 ( X ) )
and
δ ^ r ( X ) = D r ( L 1 ( X ) , L 2 ( X ) L 1 2 ( X ) )
for the parameters r and δ are strongly consistent.
Lemma 3.
For fixed parameters r, t and s of the distribution G E ( r , ν , s , t , δ ) the estimators
ν ^ ( X ) = V ( L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) )
and
δ ^ ν ( X ) = D ν ( L 1 ( X ) , L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) )
for the parameters ν and δ are strongly consistent.

3. Auxiliary Relations and Statements

In what follows, we need the derivatives of Functions (4)–(7):
R 1 ( x ) d R d x ( x ) = ν 2 2 ψ ( s ) ψ ( s ) ν 2 x ψ ( t ) ;
D r , 1 ( x , y ) D r x ( x , y ) = exp x ψ ( t ) ν + ψ ( s ) ν 2 y ψ ( t ) ν 2 ψ ( s ) ;
D r , 2 ( x , y ) D r y ( x , y ) = ψ ( s ) 2 ψ ( s ) ν 2 ψ ( s ) ν 2 y ψ ( t ) exp x ψ ( t ) ν + ψ ( s ) ν 2 y ψ ( t ) ν 2 ψ ( s ) ;
V 1 ( x ) d V d x ( x ) = ψ ( t ) r 3 ψ ( s ) 3 x 2 x 2 ( ψ ( t ) r 3 ψ ( s ) ) 2 3 ;
D ν , 1 ( x , y ) D ν x ( x , y ) = exp x ( ψ ( t ) r ψ ( s ) ) y ψ ( t ) r 3 ψ ( s ) 3 ;
D ν , 2 ( x , y ) D ν y ( x , y ) = ψ ( t ) r ψ ( s ) 3 ( ψ ( t ) r 3 ψ ( s ) ) ( ψ ( t ) r 3 ψ ( s ) ) 2 y 2 3 ×
× exp x ( ψ ( t ) r ψ ( s ) ) y ψ ( t ) r 3 ψ ( s ) 3 .
We also need some moment characteristics of the gamma-exponential distribution (1).
Consider the Mellin transform
M ζ ( z ) = 0 x z d F ζ ( x ) , z C .
We use Lemma 1 and the representation ζ = d λ / μ , where the independent random variables λ and μ have distributions G G ( ν , t , δ ) and G G ( ν / r , s , 1 ) , respectively. Since for λ G G ( ν , t , δ ) the Mellin transform has the form
M λ ( z ) = δ z Γ ( t ) Γ t + z ν , t + Re ( z ) ν > 0 ,
for the ratio of λ G G ( ν , t , δ ) to μ G G ( ν / r , s , 1 )
M λ / μ ( z ) = δ z Γ ( t ) Γ ( s ) Γ t + z ν Γ s r z ν , t + Re ( z ) ν > 0 , s r Re ( z ) ν > 0 ,
from where we get the form of the characteristic function of the logarithm of ζ :
E e i y ln ζ = δ i y Γ ( t ) Γ ( s ) Γ t + i y ν Γ s i r y ν , y R .
Differentiating the relation (18) four times, we obtain (the parameters t and s are assumed to be fixed)
μ 1 ( r , ν , δ ) E ln ζ = ν ln δ + ψ ( t ) r ψ ( s ) ν ;
μ 2 ( r , ν , δ ) E ln 2 ζ = ν ln δ + ψ ( t ) r ψ ( s ) 2 ν 2 + ψ ( t ) + r 2 ψ ( s ) ν 2 ;
μ 3 ( r , ν , δ ) E ln 3 ζ = ν ln δ + ψ ( t ) r ψ ( s ) 3 ν 3 +
+ 3 ( ψ ( t ) + r 2 ψ ( s ) ) ν ln δ + ψ ( t ) r ψ ( s ) ν 3 + ψ ( t ) r 3 ψ ( s ) ν 3 ;
μ 4 ( r , ν , δ ) E ln 4 ζ = [ ν ln δ + ψ ( t ) r ψ ( s ) ] 4 ν 4 +
+ 6 ( ψ ( t ) + r 2 ψ ( s ) ) [ ν ln δ + ψ ( t ) r ψ ( s ) ] 2 ν 4 +
+ 4 ( ψ ( t ) r 3 ψ ( s ) ) [ ν ln δ + ψ ( t ) r ψ ( s ) ] ν 4 +
+ 3 ( ψ ( t ) + r 2 ψ ( s ) ) 2 ν 4 + ψ ( t ) + r 4 ψ ( s ) ν 4 ;
σ 1 2 ( r , ν ) D ln ζ = ψ ( t ) + r 2 ψ ( s ) ν 2 ;
σ 2 2 ( r , ν , δ ) D ln 2 ζ = 4 ( ψ ( t ) + r 2 ψ ( s ) ) [ ν ln δ + ψ ( t ) r ψ ( s ) ] 2 ν 4 +
+ 4 ( ψ ( t ) r 3 ψ ( s ) ) [ ν ln δ + ψ ( t ) r ψ ( s ) ] ν 4 + 2 ( ψ ( t ) + r 2 ψ ( s ) ) 2 ν 4 + ψ ( t ) + r 4 ψ ( s ) ν 4 .
Further arguments are based on the following statements [19].
Lemma 4.
In R n , the random vector X n converges in distribution to the random vector X if and only if each linear combination of the components of X n converges in distribution to the same linear combination of the components of X.
Lemma 5.
Suppose that in R k
n ( T n 1 , , T n k ) N μ , Σ , n ,
with Σ a covariance matrix. Let g ( t ) = g ( t 1 , , t k ) be a real-valued function with a nonzero differential at t = μ . Put
d = g t 1 | t = μ , , g t k | t = μ .
Then n g ( T n 1 , , T n k ) N ( g ( μ ) , d Σ d T ) .

4. Asymptotic Normality of the Estimators for the Parameters of the Gamma-Exponential Distribution

Let us formulate the statements about the asymptotic normality of the estimators (8)–(11). Let us fix the parameters of shape ν and concentration t and s. The following statements hold.
Theorem 1.
The estimator (8) for the unknown parameter r is asymptotically normal:
n r ^ ( X ) r ν 4 σ 2 2 ( r , ν , δ ) 4 r 2 ( ψ ( s ) ) 2 N ( 0 , 1 ) , n ,
where σ 2 2 ( r , ν , δ ) is given by (24).
Proof of Theorem 1.
The statistic L 2 ( X ) L 1 2 ( X ) is a sample logarithmic variance that is representable as a sum of independent identically distributed random variables and has the mean σ 1 2 ( r , ν ) and the variance σ 2 2 ( r , ν , δ ) / n . Thus, when n
n L 2 ( X ) L 1 2 ( X ) σ 1 2 ( r , ν ) σ 2 2 ( r , ν , δ ) N 0 , 1 ,
where σ 1 2 ( r , ν ) is defined in (23). In addition, at the point x = σ 1 2 ( r , ν ) the function R ( x ) , defined in (4), has a nonzero derivative R 1 ( x ) , defined in (12).
Next, we use Lemma 5. Since in this case we consider a one-dimensional space ( k = 1 ), in terms of the notation of Lemma 5
d Σ d T = R 1 σ 1 2 ( r , ν ) 2 σ 2 2 ( r , ν , δ ) ,
and it follows that
n R ( L 2 ( X ) L 1 2 ( X ) ) R ( σ 1 2 ( r , ν ) ) R 1 σ 1 2 ( r , ν ) 2 σ 2 2 ( r , ν , δ ) N ( 0 , 1 ) ,
which concludes the proof. □
Repeating the reasoning from [19], it can be shown that Lemma 2 and the continuity of the function (24) in r and δ and the function (12) at x = σ 1 2 ( r , ν ) imply the following statement.
Corollary 1.
When n
n r ^ ( X ) r R 1 ( L 2 ( X ) L 1 2 ( X ) ) 2 σ 2 2 ( r ^ ( X ) , ν , δ ^ r ( X ) ) N ( 0 , 1 ) ,
where the functions R 1 and σ 2 2 are given by the relations (12) and (24), and the statistics r ^ ( X ) and δ ^ r ( X ) are defined in (8) and (9).
Consider the estimator for the parameter δ with fixed parameters ν , t and s. Let us introduce the notation
σ 12 ( r , ν ) = ψ ( t ) r 3 ψ ( s ) ν 3 ;
s 1 2 ( r , ν , δ ) = δ 2 σ 1 2 ( r , ν ) + ν δ 2 ψ ( s ) σ 12 ( r , ν ) r ψ ( s ) + ν 2 δ 2 ψ 2 ( s ) σ 2 2 ( r , ν , δ ) 4 r 2 ( ψ ( s ) ) 2 ,
where σ 1 2 ( r , ν ) , σ 2 2 ( r , ν , δ ) and σ 12 ( r , ν ) are given by the relations (23), (24) and (26), respectively.
Theorem 2.
The estimator (9) for the unknown parameter δ is asymptotically normal:
n δ ^ r ( X ) δ s 1 2 ( r , ν , δ ) N ( 0 , 1 ) , n ,
where s 1 2 ( r , ν , δ ) is given by (27).
Proof of Theorem 2.
Note that L 1 ( x ) is the sum of independent identically distributed random variables, therefore
n L 1 ( X ) μ 1 ( r , ν , δ ) σ 1 2 ( r , ν ) N 0 , 1 , n ,
where μ 1 ( r , ν , δ ) and σ 1 2 ( r , ν ) are given by (19) and (23). The statisitcs L 2 ( X ) L 1 2 ( X ) is also asymptotically normal, and (25) holds. Consider the covariance matrix
Σ = σ 1 2 ( r , ν ) σ 12 ( r , ν ) σ 12 ( r , ν ) σ 2 2 ( r , ν , δ ) .
Note that
cov ( L 1 ( X ) , L 2 ( X ) L 1 2 ( X ) ) = n 1 n 2 σ 12 ( r , ν ) = 1 n σ 12 ( r , ν ) + o ( n 1 ) , n .
The statistics L 1 ( X ) and L 2 ( X ) L 1 2 ( X ) together with any of their linear combinations have the property of asymptotic normality with corresponding limit means depending on ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) , and variances defined by the matrix Σ . Therefore, Lemma 4 implies the convergence of vectors
n ( L 1 ( X ) , L 2 ( X ) L 1 2 ( X ) ) N ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) , Σ .
In addition, the partial derivatives D r , 1 ( x , y ) and D r , 2 ( x , y ) , defined in (13) and (14), of the function D r ( x , y ) , defined in (5), are nonzero at ( x , y ) = ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) . Hence, Lemma 5 implies the convergence
n D r ( L 1 ( X ) , L 2 ( X ) L 1 2 ( X ) ) N ( D r ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) , d Σ d T ) ,
where
d = D r , 1 ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) , D r , 2 ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) .
The relation
d Σ d T = δ 2 σ 1 2 ( r , ν ) + ν δ 2 ψ ( s ) r ψ ( s ) σ 12 ( r , ν ) + ν 2 δ 2 ψ 2 ( s ) 4 r 2 ( ψ ( s ) ) 2 σ 2 2 ( r , ν , δ )
concludes the proof. □
Lemma 2 and the continuity of the function (27) in r and δ and functions (13) and (14) at ( x , y ) = ( μ 1 ( r , ν , δ ) , σ 1 2 ( r , ν ) ) imply the following statement.
Corollary 2.
When n
n δ ^ r ( X ) δ s 1 2 ( r ^ ( X ) , ν , δ ^ r ( X ) ) N ( 0 , 1 ) ,
where the function s 1 2 is given by (27), and the statistics r ^ ( X ) and δ ^ r ( X ) are defined in (8) and (9).
Let us fix the parameters r, t and s, and introduce the notation
σ 3 2 ( r , ν , δ ) D ln 3 ζ .
Remark 1.
The analytical form of the variance in (29) is obtained similarly to (23) and (24) by differentiating the characteristic function of the random variable ln ζ . The explicit form of this expression is not given due to its cumbersomeness.
Theorem 3.
The estimator (10) for the unknown parameter ν is asymptotically normal:
n ν ^ ( X ) ν ν 8 σ 3 2 ( r , ν , δ ) 9 ( ψ ( t ) r 3 ψ ( s ) ) 2 N ( 0 , 1 ) , n ,
where σ 3 2 ( r , ν , δ ) is given by (29).
Proof of Theorem 3.
Based on the form of statistics L k ( X ) as sums of random variables, we obtain
n L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) σ 12 ( r , ν ) σ 3 2 ( r , ν , δ ) N 0 , 1 , n ,
since
E ( L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) = n 2 3 n + 2 n 2 · ψ ( t ) r 3 ψ ( s ) ν 3 = σ 12 ( r , ν ) + o ( 1 )
and
D ln 3 X 1 3 ln X 1 L 2 ( X ) + 2 ln X 1 L 1 2 ( X ) = σ 3 2 ( r , ν , δ ) + o ( 1 ) ,
which obviously follows from the inequality
( D ξ D η ) 2 D ( ξ + η ) ( D ξ + D η ) 2 ,
valid for any random variables ξ and η with finite variances. In addition, at the point x = σ 12 ( r , ν ) the function V ( x ) , defined in (6), has a nonzero derivative V 1 ( x ) , defined in (15).
In terms of Lemma 5 for the one-dimensional case we obtain the relation
d Σ d T = V 1 σ 12 ( r , ν ) 2 σ 3 2 ( r , ν , δ ) ,
from which it follows that
n V ( L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) V ( σ 12 ( r , ν ) ) V 1 σ 12 ( r , ν ) 2 σ 3 2 ( r , ν , δ ) N ( 0 , 1 ) .
This concludes the proof. □
Lemma 3 and the continuity of the function (29) in ν and δ and the function (15) at x = σ 12 ( r , ν ) imply the following statement.
Corollary 3.
When n
n ν ^ ( X ) ν V 1 ( L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) 2 σ 3 2 ( r , ν ^ ( X ) , δ ^ ν ( X ) ) N ( 0 , 1 ) ,
where the functions V 1 and σ 3 2 are given by (15) and (29), and the statistics ν ^ ( X ) and δ ^ ν ( X ) are defined in (10) and (11).
Consider the estimator for the parameter δ with fixed parameters r, t and s. Let us introduce the notation
s 12 ( r , ν , δ ) = ψ ( t ) 6 ψ ( t ) [ ν ln δ + ψ ( t ) r ψ ( s ) ] 2 ν 4 +
+ r 4 ψ ( s ) 6 r 2 ψ ( s ) [ ν ln δ + ψ ( t ) r ψ ( s ) ] 2 ν 4 ;
s 2 2 ( r , ν , δ ) = δ 2 σ 1 2 ( r , ν ) 2 ν 2 δ 2 ( ψ ( t ) r ψ ( s ) ) s 12 ( r , ν , δ ) 3 ( ψ ( t ) r 3 ψ ( s ) ) + ν 4 δ 2 ( ψ ( t ) r ψ ( s ) ) 2 σ 3 2 ( r , ν , δ ) 9 ( ψ ( t ) r 3 ψ ( s ) ) 2 ,
where σ 1 2 ( r , ν ) , σ 3 2 ( r , ν , δ ) and s 12 ( r , ν ) are given by (23), (29) and (32), respectively.
Theorem 4.
The estimator (11) for the unknown parameter δ is asymptotically normal:
n δ ^ ν ( X ) δ s 2 2 ( r , ν , δ ) N ( 0 , 1 ) , n ,
where s 2 2 ( r , ν , δ ) is given by (33).
Proof of Theorem 4.
Note that when n the relations (28) and (31) hold. Consider the covariance matrix
Σ = σ 1 2 ( r , ν ) s 12 ( r , ν , δ ) s 12 ( r , ν , δ ) σ 3 2 ( r , ν , δ ) .
Since
cov ( L 1 ( X ) , L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) =
= μ 4 ( r , ν , δ ) n 4 μ 3 ( r , ν , δ ) μ 1 ( r , ν , δ ) n 3 μ 2 2 ( r , ν , δ ) n +
+ 4 μ 2 ( r , ν , δ ) μ 1 2 ( r , ν , δ ) n + 2 μ 1 4 ( r , ν , δ ) n + 2 D L 1 2 ( X ) + o n 1 , n ,
and
n 4 E L 1 4 ( X ) = n μ 4 ( r , ν , δ ) + 4 n ( n 1 ) μ 3 ( r , ν , δ ) μ 1 ( r , ν , δ ) +
+ 3 n ( n 1 ) μ 2 2 ( r , ν , δ ) + 3 n ( n 1 ) ( n 2 ) μ 2 ( r , ν , δ ) μ 1 2 ( r , ν , δ ) +
+ ( n 4 3 n ( n 1 ) ( n 2 ) 7 n ( n 1 ) n ) μ 1 4 ( r , ν , δ ) ,
where the moments μ k ( r , ν , δ ) , k = 1 , , 4 , are given by (19)–(22),
cov ( L 1 ( X ) , L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) = 1 n s 12 ( r , ν , δ ) + o n 1 .
The statistics L 1 ( X ) and L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) together with any of their linear combinations have the property of asymptotic normality with corresponding limit means depending on ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) , and variances defined by the matrix Σ . Therefore, Lemma 4 implies the convergence of vectors
n ( L 1 ( X ) , L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) N ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) , Σ ,
where μ 1 ( r , ν , δ ) and σ 12 ( r , ν ) are defined in (19) and (26). In addition, the partial derivatives D ν , 1 ( x , y ) and D ν , 2 ( x , y ) , defined in (16) and (17) of the function D ν ( x , y ) , defined in (7), are nonzero at ( x , y ) = ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) . Hence, by Lemma 5
n D ν ( L 1 ( X ) , L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) N ( D ν ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) , d Σ d T ) ,
where
d = D ν , 1 ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) , D ν , 2 ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) .
The realtion
d Σ d T = δ 2 σ 1 2 ( r , ν ) 2 ν 2 δ 2 ( ψ ( t ) r ψ ( s ) ) s 12 ( r , ν , δ ) 3 ( ψ ( t ) r 3 ψ ( s ) ) + ν 4 δ 2 ( ψ ( t ) r ψ ( s ) ) 2 σ 3 2 ( r , ν , δ ) 9 ( ψ ( t ) r 3 ψ ( s ) ) 2
concludes the proof. □
Lemma 3 and the continuity of the function (33) in ν and δ and the functions (16) and (17) at ( x , y ) = ( μ 1 ( r , ν , δ ) , σ 12 ( r , ν ) ) imply the following statement.
Corollary 4.
When n
n δ ^ ν ( X ) δ s 2 2 ( r , ν ^ ( X ) , δ ^ ν ( X ) ) N ( 0 , 1 ) ,
where the function s 2 2 is given by (33), and the statistics ν ^ ( X ) and δ ^ ν ( X ) are defined in (10) and (11).

5. Confidence Intervals

On the basis of Corollaries 1–4, asymptotic confidence intervals for unknown parameters of the gamma-exponential distribution can be constructed.
By u γ we denote the ( 1 + γ ) / 2 -quantile of the standard normal distribution.
Corollary 5.
The asymptotic confidence interval with the confidence level γ based on the estimator (8) for the unknown parameter r has the form
( S r ( X ) , T r ( X ) ) = r ^ ( X ) u γ n A r ( X ) , r ^ ( X ) + u γ n A r ( X ) ,
where
A r ( X ) = R 1 ( L 2 ( X ) L 1 2 ( X ) ) 2 σ 2 2 ( r ^ ( X ) , ν , δ ^ r ( X ) ) ,
the functions R 1 and σ 2 2 are given by (12) and (24), and the statistics r ^ ( X ) and δ ^ r ( X ) are defined in (8) and (9).
Corollary 6.
The asymptotic confidence interval with the confidence level γ based on the estimator (9) for the unknown parameter δ has the form
( S δ r ( X ) , T δ r ( X ) ) = δ ^ r ( X ) u γ n A δ r ( X ) , δ ^ r ( X ) + u γ n A δ r ( X ) ,
where
A δ r ( X ) = s 1 2 ( r ^ ( X ) , ν , δ ^ r ( X ) ) ,
the function s 1 2 is given by (27), and the statistics r ^ ( X ) and δ ^ r ( X ) are defined in (8) and (9).
Corollary 7.
The asymptotic confidence interval with the confidence level γ based on the estimator (10) for the unknown parameter ν has the form
( S ν ( X ) , T ν ( X ) ) = ν ^ ( X ) u γ n A ν ( X ) , ν ^ ( X ) + u γ n A ν ( X ) ,
where
A ν ( X ) = V 1 ( L 3 ( X ) 3 L 1 ( X ) L 2 ( X ) + 2 L 1 3 ( X ) ) 2 σ 3 2 ( r , ν ^ ( X ) , δ ^ ν ( X ) ) ,
the functions V 1 and σ 3 2 are given by (15) and (29), and the statistics ν ^ ( X ) and δ ^ ν ( X ) are defined in (10) and (11).
Corollary 8.
The asymptotic confidence interval with the confidence level γ based on the estimator (11) for the unknown parameter δ has the form
( S δ ν ( X ) , T δ ν ( X ) ) = δ ^ ν ( X ) u γ n A δ ν ( X ) , δ ^ ν ( X ) + u γ n A δ ν ( X ) ,
where
A δ ν ( X ) = s 2 2 ( r , ν ^ ( X ) , δ ^ ν ( X ) ) ,
the function s 2 2 is given by (33), and the statistics ν ^ ( X ) and δ ^ ν ( X ) are defined in (10) and (11).
Let us illustrate the results of Corollaries 5–8, using an example of a model sample from a gamma-exponential distribution with the given parameters r, ν , s, t and δ . The confidence level is γ = 0.95 .
Table 1, Table 2, Table 3 and Table 4 show the values of the estimates (8) and (9) for the parameters r and δ and the estimates (10) and (11) of the parameters ν and δ with the corresponding asymptotic confidence intervals obtained from the sample size n.
Table 5 and Table 6 show the proportions of the true values of parameters that fall into the asymptotic confidence intervals for 1000 runs for the sample size n.

6. Discussion

The paper considers estimates and asymptotic confidence intervals for the parameters of the gamma-exponential distribution, represented as a scale mixture of generalized gamma distributions. Distributions from the gamma class have a rich history of modeling applications in many areas of knowledge. Back in the second half of the 19th century, the Rayleigh distribution was used to describe the resulting amplitude in the problem of summation of oscillations with random phases, and the Maxwell–Boltzmann distribution was used to statistically describe the behavior of the parameters of ideal gas particles. In the 1920s, the Italian economist L. Amoroso in the study of dynamical equilibrium theory proposed a type of distribution, which is now known as the generalized gamma distribution. The use of the generalized gamma distribution and its special cases for describing various models is still relevant. For example, in [9,10,11,12], it is proposed to use distributions from the gamma class in problems of processing radar signals and images, for studying the strength of materials and reliability of equipment, as well as for estimating the concentration of NO2 in industrial areas and studying the remission period of cancer patients. A number of other examples can be cited. The gamma-exponential distribution considered in the paper generalizes the Amoroso distribution. Therefore, it can be argued that the results of the article will be in demand when studying various models that describe real processes using distributions with a non-negative unbounded support.

Author Contributions

Conceptualization, A.K. and O.S.; methodology, A.K. and O.S.; formal analysis, A.K. and O.S.; investigation, A.K. and O.S.; writing—original draft preparation, A.K. and O.S.; writing—review and editing, A.K. and O.S.; supervision, A.K. and O.S.; funding acquisition, O.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Ministry of Science and Higher Education of the Russian Federation, project No. 075-15-2020-799.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kudryavtsev, A.A. On the representation of gamma-exponential and generalized negative binomial distributions. Inform. Appl. 2019, 13, 78–82. [Google Scholar]
  2. Kudryavtsev, A.A.; Titova, A.I. Gamma-exponential function in Bayesian queueing models. Inform. Appl. 2017, 11, 104–108. [Google Scholar]
  3. Le Roy, É. Sur les séries divergentes et les fonctions définies par un développement de Taylor. Ann. Fac. Sci. Toulouse 1900, 2, 317–384. [Google Scholar] [CrossRef] [Green Version]
  4. Srivastava, H.M.; Tomovski, Ž. Fractional calculus with an integral operator containing a generalized Mittag-Leffler function in the kernel. Appl. Math. Comput. 2009, 211, 198–210. [Google Scholar] [CrossRef]
  5. Gorenlo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014. [Google Scholar]
  6. Kudryavtsev, A.A. Bayesian balance models. Inform. Appl. 2018, 12, 18–27. [Google Scholar]
  7. Amoroso, L. Ricerche intorno alla curva dei redditi. Ann. Mat. Pura Appl. 1925, 21, 123–159. [Google Scholar] [CrossRef]
  8. Stacy, E.W. A Generalization of the Gamma Distribution. Ann. Math. Stat. 1962, 33, 1187–1192. [Google Scholar] [CrossRef]
  9. Gao, G.; Ouyang, K.; Luo, Y.; Liang, S.; Zhou, S. Scheme of Parameter Estimation for Generalized Gamma Distribution and Its Application to Ship Detection in SAR Images. IEEE Trans. Geosci. Remote Sens. 2017, 55, 1812–1832. [Google Scholar] [CrossRef]
  10. Zhou, Y.; Zhu, H. Image Segmentation Using a Trimmed Likelihood Estimator in the Asymmetric Mixture Model Based on Generalized Gamma and Gaussian Distributions. Math. Probl. Eng. 2018, 2018, 3468967. [Google Scholar] [CrossRef] [Green Version]
  11. Iriarte, Y.A.; Varela, H.; Gómez, H.J.; Gómez, H.W. A Gamma-Type Distribution with Applications. Symmetry 2020, 12, 870. [Google Scholar] [CrossRef]
  12. Rivera, P.A.; Barranco-Chamorro, I.; Gallardo, D.I.; Gómez, H.W. Scale Mixture of Rayleigh Distribution. Mathematics 2020, 8, 1842. [Google Scholar] [CrossRef]
  13. Lancaster, H.O. Expectations of Life: A Study in the Demography, Statistics, and History of World Mortality; Springer: New York, NY, USA, 1990. [Google Scholar]
  14. Bakshi, V.U.; Bakshi, U.A. Basics of Electrical Engineering; Technical Publications: Pune, India, 2008. [Google Scholar]
  15. Bocharov, P.P.; DApice, C.; Pechinkin, A.V.; Salerno, S. Queueing Theory. Series “Modern Probability and Statistics”; VSP Publishing: Utrecht, The Netherlands, 2003. [Google Scholar]
  16. Bostrom, A.; French, S.P.; Gottlieb, S.J. Risk Assessment, Modeling and Decision Support; Springer: Berlin, Germany, 2008. [Google Scholar]
  17. Rausand, M.; Hoyland, A. System Reliability Theory: Models, Statistical Methods, and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  18. Kudryavtsev, A.A.; Shestakov, O.V. Method of logarithmic moments for estimating the gamma-exponential distribution parameters. Inform. Appl. 2020, 14, 49–54. [Google Scholar]
  19. Serfling, R.J. Approximation Theorems of Mathematical Statistics; John Wiley & Sons, Inc.: New York, NY, USA, 2002. [Google Scholar]
Figure 1. Gamma-exponential distribution densities for different values of parameters.
Figure 1. Gamma-exponential distribution densities for different values of parameters.
Mathematics 09 00273 g001
Table 1. Values of the estimates of the parameters r and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.4 , ν = 1.7 , s = 0.7 , t = 1.5 , δ = 1 .
Table 1. Values of the estimates of the parameters r and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.4 , ν = 1.7 , s = 0.7 , t = 1.5 , δ = 1 .
r ^ ( X ) S r ( X ) T r ( X ) δ ^ r ( X ) S δ r ( X ) T δ r ( X )
n = 100 0.5170.3440.6900.7820.6450.919
n = 1000 0.3910.3320.4490.9820.9241.041
n = 10 , 000 0.3990.3810.4180.9940.9751.012
Table 2. Values of the estimates of the parameters ν and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.4 , ν = 1.7 , s = 0.7 , t = 1.5 , δ = 1 .
Table 2. Values of the estimates of the parameters ν and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.4 , ν = 1.7 , s = 0.7 , t = 1.5 , δ = 1 .
ν ^ ( X ) S ν ( X ) T ν ( X ) δ ^ ν ( X ) S δ ν ( X ) T δ ν ( X )
n = 100 1.048 0.313 2.4090.7020.2071.197
n = 1000 1.7250.9772.4740.9810.8471.115
n = 10 , 000 1.7181.4801.9560.9970.9531.040
Table 3. Values of the estimates of the parameters r and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.5 , ν = 0.7 , s = 1.5 , t = 0.6 , δ = 1 .
Table 3. Values of the estimates of the parameters r and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.5 , ν = 0.7 , s = 1.5 , t = 0.6 , δ = 1 .
r ^ ( X ) S r ( X ) T r ( X ) δ ^ r ( X ) S δ r ( X ) T δ r ( X )
n = 100 0.840 0.637 2.3180.8270.4401.214
n = 1000 0.445 0.373 1.2641.0570.9141.199
n = 10 , 000 0.5070.2790.7351.0641.0181.110
Table 4. Values of the estimates of the parameters ν and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.5 , ν = 0.7 , s = 1.5 , t = 0.6 , δ = 1 .
Table 4. Values of the estimates of the parameters ν and δ with the corresponding asymptotic confidence intervals for the true values of the parameters r = 0.5 , ν = 0.7 , s = 1.5 , t = 0.6 , δ = 1 .
ν ^ ( X ) S ν ( X ) T ν ( X ) δ ^ ν ( X ) S δ ν ( X ) T δ ν ( X )
n = 100 0.6680.1901.1470.902 0.688 2.493
n = 1000 0.7130.5550.8711.0160.4991.533
n = 10 , 000 0.7000.6520.7491.0620.8891.234
Table 5. Proportions of true parameter values falling into asymptotic confidence intervals for 1000 runs for the sample size n with true parameter values r = 0.4 , ν = 1.7 , s = 0.7 , t = 1.5 , δ = 1 .
Table 5. Proportions of true parameter values falling into asymptotic confidence intervals for 1000 runs for the sample size n with true parameter values r = 0.4 , ν = 1.7 , s = 0.7 , t = 1.5 , δ = 1 .
r δ r ν δ ν
n = 100 0.9540.9260.8490.878
n = 1000 0.9260.9180.9690.897
n = 10 , 000 0.9200.9050.9770.946
Table 6. Proportions of true parameter values falling into asymptotic confidence intervals for 1000 runs for the sample size n with true parameter values r = 0.5 , ν = 0.7 , s = 1.5 , t = 0.6 , δ = 1 .
Table 6. Proportions of true parameter values falling into asymptotic confidence intervals for 1000 runs for the sample size n with true parameter values r = 0.5 , ν = 0.7 , s = 1.5 , t = 0.6 , δ = 1 .
r δ r ν δ ν
n = 100 0.5100.4800.9930.976
n = 1000 0.7720.7380.9900.996
n = 10 , 000 0.9730.9160.9930.999
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kudryavtsev, A.; Shestakov, O. Asymptotically Normal Estimators for the Parameters of the Gamma-Exponential Distribution. Mathematics 2021, 9, 273. https://0-doi-org.brum.beds.ac.uk/10.3390/math9030273

AMA Style

Kudryavtsev A, Shestakov O. Asymptotically Normal Estimators for the Parameters of the Gamma-Exponential Distribution. Mathematics. 2021; 9(3):273. https://0-doi-org.brum.beds.ac.uk/10.3390/math9030273

Chicago/Turabian Style

Kudryavtsev, Alexey, and Oleg Shestakov. 2021. "Asymptotically Normal Estimators for the Parameters of the Gamma-Exponential Distribution" Mathematics 9, no. 3: 273. https://0-doi-org.brum.beds.ac.uk/10.3390/math9030273

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