Next Article in Journal
A Non-Equilibrium Interpolation Scheme for IB-LBM Optimized by Approximate Force
Next Article in Special Issue
Consistency of the Estimator for the Common Mean in Fixed-Effect Meta-Analyses
Previous Article in Journal
Cooperative Games Based on Coalition Functions in Biform Games
Previous Article in Special Issue
Different Estimation Methods for New Probability Distribution Approach Based on Environmental and Medical Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bayesian Inference and Data Analysis of the Unit–Power Burr X Distribution

1
Department of Statistics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Faculty of Graduate Studies for Statistical Research (FGSSR), Cairo University, Giza 12613, Egypt
3
Faculty of Business Administration, Delta University for Science and Technology, Gamasa 11152, Egypt
*
Authors to whom correspondence should be addressed.
Axioms 2023, 12(3), 297; https://doi.org/10.3390/axioms12030297
Submission received: 8 February 2023 / Revised: 11 March 2023 / Accepted: 12 March 2023 / Published: 14 March 2023

Abstract

:
The unit–power Burr X distribution (UPBXD), a bounded version of the power Burr X distribution, is presented. The UPBXD is produced through the inverse exponential transformation of the power Burr X distribution, which is also beneficial for modelling data on the unit interval. Comprehensive analysis of its key characteristics is performed, including shape analysis of the primary functions, analytical expression for moments, quantile function, incomplete moments, stochastic ordering, and stress–strength reliability. Rényi, Havrda and Charvat, and d-generalized entropies, which are measures of uncertainty, are also obtained. The model’s parameters are estimated using a Bayesian estimation approach via symmetric and asymmetric loss functions. The Bayesian credible intervals are constructed based on the marginal posterior distribution. Monte Carlo simulation research is intended to test the accuracy of various estimators based on certain measures, in accordance with the complex forms of Bayesian estimators. Finally, we show that the new distribution is more appropriate than certain other competing models, according to their application for COVID-19 in Saudi Arabia and the United Kingdom.

1. Introduction

Utilizing differential equations, Burr [1] introduced twelve distributions. In the literature, Burr type XII distributions and single-parameter Burr type X, have drawn much interest. Surles and Padgett [2] have proposed the two-parameter Burr X distribution (BXD), often known as the generalized Rayleigh distribution. For data modelling, the BXD can be used as an alternative to the Weibull and Rayleigh distributions. However, the model has a considerable impact on the prediction of failure rates and has generated a lot of interest in modelling across a wide range of disciplines, including hydrology, medicine and reliability analysis. The cumulative distribution function (CDF) of the BXD is given by:
G ( x ) = 1 exp ( a x b ) 2 ;     x > 0 ,
where, b > 0, and a > 0 are the shape and scale parameters, respectively. The probability density function (PDF) associated with (1) is given by:
g ( x ) = 2 a 2 b x 2 b 1 exp ( a x b ) 2 ;     x > 0 .
According to Raqab and Kundu [3], the shape parameter (b) determines whether the hazard rate function (HF) of the BXD is a bathtub or an increasing function. The HF is bathtub for b 1 , and is an increasing function for b > 1 . In the literature, numerous studies have been undertaken in recent years to create modified or generalized forms of the BXD in order to increase the viability of BXDs, see, for example, [4,5,6,7,8,9]. Our focus is on the recently established power BXD (PBXD) by Usman and Ilyas [10], with an additional shape parameter that depends on the transformation Y = 1 / X 1 / δ , δ > 0 . The CDF and PDF of the PBXD are, respectively, given by:
G ( y ) = 1 exp ( a y b ) 2   δ ;   y > 0 ,
g ( y ) = 2 δ a 2 b y 2 b 1 e ( a y b ) 2 1 exp ( a y b ) 2 δ 1 ;     y > 0 .
For δ = 1 , the CDF (3) reduces to BXD. Usman and Ilyas [10] mentioned that, subject to certain restrictions, their model can handle both symmetrical and heavy-tailed skewed data sets.
A significant challenge in data modelling is the selection of an adequate lifetime probability. However, over time, a variety of probability models have been widely proposed for the analysis of data sets in a variety of fields, including the medical sciences, actuarial sciences, engineering, finance and insurance, demography, biological sciences, and economics. In many practical scenarios, we are required to deal with the uncertainty of bounded situations. We commonly encounter variables that fall within the range of (0, 1), such as the percentage of a particular trademark, the results of some capacity tests, different lists, and rates. In order to model these variables effectively, continuous unit distributions, or probability distributions with support for (0, 1), are crucial. Due to this, some authors have recently concentrated on the creation of distributions that are specified on the bounded interval using any one of the parent distribution modification strategies. Among distributions that are specified in the (0, 1) interval, the beta distribution is obviously the most well-known. The beta distribution is helpful for simulating data on the unit interval, but different distributions have also been proposed and researched over time. The Topp–Leone distribution (see [11]) and the Kumaraswamy distribution (see [12]) can all be used as examples by the reader. The idea of offering distributions defined by the unit interval corresponding to any continuous distribution, however, has recently attracted the interest of statisticians. The following are a few of the most practical unit–interval distributions: the log–Lindley (Gómez-Déniz et al. [13]), unit–Birnbaum–Saunders (Mazucheli et al. [14]), unit–inverse Gaussian (Ghitany et al. [15]), unit–Lindley (Mazucheli et al. [16]), unit–BurrIII (Modi and Gill [17]), unit–Weibull (Mazucheli et al. [18]), unit–Burr XII (Korkmaz and Chesneau [19]), unit–odd Fréchet power function (Haq et al. [20]), unit–Teissier (Krishna et al. [21]), unit–exponentiated exponential (Jha et al. [22]) and unit–exponentiated half-logistic (Hassan et al. [23]) among others.
In this study, we propose a new unit probability distribution, based on the PBXD, that has three parameters. A new unit-PBXD (UPBXD) is provided based on the transformation W = e Y , where Y represents the PBXD. The UPBXD has the following desirable characteristics:
The UPBXD is a flexible model and can be used to describe a variety of datasets with a range between zero and one.
The new density function of the UBBXD takes several shapes, including unimodal, reversed J-shaped, U-shaped, left-skewed, and symmetric (see Section 2).
The HF shapes of the UPBXD can be increasing, J-shaped, or bathtub (U-HF) (see Section 2).
We derive some of the most important statistical characteristics of the UPBXD, such as the analytical expression for moments, the quantile function, incomplete moments, stochastic ordering, some uncertainty measures, and stress–strength reliability.
The parameter estimators of the UPBXD are explored using a Bayesian technique. The Bayesian credible intervals are also created.
To examine the effectiveness of estimators based on accuracy criteria, an exclusive simulation study was conducted.
Application to COVID-19 datasets from Saudi Arabia and the United Kingdom are used to show the superiority of the proposed model over other well-known models.
An outline of the paper’s structure is provided. Section 2 provides a definition of the suggested distribution. The distributional characteristics of the UPBXD are covered in Section 3. The maximum likelihood (ML) and Bayesian estimators utilizing various loss functions are covered in Section 4. The effectiveness of the suggested point and interval estimators is assessed using a Monte Carlo simulation in Section 5. Section 6 shows that the UPBXD outperforms the other unit distributions when employed with COVID-19 data. The paper conclusion is completed in Section 7.

2. Unit Power Burr X Distribution

In this section, we present the UPBXD, which results from the transformation of the type W = e Y , where Y is the PBXD and is a new bounded distribution with support on (0, 1). Thus, the following is how the CDF of the PBXD can be obtained:
F ( w ) = P ( W w ) = P ( e Y w ) = P ( Y ln ( w ) ) = 1 P ( Y ln ( w ) ) = 1 F Y ( ln ( w ) ) ,
which gives
F ( w ) = 1 1 exp a ( ln w ) b 2 δ ;     0 < w < 1 , a , b , δ > 0 .
Based on (5), we have F ( w ) = 0 , for w ≤ 0, and F ( w ) = 1 , for w ≤ 1. The PDF of the UPBXD related to (5) can be acquired as follows:
f ( w ) = 2 a 2 b δ w 1 ( ln w ) 2 b 1 e a ( ln w ) b 2 1 exp a ( ln w ) b 2 δ 1 ;     0 < w < 1 .
A random variable with PDF (6) is represented by UPBXD ( a , b , δ ) . For b = 1, the PDF (6) gives UBXD as a new sub-model. The following is the HF of the UPBXD:
h ( w ) = 2 a 2 b δ w 1 ( ln w ) 2 b 1 e a ( ln w ) b 2 1 exp a ( ln w ) b 2 1 .
The related plots for various selections of the parameters a , b , and δ are shown in Figure 1 and Figure 2 to provide a general overview of the shapes of the PDF (2) and HF (7).
In Figure 1, the PDF graphs for various parameter combinations display a variety of shapes, such as (a = 2, b = 2) symmetric normal, (a = 0.5, b = 0.3) U-shaped, (a = 0.5, b = 2) right-skewed, (a = 2, b = 0.3) J-shaped, and (a = 2, b = 2) normal tapered. In Figure 2, the UPBXD’s HF shapes in (a = 0.5, b = 2), (a = 2, b = 0.3), and (a = 2, b = 2) have increasing and J shapes, while (a = 0.5, b = 0.3) has a bathtub shape.
The parameter δ is responsible for the bathtub shapes given that the other two parameters (a and b) are less than one. The δ parameter is responsible for the J shapes where a > 1 and b < 1.
By inverting (5), we can get the quantile function (QF) of the UPBXD, which looks like this:
w q = exp 1 a ln 1 1 q 1 / δ 1 2 1 / b , 0 < q < 1 ,
where q is the uniform random variables. The first, median, and third quantiles are produced by setting q = 0.25, 0.5, and 0.75 in (8). It is simple to simulate the random variable of the UPBXD from (8).

3. The UPBXD’s Properties

In this section, we examine aspects of the UPBXD’s structural characteristics, such as some moment’s measures, information measures, stochastic ordering (SO), and stress–strength (SS) reliability.

3.1. Some Moments Measures

The mth moment for W~UPBXD ( a , b , δ ) , is determined as follows:
μ m = 2 a 2 b δ 0 1 w m 1 ( ln w ) 2 b 1 e a ( ln w ) b 2 1 exp a ( ln w ) b 2 δ 1 d w .
Using the binomial expansion in (9) provides
μ m = 2 a 2 b δ j = 0 ( 1 ) j δ 1 j 0 1 w m 1 ( ln w ) 2 b 1 e ( j + 1 ) a ( ln w ) b 2 d w .
Let y = a ( ln w ) b 2 , then the mth moment of W, is given by
μ m = δ j = 0 ( 1 ) j δ 1 j 0 e m a 1 b y 1 / 2 b e ( j + 1 ) y d y .
Use the exponential expansion then μ m , obtains the following form:
μ m = j , k = 0 ( 1 ) j + k m k δ k ! a k b δ 1 j 0 y k 2 b e ( j + 1 ) y d y = j , k = 0 ( 1 ) j + k δ m k k ! ( j + 1 ) ( k / 2 b ) + 1 a k b δ 1 j Γ k 2 b + 1 ,
where, Γ ( . ) is a gamma function. Furthermore, the mth central moment of W, is defined by
μ m = E ( W μ 1 ) m = i = 0 m ( 1 ) i m i ( μ 1 ) i μ m i .
Some moments measures including, first four moments, variance ( σ 2 ), coefficient of skewness ( α 3 ) and coefficient of kurtosis ( α 4 ) for the UPBXD are calculated for specific parameter values. Table 1 provides these measures considering parameter values as: (i) ( a = 1.5 , b = 1.3 , δ = 1.4 ) , (ii) ( a = 0.5 , b = 0.4 , δ = 0.4 ) , (iii) ( a = 4 , b = 2 , δ = 0.7 ) , (iv) ( a = 0.7 , b = 0.4 , δ = 2 ) , (v) ( a = 1.5 , b = 0.5 , δ = 0.5 ) , (vi) ( a = 5 , b = 1.6 , δ = 0.4 ) , and (vii) ( a = 1.5 , b = 1.5 , δ = 3 ) .
Table 1 displays that the UPBXD is right- and left- skewed in accordance with the values of α 3 . Additionally, the distribution is leptokurtic and platykurtic according to the values of α 4 .  Figure 3 shows the 3-dimensional plots for coefficient of skewness and kurtosis for UPBXD with different values of parameters. Looking at Figure 3, we can see that the coefficient of skewness and kurtosis increases when b and δ increases, while a increases then the coefficient of skewness decreases and coefficient of kurtosis increases.
Furthermore, the mth lower incomplete moment, say υ m ( x ) , of the UPBXD is given by:
υ m ( x ) = 2 a 2 b δ 0 x w m 1 ( ln w ) 2 b 1 e a ( ln w ) b 2 1 exp a ( ln w ) b 2 δ 1 d w .
Let z = a ( ln w ) b 2 , and using the binomial expansion, then the mth incomplete moment of W is
υ m ( x ) = δ j = 0 ( 1 ) j δ 1 j ( a ( ln x ) b ) 2 e m a 1 b y 1 / 2 b e ( j + 1 ) y d y .
Using exponential expansion and after simplification, the mth moment is as below:
υ m ( x ) = j , k = 0 ( 1 ) j + k δ m k k ! ( j + 1 ) ( k / 2 b ) + 1 a k b δ 1 j γ k 2 b + 1 , ( j + 1 ) ( a ( ln x ) b ) 2 ,
where γ ( . , x ) is an upper incomplete gamma function. The Lorenz and Bonferroni curves are well-known applications of the first incomplete moment. In the fields of economics, demographics, insurance, engineering, and medicine, these curves are especially helpful.

3.2. Information Measures

In this sub-section, we examine the entropies of Rényi, Havrda and Charvat, as well as d-generalized entropy as information metrics. These measures collectively provide information about the system’s overall amounts of data. The Rényi entropy presented by Rényi [24], is conceptually the quantity of information contained in a random process, it is defined by:
η ( d ) = ( 1 d ) 1 log f ( w ) d d w , d > 0 , d 1 .
Inserting (6) in (10), and using binomial expansion, then η ( d ) is as follows:
η ( d ) = 1 ( 1 d ) log i = 0 ( 1 ) i d ( δ 1 ) i ( 2 a 2 b δ ) d 0 1         w d ( ln w ) d ( 2 b 1 ) e ( d + i ) ( a ( ln w ) b ) 2 d w .
Let z = a ( ln w ) b 2 , and using exponential expansion in (11), we obtain
η ( d ) = 1 ( 1 d ) log i , m = 0 Ξ i , m ( d , δ , b , a ) Γ d ( 2 b 1 ) + 1 + m 2 b ,
where
Ξ i , m ( d , δ , b , a ) = ( 1 ) i + m d ( δ 1 ) i ( 2 b ) d 1 δ d ( d 1 ) m ( m ! ) ( d + i ) d ( 2 b 1 ) + 1 + m 2 b a d b 1 b m b
Reference [25] proposed another uncertainty measure, the Havrda and Charvat. Here we assume H ( d ) , and this is represented mathematically by:
H d = 1 2 1 d 1 f ( w ) d d w 1 d 1 , d 1 , d > 0 .
Using the same procedure above, we obtain H ( d ) , as follows:
H ( d ) = 1 2 1 d 1 i , m = 0 Ξ i , m ( d , δ , b , a ) Γ d ( 2 b 1 ) + 1 + m 2 b 1 d 1
In reference [26], a further generalized Shannon entropy form known as d-generalized entropy was developed. It is represented mathematically as below:
Κ ( d ) = ( d 1 ) 1 f ( w ) 2 d d w 1 , 0 < d < 2 , d 1 .
Using a similar way as above, we obtain the d-generalized entropy as follows,
Κ ( d ) = 1 ( d 1 ) i , m = 0 Ξ i , m * * ( d , δ , b , a ) Γ ( 2 d ) ( 2 b 1 ) + 1 + m 2 b ,
where
Ξ i , m * * ( d , δ , b , a ) = ( 1 ) i + m ( 2 d ) ( δ 1 ) i ( 2 b ) 1 d δ 2 d ( d + 1 ) m ( m ! ) ( 2 d + i ) ( 2 d ) ( 2 b 1 ) + 1 + m 2 b a 1 b d b m b
We use the following sets of parameters to provide entropy numerical values for the measurements under consideration: (i) ( a = 1.5 ,   b = 1.3 ,   δ = 1.4 ) , (ii) ( a = 0.5 ,   b = 0.4 ,   δ = 0.4 ) , (iii) ( a = 4 ,   b = 2 ,   δ = 0.7 ) , (iv) ( a = 0.7 ,   b = 0.4 ,   δ = 2 ) , (v) ( a = 5 ,   b = 1.6 ,   δ = 0.4 ) , (vi) ( a = 0.7 ,   b = 0.5 ,   δ = 2 ) , and (vii) ( a = 1.5 ,   b = 1.5 ,   δ = 3 ) .  Table 2 provides some numerical values for the provided three entropy measures.

3.3. Stochastic Ordering

The statistical literature places a great emphasis on the ordering of distributions, especially among lifetime distributions. A significant part of the ranking of various lifetime distributions is found in Johnson et al. [27]. Here, we take into account four distinct SO for two independent UPBX random variables with a restricted parameter space: the usual, the hazard rate, the mean residual life, and the likelihood ratio order. Recall that a family has the monotone likelihood ratio property if it has a likelihood ratio ordering. This suggests that, when the other parameters are known, there exists a test that is consistently the strongest for any one-sided hypothesis. According to Shaked and Shanthikumar [28], when two independent random variables, W1 and W2, have CDFs that are F W 1 ( w ) and F W 2 ( w ) , respectively, W1 is said to be smaller than W2 in the
Stochastic order (W1st (W2)) if F W 1 ( w ) F W 2 ( w ) w
Hazard rate order (W1hr (W2)) if h W 1 ( w ) h W 2 ( w ) w
Mean residual life order (W1mrl (W2)) if m W 1 ( w ) m W 2 ( w ) w
Likelihood ratio order (W1lr (W2)) if f W 1 ( w ) / f W 2 ( w ) decreases in w.
Assume that Wi, i = 1, 2 have the UPBXD with parameters ( a i , b i , δ i ) . Further, assume that F i ( w ) and f i ( w ) indicate, respectively, Wi’s CDF and PDF.
If f W 1 ( w ) / f W 2 ( w ) is a decreasing function ∀ w, then, in terms of likelihood ratio order; W1 is said to be stochastically less than W2 (W1lrW2)
Let W1~UPBXD ( a 1 , b 1 , δ 1 ) and W2~UPBXD ( a 2 , b 2 , δ 2 ) , then the likelihood ratio ordering is as follows:
f W 1 ( w ) f W 2 ( w ) = 2 a 1 2 b 1 δ 1 ( ln w ) 2 b 1 1 e a 1 ( ln w ) b 1 2 1 exp a 1 ( ln w ) b 1 2 δ 1 2 a 2 2 b 2 δ 2 ( ln w ) 2 b 2 1 e a 2 ( ln w ) b 2 2 1 exp a 2 ( ln w ) b 2 2 δ 2 , d d w log f W 1 ( w ) f W 2 ( w ) = 2 ( b 1 b 2 ) w ln w + 2 a 1 2 b 1 ( ln w ) 2 b 1 1 a 2 2 b 2 ( ln w ) 2 b 2 1 w 2 a 1 2 b 1 δ 1 exp a 1 ( ln w ) b 1 2 ( ln w ) 2 b 1 1 w 1 exp a 1 ( ln w ) b 1 2 δ 1 + 2 a 2 2 b 2 δ 2 ( ln w ) 2 b 2 1 exp a 2 ( ln w ) b 2 2 w 1 exp a 2 ( ln w ) b 2 2 δ 2 .
For a 1 > a 2 , b 1 > b 2 , δ 1 > δ 2 , we get d d w log f W 1 ( w ) f W 2 ( w ) < 0 , for all 0 w < 1 , hence f 1 ( w ) f 2 ( w ) is decreasing in w and hence W1 lr W2. Moreover, W1 is said to be smaller than W2 in other orderings such as SO (W1 stW2), HF(W1hrW2), and mean residual order (W1mrlW2).

3.4. Stress–Stress Reliability

In statistical literature, the term “SS reliability” is used to characterize the reliability of a system subjected to random stress W2 and having random strength W1, with the system failing if W2 is greater than W1, that is; R = P(W2 < W1). Let us assume that W1∼UPBXD ( a 1 , b , δ 1 ) and W2∼UPBXD ( a 2 , b , δ 2 ) are two independent random variables. The SS reliability of the UPBXD is then calculated as follows:
R = 2 a 1 b δ 1 0 1 w 1 ( ln w ) 2 b 1 e ( a 1 ( ln w ) b ) 2 1 e ( a 1 ( ln w ) b ) 2 δ 1 1 1 1 e ( a 2 ( ln w ) b ) 2 δ 2 d w .
Using the binomial expansions in (13), we get
R = 1 2 a 1 b δ 1 i 1 , i 2 = 0 ( 1 ) i 1 + i 2 δ 1 1 i 1 δ 2 i 2 0 1 w 1 ( ln w ) 2 b 1 e ( a 1 2 i 1 + i 2 a 2 2 ) ( ln w ) 2 b d w = 1 k = 0 ( 1 ) i 1 + i 2 a 1 δ 1 ( a 1 2 i 1 + i 2 a 2 2 ) δ 1 1 i 1 δ 2 i 2 .
As seen in (14) the SS reliability dependent on the parameters a 1 , a 2 , δ 1 , and δ 2 .

4. Parameter Estimation

The estimation methodologies for the parameters ( a , b , δ ) T of the UPBXD are obtained in this part using Bayesian and non-Bayesian estimation approaches. We provide classical method for the UPBXD as ML and Bayesian estimation utilizing various loss functions, including the squared error loss function (SELF), the linear exponential (LINEX) loss function and entropy loss function (ELF).

4.1. Maximum Likelihood Method

Consider a population that has a UPBXD described by PDF (6) with an unknown parameter vector ( a , b , δ ) , and that a random sample of size n is taken from that population. Following that, the likelihood of UPBXD for ( a , b , δ ) , say L w ¯ | , will be
L w ¯ | = 2 n a 2 n b n δ n e i = 1 n ( a ( ln w i ) b ) 2 i = 1 n ( ln w i ) 2 b 1 w i 1 exp a ( ln w i ) b 2 δ 1 .
The log likelihood function for ( a , b , δ ) , say l w ¯ | , will be
l w ¯ | = n ln ( 2 δ a 2 b ) i = 1 n ln w i + ( 2 b 1 ) i = 1 n ln ( ln w i ) a 2 i = 1 n ( ln w i ) 2 b + ( δ 1 ) i = 1 n ln 1 exp a ( ln w i ) b 2 .
The nonlinear equations created by differentiating (16) with respect to a , b , and δ , are solved to obtain the ML estimator for the unknown parameters. The score vector components, say U ( ) = l ( w ) = l ( w ) a , l ( w ) b , l ( w ) δ T , are given by
U ( a ) = 2 n a 2 a i = 1 n ( ln w i ) 2 b + ( δ 1 ) i = 1 n 2 a ( ln w i ) 2 b exp a ( ln w i ) b 2 1 ,
U ( b ) = n b + 2 i = 1 n ln ( ln w i ) 2 a 2 i = 1 n ( ln w i ) 2 b ln ln w i + i = 1 n 2 a 2 ( δ 1 ) ( ln w i ) 2 b ln ln w i exp a 2 ( ln w i ) 2 b 1 ,
U ( δ ) = n δ + i = 1 n ln 1 exp a ( ln w i ) b 2 .
The ML estimator of , say ^ , is achieved by solving the nonlinear system (17)–(19). These equations cannot be resolved analytically, but they can be resolved numerically by iterative statistical software techniques. We can use iterative methods, such as a Newton–Raphson algorithm, to obtain these estimates.

4.2. Bayesian Estimation

In this section, the Bayesian estimators based on different loss functions and associated highest posterior density (HPD) intervals of the UPBXD parameters are developed. The posterior distribution of is described in the following if we assume that the prior PDF of is unknown.
π ( | w ¯ ) L w ¯ | g ( ) .
The posterior density of is defined in Equation (20) as π ( | w ¯ ) , where on the right hand side L w ¯ | is the likelihood function of UPBXD ( ) and g ( ) is the prior density of .

4.2.1. Prior Information

For the purpose of discussing Bayesian estimate, we assume that the parameters a , b , and δ are independently distributed using the gamma distribution. Let q j and h j where j =1, 2, 3, be the scale and shape parameters for the gamma priors of a , b , and δ . The following is a proportionate representation of the joint density of a , b , and δ .
π a q 1 1 e h 1 a b q 2 1 e h 2 b δ q 2 1 e h 3 δ ; q j , h j > 0 , j = 1 , 2 , 3 .
The hyper-parameters will be elicited using the informative priors. When j = 1, …, L and k are the number of samples available from the UPBXD simulation, the mean and variance obtained using the ML estimates of the UPBXD a , b , and δ will be equal to the mean and variance of the considered priors (Gamma priors) q j and h j . By equating a , b , and δ with the mean and variance of gamma priors, we may determine their respective means and variances. Thus, we obtain
q 1 = 1 L j = 1 L a ^   j 2 1 L 1 j = 1 L a ^   j 1 L j = 1 L a ^   j 2 , & h 1 = 1 L j = 1 L a ^   j 1 L 1 j = 1 L a ^   j 1 L j = 1 L a ^   j 2 , q 2 = 1 L j = 1 L b ^   j 2 1 L 1 j = 1 L b ^   j 1 L j = 1 L b ^   j 2 , & h 2 = 1 L j = 1 L b ^   j 1 L 1 j = 1 L b ^   j 1 L j = 1 L b ^   j 2 , q 3 = 1 L j = 1 L δ ^   j 2 1 L 1 j = 1 L δ ^   j 1 L j = 1 L δ ^   j 2 , & h 3 = 1 L j = 1 L δ ^   j 1 L 1 j = 1 L δ ^   j 1 L j = 1 L δ ^   j 2 .
In regard to be solving the above two equations, the estimated hyper-parameters can be written as described in the following subsections.

4.2.2. Posterior Distribution

Here, the symmetric loss function (SELF), and asymmetric loss function (LINEX and ELF) are used to develop the Bayesian estimators for the same unknown parameters by utilizing independent gamma priors.
The likelihood function (15) and the joint prior function (21) are combined to form the joint posterior distribution. Hence, the joint posterior density function is
π ( | w ¯ ) a 2 n + q 1 1 e h 1 a b n + q 2 1 e h 2 b δ n + q 2 1 e h 3 δ e i = 1 n ( a ( ln w i ) b ) 2 i = 1 n ( ln w i ) 2 b 1 w i 1 exp a ( ln w i ) b 2 δ 1 .
The SELF, is defined as follows:
SELF ˜ , ˜ 2 .
The Bayesian estimator of under SELF is as follows:
˜ = 0 0 0 π | w _ d a   d b   d δ .
The LINEX, as asymmetric loss function, which is denoted by ˜ ˜ , is the derived as follows:
Linex ˜ ˜ , e c   c   1 ;         c 0 .
The Bayesian estimator of under LINEX loss function is as follows:
˜ ˜ = 0 0 0 e c   π | w _ d a   d b   d δ .
The ELF was first suggested by James and Stein [29] to estimate the Variance–Covariance (i.e., dispersion) matrix of the multivariate normal distribution. According to Calabria and Pulcini [30], the ELF is an excellent asymmetric loss function. The form’s ELF is thought of as
ELF , c 2 2 ln ln ;         c 0 .
The Bayesian estimator of under ELF is as follows:
= 0 0 0 c π | w _ d a   d b   d δ .
The Bayes estimator of a , b , and δ via different loss functions cannot be expressed in an explicit statement, as is evident from Equations (23)–(25). To do this, we suggest generating samples from conditional posterior distribution using Bayes Monte Carlo Markov chain (MCMC) techniques in order to compute the acquired Bayes estimates and create associated HPD intervals.

4.2.3. Markov Chain Monte Carlo

Since it is challenging to solve these integrals analytically, the MCMC method will be used. The most important sub-classes of MCMC algorithms are Gibbs sampling and the Metropolis–Hastings (MH) samplers. To do this, it regards a candidate value produced from a proposal distribution as normal for each iteration of the process, the MH method is comparable to acceptance–rejection sampling. From Equation (22), the full conditional density of a , b , and δ are provided, respectively, to execute the MCMC sampler as follows:
π ( a | b , δ , w ¯ ) a 2 n + q 1 1 e w 1 a e i = 1 n ( a ( ln w i ) b ) 2 i = 1 n 1 exp a ( ln w i ) b 2 δ 1 , π ( b | a , δ , w ¯ ) b n + q 2 1 e w 2 b e i = 1 n ( a ( ln w i ) b ) 2 i = 1 n ( ln w i ) 2 b 1 w i 1 exp a ( ln w i ) b 2 δ 1 ,
and
π ( δ | a , b , w ¯ ) δ n + q 2 1 e w 2 δ i = 1 n 1 exp a ( ln w i ) b 2 δ 1 .
It is thought that the MH algorithm can resolve this issue (for detail, see Alrumayh et al. [31] and Almetwally et al. [32]). The MH algorithm’s sampling procedure is carried out as follows:
Step 1: Set the initial values a ( 0 ) = a ^ , b ( 0 ) = b ^ , and δ ( 0 ) = δ ^ . Step 2: Set I = 1.
Step 3: Generate a , b and δ from N ( a ^ , V a ^ ) , N ( b ^ , V b ^ ) and N ( δ ^ , V δ ^ ) , respectively.
Step 4: Obtain a = min 1 , π a * | b I 1 , δ I 1 , w _ π a I 1 | b I 1 , δ I 1 , w _ , b = min 1 , π b * | a I 1 , δ I 1 , w _ π b I 1 | a I 1 , δ I 1 , w _ ,
δ = min 1 , π δ * | a I 1 , b I 1 , w _ π δ I 1 | a I 1 , b I 1 , w _ .
Step 5: Generate samples Uj j =1,2,3 from the uniform U(0, 1) distribution.
Step 6: If U 1 a , U 2 b , and U 3 δ , then set a ( I ) = a * , b ( I ) = b * , δ ( I ) = δ ; otherwise a ( I ) = a ( I 1 ) , b ( I ) = b ( I 1 ) , and δ I = δ I 1 .
Step 7: Set I = I+ 1.
Step 8: Repeat steps 3–7 B times and obtain a ( I ) , b ( I ) , and δ ( I ) , for I = 1, 2,..., B.

4.2.4. Highest Posterior Density Interval

Using the technique suggested by Chen and Shao [33], 100 1 α % HPD interval estimates of a , b , and δ are created. The MCMC samples of j for j = 1, …, B are first ordered. Therefore, the two-sided 100 1 α % HPD interval of is given by
a α 2 B , a 1 α 2 B ,   b α 2 B , b 1 α 2 B and δ α 2 B , δ 1 α 2 B ,
where a 1 < a 2 < < a B , b 1 < b 2 < < b B , and δ 1 < δ 2 < < δ B .

5. Simulation

A Monte Carlo simulation was run to evaluate the performance of the proposed point and interval estimators that were introduced in the previous sections. Based on various selections for sample size n as 40, 80, and 160, UPBXD was used to create a total of 5000 samples. To compare the results of Bayesian estimate based on various loss functions, the bias and mean squared errors (MSE) were calculated. The UPBXD was used to generate the data for the lifetime of various parameters a , b , and δ , as follows.
In Table 3: a = 0.5 ,   b = 0.6 and δ = 0.5 ,   1.2 and 3. In Table 4: a = 0.5 ,   b = 1.7 and δ = 0.5 ,   1.2 and 3. In Table 5: a = 2 ,   b = 0.6 and δ = 0.5 ,   1.2 and 3. In Table 6: a = 2 ,   b = 1.7 and δ = 0.5 ,   1.2 and 3.
The hybrid MCMC algorithm described in Section 4.2.3 was adopted to generate 12,000 MCMC samples, and we discarded the first 2000 values as ‘burn-in’. Accordingly, the 10,000 MCMC samples were used to produce the average Bayes MCMC estimates and 95% two-sided Bayesian credible intervals.
  • Algorithm for simulation: By establishing all simulation controls, we can build our model. The following actions must be finished in this stage in the correct order:
  • Assume different values for the UPBXD parameter vector and sample size.
  • Make the sample random values for the UPBXD using uniform and the QF in Equation (7).
  • We calculated the accuracy measures for each Bayes estimates of the UPBXD parameters using MH algorithm.
  • This experiment should be run (L-1) times.

5.1. Simulation Results

Table 3, Table 4, Table 5 and Table 6 show the results of the suggested techniques for calculating the point and interval parameter estimates. They offer the findings as well as some intriguing data. The following observations are permissible:
  • The estimates are asymptotically unbiased since they are more accurate as the sample size increases.
  • The parameter estimates come from the best unbiased estimator when the MSE value is near zero.
  • As the sample size grows, the MSE declines for each estimate, demonstrating consistency between the various estimates.
  • When the true value of δ increases, the bias, MSE, and length of the credible confidence interval (LCCI) of all estimates decrease.
  • The MSE and LCCI for the Bayesian estimates with positive weight for the asymmetric loss function are smaller than the Bayesian estimates with negative weight for asymmetric loss function.
  • The LCCI for estimates obtains its largest value, based on the suggested method, as the true values of the parameters increase.
  • An entropy loss function with positive weight is better than the other loss functions.

5.2. Represention Results

Figure 4, Figure 5, Figure 6 and Figure 7 show heatmap descriptions for the MSE results, where the bold color represents the highest values of MSE and the white color represents the lowest values of MSE.
The X-label belongs to SELFj, (j = 1, 2, 3) which are the MSE of Bayes estimates based on SELF with different parameters;
LINEXaj, (j =1, 2, 3) are the MSE of Bayes estimates based on LINEX c = 1.5 with different parameters;
LINEXbj, (j =1, 2, 3) are the MSE of Bayes estimates based on LINEX c = 1.5 with different parameters;
ELFaj, (j =1, 2, 3) are the MSE of Bayes estimates based on ELF c = 1.5 with different parameters;
ELFbj, (j =1, 2, 3) are the MSE of Bayes estimates based on ELF c = 1.5 with different parameters.
The Y-label belongs to cases and sample sizes, where C1n1 for δ = 0.5 and n = 40; C1n2 for δ = 0.5 and n = 80; C1n3 for δ = 0.5 and n = 80.

6. Application of Real Data

This section analyses two real-world datasets to show the adaptability and practical application of the UPBXD. The UPBXD is compared with the following models: unit-exponentiated half-logistic (UEHL) [23], Type II power Topp–Leone exponential (TIIPTLE) [34], Topp–Leone generalized exponential (TLGE) [35], Kumaraswamy (K), Beta, unit Weibull (UW), and Marshall–Olkin–Kumaraswamy (MOK). Two actual COVID-19 mortality rate datasets from Saudi Arabia and the United Kingdom are provided in this section to evaluate the UPBXD goodness of fit. The two real datasets were utilized to estimate the unknown parameters of the specified models using the maximum likelihood and Bayesian approaches. Kolmogorov–Smirnov statistics (KSS) with p-value, Cramer–von Mises statistics (WS), and Anderson–Darling statistics (AS) were used to compare all of the models.

6.1. Analysis for First Data

Data set I: The first set of data shows Saudi Arabia’s COVID-19 mortality rates over a 36-day period (22 July 2021 to 26 August 2021). The information is as follows: 0.1310, 0.1319, 0.1497, 0.1504, 0.1686, 0.1689, 0.1706, 0.1716, 0.1879, 0.1890, 0.1924, 0.1951, 0.2063, 0.2077, 0.2091, 0.2113, 0.2126, 0.2140, 0.2167, 0.2249, 0.2259, 0.2271, 0.2278, 0.2314, 0.2329, 0.2347, 0.2353, 0.2375, 0.2452, 0.2487, 0.2666, 0.2674, 0.2683, 0.2711, 0.2752, 0.2962. Table 7 shows the ML estimate of parameters with their standard errors (SEs) for each distribution and obtained the goodness of fit measures as KSS, WS, and AD. By the results shown in Table 7, we are able to see that the UPBXD is better than the other distributions, such as TLPTLE, TLGE, K, Beta, UW, UEHL, and MOK, for COVID-19 mortality rates in the Saudi Arabia data set.
As can be seen, the TLPTLE, TLGE, K, Beta, UW, UEHL, and MOK distributions work well for modelling the COVID-19 mortality rates indicated in the Saudi Arabia data set, but that the UPBXD is the best. This is based on a significance level of 0.05. Figure 8 illustrates the estimated CDF in the red line with empirical CDF in the black line. It also shows the probability–probability (PP) plots of the UPBXD in the red line, also known as “parametric plots”, for the COVID-19 mortality rates of the Saudi Arabia data set, which demonstrate the empirical findings, reported in Table 7 and the empirical CDF line the (black) with the estimated CDF line (red).
Figure 9 shows three plots of COVID-19 mortality rates for the Saudi Arabia data set, where the left is a boxplot of data that explains that the data have no outlier values, the center is a TTT plot of data that explains this data set is increasing, and the right is a hazard estimated plot line that indicates the HF is increasing.

6.2. Analysis for Second Data

Data set II: The second set of data shows the United Kingdom COVID-19 mortality rates over a 28-day period (1 January 2022 to 28 January 2022). The information is as follows: 0.1484, 0.1174, 0.0522, 0.0296, 0.0339, 0.2274, 0.1555, 0.1530, 0.2079, 0.0640, 0.1407, 0.2463, 0.2569, 0.2150, 0.1723, 0.1823, 0.1807, 0.1823, 0.2736, 0.2228, 0.2036, 0.1767, 0.1814, 0.1361, 0.1620, 0.2639, 0.2067, 0.2008.
Table 8 shows the ML estimate of parameters for each distribution and obtained the goodness of fit measures as KSS, WS, and AD. By the results of Table 8, we can see that the UPBXD is better than the other distributions, such as TLPTLE, TLGE, K, Beta, UW, UEHL, and MOK, for COVID-19 mortality rates in the United Kingdom data set. Additionally, we can see that the TLPTLE, TLGE, K, Beta, UW, UEHL, and MOK distributions work well for modelling the COVID-19 mortality rates of the United Kingdom data set, though the UPBXD is the best. This is based on a significance level of 0.05.
Figure 10 illustrates the PP plots for the COVID-19 mortality rates of the United Kingdom data set, which demonstrate the empirical findings reported in Table 8 and the empirical CDF line (black) with the estimated CDF line (red). Figure 11 shows three plots of COVID-19 mortality rates for the United Kingdom data set, where the left is a boxplot of data that explains that these data have no outlier values, the center is a TTT plot of data that explains that these data are increasing, and the right is a hazard estimated plot line that indicates the hazard is increasing.

6.3. Data Analysis via Bayesian Method

Here, we analyze data sets presented in previous sub-sections using the proposed Bayesian estimation method.
The Bayesian estimation parameters of UPBXD for each of the data sets, respectively are given in Table 9. The Bayesian estimates of UPBXD parameters under SELF and the corresponding SEs are calculated. The lower and upper HPD intervals are also calculated.
Figure 12 and Figure 13 display the trace plot of the UPBXD’s parameter values for the MCMC finding.
The autocorrelation function (ACF) is generated as shown in Figure 14 and Figure 15. Figure 16 and Figure 17 demonstrate the symmetric normal distribution of the posterior density for the parameters of the UPBXD.
Figure 18 and Figure 19 display the parameter convergence charts for UPBXD draws as well as the parameter random draw plot, respectively.

7. Conclusions

This article focuses on a three-parameter unit distribution created based on the power Burr X distribution and called the UPBXD. The statistical properties of the UPBXD have been derived and expressed in closed forms. The presented unit distribution can be used as a statistical tool to model different types of HFs, including those that are bathtub, increasing and unimodally shaped. Its important features are carefully studied, including the analytical expression of moments, quantile function, incomplete moments, stochastic ordering, and stress–strength reliability. Moreover, the uncertainty-measuring metrics Rényi, Havrda, and Charvat as well as d-generalized entropy were obtained. The UPBXD parameters have been estimated utilizing ML approach as well as Bayesian estimation approach with different loss functions. Additionally, Bayesian credible intervals were constructed based on the marginal posterior distribution. For some difficult calculations, the Markov chain Monte Carlo method was used. To assess how various estimates work, simulation studies based on various sample sizes have been carried out. In light of the simulation study’s findings, it was found that the Bayesian-based symmetric loss function and LINEX loss function techniques work quite effectively for estimating the UPBXD parameters. Bayesian estimates under an entropy loss function with positive weight are superior to those under other loss functions. The MSE and length of the credible confidence interval for Bayesian estimates with positive weight are smaller than the corresponding values with negative weight. Finally, two actual COVID-19 mortality rate data sets from Saudi Arabia and the United Kingdom have been analyzed and discussed to illustrate the notability of the UPBXD. The UPBXD gives superior fits over several other competing models, as shown by a real data application. Future discussions can be expanded on the use of Bayesian estimation in stress–strength reliability for the UPBXD based on some sampling techniques [36,37,38]. Furthermore, the proposed methodology can be expanded in multivariate and Bivariate case as [39,40,41].

Author Contributions

Conceptualization, A.S.H. and E.M.A.; Methodology, A.F., A.S.H., H.B. and E.M.A.; software, A.S.H. and E.M.A.; validation, A.S.H. and H.B.; formal analysis, A.F. and H.B.; investigation, A.F. and H.B.; data curation, A.S.H. and E.M.A.; writing—original draft, A.F., A.S.H., H.B. and E.M.A.; writing—review and editing, A.S.H., H.B. and E.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data sets have been provided in Section 6.

Acknowledgments

The authors would like to thank the editor and the anonymous referees for their efforts, wise observations, and constructive criticism, which greatly enhanced the manuscript’s contents.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burr, I.W. Cumulative frequency functions. Ann. Math. Stat. 1942, 13, 215–232. [Google Scholar] [CrossRef]
  2. Surles, J.; Padgett, W. Some properties of a scaled Burr type X distribution. J. Stat. Plan Inference 2005, 128, 271–280. [Google Scholar] [CrossRef]
  3. Raqab, M.Z.; Kundu, D. Burr type X distribution: Revisited. J. Probab. Stat. Sci. 2006, 4, 179–193. [Google Scholar]
  4. Merovci, F.; Khaleel, M.A.; Ibrahim, N.A.; Shitan, M. The beta Burr type X distribution properties with application. SpringerPlus 2016, 5, 697. [Google Scholar] [CrossRef] [Green Version]
  5. Yousof, H.M.; Afify, A.Z.; Hamedani, G.; Aryal, G. The Burr X generator of distributions for lifetime data. J. Stat. Theory Appl. 2017, 16, 288–305. [Google Scholar] [CrossRef] [Green Version]
  6. Ibrahim, N.A.; Khaleel, M.A.; Merovci, F.; Shitan, M. Weibull Burr Type X distribution properties with application. Pak. J. Stat. 2017, 33, 315–336. [Google Scholar]
  7. Nasiru, S.; Mwita, P.N.; Ngesa, O. Exponentiated generalized half logistic Burr X distribution. Adv. Appl. Statist. 2018, 52, 145–169. [Google Scholar] [CrossRef]
  8. Shrahili, M.; Elbatal, I.; Muhammad, M. The type I half-logistic Burr X distribution: Theory and practice. J. Nonlinear Sci. Appl. 2019, 12, 262–277. [Google Scholar] [CrossRef] [Green Version]
  9. Khan, M.S.; King, R.; Hudson, I.L. Transmuted Burr Type X distribution with covariates regression modeling to analyze reliability data. Am. J. Math. Manag. 2019, 39, 99–121. [Google Scholar] [CrossRef]
  10. Usman, R.M.; Ilyas, M. The power Burr Type X distribution: Properties, regression modeling and applications. Punjab Univ. J. Math. 2020, 52, 27–44. [Google Scholar]
  11. Topp, C.W.; Leone, F.C. A family of J-shaped frequency functions. J. Am. Stat. Assoc. 1955, 50, 209–219. [Google Scholar] [CrossRef]
  12. Kumaraswamy, P. A generalized probability density function for double-bounded random processes. J. Hydrol. 1980, 46, 79–88. [Google Scholar] [CrossRef]
  13. Gómez-Déniz, E.; Sordo, M.A.; Calderín-Ojeda, E. The log-Lindley distribution as an alternative to the beta regression model with applications in insurance. Insur. Math. Econ. 2014, 54, 49–57. [Google Scholar] [CrossRef]
  14. Mazucheli, J.; Menezes, A.F.B.; Dey, S. The unit-Birnbaum–Saunders distribution with applications. Chil. J. Stat. 2018, 9, 47–57. [Google Scholar]
  15. Ghitany, M.E.; Mazucheli, J.; Menezes, A.F.B.; Alqallaf, F. The unit-inverse Gaussian distribution: A new alternative to two-parameter distributions on the unit interval. Commun. Stat. Theory Methods 2019, 48, 3423–3438. [Google Scholar] [CrossRef]
  16. Mazucheli, J.; Menezes, A.F.B.; Chakraborty, S. On the one parameter unit-Lindley distribution and its associated regression model for proportion data. J. Appl. Stat. 2019, 46, 700–714. [Google Scholar] [CrossRef] [Green Version]
  17. Modi, K.; Gill, V. Unit Burr-III distribution with application. J. Stat. Manag. Syst. 2020, 23, 579–592. [Google Scholar] [CrossRef]
  18. Mazucheli, J.; Menezes, A.F.B.; Fernandes, L.B.; de Oliveira, R.P.; Ghitany, M.E. The unit-Weibull distribution as an alternative to the Kumaraswamy distribution for the modeling of quantiles conditional on covariates. J. Appl. Stat. 2020, 47, 954–974. [Google Scholar] [CrossRef]
  19. Korkmaz, M.C.; Chesneau, C. On the unit Burr-XII distribution with the quantile regression modeling and applications. Comput. Appl. Math. 2021, 40, 29. [Google Scholar] [CrossRef]
  20. Haq, M.A.U.; Albassam, M.; Aslam, M.; Hashmi, S. Statistical inferences on odd Fr´echet power function distribution. J. Reliab. Stat. Stud. 2021, 14, 141–172. [Google Scholar] [CrossRef]
  21. Krishna, A.; Maya, R.; Chesneau, C.; Irshad, M.R. The unit Teissier distribution and its applications. Math. Comput. Appl. 2022, 27, 12. [Google Scholar] [CrossRef]
  22. Jha, M.K.; Dey, S.; Alotaibi, R.; Alomani, G.; Tripath, Y.M. Multicomponent stress-strength reliability estimation based on unit generalized exponential distribution. Ain Shams Eng. J. 2022, 13, 101627. [Google Scholar] [CrossRef]
  23. Hassan, A.S.; Fayomi, A.; Algarni, A.; Almetwally, E.M. Bayesian and non-Bayesian inference for unit-exponentiated half-logistic distribution with data analysis. Appl. Sci. 2022, 12, 11253. [Google Scholar] [CrossRef]
  24. Rényi, A. On measures of entropy and information. In Proceedings of the Fourth Berkeley Symposium on Mathematical Statistics and Probability, Volume 1: Contributions to the Theory of Statistics; University of California Press: Berkeley, CA, USA, 1960; Volume 1, pp. 547–561. [Google Scholar]
  25. Havrda, J.; Charvát, F. Quantification method of classification processes. Concept Struct. -Entropy. Kybern. 1967, 3, 30–35. [Google Scholar]
  26. Mathai, A.M.; Haubold, H.J. On a generalized entropy measure leading to the pathway model with a preliminary application to solar neutrino data. Entropy 2013, 15, 4011–4025. [Google Scholar] [CrossRef] [Green Version]
  27. Johnson, N.L.; Kotz, S.; Balakrishnan, N. Continuous Univariate Distributions; Wiley: New York, NY, USA, 1995; Volume 2. [Google Scholar]
  28. Shaked, M.; Shanthikumar, J.G. Stochastic Orders; Wiley: New York, NY, USA, 2007. [Google Scholar]
  29. James, W.; Stein, C. Estimation with quadratic loss. In Breakthroughs in Statistics; Springer: New York, NY, USA, 1992; pp. 443–460. [Google Scholar]
  30. Calabria, R.; Pulcini, G. An engineering approach to Bayes estimation for the Weibull distribution. Microelectron. Reliab. 1994, 34, 789–802. [Google Scholar] [CrossRef]
  31. Alrumayh, A.; Weera, W.; Khogeer, H.A.; Almetwally, E.M. Optimal analysis of adaptive type-II progressive censored for new unit-lindley model. J. King Saud Univ. Sci. 2023, 35, 102462. [Google Scholar] [CrossRef]
  32. Almetwally, E.M.; Jawa, T.M.; Sayed-Ahmed, N.; Park, C.; Zakarya, M.; Dey, S. Analysis of unit-Weibull based on progressive type-II censored with optimal scheme. Alex. Eng. J. 2023, 63, 321–338. [Google Scholar] [CrossRef]
  33. Chen, M.H.; Shao, Q.M. Monte Carlo estimation of Bayesian credible and HPD intervals. J. Comput. Graph. Stat. 1999, 8, 69–92. [Google Scholar]
  34. Bantan, R.A.; Jamal, F.; Chesneau, C.; Elgarhy, M. Type II Power Topp-Leone generated family of distributions with statistical inference and applications. Symmetry 2020, 12, 75. [Google Scholar] [CrossRef] [Green Version]
  35. Sangsanit, Y.; Bodhisuwan, W. The Topp-Leone generator of distributions: Properties and inferences. Songklanakarin J. Sci. Technol. 2016, 38, 537–548. [Google Scholar]
  36. Yousef, M.M.; Hassan, A.S.; Alshanbari, H.M.; El-Bagoury, A.-A.H.; Almetwally, E.M. Bayesian and non-Bayesian analysis of exponentiated exponential stress–strength model based on generalized progressive hybrid censoring process. Axioms 2022, 11, 455. [Google Scholar] [CrossRef]
  37. Hassan, A.S.; Nagy, H.F. Reliability estimation in multicomponent stress strength for generalized inverted exponential distribution based on ranked set sampling. Gazi Univ. J. Sci. 2022, 35, 314–331. [Google Scholar]
  38. Hassan, A.S.; Almanjahie, I.M.; Al-Omari, A.I.; Alzoubi, L.; Nagy, H.F. Stress–strength modeling using median-ranked set sampling: Estimation, simulation, and application. Mathematics 2023, 11, 318. [Google Scholar] [CrossRef]
  39. Chesneau, C. On new three- and two-dimensional ratio-power copulas. Comput. J. Math. Stat. Sci. 2023, 2, 106–122. [Google Scholar] [CrossRef]
  40. El-Sherpieny, E.S.A.; Muhammed, H.Z.; Almetwally, E.M. Bivariate Chen distribution based on copula function: Properties and application of diabetic nephropathy. J. Stat. Theory Pract. 2022, 16, 54. [Google Scholar] [CrossRef]
  41. El-Sherpieny, E.S.A.; Almetwally, E.M.; Muhammed, H.Z. Bivariate Weibull-G family based on copula function: Properties, Bayesian and non-Bayesian estimation and applications. Stat. Optim. Inf. Comput. 2022, 10, 678–709. [Google Scholar] [CrossRef]
Figure 1. Plots of various PDF shapes of the UPBXD for different parameter values.
Figure 1. Plots of various PDF shapes of the UPBXD for different parameter values.
Axioms 12 00297 g001
Figure 2. Plots of various HF shapes of the UPBXD for different parameter values.
Figure 2. Plots of various HF shapes of the UPBXD for different parameter values.
Axioms 12 00297 g002
Figure 3. Coefficient of skewness and kurtosis for UPBXD.
Figure 3. Coefficient of skewness and kurtosis for UPBXD.
Axioms 12 00297 g003
Figure 4. Heatmap for MSE when a = 0.5 ,   b = 0.6 .
Figure 4. Heatmap for MSE when a = 0.5 ,   b = 0.6 .
Axioms 12 00297 g004
Figure 5. Heatmap for MSE when a = 0.5 ,   b = 1.7 .
Figure 5. Heatmap for MSE when a = 0.5 ,   b = 1.7 .
Axioms 12 00297 g005
Figure 6. Heatmap for MSE when a = 2 ,   b = 0.6 .
Figure 6. Heatmap for MSE when a = 2 ,   b = 0.6 .
Axioms 12 00297 g006
Figure 7. Heatmap for MSE when a = 2 ,   b = 1.7 .
Figure 7. Heatmap for MSE when a = 2 ,   b = 1.7 .
Axioms 12 00297 g007
Figure 8. The CDF plot with empirical line and PP plot for Saudi Arabia data set.
Figure 8. The CDF plot with empirical line and PP plot for Saudi Arabia data set.
Axioms 12 00297 g008
Figure 9. Boxplot, TTT plots and hazard line of UPBXD plot for Saudi Arabia data set.
Figure 9. Boxplot, TTT plots and hazard line of UPBXD plot for Saudi Arabia data set.
Axioms 12 00297 g009
Figure 10. The CDF plot with empirical line and PP plot for the United Kingdom data set.
Figure 10. The CDF plot with empirical line and PP plot for the United Kingdom data set.
Axioms 12 00297 g010
Figure 11. Boxplot, TTT plots and hazard line of UPBXD plot for the United Kingdom data set.
Figure 11. Boxplot, TTT plots and hazard line of UPBXD plot for the United Kingdom data set.
Axioms 12 00297 g011
Figure 12. Trace plots of MCMC results with interval limit line for Saudi Arabia data set.
Figure 12. Trace plots of MCMC results with interval limit line for Saudi Arabia data set.
Axioms 12 00297 g012
Figure 13. Trace plots of MCMC results with interval limit line for the United Kingdom data set.
Figure 13. Trace plots of MCMC results with interval limit line for the United Kingdom data set.
Axioms 12 00297 g013
Figure 14. The ACF plot of MCMC results for Saudi Arabia data set.
Figure 14. The ACF plot of MCMC results for Saudi Arabia data set.
Axioms 12 00297 g014
Figure 15. ACF plot of MCMC results for the United Kingdom data set.
Figure 15. ACF plot of MCMC results for the United Kingdom data set.
Axioms 12 00297 g015
Figure 16. Histogram plots of MCMC results for Saudi Arabia data set.
Figure 16. Histogram plots of MCMC results for Saudi Arabia data set.
Axioms 12 00297 g016
Figure 17. Histogram plots of MCMC results for the United Kingdom data set.
Figure 17. Histogram plots of MCMC results for the United Kingdom data set.
Axioms 12 00297 g017
Figure 18. Convergence lines of MCMC results for Saudi Arabia data set.
Figure 18. Convergence lines of MCMC results for Saudi Arabia data set.
Axioms 12 00297 g018
Figure 19. Convergence lines of MCMC results for the United Kingdom data set.
Figure 19. Convergence lines of MCMC results for the United Kingdom data set.
Axioms 12 00297 g019
Table 1. Several UPBXD moment values.
Table 1. Several UPBXD moment values.
μ m (i)(ii)(iii)(iv)(v)(vi)(vii)
μ 1 0.4990.3960.670.8020.8030.4770.416
μ 2 0.2640.2310.4580.6530.6880.2390.18
μ 3 0.1480.1600.320.540.6120.1250.08
μ 4 0.0870.1220.2280.4540.5560.0690.037
σ 2 0.0150.0740.0090.0110.0440.0110.0062
α 3 0.3770.5100.38−0.135−1.2370.4650.375
α 4 2.8562.1322.752.2483.7993.0703.135
Table 2. Numerical values for the UPBXD’s entropy measures.
Table 2. Numerical values for the UPBXD’s entropy measures.
dMeasures(i)(ii)(iii)(iv)(v)(vi)(vii)
0.5 η ( d ) −0.619−0.308−0.819−0.829−0.499−0.669−0.944
H ( d ) −1.115−0.64−1.35−1.36−0.948−1.177−1.475
K ( d ) −1.1010.792−1.354−1.201−2.953−1.179−1.688
1.5 η ( d ) −0.8771.01−1.034−0.94−1.814−0.927−1.224
H ( d ) −1.160.976−1.405−1.257−2.835−1.236−1.72
K ( d ) −0.533−0.286−0.672−0.678−0.442−0.568−0.753
Table 3. Bayesian inference with different loss functions when a = 0.5 ,   b = 0.6 .
Table 3. Bayesian inference with different loss functions when a = 0.5 ,   b = 0.6 .
a = 0.5 ,
b = 0.6
SELF LINEX   c = 1.5 LINEX   c = 1.5 ELF   c = 1.5 ELF   c = 1.5
δ n BiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCI
0.540 a 0.13920.02240.20750.14480.02420.21400.13350.02070.20210.14220.02330.20980.12370.01810.1985
b −0.02610.00080.0423−0.02530.00080.0420−0.02690.00070.0427−0.02560.00080.0421−0.02840.00070.0426
δ 0.61230.40950.70790.65870.47930.79900.55920.33640.59540.62670.42980.72750.52220.29320.5536
80 a 0.09000.00940.13900.09270.00990.14130.08720.00880.13650.09150.00970.14000.08180.00790.1331
b −0.01360.00020.0262−0.01330.00020.0262−0.01390.00020.0265−0.01340.00020.0263−0.01440.00030.0267
δ 0.42290.18780.35620.44790.21170.39960.39540.16320.31340.43230.19650.37210.36720.14050.2908
160 a 0.08430.00780.10170.08650.00820.10320.08210.00740.09940.08560.00800.10260.07790.00670.0967
b −0.01350.00020.0237−0.01260.00020.0237−0.01270.00020.0238−0.01360.00020.0237−0.01370.00020.0241
δ 0.39060.15670.24960.41210.17510.27320.36700.13800.21860.39910.16380.25640.34110.11900.2017
1.240 a 0.08040.00820.16040.08410.00880.16140.07670.00750.15850.08250.00850.16090.06950.00650.1600
b −0.02040.00080.0769−0.01870.00080.0775−0.02210.00090.0757−0.01950.00080.0772−0.02520.00100.0756
δ 0.69590.53210.80340.79360.69990.98370.59050.37790.65590.71380.56040.82490.59860.39130.7033
80 a 0.04170.00250.11300.04310.00270.11350.04030.00240.11120.04260.00260.11310.03720.00220.1105
b −0.00760.00020.0494−0.00690.00020.0498−0.00820.00020.0491−0.00720.00020.0496−0.00940.00030.0486
δ 0.37240.15150.43170.40440.17940.48840.33940.12520.38130.37920.15710.43870.33700.12380.3846
160 a 0.03710.00180.07800.03800.00180.07890.03630.00170.07740.03770.00180.07850.03440.00160.0769
b −0.00820.00020.0370−0.00790.00020.0371−0.00860.00020.0370−0.00800.00020.0371−0.00920.00020.0371
δ 0.34290.12420.31010.36890.14430.34530.31610.10500.27320.34860.12840.31680.31350.10340.2735
340 a 0.02240.00190.14500.02490.00200.14570.01980.00180.14380.02400.00200.14480.01430.00160.1431
b −0.00360.00110.1197−0.00130.00110.1215−0.00590.00110.1188−0.00230.00110.1199−0.01010.00110.1184
δ 0.42440.24851.01590.48970.32801.14890.35760.17930.85980.43070.25561.03060.39190.21420.9512
80 a 0.00910.00080.10130.01010.00080.10210.00800.00070.10110.00980.00080.10190.00570.00070.1005
b −0.00100.00050.09040.00000.00050.0907−0.00190.00050.0898−0.00040.00050.0906−0.00370.00050.0899
δ 0.19500.05670.52790.21130.06610.56900.17860.04810.49050.19670.05770.53220.18660.05230.5095
160 a 0.00860.00050.07670.00910.00050.07660.00800.00040.07570.00890.00050.07660.00670.00040.0753
b −0.00160.00040.0789−0.00100.00040.0791−0.00220.00040.0783−0.00130.00040.0791−0.00320.00040.0791
δ 0.18040.04090.34970.19110.04590.37340.16960.03610.32940.18150.04140.35190.17480.03840.3404
Table 4. Bayesian inference with different loss functions when a = 0.5 ,   b = 1.7 .
Table 4. Bayesian inference with different loss functions when a = 0.5 ,   b = 1.7 .
a = 0.5 ,
  b = 1.7
SELF LINEX   c = 1.5 LINEX   c = 1.5 ELF   c = 1.5 ELF   c = 1.5
δ n BiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCI
0.540 a 0.12120.01790.21550.12570.01910.21990.11670.01670.21100.12360.01850.21720.10860.01470.2073
b −0.03080.00130.0722−0.02810.00120.0716−0.03360.00150.0742−0.03030.00130.0722−0.03360.00150.0744
δ 0.59440.39010.72860.63740.45340.81820.54570.32380.62400.60780.40860.74700.51130.28440.5813
80 a 0.08450.00840.13820.08700.00890.14020.08200.00790.13650.08590.00870.13860.07710.00710.1335
b −0.02290.00090.0662−0.02150.00080.0663−0.02420.00100.0661−0.02260.00090.0662−0.02420.00100.0661
δ 0.42670.19110.37140.45110.21470.40760.39950.16660.32430.43590.19960.38100.37170.14390.2986
160 a 0.07540.00620.08970.07710.00650.09100.07360.00590.08890.07640.00640.09040.07010.00540.0877
b −0.01960.00050.0476−0.01890.00050.0477−0.02030.00060.0476−0.01950.00050.0477−0.02030.00060.0476
δ 0.38900.15500.23390.40990.17250.25730.36600.13680.21170.39730.16180.24140.34040.11820.1901
1.240 a 0.07220.00690.16230.07510.00740.16450.06930.00650.16000.07390.00720.16310.06350.00570.1576
b −0.03110.00330.1871−0.02410.00300.1885−0.03800.00370.1841−0.02970.00320.1865−0.03800.00370.1853
δ 0.73330.59230.88270.83770.77921.06650.61870.41660.70140.75230.62360.91190.62880.43420.7495
80 a 0.04230.00260.11020.04360.00270.11110.04100.00250.10940.04310.00260.11080.03830.00220.1089
b −0.01790.00200.1693−0.01520.00200.1689−0.02070.00210.1678−0.01740.00200.1694−0.02070.00220.1681
δ 0.39620.17090.46260.43100.20320.52010.35990.14030.40130.40360.17740.47220.35770.13900.4080
160 a 0.03650.00180.08390.03730.00190.08470.03570.00170.08310.03700.00180.08450.03400.00160.0822
b −0.01630.00170.1472−0.01460.00160.1465−0.01800.00180.1486−0.01600.00170.1468−0.01800.00180.1488
δ 0.34760.12710.30370.37380.14750.33670.32050.10760.27070.35330.13140.31020.31790.10600.2717
340 a 0.02020.00180.14840.02220.00190.15030.01820.00170.14830.02150.00190.14930.01380.00160.1477
b −0.00370.01000.39230.00870.01040.3968−0.01610.01000.3886−0.00130.01000.3898−0.01600.01020.3926
δ 0.45030.28041.07270.52440.37531.22600.37480.19810.92090.45750.28881.09000.41380.23981.0071
80 a 0.01060.00080.10230.01130.00080.10190.00980.00080.10180.01100.00080.10190.00810.00080.1018
b −0.00560.00520.2814−0.00150.00520.2815−0.00960.00530.2807−0.00480.00520.2808−0.00960.00530.2821
δ 0.20740.06180.52020.22510.07220.55200.18970.05220.48590.20920.06280.52250.19830.05680.5058
160 a 0.00820.00040.07080.00850.00040.07090.00780.00040.07060.00840.00040.07090.00690.00040.0702
b −0.00370.00250.1941−0.00180.00250.1921−0.00560.00250.1949−0.00330.00250.1941−0.00560.00250.1953
δ 0.18950.04570.39200.20140.05170.42200.17750.04000.36620.19080.04630.39520.18330.04280.3771
Table 5. Bayesian inference with different loss functions when a = 2 ,   b = 0.6 .
Table 5. Bayesian inference with different loss functions when a = 2 ,   b = 0.6 .
a = 2 ,
  b = 0.6
SELF LINEX   c = 1.5 LINEX   c = 1.5 ELF   c = 1.5 ELF   c = 1.5
δ n BiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCI
0.540 a 0.14420.03300.41350.16070.03930.42500.12750.02730.39410.14680.03390.41810.13130.02860.4004
b −0.03880.00210.0635−0.03800.00200.0621−0.03960.00210.0652−0.03820.00200.0625−0.04110.00230.0681
δ 0.56790.35970.73170.60740.41690.81790.52320.29980.60750.58040.37650.75450.49070.26340.5613
80 a 0.06780.00960.27570.07350.01070.28520.06220.00870.27210.06870.00980.27640.06340.00890.2724
b −0.02750.00090.0399−0.02710.00080.0395−0.02790.00090.0402−0.02720.00090.0397−0.02860.00090.0410
δ 0.39920.16760.34630.41990.18610.37520.37640.14820.31210.40730.17450.35590.35180.12930.2869
160 a 0.05820.00620.20800.06140.00670.21160.05500.00580.20650.05870.00630.20840.05560.00580.2069
b −0.02820.00080.0342−0.02790.00090.0342−0.02850.00100.0341−0.02800.00090.0334−0.02910.00100.0332
δ 0.37020.14100.23690.38860.15570.26400.35010.12570.20960.37760.14680.24730.32700.10940.1937
1.240 a 0.09950.01780.34230.11260.02100.34830.08660.01500.33730.10160.01820.34290.08920.01560.3408
b −0.03730.01040.0748−0.03530.01320.0747−0.04000.00510.0744−0.03630.01090.0746−0.04290.00620.0747
δ 0.72560.58340.91880.83050.77241.13040.61210.41040.72130.74460.61460.94360.62180.42770.7719
80 a 0.03950.00570.25320.04380.00620.25770.03520.00530.25090.04020.00580.25420.03600.00540.2511
b −0.01980.00060.0506−0.01910.00050.0511−0.02050.00060.0497−0.01940.00050.0509−0.02180.00060.0492
δ 0.38380.16120.44500.41720.19130.49720.34930.13270.39170.39090.16720.45650.34700.13130.3945
160 a 0.03520.00350.18020.03750.00370.18160.03290.00330.17800.03560.00360.18050.03340.00340.1788
b −0.01940.00060.0366−0.01910.00070.0365−0.01980.00080.0366−0.01920.00070.0365−0.02060.00090.0368
δ 0.34370.12440.29400.36930.14410.33340.31720.10550.26350.34930.12850.30260.31460.10390.2647
340 a 0.03820.00610.26460.04610.00690.26660.03040.00540.25820.03950.00620.26400.03180.00550.2594
b −0.00380.01670.0996−0.00020.03180.1001−0.00810.00360.0979−0.00250.01870.0993−0.01130.00510.0978
δ 0.46440.29521.08150.53880.39411.20130.38830.20920.91740.47160.30391.08660.42790.25301.0157
80 a 0.01030.00310.21400.01340.00320.21450.00710.00300.21100.01080.00310.21430.00770.00300.2116
b −0.00580.00050.0801−0.00470.00050.0800−0.00680.00050.0796−0.00520.00050.0802−0.00860.00050.0795
δ 0.21690.06750.55250.23620.07970.59740.19740.05640.51000.21880.06870.55770.20690.06180.5310
160 a 0.01150.00140.13560.01290.00140.13650.01000.00130.13520.01170.00140.13570.01020.00130.1353
b −0.00510.00020.0534−0.00460.00020.0532−0.00550.00020.0531−0.00480.00020.0533−0.00630.00020.0529
δ 0.18940.04610.38620.20140.05210.40960.17730.04040.36000.19070.04670.38960.18320.04310.3755
Table 6. Bayesian inference with different loss functions when a = 2 ,   b = 1.7 .
Table 6. Bayesian inference with different loss functions when a = 2 ,   b = 1.7 .
a = 2 ,
  b = 1.7
SELF LINEX   c = 1.5 LINEX   c = 1.5 ELF   c = 1.5 ELF   c = 1.5
δ n BiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCIBiasMSELCCI
0.540 a 0.12550.02320.32280.13880.02740.33370.11220.01940.31500.12760.02380.32370.11510.02030.3166
b −0.05380.00360.0995−0.05070.00320.0975−0.05680.00400.1018−0.05320.00350.0990−0.05690.00400.1021
δ 0.53810.31720.60430.57100.36030.67710.50120.27190.53600.54890.33050.62530.47230.24130.5051
80 a 0.07150.00940.25490.07740.01050.26210.06570.00840.24830.07240.00960.25670.06680.00860.2506
b −0.05050.00320.0934−0.04850.00300.0927−0.05250.00340.0944−0.05010.00320.0932−0.05260.00340.0946
δ 0.40030.16810.33320.42110.18680.36690.37750.14870.30080.40830.17510.34390.35310.12970.2722
160 a 0.05800.00540.17840.06080.00580.18050.05510.00500.17630.05840.00550.17840.05570.00510.1770
b −0.03640.00160.0611−0.03550.00150.0604−0.03720.00170.0622−0.03620.00160.0610−0.03720.00170.0623
δ 0.36240.13460.22580.37900.14750.24700.34420.12110.20640.36920.13980.23380.32250.10620.1864
1.240 a 0.10310.01800.31870.11560.02130.32850.09060.01500.31320.10500.01840.32060.09320.01570.3156
b −0.07800.00890.2048−0.06960.00770.2064−0.08630.01030.2051−0.07630.00860.2049−0.08660.01030.2065
δ 0.73170.59270.90070.83580.78211.07990.61800.41710.70110.75050.62410.93420.62820.43490.7433
80 a 0.04230.00530.22950.04640.00580.23220.03810.00480.22860.04290.00530.23030.03890.00490.2297
b −0.03470.00260.1457−0.03190.00240.1454−0.03760.00280.1452−0.03410.00260.1456−0.03760.00280.1455
δ 0.39350.16760.44090.42830.19950.48710.35730.13760.38510.40090.17410.45270.35490.13620.3910
160 a 0.03390.00320.17650.03610.00330.17960.03170.00300.17450.03420.00320.17680.03210.00300.1751
b −0.03330.00210.1238−0.03160.00200.1235−0.03500.00230.1249−0.03300.00210.1237−0.03500.00230.1250
δ 0.34080.12240.30730.36610.14170.34020.31470.10390.26880.34640.12650.31490.31210.10230.2702
340 a 0.03650.00580.24450.04400.00650.24800.02910.00520.24130.03770.00590.24440.03050.00530.2415
b −0.01650.00670.2990−0.00580.00670.3036−0.02710.00690.2941−0.01440.00660.2979−0.02710.00700.2965
δ 0.48060.32341.16720.56380.43891.32600.39580.22411.01170.48870.33341.18390.43990.27461.0958
80 a 0.01580.00250.17990.01880.00270.18250.01280.00240.17860.01630.00260.18010.01330.00240.1791
b −0.00970.00590.2519−0.00530.00530.2507−0.01400.00650.2519−0.00870.00560.2512−0.01420.00730.2523
δ 0.22420.07150.55630.24390.08420.59150.20430.05980.51350.22620.07270.55980.21410.06550.5382
160 a 0.01100.00120.12510.01240.00120.12550.00960.00120.12450.01120.00120.12510.00980.00120.1247
b −0.01070.00190.1672−0.00890.00180.1683−0.01240.00190.1684−0.01030.00190.1673−0.01240.00190.1686
δ 0.19700.04750.35680.20920.05370.38040.18470.04170.33250.19820.04810.36030.19060.04450.3454
Table 7. ML estimates with SE and goodness of fit statistics: Saudi Arabia data set.
Table 7. ML estimates with SE and goodness of fit statistics: Saudi Arabia data set.
ab δ KSSp-ValueWSAS
UPBXDEstimates2.27170.56092813.28860.07780.96930.03270.2364
SE0.07080.0692469.2354
TLPTLEEstimates693.17740.64710.64760.09380.88000.04790.3110
SE1626.84450.09390.8346
TLGEEstimates0.368220.4075179.70440.14030.43780.09650.5911
SE0.28432.7038178.4010
KEstimates3.3085125.2161 0.18210.16210.04210.2793
SE0.282149.4480
BetaEstimates20.817476.5218 0.11270.70890.06360.3992
SE4.869018.0555
UWEstimates0.02037.7557 0.16330.26240.18241.0950
SE0.01100.9132
UEHLEstimates6.06553670.3422 0.07920.96410.03300.2393
SE0.7918405.8862
MOKEstimates703.31301.309745.44760.08110.95670.03360.2554
SE4615.48971.571272.4483
Table 8. Estimates with SE and goodness of fit statistics of ML: The United Kingdom data set.
Table 8. Estimates with SE and goodness of fit statistics of ML: The United Kingdom data set.
ab δ KSSp-ValueWSAS
UPBXDEstimates2.34710.33542154.87420.11000.85120.11150.7238
SE0.07160.0478402.9847
TLPTLEEstimates3760.83720.11270.02160.25710.04040.46562.6059
SE5117.25710.02870.0151
TLGEEstimates0.296810.389613.89630.20970.14710.28041.6718
SE0.30762.341015.4726
KEstimates2.9163125.0007 0.13290.65700.14920.9456
SE0.468994.0394
BetaEstimates3.927719.1899 0.19250.22020.26091.5673
SE1.00905.1885
UWEstimates0.09043.2548 0.24680.05480.46222.5905
SE0.03860.4231
UEHLEstimates2.978969.5723 0.12940.68880.14760.9369
SE0.480853.4192
MOKEstimates0.01243.54836.70140.15060.50200.31171.8656
SE0.02990.602614.9920
Table 9. Bayesian estimation based on SELF for parameters of UPBXD.
Table 9. Bayesian estimation based on SELF for parameters of UPBXD.
Data EstimatesSELowerUpper
Saudi Arabia a 2.32350.06191.95412.7490
b 0.55560.05580.39900.7320
δ 2979.20332.49312974.21662984.1337
The United Kingdom a 2.38520.05602.01082.8278
b 0.33530.03140.20640.4830
δ 2154.87430.07872154.71692155.0299
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fayomi, A.; Hassan, A.S.; Baaqeel, H.; Almetwally, E.M. Bayesian Inference and Data Analysis of the Unit–Power Burr X Distribution. Axioms 2023, 12, 297. https://doi.org/10.3390/axioms12030297

AMA Style

Fayomi A, Hassan AS, Baaqeel H, Almetwally EM. Bayesian Inference and Data Analysis of the Unit–Power Burr X Distribution. Axioms. 2023; 12(3):297. https://doi.org/10.3390/axioms12030297

Chicago/Turabian Style

Fayomi, Aisha, Amal S. Hassan, Hanan Baaqeel, and Ehab M. Almetwally. 2023. "Bayesian Inference and Data Analysis of the Unit–Power Burr X Distribution" Axioms 12, no. 3: 297. https://doi.org/10.3390/axioms12030297

APA Style

Fayomi, A., Hassan, A. S., Baaqeel, H., & Almetwally, E. M. (2023). Bayesian Inference and Data Analysis of the Unit–Power Burr X Distribution. Axioms, 12(3), 297. https://doi.org/10.3390/axioms12030297

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