|Computer Modeling in Engineering & Sciences|
Quantile Version of Mathai-Haubold Entropy of Order Statistics
1Department of Mathematics, College of Science, King Khalid University, Abha, 62529, Saudi Arabia
2Statistical Research and Studies Support Unit, King Khalid University, Abha, 62529, Saudi Arabia
3Department of Mathematical Sciences, IUST, Kashmir, 192231, India
4Department of Mathematics, Faculty of Science, Al al-Bayt University, Mafraq, 25113, Jordan
5Department of Mathematics, Govt. Degree College Kilam, Higher Education, J & K, 192231, India
*Corresponding Author: Ibrahim M. Almanjahie. Email: firstname.lastname@example.org
Received: 06 November 2020; Accepted: 18 March 2021
Abstract: 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.
Keywords: Shannon entropy; Mathai-Haubold entropy; quantile function; residual entropy; order statistics; failure time; reliability measures
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. , Malik et al. ). 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.  and Arnold et al. .
Based on a random sample of , let the corresponding order statistics to be . Then, as in David  and Arnold et al. , the density of is
with . Eq. (1) can be used to determine the smallest (when ) and largest (when ) probability density functions; they, respectively, are and . The corresponding distribution functions are obtained, respectively, by and .
Shannon  was the first author who introduced the entropy idea for a random variable (r.v.) in the field of information theory and defined it as
Here, indicates the pdf of the r.v. . Based on the Shannon entropy measure, Mathi et al.  (now onwards M-H entropy) developed a generalized version of Eq. (2) and defined it as
When the M-H entropy measure will be reduced to the Shannon entropy measure defined in Eq. (2).
Mathai et al.  and Sebastian  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 . 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
We refer the readers to Nair et al.  and Sunoj et al.  and references therein for more details about quantile function. Recently, Sunoj et al.  studied Shannon entropy and as well as its residual and introduced quantile versions of them defined as
respectively, where denotes the density of quantile function. If we define the quantile density function by , then, we obtain
For Shannon past entropy, Sunoj et al.  also introduced its quantile version and defined it as
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 order statistics is given in Section 4. Quantile Residual Entropy of M-H for 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.
2 Quantile Based M-H Entropy of Order Statistics
Wong et al. , Park , Ebrahimi et al.  and Baratpour et al.  are the authors who discuss in detail the aspects of information-theoretic based on order statistics. Paul et al.  considered the M-H entropy and, based on record values, studied some of its essential properties. For order statistics , the M-H entropy is defined as
where is given in Eq. (1). Now, then, the pdf of order statistics becomes
where denotes beta-distribution density with and as its parameters. The quantile-based M-H entropy of is determined by
Remark 2.1: For , Eq. (6) reduces to
which is the quantile entropy of order statistics investigated by Sunoj et al. .
3 Expressions for Some Distributions
In the following, we provide expressions for Quantile-based M-H entropy of order statistics for some life time distributions:
(i) Govindarajulu’s Distribution: The quantile version and the corresponding density functions, respectively, are
Using Eq. (6), we can easily obtain quantile-based M-H Entropy of order statistics for Govindarajulu distribution as
Similarly, based on the quantile and quantile density functions, we obtain the quantile-based for the M-H Entropy ( of order statistics for the following distributions.
(ii) Uniform Distribution:
(iii) Pareto-I Distribution:
(iv) Exponential distribution:
(v) Power distribution
Figs. 1–3 give the quantile version of M-H entropy plots of smallest order statistics under exponential, Pareto-I and uniform distributions, respectively.
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.
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.
4 Quantile-Based Generalized Divergence Measure of 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 and are the density functions of the non-negative r.vs and respectively. The direct divergence of from is measured by the Kullback et al.  and is
With order , the divergence measure of the M-H or the relative-entropy of concerning is obtained by
where, for , the expression in Eq. (8) reduces to Eq. (7) (see Kullback et al. ).
Theorem 4.1: The quantile-based-generalized divergence measure between the order statistics distribution and the primary distribution is a distribution-free.
Proof: From equation Eq. (8), we have
Now, using the value of in Eq. (9), we obtain .
Using the fact that we determine the quantile-based-generalized divergence measure between the distribution of order statistics and primary distribution, as
which is a distribution-free. Hence, the theorem is proved.
5 M-H Quantile Residual Entropy for 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  describes a more practical approach that considers the use of age, defined as
with indicates the survival-function. Note that for the Eq. (5) reduces to Eq. (4). Denote by , a non-negative r.v., the unit’s life-time at time . Then, the residual function of M-H  is
The M-H residual entropy for the order statistics is given by
where refers to the survival-function of the order statistics and is the incomplete gamma function.
Considering the order statistic, the quantile residual entropy function of M-H of is given by
The following theorem will state important result.
Theorem 5.1: Considering the order statistic, the quantile residual entropy function of M-H determines the underlying distribution uniquely.
Proof: Using equation Eq. (10), we obtain
Differentiate both sides with respect to (w.r.t) to obtain
where denotes the differentiation w.r.t This equation involved an immediate connection between the and , which implies that the quantile residual entropy function of M-H of 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 order statistic for some lifetime models.
(i) Govindarajulu’s Distribution
The quantile for the Govindarajulu distribution is
and the corresponding density is
The quantile residual entropy function of M-H of 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 order statistics for the following distributions.
(ii) Uniform Distribution and
(iii) Pareto-I Distribution
(iv) Exponential distribution and
(v) Power distribution and
Based on residual M-H quantile entropy of order statistics , the following nonparametric classes of life distribution are defined.
Definition 5.1: is said to have an increasing (a decreasing) M-H quantile entropy of order statistics if is increasing (decreasing) in
The following lemma is useful in proving the results in monotonicity of .
Lemma 5.1: Let and be any two functions. If is increasing and is increasing (decreasing) in then is increasing (decreasing) in provided the existence of integrals.
Theorem 5.2: Let be a non-negative and continuous r.v. with quantile and density . Define , where is nonnegative, increasing and convex(concave) function. Then,
(i) For , increases (or decreases) in whenever increases (or decreases) in .
(ii) For , increases (or decreases) in whenever increases (or decreases) in .
Proof: (i) The quantile density of is given by .
Thus, we have
From the given condition, is increasing in , therefore,
is increasing in .
Since and is non-negative, increasing and convex (concave) function, the 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, , increases (or decreases) in , because increases and it is convex. Consequently, Eq. (11) is decreasing (increasing) in , which proves (ii) of the Theorem. The immediate application of Theorem 5.2 is given below:
Let be an r.v. following the distribution of exponential and having a failure rate . Also, let Therefore, follows Weibull distribution where . The function is a convex (concave) if Hence, based on Theorem 5.2, the Weibull distribution is increasing (decreasing) M-H quantile entropy of order statistics if
6 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 (for smallest) and (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 and .
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.  to build up the idea of past entropy on . If denotes the life-time of a component, thus, the past entropy of is obtained by
with is the cumulative distribution function. For Eq. (13) reduces to Eq. (2).
The quantile form of past M-H residual entropy of order statistic is determined by
For sample maxima , the past quantile entropy of M-H is given by
Next, we state some properties based on quantile M-H residual entropy of the smallest order statistics.
Theorem 6.1: Let represents the first order statistics with survival and hazard quantile functions and Then, the is determined by
if and only if:
(i) follows an exponential distribution, if
(ii) follows Pareto distribution where the quantile density function is , if
(iii) follows a finite range distribution with quantile density function , if
Proof: Assume that the conditions in Eq. (15) are held. Then, using Eq. (12), we have
Therefore, using in Eq. (16), then, differentiating (w.r.t) , we obtain
This implies , where is a constant. Thus, has exponential distribution, Pareto distribution and finite range distribution if , 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 and M-H residual quantile entropy of the life-time of each component is independent of and relies only on and the number of components of the system.
Proof: For the exponential distribution, we have
which complete the prove of Theorem 6.2.
Theorem 6.3: Let represents the largest order statistics with hazard and survival quantile functions, and . Then, for sample maxima , the past quantile entropy of M-H, , is
if and only if follows the power distribution.
Proof: The quantile and quantile density functions for the power distribution are, respectively,
It is simple to show that . Taking gives the if part of Theorem.
Conversely, let Eq. (17) is valid. Thus, using Eq. (14), we determine
Substituting , we have
Taken the derivative (w.r.t) yields The latter gives
where denotes a constant. Hence, the power distribution is a characterized for
7 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.
7.1 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 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:
1. The ML estimates of approaches to true value when sample size n increases.
2. When sample size n is increased, the MSE and variance of decreases.
7.2 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.  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 () in the exponential distribution case, and (ii) for comparing to the quantile-based Tsallis entropy that was proposed by Kumar .
Based on this data, we used first the maximum likelihood method to estimate the exponential distribution parameter, . Then, for different values of , varied from 0.1 to 0.9, we calculate the estimated values of and of smallest order statistics under exponential distribution. The results are displayed in Tab. 6.
It should be noted that the estimated values of is generally increased when . Also, the results in Tab. 6 indicate clearly that the estimated entropy values based on are less than those given by .
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 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. . Nisa et al.  presented a quantile version of two parametric generalized entropy of order statistics residual and past lifetimes and derived some characterization results. Moreover, Qiu  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  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  and derive a quantile-based cumulative residual M-H entropy and extend it in the context of order statistics. Finally, Krishnan  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.
Acknowledgement: 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.
|This work is licensed under a Creative Commons Attribution 4.0 International License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.|