Many researchers measure the uncertainty of a random variable using quantile-based entropy techniques. These techniques are useful in engineering applications and have some exceptional characteristics than their distribution function method. Considering order statistics, the key focus of this article is to propose new quantile-based Mathai-Haubold entropy and investigate its characteristics. The divergence measure of the Mathai-Haubold is also considered and some of its properties are established. Further, based on order statistics, we propose the residual entropy of the quantile-based Mathai-Haubold and some of its property results are proved. The performance of the proposed quantile-based Mathai-Haubold entropy is investigated by simulation studies. Finally, a real data application is used to compare our proposed quantile-based entropy to the existing quantile entropies. The results reveal the outperformance of our proposed entropy to the other entropies.
The order statistics are considered in a varied scope of complicated problems, including characterization of a probability distribution, quality control, robust statistical estimation and identifying outliers, analysis of a censored sample, the goodness of fit-tests, etc. Based on order statistics, the usage of the recurrence relationships for moments is well recognized by many researchers (see, for instance, Arnold et al. [1], Malik et al. [2]). For an enhancement, many recurrence relations and identities for the order statistics moments originating from numerous particular continuous probability distributions (i.e., gamma, Cauchy, normal, logistic, and exponential) have been reviewed by Samuel et al. [3] and Arnold et al. [1].
Based on a random sample of X1,X2,…,Xn, let the corresponding order statistics to be X1:n≤X2:n≤…≤Xn:n. Then, as in David [4] and Arnold et al. [1], the density of Xr:n,1≤r≤n, is
fr:n(x)=Cr:n{[F(x)]r−1[1−F(x)]n−rf(x)},0<x<∞,
with Cr:n=n!(r−1)!(n−r)!. Eq. (1) can be used to determine the smallest (when r=1) and largest (when r=n) probability density functions; they, respectively, are f(1,n)(x)=n[1−F(x)]n−1f(x) and f(n,n)(x)=n[F(x)]n−1f(x). The corresponding distribution functions are obtained, respectively, by F(1,n)(x)=1−[1−F(x)]n and F(n,n)(x)=[F(x)]n.
Shannon [5] was the first author who introduced the entropy idea for a random variable (r.v.) X in the field of information theory and defined it as
H(X)=−∫0∞fX(x)logfX(x)dx.
Here, fX(x) indicates the pdf of the r.v. X. Based on the Shannon entropy measure, Mathi et al. [6] (now onwards M-H entropy) developed a generalized version of Eq. (2) and defined it as
Mα(X)=1α−1(∫0∞[f(x)]2−αdx−1),0<α<2,α≠1.
When α→1 the M-H entropy measure Mα(X) will be reduced to the Shannon entropy measure defined in Eq. (2).
Mathai et al. [6] and Sebastian [7] discussed the main property allied with Eq. (3). In other words, applying the maximum entropy and using its normalization version together with energy restrictions will result in the well-recognized pathway-model as provided by Mathai [8]. However, this model contains many special cases of familiar probability distributions.
Theoretical surveys and applications employing the measurement information are distributional dependents, and they may be found to be not appropriate in circumstances once the distribution is analytically not tractable. Hence, utilizations of quantile function are considered as an alternative method, where
Q(u)=F−1(u)=inf{x/F(x)≥u},0≤u≤1.
We refer the readers to Nair et al. [9] and Sunoj et al. [10] and references therein for more details about quantile function. Recently, Sunoj et al. [11] studied Shannon entropy and as well as its residual and introduced quantile versions of them defined as
H=∫01logq(p)dp
and
H(u)=H(X;Q(u))=log(1−u)+11−u∫u1logq(p)dp
respectively, where q(u)=dQ(u)du denotes the density of quantile function. If we define the quantile density function by fQ(u)=f(Q(u)), then, we obtain
q(u)f(Q(u))=1
For Shannon past entropy, Sunoj et al. [11] also introduced its quantile version and defined it as
H¯(u)=H(X;Q(u))=log(u)+1u∫0ulogq(p)dp
In the present paper, we work with the order statistics, propose the quantile-based version of M-H entropy and discuss its properties. The M-H divergence measure is also considered and we establish some of its distribution free properties. In addition, we introduce the version of the quantile-based residual for the M-H entropy and prove some characterization results. To the best of our knowledge, the results presented here, treat a research gap that has not been addressed or studied systematically by others, which was the primary motivation of our paper.
The paper is outlined as follows. Section 2 is devoted to the construction of our quantile-based M-H entropy and its properties. Next, expressions for the quantile-based version of M-H entropy for some life-time distributions are presented in Section 3. A quantile-based generalized divergence measure of rth order statistics is given in Section 4. Quantile Residual Entropy of M-H for rth order statistics and also for some lifetime models are introduced in Section 5. Characterization theorems based on M-H Quantile Residual Entropy are presented in Section 6. In Section 7, simulation studies for investigating the performance of our proposed quantiles and real data life application are presented. Our conclusion is stated in Section 8.
Quantile Based M-H Entropy of <inline-formula id="ieqn-20"><mml:math id="mml-ieqn-20"><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mrow><mml:mi mathvariant="bold-italic">t</mml:mi><mml:mi mathvariant="bold-italic">h</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> Order Statistics
Wong et al. [12], Park [13], Ebrahimi et al. [14] and Baratpour et al. [15] are the authors who discuss in detail the aspects of information-theoretic based on order statistics. Paul et al. [16] considered the M-H entropy and, based on record values, studied some of its essential properties. For rth order statistics Xr;n, the M-H entropy is defined as
Mα(Xr;n)=1α−1(∫0∞[fr;n(x)]2−αdx−1),0<α<2,α≠1,
where fr;n(x) is given in Eq. (1). Now, FQ(u)=u, then, the pdf of rth order statistics becomes
Similarly, based on the quantile and quantile density functions, we obtain the quantile-based for the M-H Entropy (MXr:nα) of rth order statistics for the following distributions.
Figs. 1–3 give the quantile version of M-H entropy plots of smallest order statistics under exponential, Pareto-I and uniform distributions, respectively.
Quantile M-H entropy plots of smallest order statistics (Exponential distribution)Quantile M-H entropy plots of smallest order statistics (Pareto-I Distribution)
For an increasing value of parameters α and λ, the entropy plot, based on the exponential distribution, increases. In the case of entropy plot under the Pareto-I distribution, the plot has an increasing (a decreasing) behaviour for different parameter combinations. The entropy plot under uniform distribution also has increasing behaviour for different parameter values. Tabs. 1--3 give entropy values when the parameters α and λ are varied.
Quantile M-H entropy plots of smallest order statistics (Uniform Distribution)Quantile M-H Entropy value for the smallest order statistics, <inline-formula id="ieqn-51"><mml:math id="mml-ieqn-51"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mo>,</mml:mo><mml:mi mathvariant="bold-italic">n</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">10</mml:mn></mml:mrow></mml:math></inline-formula> (Exponential distribution)
Parameter values
λ
0.1
0.6
0.9
1.5
3.5
9.4
α
0.1
−97.4989
−2.1656
−0.4055
0.5365
0.9962
1.0935
0.3
−78.5064
−2.3723
−0.4792
0.6280
1.2390
1.3932
0.5
−68.5431
−2.7999
−0.6127
0.7857
1.6593
1.9226
0.7
−70.7597
−3.8807
−0.9252
1.1412
2.6047
3.1316
0.9
−130.0233
−9.5090
−2.4892
2.8797
7.1963
9.0543
1.1
78.1511
7.5748
2.2014
−2.2955
−6.4061
−8.5229
1.3
15.1448
1.9384
0.6357
−0.5575
−1.7994
−2.5651
1.5
4.9570
0.8402
0.3190
−0.2037
−0.8241
−1.2824
1.7
1.6733
0.3835
0.1760
−0.0520
−0.3610
−0.6348
1.9
0.3654
0.1232
0.0742
0.0151
−0.0764
−0.1737
Quantile M-H Entropy value for smallest order statistics, <inline-formula id="ieqn-54"><mml:math id="mml-ieqn-54"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mo>,</mml:mo><mml:mi mathvariant="bold-italic">n</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">10</mml:mn></mml:mrow></mml:math></inline-formula> (Pareto-I Distribution)
Parameters
a = 0.5, b = 0.3, u = 0.5
a = 1.4, b = 0.3, u = 0.5
a = 2.5, b = 0.3, u = 0.5
a = 0.5, b = 0.9, u = 0.9
a = 0.5, b = 2.3, u = 0.9
a = 5.3, b = 0.3, u = 0.9
α
0.1
−0.1092
−1.9714
−4.0833
0.6571
0.9160
1.0191
0.3
−0.4596
−2.4535
−4.3969
0.5535
0.9748
1.1756
0.5
−0.5820
−2.3205
−3.7735
0.5093
1.0675
1.3857
0.7
−0.7351
−2.2075
−3.2602
0.4072
1.1251
1.6143
0.9
−1.4585
−2.7012
−3.4594
−0.2664
0.6531
1.4017
1.1
0.7963
−0.2600
−0.8087
2.0500
3.2353
4.3876
1.3
0.1063
−0.8077
−1.2110
1.4491
3.0038
4.8073
1.5
0.0449
−0.7779
−1.0855
1.5419
3.6622
6.5952
1.7
0.1869
−0.6428
−0.9049
2.0572
5.2940
10.6306
1.9
2.0158
0.1268
−0.3765
7.2936
18.4441
40.3418
Quantile M-H Entropy value for smallest order statistics, <inline-formula id="ieqn-56"><mml:math id="mml-ieqn-56"><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mo>,</mml:mo><mml:mi mathvariant="bold-italic">n</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">10</mml:mn></mml:mrow></mml:math></inline-formula> (Uniform Distribution)
Parameters
a = 0.9, b = 2.5, u = 0.5
a = 1.4, b = 2.5, u = 0.5
a = 2.3, b = 2.5, u = 0.5
a = 0.5, b = 0.6, u = 0.9
a = 0.5, b = 0.9, u = 0.9
a = 0.5, b = 1.5, u = 0.9
α
0.1
0.38325
0.09134
−3.61856
−7.71476
−1.42345
−0.64852
0.3
0.40051
0.09220
−2.97881
−5.73125
−1.28449
−0.61409
0.5
0.41886
0.09307
−2.47214
−4.32456
−1.16228
−0.58199
0.7
0.43837
0.09396
−2.06886
−3.31754
−1.05461
−0.55205
0.9
0.45913
0.09486
−1.74619
−2.58925
−0.95958
−0.52410
1.1
0.48122
0.09577
−1.48660
−2.05672
−0.87556
−0.49800
1.3
0.50475
0.09669
−1.27655
−1.66271
−0.80114
−0.47361
1.5
0.52982
0.09762
−1.10557
−1.36754
−0.73509
−0.45081
1.7
0.55654
0.09856
−0.96553
−1.14353
−0.67635
−0.42947
1.9
0.58504
0.09952
−0.85008
−0.97123
−0.62402
−0.40950
Clearly, we see from Tabs. 1–3 that the entropy values under exponential, Pareto-I and uniform distributions portray the same behaviour as discussed in the graphical plots.
Quantile-Based Generalized Divergence Measure of <inline-formula id="ieqn-58"><mml:math id="mml-ieqn-58"><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mrow><mml:mi mathvariant="bold-italic">t</mml:mi><mml:mi mathvariant="bold-italic">h</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> Order Statistics
Different measures deal with the dissimilarity or the distance between two probability distributions. Certainly, these measures are essential in theory, inferential statistics, applied statistics and data processing sciences, such as comparison, classification, estimation, etc.
Assume f and g are the density functions of the non-negative r.vs X and Y, respectively. The direct divergence of f from g is measured by the Kullback et al. [17] and is
D(f/g)=∫0∞f(x)logf(x)g(x)dx.
With order α, the divergence measure of the M-H or the relative-entropy of g concerning f is obtained by
Mα(f/g)=11−α{∫0∞f(x)(f(x)g(x))1−αdx−1},
where, for α→1, the expression in Eq. (8) reduces to Eq. (7) (see Kullback et al. [17]).
Theorem 4.1: The quantile-based-generalized divergence measure between the rth order statistics distribution and the primary distribution is a distribution-free.
Proof: From equation Eq. (8), we have
Mα(Xr:n/X)=11−α{∫0∞(fr:n(x)f(x))1−αfr:n(x)dx−1}.
Now, using the value of fr:n(x) in Eq. (9), we obtain Mα(Xr:n/X)=11−α{∫0∞F(2−α)(r−1)(x)(1−F(x))(n−r)(2−α)(β(r,n−r+1))2−αf(x)dx−1}.
Using the fact that q(u)f(Q(u))=1, we determine the quantile-based-generalized divergence measure between the distribution of rth order statistics and primary distribution, as
which is a distribution-free. Hence, the theorem is proved.
M-H Quantile Residual Entropy for <inline-formula id="ieqn-74"><mml:math id="mml-ieqn-74"><mml:msup><mml:mi mathvariant="bold-italic">r</mml:mi><mml:mrow><mml:mi mathvariant="bold-italic">t</mml:mi><mml:mi mathvariant="bold-italic">h</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula> Order Statistics
Entropy functions are very popular in the applications of finance and tectonophysics, machine learning, reliability theory, etc. However, in reliability and real-life applications, the life test time is truncated at a specific time, and in such situations, Eq. (2) is not an appropriate measure. Therefore, Shannon's entropy is not an adequate measure when we have knowledge about the component's current age, which can be used when determining its uncertainty. Ebrahimi [14] describes a more practical approach that considers the use of age, defined as
H(X;t)=−∫t∞f(x)F¯(t)logf(x)F¯(t)dx,
with F¯(t) indicates the survival-function. Note that for t=0, the Eq. (5) reduces to Eq. (4). Denote by X, a non-negative r.v., the unit’s life-time at time t. Then, the residual function of M-H [18] is
Mα(X;t)=1α−1{∫t∞f2−α(x)F¯2−α(t)dx−1},0<α<2,α≠1.
The M-H residual entropy for the rth order statistics is given by
where ′ denotes the differentiation w.r.t u. This equation involved an immediate connection between the q(u) and Mα(Xr;n;u), which implies that the quantile residual entropy function of M-H of rth order statistic leads to the unicity of the underlying distribution.
Next, we make the derivation of the quantile form of M-H residual entropy of the rth order statistic for some lifetime models.
(i) Govindarajulu’s Distribution
The quantile for the Govindarajulu distribution is
Q(u)=a{(b+1)ub−bub+1},
and the corresponding density is
q(u)=ab(b+1)(1−u)ub−1,0≤u≤1;a,b>0.
The quantile residual entropy function of M-H of rth order statistics for the distribution of Govindarajulu is
Similarly, based on the quantile and quantile density functions, we obtain the quantile-based residual M-H Entropy of rth order statistics for the following distributions.
(ii) Uniform DistributionQ(u)=a+(b−a)u and q(u)=(b−a),0≤u≤1;a<b.
Based on residual M-H quantile entropy of order statistics Mα(Xr;n,u), the following nonparametric classes of life distribution are defined.
Definition 5.1:X is said to have an increasing (a decreasing) M-H quantile entropy of order statistics if Mα(Xr;n,u) is increasing (decreasing) in u≥0.
The following lemma is useful in proving the results in monotonicity of Mα(Xr;n,u).
Lemma 5.1: Let f(u,x):R+2→R+ and g:R+→R+ be any two functions. If ∫u∞f(u,x)dx is increasing and g(u) is increasing (decreasing) in u, then ∫u∞f(u,x)g(x)dx is increasing (decreasing) in u, provided the existence of integrals.
Theorem 5.2: Let X be a non-negative and continuous r.v. with quantile QX(.) and density qX(.). Define Y=∅(X), where ∅(.) is nonnegative, increasing and convex(concave) function. Then,
(i) For 1<α<2, Mα(Yr;n,u) increases (or decreases) in u whenever Mα(Xr;n,u) increases (or decreases) in u.
(ii) For 0<α<1, Mα(Yr;n,u) increases (or decreases) in u whenever Mα(Xr;n,u) increases (or decreases) in u.
Proof: (i) The quantile density of Y is given by
g(QY(u))=1qY(u)=1qX(u)∅′(QX(u)).
Since 1<α<2 and ∅ is non-negative, increasing and convex (concave) function, the (∅′(QX(p))α−1 increases (or decreases) and it is also non-negative. Consequently, using Lemma 3.1, Eq. (11) is increasing (decreasing), which gives the proof of (i) of the Theorem. Similarly, 0<α<1, (∅′(QX(p))α−1 increases (or decreases) in p, because ∅ increases and it is convex. Consequently, Eq. (11) is decreasing (increasing) in u, which proves (ii) of the Theorem. The immediate application of Theorem 5.2 is given below:
Let X be an r.v. following the distribution of exponential and having a failure rate λ. Also, let Y=X1α,α>0. Therefore, Y follows Weibull distribution where Q(u)=λ−1α(−log(1−u))1α. The function ϕ(x)=x1α,x>0,α>0 is a convex (concave) if 1<α<2,(0<α<1). Hence, based on Theorem 5.2, the Weibull distribution is increasing (decreasing) M-H quantile entropy of order statistics if 1<α<2,(0<α<1).
Characterization Theorems Based on M-H Quantile Residual Entropy
This section provides some characterizations for the quantile M-H residual entropies of the smallest and largest order statistics. The corresponding quantile M-H residual entropy can be determined by substituting r=1 (for smallest) and r=n (for largest) in Eq. (10) and are, respectively, given by
Now, we define the hazard and the reversed hazard functions for the quantile version which are, respectively, corresponding to the well-recognized hazard rate and reversed hazard rate functions, as
K(u)=h(Q(u))=fQ(u)(1−u)=1(1−u)q(u) and K¯(u)=K¯(Q(u))=fQ(u)F(Q(u))=(uq(u))−1.
In numerous practical circumstances, the uncertainty is essentially not identified with the future. Therefore, it can likewise allude to the past. This thought empowered Crescenzo et al. [19] to build up the idea of past entropy on (0,t). If X denotes the life-time of a component, thus, the past entropy of X is obtained by
H0(X;t)=−∫0tf(x)F(t)logf(x)F(t)dx,
with F(t) is the cumulative distribution function. For t=0,Eq. (13) reduces to Eq. (2).
The quantile form of past M-H residual entropy of rth order statistic is determined by
Next, we state some properties based on quantile M-H residual entropy of the smallest order statistics.
Theorem 6.1: Let X1;n represents the first order statistics with survival and hazard quantile functions F¯1;n(x) and KX1;n(u). Then, the Mα(X1;n;u) is determined by
Mα(X1;n;u)=1α−1{C(KX1;n(u))1−α−1},0<α<2,α≠1,
if and only if:
(i) X follows an exponential distribution, if C=12−α.
(ii) X follows Pareto distribution where the quantile density function is q(u)=ba{(1−u)−(1+1a)},0≤u≤1;a,b>0 if C<12−α.
(iii) X follows a finite range distribution with quantile density function q(u)=ba{(1−u)(1a−1)},0≤u≤1;b>0,a>1 if C>12−α.
Proof: Assume that the conditions in Eq. (15) are held. Then, using Eq. (12), we have
Therefore, using KX1;n(u)=n(1−u)q(u) in Eq. (16), then, differentiating (w.r.t) u, we obtain
q′(u)q(u)=(2nc−ncα+αc−c−nc−cα)(11−u).
This implies q(u)=A(1−u)n(2c−cα−1c−αc−1), where A is a constant. Thus, X has exponential distribution, Pareto distribution and finite range distribution if C=12−α,C<12−α,C>12−α, respectively.
Theorem 6.2: For the exponential distribution, the difference between quantile-based M-H residual entropy of the life-time of the series system (Mα(X1;n;u)) and M-H residual quantile entropy of the life-time of each component (Mα(X;u)) is independent of u and relies only on α and the number of components of the system.
Proof: For the exponential distribution, we have
Mα(X1;n;u)=1α−1((nλ)1−αα−1) and Mα(X;u)=1α−1((λ)1−αα−1).
Therefore,
Mα(X1;n;u)−Mα(X;u)=(λ1−αα(α−1))(n1−α−1),
which complete the prove of Theorem 6.2.
Theorem 6.3: Let Xn;n represents the largest order statistics with hazard and survival quantile functions, K¯Xn;n(u) and F¯n;n(x). Then, for sample maxima Xn;n, the past quantile entropy of M-H, M¯α(Xn;n;u), is
M¯α(Xn;n;u)=1α−1{C(K¯Xn;n(u))1−α−1},0<α<2,
if and only if X follows the power distribution.
Proof: The quantile and quantile density functions for the power distribution are, respectively,
Q(u)=au1b and q(u)=abu1b−1,0≤u≤1;a,b>0.
It is simple to show that K¯Xn;n(u)=nbu1ba. Taking C=bn3−α2nb−nbα+α−1 gives the if part of Theorem.
Conversely, let Eq. (17) is valid. Thus, using Eq. (14), we determine
Substituting K¯Xn;n(u)=nuq(u), we have {n3−α∫0up(2−α)(n−1)qα−1(p)dp}=Cu2n−nα+α−1(q(u))α−1.
Taken the derivative (w.r.t) u yields q′(u)q(u)=(n3−α−C(2n−nα+α−1)Cα−C)(1u). The latter gives
q(u)=Au(n3−α−C(2n−nα+α−1)C(α−1)),
where A denotes a constant. Hence, the power distribution is a characterized for C=bn3−α2nb−nbα+α−1.
Simulation Study and Application to Real Life Data
In this paper, the quantile-based M-H entropy is proposed for some distributions. However, based on the available real data and to keep the simulation study related to the application part, we investigate the performance of the quantile-based M-H entropy for the exponential distribution.
Simulation Study
We conducted simulation studies to investigate the efficiency of the quantile-based M-H entropy estimators of smallest order statistics for exponential distribution (MX1:nα) in terms of the average bias (Bias), variance and mean squared error (MSE), based on sample sizes 10, 25, 100, 200 and 500 for different parameter combinations. The estimation of parameter λ was achieved using ML estimation and the process was repeated 2000 times.
From the results of the simulation study (see Tabs. 4 and 5), conclusions are drawn regarding the behaviour of the entropy estimator in general, which are summarized below:
The ML estimates of MX1:nα approaches to true value when sample size n increases.
When sample size n is increased, the MSE and variance of MX1:nα decreases.
Average estimates, Bias, Variance and Mean Squared Error for <inline-formula id="ieqn-210"><mml:math id="mml-ieqn-210"><mml:msubsup><mml:mi mathvariant="bold-italic">M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">X</mml:mi><mml:mrow><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mrow><mml:mo>:</mml:mo></mml:mrow><mml:mi mathvariant="bold-italic">n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>∝</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula>under exponential distribution for different values of <inline-formula id="ieqn-211"><mml:math id="mml-ieqn-211"><mml:mi mathvariant="bold-italic">λ</mml:mi></mml:math></inline-formula> and fixed value of <inline-formula id="ieqn-212"><mml:math id="mml-ieqn-212"><mml:mi mathvariant="bold-italic">α</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">0.2</mml:mn></mml:mrow></mml:math></inline-formula>
n
Criterion
λ=0.5,α=0.2
λ=0.8,α=0.2
λ=1.5,α=0.2
10
E(MX1:nα)
−4.071545
0.529750
0.513037
Bias
−0.460029
−0.047345
−0.064059
Variance
9.852729
0.175090
0.183407
MSE
10.064355
0.177332
0.187510
25
E(MX1:nα)
−3.477109
0.598628
0.598179
Bias
−0.189569
−0.023310
−0.023759
Variance
2.910436
0.050590
0.053469
MSE
2.946372
0.051133
0.054034
100
E(MX1:nα)
−3.202870
0.632070
0.635914
Bias
−0.055952
−0.009332
−0.005489
Variance
0.659129
0.012280
0.012697
MSE
0.662260
0.012367
0.012727
200
E(MX1:nα)
−3.135189
0.645690
0.639338
Bias
−0.010514
0.001209
−0.005143
Variance
0.301033
0.006124
0.006095
MSE
0.301144
0.006125
0.006122
500
E(MX1:nα)
−3.101819
0.647859
0.647135
Bias
0.009665
0.001552
0.000827
Variance
0.124756
0.002350
0.002462
MSE
0.124849
0.002353
0.002463
Average estimates, Bias, Variance and Mean Squared Error for <inline-formula id="ieqn-221"><mml:math id="mml-ieqn-221"><mml:msubsup><mml:mi mathvariant="bold-italic">M</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">X</mml:mi><mml:mrow><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mrow><mml:mo>:</mml:mo></mml:mrow><mml:mi mathvariant="bold-italic">n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mo>∝</mml:mo></mml:mrow></mml:msubsup></mml:math></inline-formula> under exponential distribution for different values of <inline-formula id="ieqn-222"><mml:math id="mml-ieqn-222"><mml:mi mathvariant="bold-italic">α</mml:mi></mml:math></inline-formula> and fixed value of <inline-formula id="ieqn-223"><mml:math id="mml-ieqn-223"><mml:mi mathvariant="bold-italic">λ</mml:mi><mml:mo>=</mml:mo><mml:mrow><mml:mn mathvariant="bold">1.3</mml:mn></mml:mrow></mml:math></inline-formula>
n
Criterion
λ=1.3,α=0.1
λ=1.3,α=0.9
λ=1.3,α=1.5
10
E(MX1:nα)
0.322104
1.655315
−0.099460
Bias
−0.034907
−0.010507
−0.028988
Variance
0.231985
7.787576
0.091822
MSE
0.233204
7.787687
0.092662
25
E(MX1:nα)
0.387186
2.185596
−0.189670
Bias
−0.020296
−0.007890
−0.013950
Variance
0.082536
3.056875
0.033740
MSE
0.082948
3.056938
0.033935
100
E(MX1:nα)
0.425782
2.430850
−0.232748
Bias
−0.003539
−0.000279
−0.004406
Variance
0.017327
0.737979
0.007372
MSE
0.017340
0.737979
0.007392
200
E(MX1:nα)
0.434903
2.468054
−0.238231
Bias
0.002132
−0.001162
−0.001118
Variance
0.008627
0.339305
0.003970
MSE
0.008631
0.339307
0.003971
500
E(MX1:nα)
0.431925
2.476880
−0.242884
Bias
−0.002892
−0.014988
−0.000509
Variance
0.003127
0.134222
0.001495
MSE
0.003136
0.134447
0.001496
Application to Real Life Data
The real data in this section represents the failure times of 20 mechanical components that were used previously by Murthy et al. [20] for investigating some of Weibull models. The data values are: 0.067, 0.068, 0.076, 0.081, 0.084, 0.085, 0.085, 0.086, 0.089, 0.098, 0.098, 0.114, 0.114, 0.115, 0.121, 0.125, 0.131, 0.149, 0.160, 0.485. We use this data for two main purposes: (i) for investigating the performance of our quantile-base M-H entropy (MX1:nα) in the exponential distribution case, and (ii) for comparing MX1:nαto the quantile-based Tsallis entropy (HX1:nα) that was proposed by Kumar [21].
Based on this data, we used first the maximum likelihood method to estimate the exponential distribution parameter, λ^=0.122. Then, for different values of α, varied from 0.1 to 0.9, we calculate the estimated values of MX1:nα and HX1:nαof smallest order statistics under exponential distribution. The results are displayed in Tab. 6.
Estimates of <inline-formula id="ieqn-239"><mml:math id="mml-ieqn-239"><mml:msubsup><mml:mrow><mml:mover><mml:mi mathvariant="bold-italic">M</mml:mi><mml:mo stretchy="false">^</mml:mo></mml:mover></mml:mrow><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">X</mml:mi><mml:mrow><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mrow><mml:mo>:</mml:mo></mml:mrow><mml:mi mathvariant="bold-italic">n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">α</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula> and <inline-formula id="ieqn-240"><mml:math id="mml-ieqn-240"><mml:msubsup><mml:mi mathvariant="bold-italic">H</mml:mi><mml:mrow><mml:msub><mml:mi mathvariant="bold-italic">X</mml:mi><mml:mrow><mml:mrow><mml:mn mathvariant="bold">1</mml:mn></mml:mrow><mml:mrow><mml:mo>:</mml:mo></mml:mrow><mml:mi mathvariant="bold-italic">n</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mrow><mml:mi mathvariant="bold-italic">α</mml:mi></mml:mrow></mml:msubsup></mml:math></inline-formula>for different values of <inline-formula id="ieqn-241"><mml:math id="mml-ieqn-241"><mml:mi mathvariant="bold-italic">α</mml:mi></mml:math></inline-formula>
α
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
M^X1:nα
0.948
1.024
1.209
1.205
1.314
1.439
1.582
1.747
1.937
HX1:nα
9.625
6.846
5.546
4.678
4.023
3.502
3.074
2.716
2.415
It should be noted that the estimated values of M^X1:nαis generally increased when 0<α<1. Also, the results in Tab. 6 indicate clearly that the estimated entropy values based on M^X1:nαare less than those given by HX1:nα.
Conclusion
The key focus of this article is to propose new quantile-based Mathai-Haubold entropy and investigate its characteristics. We also considered the divergence measure of the Mathai-Haubold and established some of its properties. Further, based on order statistics, we propose the residual entropy of the quantile-based Mathai-Haubold and some of its property results are proved. The performance of the proposed quantile-based Mathai-Haubold entropy is investigated by simulation studies and also by using real data application example. The proposed quantile-based Mathai-Haubold entropy's performance is investigated by simulation studies and by using a real data application example. We found that the ML estimates of MX1:nα approach true value when sample size n increases for the simulation part. For the application part, we compared our proposed quantile-based entropy to the existing quantile entropies and the results showed the outperformance of our proposed entropy to the other entropies. Our proposed quantile-based Mathai-Haubold entropy is useful for many future engineering applications such as reliability and mechanical components analysis.
Quantile functions are efficient and equivalent alternatives to distribution functions in modeling and analysis of statistical data. The scope of these functions and the probability distributions are essential in studying and analyzing real lifetime data. One reason is that they convey the same information about the distribution of the underlying random variable X. However, even if sufficient literature is available on probability distributions' characterizations employing different statistical measures, little works have been observed for modeling lifetime data using quantile versions of order statistics. Therefore, future work is necessary for enriching this area, and for this reason, we give precise recommendations for future research. First, the results obtained in this article are general because they can be reduced to some of the results for quantile based Shannon entropy for order statistics once parameter approaches unity. Recently, a quantile version of generalized entropy measure for order statistics for residual and past lifetimes was proposed by Kumar et al. [22]. Nisa et al. [23] presented a quantile version of two parametric generalized entropy of order statistics residual and past lifetimes and derived some characterization results. Moreover, Qiu [24] studied further results on quantile entropy in the past lifetime and gave the quantile entropy bounds in the past lifetime for some ageing classes. The ideas presented by these mentioned papers can be somehow combined/merged with our results in this paper to produce more results and properties for the quantile-based Mathai-Haubold. Second, Krishnan [25] recently introduced a quantile-based cumulative residual Tsallis entropy (CRTE) and extended the quantile-based CRTE in the context of order statistics. Based on these new results and our proposed quantile M-H entropy in this paper, one can follow Krishnan [25] and derive a quantile-based cumulative residual M-H entropy and extend it in the context of order statistics. Finally, Krishnan [25] also proposed a cumulative Tsallis entropy in a past lifetime based on quantile function. As an extension, the cumulative M-H entropy in a past lifetime based on quantile function can also be derived.
The authors would like to thank the anonymous reviewers and the editor for their useful suggestions and comments which increase and improve the quality of this paper.
Funding Statement: Authors thank and appreciate funding this work by the Deanship of Scientific Research at King Khalid University through the Research Groups Program under the Grant No. (R.G.P. 2/82/42).
Conflicts of Interest: The authors declare that they have no conflicts of interest to report regarding the present study.
ReferencesArnold, B. C., Blakrishan, N., Nagraja, N. H. (1992). Malik, J. H., Balakrishnan, N., Ahmad, S. E. (1998). Recurrence relations and identities for moments of order statistics I: Arbitrary continuous distributions. Samuel, P., Thomas, P. Y. (2000). An improved form of a recurrence relation on the product moment of order statistics. David, H. A. (1981). Shannon, C. E. (1948). A mathematical theory of communication. Mathai, A., Haubold, H. (2007). On generalized entropy measures and pathways. Sebastian, N. (2015). Generalized pathway entropy and its applications in diffiusion entropy analysis and fractional calculus. Mathai, A. (2005). A pathway to matrix-variate gamma and normal densities. Nair, N. U., Sankaran, P. G. (2009). Quantile-based reliability analysis. Sunoj, S. M., Krishnan, A. S., Sankaran, P. G. (2017). Quantile-based entropy of order statistics. Sunoj, S. M., Sankaran, P. G. (2012). Quantile based entropy function. Wong, K. M., Chen, S. (1990). The entropy of ordered sequences and order statistics. Park, S. (2007). The entropy of consecutive order statistics. Ebrahimi, N. (1996). How to measure uncertainty in the residual lifetime distribution. Baratpour, S., Ahmadi, J., Arghami, N. R. (2007). Some characterizations based on entropy of order statistics and record values. Paul, J., Yageen, P. (2019). On some properties of Mathai-Haubold entropy of record values. Kullback, S., Leibler, R. A. (1951). On information and sufficiency. Dar, J. G., Al-Zahrani, B. (2013). On some characterization results of life time distributions using Mathai-Haubold residual entropy. Crescenzo, A. D., Longobardi, M. (2002). Entropy based measure of uncertainty in past lifetime distribution. Murthy, D. N. P., Xie, M., Jiang, R. Y. (2004). Kumar, V. R. (2018). A quantile approach of Tsallis entropy for order statistics. Kumar, V., Singh, N. (2019). Quantile-based generalized entropy of order (α, β) for order statistics. Nisa, R., Baig, M. A. K. (2019). Characterization results based on quantile version of two parametric generalized entropy of order statistics. Qiu, G. (2019). Further results on quantile entropy in the past lifetime. Krishnan, A. S. (2019).