Next Article in Journal
Multi-Objective Optimal Design of a Hydrogen Supply Chain Powered with Agro-Industrial Wastes from the Sugarcane Industry: A Mexican Case Study
Next Article in Special Issue
Optimal Constant-Stress Accelerated Life Test Plans for One-Shot Devices with Components Having Exponential Lifetimes under Gamma Frailty Models
Previous Article in Journal
Almost Automorphic Solutions in the Sense of Besicovitch to Nonautonomous Semilinear Evolution Equations
Previous Article in Special Issue
A Study on Computational Algorithms in the Estimation of Parameters for a Class of Beta Regression Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Estimation Using Expected LINEX Loss Function: A Novel Approach with Applications

1
Department of Statistics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Statistics, Faculty of Commerce, Zagazig University, Zagazig 44519, Egypt
3
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
4
Department of Mathematics, Faculty of Science, Al-Azhar University, Cairo 11884, Egypt
5
School of Mathematics, Yunnan Normal University, Kunming 650092, China
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(3), 436; https://doi.org/10.3390/math10030436
Submission received: 18 January 2022 / Revised: 24 January 2022 / Accepted: 25 January 2022 / Published: 29 January 2022
(This article belongs to the Special Issue Statistical Simulation and Computation II)

Abstract

:
The loss function plays an important role in Bayesian analysis and decision theory. In this paper, a new Bayesian approach is introduced for parameter estimation under the asymmetric linear-exponential (LINEX) loss function. In order to provide a robust estimation and avoid making subjective choices, the proposed method assumes that the parameter of the LINEX loss function has a probability distribution. The Bayesian estimator is then obtained by taking the expectation of the common LINEX-based Bayesian estimator over the probability distribution. This alternative proposed method is applied to estimate the exponential parameter by considering three different distributions of the LINEX parameter, and the associated Bayes risks are also obtained in consequence. Extensive simulation studies are conducted in order to compare the performance of the proposed new estimators. In addition, three real data sets are analyzed to investigate the applicability of the proposed results. The results of the simulation and real data analysis show that the proposed estimation works satisfactorily and performs better than the conventional standard Bayesian approach in terms of minimum mean square error and Bayes risk.

1. Introduction

In Bayesian analysis and decision theory, the loss function plays an important role as it can be used to describe the overestimation and underestimation in analysis. Conventionally used loss functions include symmetric loss and asymmetric loss. In the symmetric loss function, the estimation treats overestimation and underestimation equally, whereas the asymmetric loss function gives different weights to overestimation and underestimation. Compared with the symmetric loss, the asymmetric loss is more realistic and useful in practical applications. For example, in reliability and survival analysis, overestimation of the reliability function or average failure time is usually much more serious than underestimation of the reliability function or mean failure time. Similarly, an underestimate of the failure rate results in more serious consequences than an overestimate of the failure rate. Moreover, many other authors, including Zellner [1], Chang and Huang [2] and Khatun and Matin [3], also pointed out that commonly used symmetric loss functions such as squared error (SE) loss may be inappropriate in practical applications. In consequence, asymmetric loss functions are widely applied in statistical inference; see, for example, Ali [4,5]. One of the most useful asymmetric loss functions is called LINEX loss, first introduced by Klebanov [6] and used by Varian [7] in his study of real estate assessment. The LINEX loss function rises approximately exponentially on one side of zero and approximately linearly on the other side.
Let θ be the unknown parameter to be estimated and θ ˜ be an estimate of θ ; in this case, the LINEX loss function can be written as
L ( θ ˜ , θ ) = e φ ( θ ˜ θ ) φ ( θ ˜ θ ) 1 , φ 0 .
From (1), the sign and magnitude of the shape parameter φ represent the direction and degree of symmetry of LINEX loss. Several comments are necessary regarding (1). First, for φ > 0 , overestimation is more serious than underestimation, and vice versa. When φ is close to zero, the LINEX loss is approximately the squared error loss and therefore almost symmetric. For clarity of illustration, various plots of the LINEX loss function L ( Δ ) with Δ = θ ˜ θ are presented in Figure 1 under different choices of parameter φ . The plots also demonstrate that positive, negative and small values of φ give more weight for overestimating, underestimating and almost equal estimation, respectively. It is observed from Figure 1a that, for positive values of φ , the function is completely asymmetric, with overestimation being more harmful than underestimation. For negative values of φ , as in Figure 1b, the function grows approximately exponentially when Δ < 0 and almost linearly when Δ > 0 . Moreover, when φ tends towards zero, as in Figure 1c, the function is approximately symmetric as the squared error loss function.
Let g ( θ | d a t a ) be the posterior distribution of the parameter θ ; then, the Bayesian estimator of θ using the LINEX loss function, denoted by θ ˜ B L , is given by
θ ˜ B L = 1 φ log E θ e φ θ , φ 0 ,
where E θ [ · ] is the expectation with respect to g ( θ | d a t a ) .
Due to the practicality of the asymmetric LINEX loss in many applications, including reliability, economic and social studies, among others, many authors have studied the Bayesian inference and decision analysis problems using the LINEX loss function. Dey et al. [8] explored the robustness properties of classes of LINEX losses under exponential and discrete power series families of distributions. Soliman [9] compared the Bayesian estimators of the parameters and reliability characteristics for Rayleigh distribution using quadratic and LINEX loss functions. Misra and Meulen [10] investigated the estimation of the mean of two normal populations with unknown means and common known variance using the LINEX loss function. Micheas [11] investigated robustness to loss functions employing the LINEX loss function and applied the methods to standard examples. Hoque et al. [12] analyzed the performance of the unrestricted estimator and preliminary test estimator of the linear regression model based on the LINEX loss function. Kazmi et al. [13] studied the properties of Bayesian estimators of the parameter of a class of lifetime distributions, employing various loss functions, including the LINEX loss, through simulation and real data. Theoretical results were obtained for various lifetime distributions under different data types by many authors; see, for example, the works of Pandey et al. [14], Jaheen [15], Azimi et al. [16], Ahmed [17,18], Ashour and Nassar [19], Nassar et al. [20], Kotb and Raqab [21], Chen and Gui [22] and Nassar et al. [23], among others. From practical perspectives, due to its flexibility of measuring the overestimation and underestimation and its more appropriate data fitting feature in data analysis, the LINEX loss is also widely used in many practical discussions. For instance, in order to estimate the mean of the selected normal population from two normal populations with unknown means and common known variance, Misra and Meulen [10] proposed some admissibility results for a subclass of equivalent estimators, and a sufficient condition for the inadmissibility of an arbitrary equivalent estimator was also provided. In the classical Taguchi quality model, Chang and Hung [2] used the LINEX loss to measure the loss of quality to determine optimum process parameters for the product quality. By inheriting the asymmetric merit of LINEX loss, Tang et al. [24] proposed a general multi-view LINEX SVM framework to improve the performance of multi-view learning. Estimation of the costs of non-conforming products is one of the most important inputs in the economic design of control charts. Abolmohammadi et al. [25] proposed an extending control chart by considering the LINEX loss on the basis of the Markov chain approach, and the results showed that the practical VSIX control chart based on the LINEX loss was more efficient in terms of the expected cost per hour. However, it should be mentioned here that, after extensive documentation retrieval, the statistical inferential results were obtained under LINEX loss with fixed parameter φ , where a random chosen positive or negative φ was taken for indicating overestimation and underestimation, respectively, but there was no discussion on how to choose a proper value for the LINEX parameter φ .
In the practical Bayesian inferential approach, with respect to LINEX loss and even other asymmetric losses, it is an inevitable challenge issue to choose the values of the LINEX parameter in practice. To balance the overestimation and underestimation for practical data analysis and to give relatively explicit weights between over- and underestimation, the value(s) of the LINEX parameter φ are usually determined through historical information, expert experience, as well as some other strategies. Obviously, inference using a fixed parametrization has limitations and the randomly given value of the parameter φ merely provides overestimation or underestimation in general and cannot provide a proper balance between overestimation and underestimation in practice. To overcome this drawback, it is more reasonable to take the idea that there are many values of the parameter φ in a given issue, and these potential values feature different weights. Therefore, it is assumed that the LINEX parameter φ is a random variable. In Bayesian inference, the main advantage of such an assumption is that, if the number of choices for the parameter φ is sufficient, the subjective prefixed value of parameter φ can be viewed as a good choice in analysis. However, within the extensive literature research and empirical illustrations, such an arbitrary subjective approach may lead to poor inferential results if there is not enough information available. The value of the parameter φ is given randomly in many practical studies, and the final “proper value” may be determined by empirical comparison and cross-validation approaches. Therefore, compared with the conventional subjective procedure, it is necessary and important to adopt an alternative method for the determination of the parameter φ , where relatively objective value(s) of φ will provide a more effective balance between overestimation and underestimation in Bayesian inference. It should be mentioned here that, from the perspective of physical application, for complex macroscopical or microcosmic systems, the system has multiple normal states in its operating cycle, and the observed condition is the superposition across different detailed state combinations. Therefore, the presented approach may have potential benefits in studies of statistical physics.
Motivated by such reasons, the main aim of this paper is to introduce an alternative Bayesian analysis for parameter estimation with respect to the LINEX loss function. To overcome the subjective method for the choice of the parameter φ , it is assumed that the parameter φ is denoted by a random variable with a proper probability distribution. By taking the expectation of the Bayesian estimator using the usual LINEX loss function over the distribution of φ , the proposed Bayesian estimator can be obtained in consequence. The new Bayesian estimator is called the Bayesian estimator using the expected LINEX loss function. When the lifetime data come from an exponential distribution, Bayesian estimation using the expected LINEX loss is discussed and empirical illustrations are carried out to show the associated effect in the estimator(s) behavior compared with some arbitrary considered value(s) of the parameter φ .
The concept of the loss (or cost) function is employed in different optimization procedures in statistical theory and the decision-making process to scheme an event in a real number and the associated loss or cost. An optimization problem always aims to minimize a loss function. As shown by the empirical results in Section 4 and Section 5, the expected LINEX loss function provides Bayesian estimators with minimum loss as measured through the Bayes risks. Therefore, we can conclude that the Bayesian estimators obtained based on the suggested technique are good estimators and can be used effectively in various fields, including reliability analysis, estimation processes, financial investment, forecasting and decision-making processes.
The rest of this paper is organized as follows. The expected LINEX and model description are discussed in Section 2. In Section 3, we obtain the Bayesian estimator of the exponential distribution using the expected LINEX loss function as well as the Bayes risk. Extensive simulation studies are presented in Section 4. Section 5 demonstrates the usefulness of the new estimators by analyzing three real data sets. Finally, the paper is concluded in Section 6.

2. Expected LINEX and Model Description

Let π ( φ ) be the probability density function (PDF) of φ with range ( , ) ; then, from (2), one can obtain the Bayesian estimator of the unknown parameter θ using the expected LINEX loss function as follows:
θ ˜ B E L = θ ˜ B L ( φ ) π ( φ ) d φ ,
where θ ˜ B L ( φ ) θ ˜ B L is the Bayesian estimator of the unknown parameter θ given by (2). One can choose a proper probability model π ( φ ) to reflect being symmetric, positive or negatively skewed, which will, in turn, give proper weights to balance overestimation and underestimation. It is noted that, compared with the conventional prefixed choice of the parameter φ , the Bayesian estimation (3) under expected LINEX loss takes all possible values of the parameter φ into account, and gives different probabilistic weights to each value. Then, the overall balance between overestimation and underestimation is characterized by the weighted sum in the full range ( , ) of the parameter φ , which may be more objective and robust than the conventional approach in Bayesian analysis.
Let X be a random variable following the exponential distribution with parameter β with the PDF
f ( x ; β ) = β e β x , x > 0 , β > 0 ,
where β is the scale parameter. Suppose that X 1 , X 2 , , X n is a random sample of size n from the exponential distribution (4), and x 1 , x 2 , , x n is the associated observation of X 1 , X 2 , , X n . Then, the likelihood function of β can be written as
L ( β | d a t a ) = β n e β i = 1 n x i .
By direct computation, it is known that the maximum likelihood estimator (MLE) of the parameter β is
β ^ = n i = 1 n x i .
To obtain the Bayesian estimator of the parameter β , we need to determine the prior distribution, say g ( β ) , of β . In this paper, a conjugate prior distribution of β is used as the gamma distribution with density function
π ( β ) β a 1 e b β , β > 0 .
It is seen that the prior distribution in (7) is a natural conjugate prior because it generates a posterior distribution that relates to the same family. Furthermore, the posterior distribution generated from the natural conjugate prior can be easily mathematically manipulated. The researcher can select the values of the hyper-parameters a and b to accommodate the prior knowledge. Moreover, the Jeffreys prior can be derived as a special case of (7) by replacing a = b = 0 ; for more details, see Kundu and Howlader [26]. Then, the posterior distribution of β is given by
g ( β | d a t a ) β n + a 1 e β b + i = 1 n x i .
Under the symmetric squared error (SE) loss function, the Bayesian estimator of β is the posterior mean, i.e.,
β ˜ B S = n + a b + i = 1 n x i .
Similarly, the Bayesian estimator of β using the LINEX loss function can be written from (2) as
β ˜ B L = n + a φ log 1 + φ b + i = 1 n x i .
Further, to obtain the Bayesian estimator of β using the expected LINEX loss function, we propose the following three probabilistic distributions for the parameter φ ( , )  as
π 1 ( φ ) = c e φ 1 + e φ c + 1 , c > 0 , π 2 ( φ ) = 1 2 π e ( φ μ ) 2 2 , < μ < , π 3 ( φ ) = exp e φ d φ + d , < d < .
From (11), it is noted that π 1 ( φ ) is the generalized logistic type-I (this type has also been called the “skew-logistic”) distribution with shape parameter c, π 2 ( φ ) is the normal distribution with scale parameter 1 and location parameter μ , and π 3 ( φ ) is the standard Gumbel or extreme value distribution with location parameter d. It is seen that, in different distributions π 1 ( φ ) , π 2 ( φ ) and π 3 ( φ ) , the parameter φ features different probability weights at different values in full range ( , ) , and that overestimation and underestimation could be well modeled by choosing constants c , μ and d properly. For illustration, the cumulative probability weight for the parameter φ is calculated and the associated plots are presented in Figure 2 with a change in parameters c , μ and d, respectively. One could see that when parameters c > 1   μ > 0 and d > 0.36651 , the weights of overestimation are more serious than underestimation using distributions π 1 ( φ ) , π 2 ( φ ) and π 3 ( φ ) , respectively. The distributions in (11) will be used to obtain the Bayesian estimator of the parameter β using the expected LINEX loss function in the next section.
Remark 1.
It should be mentioned here that the choice of the distribution for the parameter φ is a tricky issue in practice. The main principle is that the model of φ should be a distribution with full range ( , ) . On this basis, due to historical information and expert experience, among others, if there are sufficient extra data used regarding φ, one could use a goodness-of-fit technique to find a proper model from among different candidates, whereas, if there is scarce information about φ’s distributional characteristics, one could use a skewed distribution with concise expression from an ease of computation perspective. Further, one could also refer to the maximum likelihood and moment methods to give proper estimates for parameters c , μ and d, respectively, by using the information of φ in consequence.
Remark 2.
A proper probability model π ( ϕ ) provides weights to balance between overestimation and underestimation compared with fixed ϕ. In general, when the random model π ( ϕ ) is available, the optimal Bayesian estimator of θ from a decision-making perspective can be obtained as follows:
θ ˜ = min δ E θ E ϕ L ( θ , δ ) ,
where L ( · , · ) is the loss function, δ is the estimator of parameter θ, E θ ( · ) denotes the expectation with the posterior distribution of θ, and E ϕ denotes the expectation with the random model of ϕ. Under this approach, it is clear that the proposed expected Bayesian estimator can be viewed as the mixture of conventional Bayesian estimation under the weight π ( ϕ ) at each value of ϕ, and that the estimator refers to (3) under the expected LINEX loss from the decision-making procedure.

3. Bayesian Estimation Using Expected LINEX

In this section, the alternative Bayesian estimators of the exponential parameter β are proposed under the expected LINEX approach with the different probability distributions in (11), respectively. In addition, in order to investigate the behavior of the different estimators, the Bayesian risks are also presented in this section.

3.1. Bayesian Estimation Using Expected LINEX

Using the following series (e.g., formula 1.511 of Gradshteyn and Ryzhik [27]),
log ( 1 + z ) = j = 1 ( 1 ) j + 1 z j j , | z | < 1 ,
the β ˜ B L in (10) with respect to the LINEX loss function can be rewritten as
β ˜ B L = ( n + a ) j = 1 v j φ j 1 with v j = ( 1 ) j + 1 j b + i = 1 n x i j .
From (3), (11) and (12), the Bayesian estimator β ˜ 1 B E L of the parameter β using the expected LINEX loss with respect to π 1 ( φ ) can be obtained as follows:
β ˜ 1 B E L = ( n + a ) j = 1 v j φ j 1 c e φ 1 + e φ c + 1 d φ .
We use the generalized Hurwitz–Lerch Zeta function introduced by Goyal and Laddha [28], defined by
ψ M * ( u , v , r ) = 1 Γ ( v ) 0 t v 1 e r t 1 u e t M d t = k = 0 ( M ) k u k ( r + k ) v k ! ,
where M C , r C Z 0 , | u | < 1 , R { v M } > 1 and ( M ) k = Γ ( M + k ) Γ ( M ) denotes the Pochhammer symbol, with the convention ( M ) 0 = 1 ; see Srivastava et al. [29]. Hence, the Bayes estimator β ˜ 1 B E L takes the form
β ˜ 1 B E L = c ( n + a ) j = 1 v j Γ ( j ) ψ c + 1 * ( 1 , j , 1 ) + ( 1 ) j 1 ψ c + 1 * ( 1 , j , c ) , j > c .
Following the same approach, the Bayesian estimator β ˜ 2 B E L of the parameter β using the expected LINEX loss with respect to π 2 ( φ ) is given by
β ˜ 2 B E L = ( n + a ) j = 1 v j φ j 1 1 2 π e ( φ μ ) 2 2 d φ = ( n + a ) e μ 2 2 2 π j = 1 v j φ j 1 e φ 2 2 + μ φ d φ = ( n + a ) e 3 4 μ 2 j = 1 v j e π 2 ( j 1 ) D j 1 ( μ ) ,
where, D u ( y ) is the parabolic cylinder function (see Gradshteyn and Ryzhik [27], formula 9.24(0,1), page 1028).
Further, the Bayesian estimator β ˜ 3 B E L of the parameter β using π 3 ( φ ) under the expected LINEX loss can be expressed as
β ˜ 3 B E L = ( n + a ) j = 1 v j φ j 1 exp e φ d φ + d d φ = ( n + a ) j = 1 v j 0 d log ( y ) j 1 e y d y , ( with y = e φ + d ) .
Using the binomial expansion for u + v m = n = 0 m m n u m n v n in the last integral (17), one further has
β ˜ 3 B E L = ( n + a ) j = 1 m = 0 j 1 ( 1 ) m j 1 m v j d j m 1 0 log ( y ) m e y d y ,
Using the result by Nadarajah [30],
I ( m , ω ) = 0 [ log ( y ) ] m e y ( ω + 1 ) d y = s m ( ω + 1 ) s Γ ( s ) s = 1 .
Then, by combining (18) and (19), the estimator β ˜ 3 B E L reduces to
β ˜ 3 B E L = ( n + a ) j = 1 m = 0 j 1 ( 1 ) m j 1 m v j d j m 1 s m Γ ( s ) s = 1 .
Clearly, the Monte Carlo integration technique can be used to approximate the Bayesian estimators (15), (16) and (20) of the parameter β under the expected LINEX loss with respect to the different probability distributions in (11). In this case, we can obtain β ˜ k B E L , k = 1 , 2 , 3 based on β ˜ B L in (10) as follows:
β ˜ k B E L n + a N t = 1 N y k t 1 log 1 + y k t b + i = 1 n x i , k = 1 , 2 , 3 ,
where y k 1 , , y k N are the generated random numbers from π k ( φ ) , k = 1 , 2 , 3 , respectively.

3.2. Bayes Risks

In order to compare the behavior of different Bayesian estimators using the expected LINEX loss function, the Bayes risks are also provided in this section. We denote by β ˜ the estimator of β , and then the associated Bayes risk of β ˜ is defined as E β [ L ( β ˜ , β ) ] , where L ( β ˜ , β ) is any loss function and the expectation is taken over the posterior distribution of β . Therefore, the Bayes risk using the LINEX loss function can be obtained from (1) and (8) as
B R L = 0 e φ ( β ˜ B L β ) φ ( β ˜ B L β ) 1 β n + a 1 e β b + i = 1 n x i d β = φ n + a b + i = 1 n x i n + a φ log 1 + φ b + i = 1 n x i .
Correspondingly, the Bayes risk using the expected LINEX loss function can be defined as follows:
B R E L = B R L ( φ ) π ( φ ) d φ ,
where B R L ( φ ) B R L . Therefore, from (11) and (22), the Bayes risk of β ˜ 1 B E L using the expected LINEX loss function with respect to π 1 ( φ ) can be expressed as
B R 1 E L = φ n + a b + i = 1 n x i n + a φ log 1 + φ b + i = 1 n x i c e φ 1 + e φ c + 1 d φ = φ n + a b + i = 1 n x i ( n + a ) j = 1 v j φ j 1 c e φ 1 + e φ c + 1 d φ = c ( n + a ) b + i = 1 n x i ψ c + 1 * ( 1 , 2 , 1 ) ψ c + 1 * ( 1 , 2 , c ) c ( n + a ) j = 1 v j Γ ( j + 1 ) ψ c + 1 * ( 1 , j + 1 , 1 ) + ( 1 ) j ψ c + 1 * ( 1 , j + 1 , c ) , j > c 1 .
Similarly, using π 2 ( φ ) in (11) and (22), the Bayes risk of β ˜ 2 B E L using the expected LINEX loss function can be obtained as
B R 2 E L = φ n + a b + i = 1 n x i n + a φ log 1 + φ b + i = 1 n x i 1 2 π e ( φ μ ) 2 2 d φ = φ n + a b + i = 1 n x i ( n + a ) j = 1 v j φ j 1 1 2 π e ( φ μ ) 2 2 d φ = μ ( n + a ) b + i = 1 n x i ( n + a ) e 3 4 μ 2 j = 1 v j e π 2 j D j ( μ ) .
Finally, the Bayes risk of β ˜ 3 B E L using the expected LINEX loss function based on π 3 ( φ ) can be expressed as
B R 3 E L = φ n + a b + i = 1 n x i ( n + a ) j = 1 v j φ j 1 e e φ d φ + d d φ = ( n + a ) d + γ b + i = 1 n x i ( n + a ) j = 1 m = 0 j 1 ( 1 ) m j 1 m v j d j m 1 s m Γ ( s ) | s = 1 ,
where γ 0.5772156649 is known as Euler’s constant.
Another simple way to obtain the Bayes risk of β ˜ k B E L , k = 1 , 2 , 3 is to use the Monte Carlo integration routine. In this case, we have from (22) the following:
B R k E L 1 N t = 1 N y k t ( n + a ) b + i = 1 n x i ( n + a ) log 1 + y k t b + i = 1 n x i , k = 1 , 2 , 3 ,
where y k 1 , , y k N are the generated random numbers from π k ( φ ) , k = 1 , 2 , 3 , respectively.

4. Simulation Study

In this section, in order to investigate and compare the performance of different Bayesian estimators using the expected LINEX loss function, along with those based on SE and LINEX loss functions, extensive Monte Carlo simulations are carried out for illustration. The simulation procedure can be described by the following steps.
Step 1
Determine the value of the parameter β and the sample size n.
Step 2
Generate a random sample from the exponential distribution, where x i = log ( 1 u i ) β , i = 1 , , n where U u n i f o r m ( 0 , 1 ) .
Step 3
For given values of gamma hyper-parameters a and b, generate β from the gamma prior density (7).
Step 4
The Bayesian estimates β ˜ B S , β ˜ B L and the associated Bayesian estimates β ˜ 1 B E L , β ˜ 2 B E L , β ˜ 3 B E L using the expected LINEX loss of the parameter β are computed.
Step 5
Based on multiple replications, the performance of the different Bayesian estimates is evaluated in terms of criteria quantities of the mean relative estimate (MRE) and mean square error (MSE) as well as the Bayes risks.
In this procedure, the simulation is performed by considering different values of β = 0.5 , 2 and sample sizes n = 20 , 40 , 60 , 80 , 100 , respectively. Meanwhile, the gamma hyper-parameters are chosen in such a way that the prior mean becomes the expected value of the corresponding parameters. We assume a = 1 , b = 2 for β = 0.5 and a = 2 , b = 1 for β = 2 in this simulation. Further, the Bayesian estimates with respect to the LINEX loss function are calculated with selected φ = 1.5 in all the settings, whereas the Bayesian estimates using the expected LINEX loss function are obtained under different parameter values c , μ and d of distributions π j ( φ ) , j = 1 , 2 , 3 , respectively. For each π j ( φ ) , five configurations (Conf) are proposed, as shown in Table 1. As mentioned before, when c > 1 μ > 0 and d > 0.36651 , the weights of overestimation are more serious than underestimation employing the distributions π 1 ( φ ) , π 2 ( φ ) and π 3 ( φ ) , respectively. Therefore, the values in Table 1 are selected to reflect the impact of these parameters on the performance of the different estimates using the expected LINEX loss function. Based on 1000 replications, the criteria quantities MRE, MSE and the Bayes risks are displayed in Table 2, Table 3 and Table 4 for β = 0.5 and β = 2 , respectively.
From the listed results in Table 2, Table 3 and Table 4, it is noted that
  • With the increase in sample size n, the MREs tend towards one and the criteria MSEs and Bayesian risks of both Bayesian estimates under common and expected losses decrease, which indicates the consistency of the associated estimates when the sample sizes increase;
  • For a fixed sample size, the Bayesian estimates obtained under common LINEX loss have relatively smaller MSE than that from SE loss, whereas the associated Bayesian risks perform similarly under these two cases in general;
  • In terms of minimum MSE, we can conclude that the Bayesian estimates using the expected LINEX loss function perform better than all other estimates for π 3 ( φ ) with d = 5 when β = 0.5 ;
  • The Bayesian estimates using π 3 ( φ ) with d = 2 have the smallest MSE among all other estimates when β = 2 ;
  • The Bayesian estimates using the expected LINEX loss function using π 2 ( φ ) with μ = 0 have the minimum Bayes risk when compared with all other estimates.
It can be observed that when μ = 0 , π 2 ( φ ) reduces to the standard normal distribution. In this case, we can see that, although the distribution does not contain any parameter, it gives the best result based on the minimum Bayes risks in all cases. This is because, when we obtain the Bayesian estimate using the expected LINEX loss function, we take the expectation over all possible values for the parameter φ . The same result can also be noted when d = 0 in π 3 ( φ ) , where the Bayes risk obtained based on the expected LINEX loss function is smaller than that based on the common LINEX loss. Therefore, when the researcher does not have any information about the parameter φ , one could take the expectation over all possible values for the parameter φ to obtain the Bayesian estimate using the expected LINEX loss.

5. Data Analysis

In this section, three real life data sets are used to show the applicability of the proposed estimators based on the expected LINEX loss function.
Example One. (Software reliability data) The first data set represents the time-between-failures (time unit in milliseconds) of software reliability. The original data were firstly presented by Lyu [31] and analyzed recently by Alotaibi et al. [32] and are also available in the reliaR package. The complete data set that is transformed from time-between-failures to failure times is presented in Table A1. Before proceeding further, we first check the validity of the exponential distribution to model these data by using the Kolmogorov–Smirnov (K-S) distance and the associated p-value based on the MLE of the parameter β . By direct computation, it is noted that the MLE of β is 2.6240 × 10 3 , the K-S distance is 0.0788 and the corresponding p-value is 0.6591. Therefore, it is found that the exponential distribution could be used as an appropriate model to fit the mentioned software reliability data. For further illustration, the empirical cumulative distribution function (CDF) plot with overlaid theoretical exponential distribution and the probability–probability (P-P) plot are also shown in Figure 3 to check the goodness of fit for the software reliability data. It is also shown that the exponential distribution can be viewed as a proper model.
Based on the software reliability data set, to find the different Bayesian estimates, we select the hyper-parameters to be a = 0.011 and b = 2.25 in all the cases, where such choices are obtained by equating the first two moments of the gamma distribution obtained from the data with the corresponding theoretical moments. Using similarly designed scenarios given in Table 1, various Bayesian estimates are presented in Table 5 under common SE, LINEX and the proposed expected LINEX losses, respectively. Comparing the performance of the different estimates, we can conclude that the Bayesian estimate using the expected LINEX loss function using π 2 ( φ ) with μ = 0 performs better than other estimates in terms of minimum Bayes risk.
Example Two. (Breakdown of insulating fluid data) The second data set is taken from Nelson [33] and has been analyzed by many authors (see, e.g., Balakrishnan and Cramer [34], Dey and Nassar [35]). The data describe the times to breakdown of an insulating fluid subjected to a voltage of 34 kV, and the detailed lifetime data are presented in Table A2. By direct calculation, the MLE of the parameter β for these data is 0.06964 . The K-S distance and the corresponding p-value are obtained as 0.2464 and 0.1990, respectively. Therefore, the exponential distribution is a suitable model to fit these data. For comparison, the associated empirical CDF plot with overlaid theoretical exponential distribution and the P-P plot are provided in Figure 4 as well. They also indicate that the exponential distribution is an acceptable model to fit the data.
Using the breakdown data of an insulating fluid, and following a similar approach as for example one, different Bayesian estimates are obtained under various designed scenarios. The results are presented in Table 6, where the hyper-parameters are chosen as a = 0.543 and b = 0.035 by equating the first two moments of the gamma distribution obtained from the data with the corresponding theoretical moments. Based on the listed results in Table 6, one could also observe that the Bayesian estimate using the expected LINEX loss function using π 2 ( φ ) with μ = 0 has the smallest Bayes risk among all other estimates.
Example Three. (COVID-19 data) The third data set consists of the COVID-19 recovery cases in Pakistan from 24 March 2020 to 1 May 2020. This data set was originally analyzed by Ali et al. [36]. The complete data set is presented in Table A3. From the data, the MLE of β is 0.008282 . Based on this result, the K-S distance and the corresponding p-value are 0.1213 and 0.6149, respectively. Hence, the exponential distribution is an appropriate model to fit these data. Moreover, the empirical CDF plot with overlaid theoretical exponential distribution and the P-P plot are presented in Figure 5. Figure 5 shows that the exponential distribution is an adequate model to fit the data.
To acquire the different Bayesian estimates, we set the hyper-parameter values to be a = 0.0016 and b = 0.1230 , which are obtained by equating the first two moments of the gamma distribution got from the data with the associated theoretical moments. The different Bayesian estimates as well as the corresponding Bayes risks are tabulated in Table 7. From the detailed outcomes in Table 7, one could also observe that the Bayesian estimate using the expected LINEX loss function employing π 2 ( φ ) with μ = 0 has the smallest Bayes risk among all other estimates.
The simulation outcomes as well as the analysis of the three real data sets show that the performance of the Bayesian estimation using the expected LINEX loss function is superior to that of the common Bayesian approach in terms of MSE and Bayes risks. The proposed expected Bayesian estimation is a robust estimation approach that provides a weight estimation for the common Bayesian estimates. The results also show that if we select an appropriate probability distribution for φ , we can reduce the MSE and Bayes risk of an estimator. The superiority of the Bayesian estimators using the expected LINEX loss function can be interpreted based on the concept of the expected value, because this estimator is obtained as the arithmetic mean of a large number of independent realizations of φ . On the other hand, the usual Bayesian estimator using the LINEX loss function is obtained by choosing an arbitrary value for φ . The new offered technique can be utilized significantly in different fields; for instance, one can obtain accurate estimators for some reliability indices, including reliability and hazard rate function. Moreover, in regression analysis, new estimators of the regression models can be derived using the proposed technique. Furthermore, the expected LINEX loss function can be employed in volatility forecasting.

6. Conclusions

In this paper, a new approach to obtain the Bayesian estimator using the LINEX loss function has been introduced. For avoiding the subjective choice of the LINEX parameter and giving a relatively objective balance between overestimation and underestimation, the proposed alternative method assumes that the parameter φ of the LINEX loss function has a probability distribution in the range of ( , ) . The Bayesian estimator in this case is obtained by taking the expected value of the usual Bayesian estimator using the LINEX loss function over the proposed distribution of φ . When an exponential sample is available, the Bayesian estimator using the expected LINEX loss is presented for the unknown exponential parameter by considering three different distributions to the LINEX parameter φ . In addition, the Bayes risks are also obtained using the different distributions. To compare the performance of the proposed Bayesian estimators with those based on the squared error and LINEX loss functions, extensive simulation studies are carried out in terms of different criteria quantities. Moreover, three real-life data sets are analyzed to show the applicability of the new method. Based on the simulation and real data analysis results, we can conclude that the Bayesian estimates using the expected LINEX loss function perform better than those based on the squared error and LINEX loss functions in terms of minimum mean square error and Bayes risk. Finally, it is of interest to employ the same approach discussed in this paper to estimate the parameters of two-parameter distributions (or more) such as Weibull and gamma distributions. This work is in progress and will be reported later.

Author Contributions

Investigation, M.N., H.O., R.A. and L.W.; methodology, M.N. and L.W.; software, M.N., H.O. and L.W.; validation, M.N., R.A., H.O. and L.W.; writing, M.N., H.O. and L.W.; funding acquisition, R.A. and L.W. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R50), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. The work of Liang Wang was supported by the National Natural Science Foundation of China (No. 12061091), the Yunnan Fundamental Research Projects (No. 202101AT070103) and the Doctoral Research Foundation of Yunnan Normal University (No. 00800205020503129).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Three data sets are contained within the article.

Acknowledgments

The authors would like to express their gratitude to the editor and anonymous referees for their helpful recommendations and beneficial comments. Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R50), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. The work of Liang Wang was supported by the National Natural Science Foundation of China (No. 12061091), the Yunnan Fundamental Research Projects (No. 202101AT070103) and the Doctoral Research Foundation of Yunnan Normal University (No. 00800205020503129).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Software reliability data set.
Table A1. Software reliability data set.
4.797.4510.2215.7626.1028.5935.5241.4942.6644.3645.53
58.2762.9674.7081.63100.71102.06104.83110.79118.36122.73145.03
149.40152.80156.85162.20164.97168.60173.82179.95182.72195.72203.93
206.06222.26238.27241.25249.99256.17282.57282.62284.11294.45318.86
323.46329.11340.30344.67353.94398.56405.70407.51422.36429.93461.47
482.62491.46511.83526.64532.23537.13543.06560.75561.60589.96592.09
610.75615.65630.52673.74687.92698.15753.05768.25801.06828.22849.97
885.02892.27911.90951.69962.59965.04976.98986.921025.94
Table A2. Breakdown of insulating fluid data set.
Table A2. Breakdown of insulating fluid data set.
0.190.780.961.312.783.164.154.674.856.50
7.358.018.2712.0631.7532.5233.9136.7172.89
Table A3. COVID-19 data set.
Table A3. COVID-19 data set.
1322234262425194
408717238105155352646928368
19912067361029690181190228111
163204192627263400

References

  1. Zellner, A. Bayesian estimation and prediction using asymmetric loss functions. J. Am. Stat. Assoc. 1986, 81, 446–451. [Google Scholar] [CrossRef]
  2. Chang, Y.C.; Hung, W.L. LINEX loss functions with applications to determining the optimum process parameters. Qual. Quant. 2007, 41, 291–301. [Google Scholar] [CrossRef]
  3. Khatun, N.; Matin, M.A. A Study on LINEX Loss Function with Different Estimating Methods. Open J. Stat. 2020, 10, 52. [Google Scholar] [CrossRef] [Green Version]
  4. Ali, S.; Aslam, M.; Abbas, N.; Kazmi, S.M.A. Scale parameter estimation of the Laplace model using different asymmetric loss functions. Int. J. Stat. Probab. 2012, 1, 105–127. [Google Scholar] [CrossRef] [Green Version]
  5. Ali, S. On the Bayesian estimation of the weighted Lindley distribution. J. Stat. Comput. Simul. 2015, 85, 855–880. [Google Scholar] [CrossRef]
  6. Klebanov, L.B. Universal loss function and unbiased estimation. Dokl. Akad. Nank SSSR Soviet Math. Dokl. 1972, 203, 1249–1251. [Google Scholar]
  7. Varian, H.R. A Bayesian Approach to Real Estate Assessment. In Studies in Bayesian Econometrics and Statistics: In Honor of L. J. Savage; Savage, L.J., Feinberg, S.E., Zellner, A., Eds.; North-Holland Pub. Co.: Amsterdam, The Netherlands, 1975; pp. 195–208. [Google Scholar]
  8. Dey, D.; Lou, K.; Bose, S. A Bayesian approach to loss robustness. Stat. Decis. 1998, 16, 65–87. [Google Scholar] [CrossRef]
  9. Soliman, A.A. Comparison of linex and quadratic bayes estimators foe the rayleigh distribution. Commun. Stat.-Theory Methods 2000, 29, 95–107. [Google Scholar] [CrossRef]
  10. Misra, N.; van der Meulen, E.C. On estimating the mean of the selected normal population under the LINEX loss function. Metrika 2003, 58, 173–183. [Google Scholar] [CrossRef]
  11. Micheas, A.C. A unified approach to Prior and Loss Robustness. Commun. Stat.-Theory Methods 2006, 35, 309–323. [Google Scholar] [CrossRef]
  12. Hoque, Z.; Khan, S.; Wesolowski, J. Performance of preliminary test estimator under linex loss function. Commun. Stat.-Theory Methods 2009, 38, 252–261. [Google Scholar] [CrossRef] [Green Version]
  13. Kazmi, S.M.A.; Aslam, M.; Ali, S. Preference of prior for the class of life-time distributions under different loss functions. Pak. J. Stat. 2012, 28, 467–487. [Google Scholar]
  14. Pandey, M.; Singh, V.P.; Srivastava, C.P.L. A Bayesian estimation of reliability model using the LINEX loss function. Microelectron. Reliab. 1994, 34, 1519–1523. [Google Scholar] [CrossRef]
  15. Jaheen, Z.F. A Bayesian analysis of record statistics from the Gompertz model. Appl. Math. Comput. 2003, 145, 307–320. [Google Scholar] [CrossRef]
  16. Azimi, R.; Yaghmaei, F.; Azimi, D. Comparison of Bayesian estimation methods for Rayleigh progressive censored data under the different asymmetric loss function. Int. J. Appl. Math. Res. 2012, 1, 452–461. [Google Scholar] [CrossRef] [Green Version]
  17. Ahmed, E.A. Bayesian estimation based on progressive Type-II censoring from two parameter bathtub-shaped lifetime model: An Markov chain Monte Carlo approach. J. Appl. Stat. 2014, 41, 752–768. [Google Scholar] [CrossRef]
  18. Ahmed, E.A. Estimation of some lifetime parameters of generalized Gompertz distribution under progressively type-II censored data. Appl. Math. Model. 2015, 39, 5567–5578. [Google Scholar] [CrossRef]
  19. Ashour, S.K.; Nassar, M. Inference for Weibull distribution under adaptive type-I progressive hybrid censored competing risks data. Commun. Stat.-Theory Methods 2017, 46, 4756–4773. [Google Scholar] [CrossRef]
  20. Nassar, M.; Abo-Kasem, O.; Zhang, C.; Dey, S. Analysis of Weibull distribution under adaptive type-II progressive hybrid censoring scheme. J. Indian Soc. Probab. Stat. 2018, 19, 25–65. [Google Scholar] [CrossRef]
  21. Kotb, M.S.; Raqab, M.Z. Statistical inference for modified Weibull distribution based on progressively type-II censored data. Math. Comput. Simul. 2019, 162, 233–248. [Google Scholar] [CrossRef]
  22. Chen, S.; Gui, W. Statistical analysis of a lifetime distribution with a bathtub-shaped failure rate function under adaptive progressive type-II censoring. Mathematics 2020, 8, 670. [Google Scholar] [CrossRef]
  23. Nassar, M.; Okasha, H.; Albassam, M. E-Bayesian estimation and associated properties of simple stepstress model for exponential distribution based on type-II censoring. Qual. Reliab. Eng. Int. 2020, 37, 997–1016. [Google Scholar] [CrossRef]
  24. Tang, J.; Xu, W.; Li, J.; Tian, Y.; Xu, S. Multi-view learning methods with the LINEX loss for pattern classification. Multi-view learning methods with the LINEX loss for pattern classification. Knowl.-Based Syst. 2021, 228, 107285. [Google Scholar] [CrossRef]
  25. Abolmohammadi, M.; Seif, A.; Behzadi, M.H.; Moghadam, M.B. Effect of Linex loss function on the VSIX control chart. J. Stat. Comput. Simul. 2019, 89, 1674–1693. [Google Scholar] [CrossRef]
  26. Kundu, D.; Howlader, H. Bayesian inference and prediction of the inverse Weibull distribution for Type-II censored data. Comput. Stat. Data Anal. 2010, 54, 1547–1558. [Google Scholar] [CrossRef]
  27. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, 7th ed.; Academic Press: London, UK, 2007. [Google Scholar]
  28. Goyal, S.P.; Laddha, R.K. On the generalized Riemann zeta functions and the generalized Lambert transform. Ganita Sandesh 1997, 11, 99–108. [Google Scholar]
  29. Srivastava, H.M.; Saxena, R.K.; Pogány, T.K.; Saxena, R. Integral and computational representations of the extended Hurwitzerch Zeta function. Integral Transform. Spec. Funct. 2011, 22, 487–506. [Google Scholar] [CrossRef]
  30. Nadarajah, S. The exponentiated Gumbel distribution with climate application. Environ. Off. J. Int. Environ. Soc. 2006, 17, 13–23. [Google Scholar] [CrossRef]
  31. Lyu, M.R. Handbook of Software Reliability Engineering; IEEE Computer Society Press: Los Alamitos, CA, USA, 1996. [Google Scholar]
  32. Alotaibi, R.; Okasha, H.; Rezk, H.; Nassar, M. A new weighted version of alpha power transformation method: Properties and applications to COVID-19 and software reliability data. Phys. Scr. 2021, 96, 125221. [Google Scholar] [CrossRef]
  33. Nelson, W.B. Accelerated Testing: Statistical Model, Test Plan and Data Analysis; Wiley: New York, NY, USA, 2005. [Google Scholar]
  34. Balakrishnan, N.; Cramer, E. The Art of Progressive Censoring: Applications to Reliability and Quality; Statistics for Industry and Technology; Springer: New York, NY, USA; Heidelberg, Germany; Dordrecht, The Netherlands; London, UK, 2014. [Google Scholar]
  35. Dey, S.; Nassar, M. Classical methods of estimation on constant stress accelerated life tests under exponentiated Lindley distribution. J. Appl. Stat. 2020, 47, 975–996. [Google Scholar] [CrossRef]
  36. Ali, A.; Mashwani, W.K.; Naeem, S.; Uddin, M.I.; Kumam, W.; Kumam, P.; Alrabaiah, H.; Jamal, F.; Chesneau, C. COVID-19 infected lung computed tomography segmentation and supervised classification approach. Comput. Mater. Contin. 2021, 68, 391–407. [Google Scholar] [CrossRef]
Figure 1. Plots of LINEX loss function for: (a) positive values of φ , (b) negative values of φ and (c) small values of φ .
Figure 1. Plots of LINEX loss function for: (a) positive values of φ , (b) negative values of φ and (c) small values of φ .
Mathematics 10 00436 g001
Figure 2. Plots of cumulative probability for π 1 ( φ ) , π 2 ( φ ) and π 3 ( φ ) .
Figure 2. Plots of cumulative probability for π 1 ( φ ) , π 2 ( φ ) and π 3 ( φ ) .
Mathematics 10 00436 g002
Figure 3. Empirical distribution and P-P plots of the exponential distribution under software reliability data.
Figure 3. Empirical distribution and P-P plots of the exponential distribution under software reliability data.
Mathematics 10 00436 g003
Figure 4. Empirical distribution and P-P plots of the exponential distribution under insulating fluid data.
Figure 4. Empirical distribution and P-P plots of the exponential distribution under insulating fluid data.
Mathematics 10 00436 g004
Figure 5. Empirical distribution and P-P plots of the exponential distribution under COVID-19 data.
Figure 5. Empirical distribution and P-P plots of the exponential distribution under COVID-19 data.
Mathematics 10 00436 g005
Table 1. Selected parameter values of π j ( φ ) , j = 1 , 2 , 3 .
Table 1. Selected parameter values of π j ( φ ) , j = 1 , 2 , 3 .
DistributionParameterConf
12345
π 1 ( φ ) c0.20.5125
π 2 ( φ ) μ −2−0.5025
π 3 ( φ ) d−2−0.5025
Table 2. The MREs of the different estimators for β = 0.5 and 2.
Table 2. The MREs of the different estimators for β = 0.5 and 2.
MethodConf n = 20 n = 40 n = 60 n = 80 n = 100
β = 0.5
β ˜ B S -1.04231.02181.00831.01321.0087
β ˜ B L -1.02251.01211.00201.00841.0050
β ˜ 1 B E L 11.11211.04721.02181.02151.0138
21.06341.03131.01441.01781.0123
31.04311.02201.00841.01331.0088
41.02951.01541.00411.01011.0062
51.01551.00850.99961.00661.0035
β ˜ 2 B E L 11.07071.03511.01691.01971.0139
21.04941.02511.01041.01481.0100
31.04641.01891.01491.01181.0044
41.01641.00900.99991.00691.0037
50.98030.99060.98770.99750.9963
β ˜ 3 B E L 11.06251.03131.01441.01781.0124
21.04161.02141.00801.01301.0086
31.03871.01521.01241.01001.0029
41.00931.00540.99761.00511.0023
50.97390.98710.98540.99580.9948
β = 2
β ˜ B S -1.04101.02241.01781.00991.0090
β ˜ B L -0.97100.98600.99320.99150.9942
β ˜ 1 B E L 11.14520.97120.93310.91090.9019
21.14211.05671.03631.02121.0165
31.05351.02521.01901.01051.0094
40.99950.99941.00190.99790.9993
50.95180.97390.98450.98470.9886
β ˜ 2 B E L 11.16691.07831.05381.03621.0299
21.07161.03631.02681.01651.0142
31.04021.03191.02001.00401.0099
40.95280.97530.98570.98580.9895
50.84770.91370.94170.95200.9619
β ˜ 3 B E L 11.12991.06211.04341.02861.0239
21.04221.02181.01711.00931.0085
31.01291.01761.01050.99701.0041
40.93050.96300.97700.97910.9841
50.81970.90320.93400.94600.9569
Table 3. The MSEs of the different estimators for β = 0.5 and 2.
Table 3. The MSEs of the different estimators for β = 0.5 and 2.
MethodConf n = 20 n = 40 n = 60 n = 80 n = 100
β = 0.5
β ˜ B S -0.01320.00660.00420.00350.0023
β ˜ B L -0.01190.00620.00410.00340.0022
β ˜ 1 B E L 10.02070.00780.00470.00370.0024
20.01490.00690.00440.00360.0023
30.01320.00660.00420.00350.0023
40.01230.00630.00420.00340.0022
50.01150.00610.00410.00330.0022
β ˜ 2 B E L 10.01550.00710.00450.00360.0023
20.01370.00670.00430.00350.0023
30.01320.00660.00410.00310.0023
40.01150.00610.00410.00330.0022
50.01000.00570.00390.00320.0021
β ˜ 3 B E L 10.01480.00690.00440.00360.0023
20.01310.00660.00420.00350.0023
30.01260.00650.00400.00310.0023
40.01120.00600.00400.00330.0022
50.00980.00560.00390.00320.0021
β = 2
β ˜ B S -0.1890.102510.06510.04650.0389
β ˜ B L -0.14130.08760.05800.04310.0365
β ˜ 1 B E L 10.53790.11720.07930.07290.0716
20.38100.13120.07540.05090.0415
30.21000.10490.06580.04670.0390
40.15760.09200.06000.99790.0371
50.13890.08570.05690.04260.0361
β ˜ 2 B E L 10.41460.14990.08530.05650.0455
50.22890.11170.06910.04840.0402
30.20270.09630.06630.04700.0382
40.13740.08570.05690.04270.0361
50.17110.09420.06030.04560.0376
β ˜ 3 B E L 10.33060.13340.07840.05300.0432
20.19280.10250.06490.04640.0388
30.17650.08850.06280.04570.0371
40.13430.08470.05630.04250.0359
50.20650.09900.06270.04720.0386
Table 4. The Bayes risks of the different estimators for β = 0.5 and 2.
Table 4. The Bayes risks of the different estimators for β = 0.5 and 2.
MethodConf n = 20 n = 40 n = 60 n = 80 n = 100
β = 0.5
β ˜ B S -0.01350.00650.00420.00320.0025
β ˜ B L -0.01480.00720.00470.00360.0028
β ˜ 1 B E L 10.40610.16170.10030.07450.0583
20.06440.02910.01860.01400.0110
30.02240.01070.00700.00530.0042
40.02120.01050.00690.00520.0041
50.03940.01960.01290.00980.0078
β ˜ 2 B E L 10.03570.01670.01080.00810.0064
20.00870.00410.00270.00200.0016
30.00680.00320.00210.00160.0013
40.03220.01590.01040.00790.0063
50.16060.08100.05340.04080.0324
β ˜ 3 B E L 10.02580.01220.00790.00590.0047
20.01090.00530.00350.00260.0021
30.01290.00630.00420.00310.0025
40.05250.02620.01720.01310.0104
50.19970.10140.06700.05120.0408
β = 2
β ˜ B S -0.20530.10190.06790.05030.0403
β ˜ B L -0.20980.10920.07390.05520.0444
β ˜ 1 B E L 13.17841.08460.65220.46150.3603
21.19290.45690.28410.20420.1608
30.37390.17010.11210.08280.0663
40.28750.15260.10450.07870.0636
50.51770.28080.19390.14660.1187
β ˜ 2 B E L 10.64780.28260.18130.13210.1047
20.14300.06690.04380.03220.0257
30.10380.05190.03410.02490.0202
40.43380.23310.15970.12020.0971
51.99621.12510.78830.60060.4890
β ˜ 3 B E L 10.44780.20190.13080.09570.0761
20.15760.08020.05420.04050.0326
30.17490.09400.06360.04710.0385
40.69670.37480.25910.19590.1588
52.69161.38730.97770.74740.6099
Table 5. Different estimates of β using the software reliability data.
Table 5. Different estimates of β using the software reliability data.
MethodConfEstimateBayes Risk
β ˜ B S - 2.62418 × 10 3 8.00631 × 10 8
β ˜ B L - 2.62412 × 10 3 9.00682 × 10 8
β ˜ 1 B E L 1 2.26912 × 10 3 6.51752 × 10 7
2 2.60649 × 10 3 3.00326 × 10 7
3 2.62394 × 10 3 1.31255 × 10 7
4 2.62390 × 10 3 1.31247 × 10 7
5 2.62350 × 10 3 2.47330 × 10 7
β ˜ 2 B E L 1 2.62426 × 10 3 2.00169 × 10 7
2 2.62420 × 10 3 5.00407 × 10 8
3 2.62418 × 10 3 4.00315 × 10 8
4 2.62410 × 10 3 2.00146 × 10 7
5 2.62398 × 10 3 1.04071 × 10 6
β ˜ 3 B E L 1 2.62423 × 10 3 1.46892 × 10 7
2 2.62417 × 10 3 6.60856 × 10 8
3 2.62415 × 10 3 7.91824 × 10 8
4 2.62407 × 10 3 3.31714 × 10 7
5 2.62395 × 10 3 1.31088 × 10 6
Table 6. Different estimates of β using the breakdown of insulating fluid data.
Table 6. Different estimates of β using the breakdown of insulating fluid data.
MethodConfEstimateBayes Risk
β ˜ B S - 0.07163 2.62502 × 10 4
β ˜ B L - 0.07143 2.94237 × 10 4
β ˜ 1 B E L 1 0.06229 2.17113 × 10 3
2 0.07131 9.94977 × 10 4
3 0.07162 4.30383 × 10 4
4 0.07149 4.27258 × 10 4
5 0.07134 8.03771 × 10 4
β ˜ 2 B E L 1 0.07189 6.60783 × 10 4
2 0.07169 1.64589 × 10 4
3 0.07163 1.31254 × 10 4
4 0.07136 6.51803 × 10 4
5 0.07098 3.36830 × 10 3
β ˜ 3 B E L 1 0.07181 4.84019 × 10 4
2 0.07162 2.15802 × 10 4
3 0.07155 2.57904 × 10 4
4 0.07129 1.07746 × 10 3
5 0.07090 4.23444 × 10 3
Table 7. Different estimates of β using COVID-19 data.
Table 7. Different estimates of β using COVID-19 data.
MethodConfEstimateBayes Risk
β ˜ B S - 8.28214 × 10 3 1.75874 × 10 6
β ˜ B L - 8.28082 × 10 3 1.97817 × 10 6
β ˜ 1 B E L 1 7.16357 × 10 3 1.43283 × 10 5
2 8.22729 × 10 3 6.60063 × 10 6
3 8.28139 × 10 3 2.88328 × 10 6
4 8.28051 × 10 3 2.88206 × 10 6
5 8.27843 × 10 3 5.43067 × 10 6
β ˜ 2 B E L 1 8.28390 × 10 3 4.39860 × 10 6
2 8.28258 × 10 3 1.09942 × 10 6
3 8.28214 × 10 3 8.79372 × 10 7
4 8.28038 × 10 3 4.39512 × 10 6
5 8.27774 × 10 3 2.28462 × 10 5
β ˜ 3 B E L 1 8.28339 × 10 3 3.22757 × 10 6
2 8.28207 × 10 3 1.45140 × 10 6
3 8.28163 × 10 3 1.73882 × 10 6
4 8.27987 × 10 3 7.28332 × 10 6
5 8.27724 × 10 3 2.87744 × 10 5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nassar, M.; Alotaibi, R.; Okasha, H.; Wang, L. Bayesian Estimation Using Expected LINEX Loss Function: A Novel Approach with Applications. Mathematics 2022, 10, 436. https://doi.org/10.3390/math10030436

AMA Style

Nassar M, Alotaibi R, Okasha H, Wang L. Bayesian Estimation Using Expected LINEX Loss Function: A Novel Approach with Applications. Mathematics. 2022; 10(3):436. https://doi.org/10.3390/math10030436

Chicago/Turabian Style

Nassar, Mazen, Refah Alotaibi, Hassan Okasha, and Liang Wang. 2022. "Bayesian Estimation Using Expected LINEX Loss Function: A Novel Approach with Applications" Mathematics 10, no. 3: 436. https://doi.org/10.3390/math10030436

APA Style

Nassar, M., Alotaibi, R., Okasha, H., & Wang, L. (2022). Bayesian Estimation Using Expected LINEX Loss Function: A Novel Approach with Applications. Mathematics, 10(3), 436. https://doi.org/10.3390/math10030436

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