Dynamic cumulative residual Rényi entropy for Lomax distribution: Bayesian and non-Bayesian methods

1.   Introduction

Entropy is a fundamental uncertainty measure of a random variable. Shannon ^[1] provided a quantitative measure of uncertainty as a measure of information. A flexible extension of Shannon entropy was introduced by Rényi ^[2]. The Rényi entropy is one parameter generalization of the Shannon entropy and can be used as a measure of randomness. Applications of Rényi entropy can be found in many fields, such as biology, genetics, electrical engineering, computer science, economics, chemistry and physics.

The Rényi entropy has been widely used in analysis of the quantum systems. It can be employed in the analysis of quantum communication protocols (Renner et al. ^[3]), quantum correlations (Lévay et al. ^[4]). Also it can be used to characterize the number of reads needed to reconstruct a DNA sequences (see, Motahari et al. ^[5]). Rényi entropy has numerous issues as a data analysis tools in many practical applications, especially in data that dealing with time-frequency representations. Gabarda and Cristobal ^[6] applied the Rényi entropy to an image fusion method and identified which pixels have a higher amount of information among the given input images. Rényi entropy can be used to calculate the risk measure of portfolio.

Entropy concepts can be related to informational market efficiency. Martina et al. ^[7] studied the complexity of crude oil prices using entropy measures to monitor the evolution of crude oil price movements. They found that the methods based on entropy concepts can shed light on the structure of crude oil markets as well as on its link to macroeconomic conditions and socio-political extreme events. Resconi et al. ^[8] provided a new suggestive reading of quantum mechanics by starting from the superposition of different Boltzmann entropies. They suggested that a quantum entropy space can be used as the fundamental arena that describing the quantum effects. Another measure of uncertainty deals with residual lifetime function is dynamic cumulative residual entropy (DCRE) which can be attractive in many fields like reliability and survival analysis. Recently, there are considerable literatures assigned to the applications, generalizations and properties of Rényi's measure of entropy.

Recently, there are many literatures studied the estimation procedure for the entropy measures. Kayal and Kumar ^[9] considered Bayesian estimation of entropy for exponential distribution under LINEX loss function. Seo et al. ^[10] produced an entropy estimate using upper record values from the generalized half-logistic distribution. Based on generalized progressive hybrid censoring scheme, Bayesian estimates of entropy for Weibull distribution were discussed by Cho et al. ^[11]. Based on record values, Chacko and Asha ^[12] studied the estimation of the entropy for generalized exponential distribution. Patra et al. ^[13] estimated a function of scale parameter of an exponential population under general loss function. Maximum likelihood (ML) estimate of Shannon entropy for inverse Weibull distribution under multiple censored data was discussed by Hassan and Zaky ^[14]. Petropoulos et al. ^[15] provided improved estimators of the entropy for the scale mixture of exponential distributions. Estimation of entropy for inverse Lomax distribution under multiple censored has been studied by Bantan et al. ^[16]. Bayesian estimate of the entropy for Lomax distribution based on upper record values has been studied by Hassan and Zaky ^[17].

If X is an absolutely continuous random variable with probability density function (PDF) f (x), then the corresponding Rényi entropy of order $\beta$ is defined as:

Recently, measurements of uncertainty for probability distributions became more interested. Sunoj and Linu ^[18] defined the cumulative residual Rényi entropy (CRRE) for residual lifetime X_t = (X −t | X > t) based on survival function instead of using PDF as follows:

where, $\bar F(t) = \, P\left( {X > {\rm{t}}} \right) = 1 - F\left( t \right)$ is the survival (reliability) function. The main features of the CRRE are always non-negative as well as it has consistent definitions in both the continuous and discrete domains. Furthermore it can be easily computed from sample data and these computations asymptotically converge to the true values (see, Rao et al. ^[19]). Sunoj and Linu ^[18] defined the DCRRE as a dynamic form of uncertainty as follows:

Therefore, when t = 0, the DCRRE tends to CRRE. In the literature, few studies had been done concerning the inferential procedures of the entropy measures incorporating DCRRE for lifetime distributions. Kamari ^[20] produced some properties of the DCRRE based on order statistics. Kundu et al. ^[21] presented the cumulative residual and past inaccuracy measures which are extensions of the corresponding cumulative entropies for truncated random variables. Renjini et al. ^[22] discussed the Bayesian estimates of the DCRRE for the Pareto distribution using Type II right censored data. Renjini et al. ^[23] discussed the Bayesian estimates of the DCRRE for Pareto distribution of the first kind under upper record values. Renjini et al. ^[24] obtained Bayesian estimate of the DCRRE for Pareto distribution based on complete data. Bayesian estimate of dynamic cumulative residual Shannon entropy for Pareto II distribution has been studied by Ahmadini et al. ^[25].

Lomax ^[26] proposed a very good alternative to the common lifetime distributions such as exponential, Weibull, or gamma when the experimenter presumes that the population distribution may be heavy-tailed. Lomax distribution can be applied in a variety of fields such as a business failure data model, income and wealth inequality, computer science, risk analysis and economics, actuarial science and reliability. Corbellini et al. ^[27] used it to model firm size and queuing problems. As an important lifetime distribution, the Lomax distribution with shape parameter $\alpha$ and scale parameter $\lambda$ has the following PDF

The cumulative distribution function (CDF) corresponding to (4) is given by,

At $\lambda = 1, \alpha \ne 1$, Lomax distribution reduces to the beta prime distribution (also known as inverted beta distribution or beta distribution of the second kind). At $\alpha = 1, \, \lambda \ne 1,$ Lomax distribution reduces to a log-logistic distribution. At $\alpha = \lambda = 1$, the Lomax distribution reduces to F-distribution, F (2, 2). The reliability and hazard rate functions of Lomax distribution are given, respectively, by

and,

The hazard rate function of Lomax distribution takes different shapes according to values of shape parameter $\alpha .$ It is a constant at $\alpha = 1$, and it is continuously decreasing at $\alpha > 1,$ which represents early failures.

Studies about the Lomax distribution have been provided by many authors. For instance, Abd-Elfattah et al. ^[28] obtained the Bayesian and non-Bayesian estimates of the sample size for the Lomax distribution in case of Type-I censored samples. Record values of Lomax distribution were discussed by Ahsanullah ^[29]. Balakrishnan and Ahsanullah ^[30] introduced some recurrence relations between the moments of record values from Lomax distribution. Hassan and Al-Ghamdi ^[31] determined the optimum test plan for simple step stress accelerated life testing. Hassan et al. ^[32] discussed the optimal times of changing stress level for k-level step stress accelerated life tests based on adaptive Type-II progressive hybrid censoring with product's life time following Lomax distribution. For more application about Lomax distribution, the reader can refer to ^{[33,34,35,36,37,38]}.

From the above, it is evident that the Lomax distribution has been received greatest attention from theoretical and statisticians primarily due to its use in various fields. In the literature, there are no reports about the statistical inference of the DCRRE for Lomax distribution. So we need to estimate the DCRRE for Lomax distribution in view of Bayesian and non-Bayesian procedures based on complete samples. The ML estimates and approximate confidence intervals of the DCRRE are obtained. The Bayesian estimate is calculated using gamma priors under squared error (SE), LINEX and precautionary (PRE) loss functions. According to the complicated forms of the DCRRE Bayesian estimate, we employ the Markov Chain Monte Carlo (MCMC) simulation technique for numerical study.

This article can be organized as follows. Section 2 gives the DCRRE of Lomax distribution as well as obtains the ML and approximate confidence intervals of DCRRE. Section 3 presents Bayesian estimate of the DCRRE for Lomax distribution under symmetric and asymmetric loss functions. The MCMC technique is given in Section 4. A real data application is illustrated in Section 5. The paper ends with some conclusions based on the results of the numerical studies.

2.   Maximum likelihood estimate of the DCRRE

This section provides an explicit expression of the DCRRE for Lomax distribution. Then, the ML estimate of the DCRRE is derived based on complete random sample.

The DCRRE of Lomax distribution is obtained by substituting (5) in (3) as follows:

Hence, after simplification the DCRRE takes the form

This is the required expression of the DCRRE for Lomax distribution. To obtain the ML estimate of ${\gamma _R}(\beta ),$ we must obtain the ML estimate of the model parameters $\alpha$ and $\lambda .$

Consider a simple random sample of size n drawn from PDF (4) and CDF (5), where $\alpha$ and $\lambda$ are unknown. Then, given the sample $\underline x = ({x_1}, {x_2}, ..., {x_n})$, the likelihood function of Lomax distribution can be written as follows

It is usually easier to maximize the natural logarithm of the likelihood function rather than the likelihood function itself. Therefore, the logarithm of likelihood function, say $\ln \ell ,$ is

The partial derivatives of the log-likelihood function with respect to $\alpha$ and $\lambda$, are obtained as follows:

and

The ML estimates of $\alpha$ and $\lambda$ are determined by solving the equations ${{\partial \ln \ell } \mathord{\left/ {\vphantom {{\partial \ln \ell } {\partial \alpha }}} \right. } {\partial \alpha }} = 0,$ and ${{\partial \ln \ell } \mathord{\left/ {\vphantom {{\partial \ln \ell } {\partial \lambda }}} \right. } {\partial \lambda }} = 0,$ simultaneously. Further the resulting equations cannot be solved analytically, so, numerical technique must be applied to solve these equations, simultaneously, to obtain ${\hat \alpha _{ML}}\,$ and ${\hat \lambda _{ML}}.$ The existence and uniqueness of ${\hat \alpha _{ML}}\,$ and ${\hat \lambda _{ML}}$ can be checked by solving, algebraically, ${{\partial \ln \ell } \mathord{\left/ {\vphantom {{\partial \ln \ell } {\partial \alpha }}} \right. } {\partial \alpha }} = 0$ and ${{\partial \ln \ell } \mathord{\left/ {\vphantom {{\partial \ln \ell } {\partial \lambda }}} \right. } {\partial \lambda }} = 0,$ simultaneously and checking that, the second derivative is negative which is a maximum point.

Hence, by using the invariance property of ML estimate, then ML estimates of the DCRRE denoted by ${\hat \gamma _{R, ML}}(\beta )$ becomes

Furthermore, a confidence interval of the DCRRE is the probability that a real value of the entropy will fall between an upper and lower bounds of a probability distribution. For large sample size, the ML estimates, under appropriate regularity conditions, are consistent and asymptotically normally distributed. Therefore, the two-sided approximate confidence limits for ${\hat \gamma _{R, ML}}(\beta )$ can be constructed, such that

where, ${z_{{\eta \mathord{\left/ {\vphantom {\eta 2}} \right. } 2}}}$ is 100 $\left( {1 - \eta } \right)$ the standard normal percentile, and $\eta$ is the significant level. Therefore, the lower (L) and upper (U) confidence limits for ${\hat \gamma _{R, \, ML}}(\beta )$ can be obtained as

where ${\sigma _{({{\hat \gamma }_{R, \, ML}}(\beta ))}}$ is the standard deviation. The two-sided approximate confidence limits for DCRRE will be constructed with confidence level 95 %.

3.   Bayesian estimate of the DCRRE

In this section, the Bayesian estimates of ${\gamma _R}(\beta )$ under symmetric (SE) and asymmetric (LINEX and PRE) loss functions are obtained by assuming the priors of parameters $\alpha$ and $\lambda$ have gamma distributions. As reported by Kayal and Kumar ^[9], the SE loss function is used when weights of error caused by observation and underestimation are equal. But asymmetric loss function is used when the overestimation may have more serious than underestimation or vice versa. To compute the Bayesian estimate of ${\gamma _R}(\beta )$, we proceed as follows. Firstly, the Bayesian estimates of $\alpha$ and $\lambda$ are obtained and then the DCRRE estimate will be calculated by substituting the parameter estimates in (7). Additionally, the width of the Bayesian credible interval (BCI) estimates is obtained.

In Bayesian method, there are many methods that can be used in selection of prior. There are two types of the prior distributions, non-informative and informative prior (IP) distributions. The prior of parameters may be selected as independent or dependent priors. In this paper, we assume the prior of parameters $\alpha$ and $\lambda ,$ denoted by ${\pi _1}(\alpha )$ and ${\pi _2}(\lambda )$ are independent. Following Pak and Mahmoudi ^[39], we take ${\pi _1}(\alpha )$ and ${\pi _2}(\lambda )$ to be gamma distributions with parameters (a, b) and (c, d) as follows

where a, b, c and d are known and non-negative. So, the joint posterior for parameters, denoted by $\pi _{1, 2}^*(\alpha , \lambda ),$ is

So, the marginal posterior PDF of parameters $\alpha$ and $\lambda$ is given respectively by:

and

Therefore, the Bayesian estimates of $\alpha$ and $\lambda$ under SE loss function, say ${\hat \alpha _{SE}}$, and ${\hat \lambda _{SE}}$ can be obtained as posterior mean as follows:

and,

While, the Bayesian estimates of parameters $\alpha$ and $\lambda$ under LINEX loss function, say ${\hat \alpha _{LINEX}}$, and $\, \, {\hat \lambda _{LINEX}}$ are obtained as follows:

and

where, $\nu$ is a real number. Additionally, the Bayesian estimates of parameters $\alpha$ and $\lambda$ under PRE loss function say ${\hat \alpha _{PRE}}\,$ and $\, {\hat \lambda _{PRE}}$ are given as follows:

and

Integrals (18)–(23) do not take a closed form, so the Metropolis-Hastings (M-H) and random-walk Metropolis algorithms are employed to generate the MCMC samples from posterior density functions (16) and (17), respectively. Hence, the Bayesian estimates of $\alpha$ and $\lambda$ under SE, LINEX and PRE loss functions are obtained as the mean of the simulated samples from their posteriors. Further, once the Bayes estimates of $\alpha$ and $\lambda$ are obtained. Therefore, the Bayesian estimate of the DCRRE, denoted by ${\hat \gamma _{R, \, BE}}{(\beta )_{SE}},$ under SE loss function is obtained, by using (7), as follows:

By similar way, we obtain the Bayesian estimate of ${\gamma _R}(\beta )$ for Lomax distribution under LINEX, and PRE loss functions. Furthermore, the BCI is a useful summary of the posterior distribution which reflects its variation that is used to quantify the statistical uncertainty. An approximate highest posterior density interval of ${\gamma _R}(\beta )$ is obtained by using the same algorithm of Chen and Shao ^[40].

4.   Monte Carlo simulation

This section evaluates the performance of the ML and Bayes estimates (BEs) based on simulation study.

The maximum likelihood estimates (MLEs) of the DCRRE for various sample sizes are computed by using Equations (9) and (10). Criteria measures include mean square error (MSE), average length (AL) and coverage probability (CP) are computed to investigate the performance of MLEs. Simulation procedures can be described as follows:

●  1000 random samples of sizes n = 30, 50, 70 and 100 are generated from Lomax distribution.

●  The values of parameters are selected as $\alpha = (1.5, 2.5),$ $\lambda = (0.5, 1.5, 4)$ and $\beta = (3, 5).$

●  The true values of the DCRRE measure are selected as ${\gamma _R}(\beta ) =$ 0.2798, 0.6264 and 0.9359 at t = 0.5, while ${\gamma _R}(\beta ) =$ 0.077075, 0.2798 and 0.58932 at t = 1.5.

●  The MLEs of $\mathrm{ }$ $\alpha$ and $\lambda$ are obtained from (9) and (10), then the MLE of ${\gamma _R}(\beta )$ is obtained by substituting ${\hat \alpha _{ML}}\,$ and ${\hat \lambda _{ML}}$ in (11). The approximate confidence interval of ${\hat \gamma _{R, \, ML}}(\beta )$ using (12) are constructed with confidence level at $\eta = 0.05.$

●  Compute the average MLEs of the DCRRE, MSEs, CPs, and ALs.

The M-H algorithm is one of the most famous subclasses of the MCMC method in Bayesian literature used to simulate the deviates from the posterior density and produce the good approximate results. The relative absolute biases (RABs), estimated risks (ERs) and the width of the BCI are computed to assess the behavior of the Bayesian estimates. A simulation study is done via R 3.1.2. Additionally, tables of simulation results are given. To compare the DCRRE Bayes estimates, the MCMC simulations are designed for different sample sizes under SE, LINEX and PRE loss functions. The values of parameters are selected as $\alpha = (1.5, 2.5),$ $\lambda = (0.5, 1.5, 4)$ and $\beta = (3, 5).$ The true values of the DCRRE measure are selected as ${\gamma _R}(\beta ) =$ 0.2798, 0.6264 and 0.9359 at t = 0.5, while ${\gamma _R}(\beta ) =$ 0.077075, 0.2798 and 0.58932 at t = 1.5. The hyper-parameters for gamma distribution are selected as a = c = 1 and b = d = 4. Also, we take (v = −2, 2) for LINEX loss function. Take N = 5000 random samples of sizes n, where each sample size (n) = 10, 30, 50, 70 and 100 are generated from Lomax distribution.

The M-H algorithm is described as follows:

Step 1: Let g(.) be the PDF of subject distribution;

Step 2: Initialize a starting value ${x_0}$ and determine the number of samples N;

Step 3: for i = 2 to N set $x = {x_{i - 1}}$;

Step 4: Generate u from uniform (0, 1) and generate y from g(.);

Step 5: if $u \leqslant \frac{{{\pi _\alpha }(y)\, g(x)}}{{{\pi _\alpha }(x)\, g(y)}}$ then set ${x_i} = y$ else set ${x_i} = x$;

Step 6: Set $i = i + 1$ and return to step 2 and repeat the previous steps N times.

4.1.   Numerical results based on ML method

Numerical observations for DCRRE estimates using ML method at $\beta = 3$ and $\beta = 5$ are presented in Tables 1–4 and represented in Figures 1–4. Some observations related to the behavior of the DCRRE estimates are outlined as follows:

●  The value of MLE, ${\hat \gamma _{R, \, ML}}(\beta )\, ,$ decreases as the value of $\lambda$ increases for fixed value of $\alpha$ (see Tables 1 and 2)

●  The value of MLE, ${\hat \gamma _{R, \, ML}}(\beta )\, ,$ decreases as $\alpha$ decreases for fixed value of $\lambda$ (see Tables 1 and 2).

●  The value of MLE, ${\hat \gamma _{R, \, ML}}(\beta )\, ,$ at $\alpha = 1.5, \, \lambda = 4$ has the smallest value compared with the corresponding estimates of the other sets of parameters (see Tables 3 and 4).

●  The MSE of ${\hat \gamma _{R, \, ML}}(\beta )$ decreases as the sample size increases.

●  For fixed value of $\lambda ,$ the MSE of ${\hat \gamma _{R, \, ML}}(\beta )\,$ decreases as the value of $\alpha$ decreases.

●  For fixed value of $\alpha ,$ the MSE of ${\hat \gamma _{R, \, ML}}(\beta )\,$ decreases as the value of $\lambda$ increases.

●  The MSE of the DCRRE estimates gets the smallest value at ${\gamma _R}(\beta )\, = 0.2798,$ where $(\alpha = 1.5, \, \lambda = 1.5),$ t = 0.5 and $\beta = 3,$ compared to the MSE of MLE of DCRRE measure for the corresponding other sets of parameters (see Figure 1).

●  The MSE of the DCRRE estimates takes the smallest value at ${\gamma _R}(\beta )\, = 0.077075,$ where $(\alpha = 1.5, \, \lambda = 1.5),$ t = 1.5 and $\beta = 3,$ compared to MSE of MLE of DCRRE for the corresponding other sets of parameters (see Figure 2).

●  At $\beta = 5,$ $(\alpha = 1.5, \, \lambda = 1.5),$ where the true value of ${\gamma _R}(\beta )\, =$ $0.2946637$ the MSE of ${\hat \gamma _{R, \, ML}}(\beta )\,$ has the smallest value compared to the MSE of MLEs for the corresponding other sets of parameters at t = 0.5 (see Figure 3).

●  At t = 1.5, the MSE of ${\hat \gamma _{R, \, ML}}(\beta )\,$ at $(\alpha = 1.5, \, \lambda = 4.0),$ has the smallest value compared to the MSE of MLEs for the corresponding other sets of parameters (see Figure 4).

Generally as seen from Tables 1–4 that the CP is very close to their corresponding nominal levels in approximately most of the cases. Also, the CP increases as the sample size increases for all values of DCRRE.

●  The MLE of the DCRRE decreases as the value of $\lambda$ increases.

●  Regarding the AL of estimates, it can be observed that, as n increases the AL of DCRRE estimates decreases.

4.2.   Numerical results based on Bayesian method

Tables 5–10 and Figures 5–8, list the numerical outcomes for different estimates of the DCRRE under different loss functions at $\beta = 3.$ Tables 11–14 and Figures 9–12, list and describe the simulation outcomes for different estimates of the DCRRE under different loss functions at $\beta = 5$. Some observations related to the behavior of the DCRRE estimates are provided as follows:

●  As anticipated, the performance of all DCRRE estimates become improved as the sample size increases.

●  The estimated value of the DCRRE increases as the value of $\alpha$ increases for fixed $\lambda .$ The estimated value of the DCRRE decreases as the value of scale parameter $\lambda$ increases for fixed $\alpha .$

●  The estimated value of the DCRRE decreases as the value of t increases for fixed values of $\alpha$ and $\lambda .$

●  Under SE loss function, the ERs for DCRRE estimates get the smallest values at true value ${\gamma _R}(\beta ) =$ 0.6264, for most values of n (see for example Figures 5 and 6).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 gets the smallest values for all values of n. Also, the width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 is the shortest compared to the width of the BCI in case of the SE and PRE loss functions for all n (see Table 5).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{SE}}$ takes the smallest values for all n except n = 100. The width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 is the shortest compared to the width of the BCI in case of the PRE and SE loss functions for all n (see Table 6).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 takes the smallest values for all values of n. The width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 is the shortest compared to the width of the BCI in case of PRE and SE loss functions at n = 10 and 30 (see Table 7).

●   At t = 1.5, the ERs for DCRRE estimates get the smallest values at ${\gamma _R}(\beta ) =$ 0.2798, for all n under SE and PRE loss functions (see for example Figures 7 and 8).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = 2 takes the smallest values at n = 50, 70 and 100, while the ER of ${\hat \gamma _{R, BE}}{(\beta )_{PRE}}$ takes the smallest values at n = 10 and 30. The width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{SE}}$ is the shortest compared to the width of the BCI in case of LINEX and PRE loss functions for all values of n (see Table 8).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 takes the smallest values at n = 50, 70 and 100, while the ER of ${\hat \gamma _{R, BE}}{(\beta )_{PRE}}$ takes the smallest values at n = 10 and 30. The width of the BCI of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 is the shortest compared to the width of the BCI in case of the SE and PRE loss functions for all values of n (see Table 9).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 takes the smallest values for n = 10, 70 and 100. The width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 is the shortest compared to the width of the BCI in case of the SE and PRE loss functions for most values of n (see Table 10).

●  At $\beta = 5$, the ERs for DCRRE estimates get the smallest values at ${\gamma _R}(\beta ) =$ 0.0919312, for most values of n under most selected loss functions (see for example Figures 9 and 10).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 gets the smallest values for all values of n. While, the width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = 2 is the shortest compared to the width of the BCI in case of the SE and PRE loss functions for all n (see Table 11).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 gets the smallest values for all values of n. While, the width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{SE}}$ is the shortest compared to the width of the BCI in case of the LINEX and PRE loss functions for all n expect at n = 10 (see Table 12).

●  Under SE loss function, the ERs for DCRRE estimates get the smallest values at ${\gamma _R}(\beta ) =$ 0.04176352, for most values of n (see for example Figures 11 and 12).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = 2 gets the smallest values for n = 10, 70 and 100. While, the width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = 2 is the shortest compared to the width of the BCI in case of the SE and PRE loss functions for all n (see Table 13).

●  The ER of ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = 2 gets the smallest values for all values of n. While, the width of the BCI for ${\hat \gamma _{R, BE}}{(\beta )_{LINEX}}$ at v = −2 is the shortest compared to the width of the BCI in case of the SE and PRE loss functions for all n (see Table 14).

5.   Application to real data

In this section, we provide an application to real data set to prove the importance and flexibility of the Lomax distribution compared with some other models. Also, the real data set can be used to illustrate the proposed methods in Sections 2 and 3.

The real data set was used by Jorgensen ^[41] which represent the active repair times (in hours) for airborne communication transceiver. The data are recorded as follows:

0.50, 0.60, 0.60, 0.70, 0.70, 0.70, 0.80, 0.80, 1.00, 1.00, 1.00, 1.00, 1.10, 1.30, 1.50, 1.50, 1.50, 1.50, 2.00, 2.00, 2.20, 2.50, 2.70, 3.00, 3.00, 3.30, 4.00, 4.00, 4.50, 4.70, 5.00, 5.40, 5.40, 7.00, 7.50, 8.80, 9.00, 10.20, 22.00, 24.50.

Firstly, some preliminary data analysis is performed. The histogram, scaled total time on test (TTT) plots and the log likelihood for the real dataset are presented in Figure 13. We can observe that the shape of TTT plots is convex curve for dataset, which demonstrates decreasing failure rate.

Secondly, an application to a real data set is given to demonstrate the importance and flexibility of the proposed distribution compared with other one, two and three parameters fitted distributions. The suggested distributions are the alpha power Lindley (APLi) (Dey et al. ^[42]), Lindley (Li) and alpha power Weibull (APW) (Nassar et al. ^[43]) distributions. The MLE of the unknown parameters and the corresponding standard error (SE) of Lomax distribution and the other competitive models are calculated. In order to compare the different models with the Lomax model, we use measures like, Akaike information criterion (AIC), and Bayesian information criterion (BIC), consistent AIC (CAIC), Hannan-Quinn information criterion (HQIC). Also, the Kolmogorov- Smirnov (K-S) goodness of fit test and its P value (PV) for all models are calculated. Table 15 reports the MLEs and the corresponding SEs for all models. Table 16 reports numerical values of AIC, BIC, CAIC, HQIC, K-S and PVs for the different competitive distributions. From the results, we found that the Lomax distribution takes the smallest values of the considered measures and the largest value of PV among all other competitive models, and so it could be chosen as the best model.

Thirdly, regarding this data, the ML and Bayes estimates of the DCRRE under SE, LINEX and PRE loss functions are obtained and listed in Table 17.

As anticipated, we conclude from Table 17 that the DCRRE estimates are decreasing function on time (t), that is the estimated values of the DCRRE decrease as the time t increases.

6.   Conclusions

In this paper, the maximum likelihood and Bayesian methods of estimation for dynamic cumulative residual Rényi entropy of Lomax distribution are considered. The maximum likelihood estimates and approximate confidence intervals of the DCRRE are obtained. The performance of the DCRRE estimates for Lomax distribution is investigated via simulation in terms of their mean square error, average length and coverage probability. The Bayesian estimate of the DCRRE for Lomax model is obtained by considering the gamma priors under symmetric and asymmetric loss functions. The MCMC procedure is employed to compute the Bayes estimates and the Bayesian credible intervals. The performance of the DCRRE estimates for Lomax distribution is inspected through their relative absolute bias, estimated risk and the width of the credible intervals. Application to real data is provided.

Regarding the simulation results, we conclude that the mean square error of maximum likelihood and Bayesian estimates of the DCRRE decreases as the sample size increases. Also, the average length of the DCRRE estimates decreases and coverage probability increases as the sample size increases. The ER and the width of credible intervals of the DCRRE Bayes estimates decrease. For small true values of DCRRE, the width of DCRRE of BCIs under LINEX loss function is smaller than the corresponding based on SE and PRE loss functions for all selected sample size at t = 0.5 and 1.5. For large true values of DCRRE, the BCI width of the DCRRE under LINEX loss function is smaller than the corresponding based on SE and PRE loss functions for large sample size at t = 0.5 and 1.5.

Generally, the DCRRE estimates for both ML and Bayesian methods approach the true value as the sample size increases. As the time (t) increases, both ML and Bayesian estimates of the DCRRE decrease. Bayesian estimates under LINEX loss function at v = −2 are more suitable than other estimates under SE and PRC loss functions in most of the situations.

Acknowledgments

The authors thank the reviewers for their detailed and constructive comments. This project is supported by Researchers Supporting Project number (RSP-2020/156) King Saud University, Riyadh, Saudi Arabia. The first author, therefore, gratefully acknowledges the KSU for technical and financial support.

Conflict of interest

The authors declare no conflict of interest.