Next Article in Journal
Pattern Dynamics of Cross Diffusion Predator–Prey System with Strong Allee Effect and Hunting Cooperation
Next Article in Special Issue
On Some Improved Class of Estimators by Using Stratified Ranked Set Sampling
Previous Article in Journal
Multi-Objective Instantaneous Center of Rotation Optimization Using Sensors Feedback for Navigation in Self-Reconfigurable Pavement Sweeping Robot
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient Estimation of Two-Parameter Xgamma Distribution Parameters Using Ranked Set Sampling Design

by
Amer Ibrahim Al-Omari
1,*,
SidAhmed Benchiha
2 and
Ibrahim M. Almanjahie
3,4
1
Department of Mathematics, Faculty of Science, Al al-Bayt University, Mafraq 25113, Jordan
2
Department of Mathematics, University of Djillali Liabes, BP 89, Sidi Bel Abbes 22000, Algeria
3
Department of Mathematics, College of Science, King Khalid University, Abha 62529, Saudi Arabia
4
Statistical Research and Studies Support Unit, King Khalid University, Abha 62529, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(17), 3170; https://doi.org/10.3390/math10173170
Submission received: 12 August 2022 / Revised: 25 August 2022 / Accepted: 30 August 2022 / Published: 2 September 2022

Abstract

:
An efficient method such as ranked set sampling is used for estimating the population parameters when the actual observation measurement is expensive and complicated. In this paper, we consider the problem of estimating the two-parameter xgamma (TPXG) distribution parameters under the ranked set sampling as well as the simple random sampling design. Various estimation methods, including the weighted least-square estimator, maximum likelihood estimators, least-square estimator, Cramer–von Mises, the maximum product of spacings estimators, and Anderson–Darling estimators, are considered. A comparison between the ranked set sampling and simple random sampling estimators, with the same number of measurement units, is conducted using a simulation study in terms of the bias, mean squared errors, and efficiency of estimators. The merit of the ranked set sampling estimators is examined using real data of bank customers. The results indicate that estimations using the ranked set sampling method are more efficient than the simple random sampling competitor considered in this study.

1. Introduction

A random variable X follows the xgamma distribution if its probability density function (pdf) is given by
f ( x , λ ) = λ 2 λ x 2 + 2 2 λ + 1 e λ x ; x > 0 , λ > 0 ,
and its cumulative distribution function (CDF) is
F ( x , λ ) = 1 x 2 λ 2 2 + x λ + λ + 1 λ + 1 e x λ ; x > 0 , λ > 0 .
Plots of the pdf of the xgamma distribution are presented in Figure 1 for some values of λ .
As an extension to the xgamma distribution, the two-parameter xgamma distribution is proposed as a new distribution by [1] by using an additional parameter to the xgamma distribution to obtain a more flexible distribution in modeling real data sets due to the wide use of the xgamma distribution in several survival analyses. When a random variable, X, follows the TPXG distribution, the probability density function and the cumulative distribution function are, respectively, given by
f ( x , β , λ ) = λ 2 β λ x 2 + 2 e λ x 2 λ + β ; x > 0 , λ > 0 , β > 0 ,
F ( x , β , λ ) = 1 β x 2 λ 2 2 + β x λ + λ + β e x λ λ + β ; x > 0 , λ > 0 , β > 0 .
For β = 1 in (3), we obtain the xgamma distribution with parameter λ as a special case of the TPXG distribution.
The rth moment for the distribution is obtained by
E ( X r ) = r ! 2 λ r ( β + λ ) [ 2 λ + β ( 1 + r ) ( 2 + r ) ] ; r = 1 , 2
The characteristic function (CF) and hazard function H ( x , β , λ ) of the model are, respectively, given by
ϕ X ( t ) = E [ e i t x ] = λ 2 β + λ ( λ i t ) 1 + β λ ( λ i t ) 3 ; t R , i = 1 .
H ( x , β , λ ) = λ 2 β λ x 2 + 2 2 λ + β ( λ x ( λ x + 2 ) + 2 ) .
Figure 2 represents some possible pdf shapes of the TPXG distribution for selected values of β and λ , which reveals the flexibility of the distribution in modeling right-skewed observations. Further, Figure 3 indicates the possible shapes of the function H ( x , β , λ ) . They are bathtub, increasing, decreasing, and decreasing-increasing shapes. For more explanations regarding the TPXG distribution, see [1].
When the variable of interest is expensive to measure or difficult to obtain, but cheap and simple to rank, ranked set sampling is recognized as an effective sampling strategy for enhancing the accuracy and efficiency of parameters estimation. McIntyre [2] proposed the ranked set sampling scheme for estimating the pasture and forage yields.
Let X TPXG distribution, with the pdfs f ( x ) and CDF F ( x ) , where μ and σ 2 represent, respectively, the population mean and variance. Let the random sample X 1 , X 2 , ⋯, X k ( i = 1 , 2 , , k ) with the same pdf f ( x ) . The method of the ranked set sampling (RSS) can be described as follows:
  • Choose k simple random sampling (SRS) each of size k (set size) from the underlying population as
    X 1 ( 1 k ) , X 1 ( 2 k ) , , X 1 ( k k ) .
    X 2 ( 1 k ) , X 2 ( 2 k ) , , X 2 ( k k ) .
    X k ( 1 k ) , X k ( 2 k ) , , X k ( k k ) .
  • Rank the units in each set of size k from lowest to the largest visually or based on any cost-free method as
    X 1 ( 1 : k ) , X 1 ( 2 : k ) , , X 1 ( k : k ) .
    X 2 ( 1 : k ) , X 2 ( 2 : k ) , , X 2 ( k : k ) .
    X k ( 1 : k ) , X k ( 2 : k ) , , X k ( k : k ) .
  • Select the ith order statistic (in bold) from the ith set ( i = 1 , 2, , k) as
    X 1 ( 1 : k ) , X 1 ( 2 : k ) , , X 1 ( k : k ) .
    X 2 ( 1 : k ) , X 2 ( 2 : k ) , , X 2 ( k : k ) .
    X k ( 1 : k ) , X k ( 2 : k ) , , X k ( k : k ) .
  • Repeat the Steps (1)–(3), n times (cycles) to obtain an RSS of size N = k n .
The selected RSS units are denoted by X 1 ( 1 : k ) , X 2 ( 2 : k ) , , X k ( k : k ) , where X j ( i : k ) is the ith largest unit in a set of size k in the jth cycle. Notice that even we selected k 2 units, we only measured k of them; these units are not identically distributed, but they are independent because they are selected from different sets.
Takahasi and Wakimoto [3] delivered the mathematical theory of the RSS, and showed that the RSS estimator of mean with the perfect ranking is unbiased and better than the SRS estimator due to its smaller variance. The SRS mean estimator is given by
X ¯ S R S = 1 k i = 1 k X i ,   with variance   Var X ¯ S R S = σ 2 k .
The RSS estimator of the population mean with its variance are given by
X ¯ R S S = 1 n k j = 1 n i = 1 k X j ( i : k ) ,
Var X ¯ R S S = 1 n k 2 i = 1 k Var X ( i : k ) = σ 2 n k 1 n k 2 i = 1 k μ ( i : k ) μ 2 .
Note that since 1 k 2 i = 1 k ( μ ( i : k ) μ ) 2 0 , we have
Var X ¯ R S S Var X ¯ S R S .
They also showed that
μ = 1 k i = 1 k μ ( i : k )   and   f ( x ) = 1 k i = 1 k f ( i : k ) ( x ) ,
where
μ ( i : k ) = x f ( i : k ) ( x ) d x ,   and   σ ( i : k ) 2 = x μ ( i : k ) 2 f ( i : k ) ( x ) d x .
Under perfect rankings, this relation emphasizes the efficiency of the X ¯ R S S mean estimator due to its variance compared to X ¯ S R S for the SRS estimator for the same number of quantified observations regardless of the distribution of nature. Even with a ranking error, Dell and Clutter [4] demonstrated that RSS is more effective than simple random sampling.
Some further modifications of RSS are suggested in the literature, such as extreme RSS by Samawi et al. [5], Mutllak [6] introduced a modification of RSS called median ranked set sampling; another scheme of RSS is proposed by Al-Saleh and Al-Kadiri [7] which is the double RSS, percentile RSS by Mutllak [8], L RSS by Al-Nasser [9], Haq et al. [10] suggested partial RSS design, and neoteric RSS by Zamanzade and Al-Omari [11]. In addition to these modifications, many authors investigated the parameter estimation of some distributions using RSS or its modifications. For example, the logistic model parameters are estimated based on SRS and RSS by Abu-Dayyeh et al. [12]. The generalized quasi-Lindley distribution parameter estimation is considered by Al-Omari et al. [13]. Yousef and Al-Subh [14] used maximum likelihood methods to estimate Gumbel parameters under RSS. Akgul and Şenoglu [15] investigated some modifications of the RSS in estimating the Weibull distribution parameters. The Bayesian and maximum likelihood estimation approaches are considered by Hussian [16] to determine parameter estimates for the Kumaraswamy distribution under RSS. Chen et al. [17] used moving extremes RSS to estimate the scale parameter for the scale distribution. Al-Omari and Bouza [18] considered ratio estimators of the population mean with missing values using RSS. Later, Al-Omari [19] considered the varied L RSS and used the MLE in location-scale families. Hassan et al. [20] used median RSS and estimated the stress–strength reliability for the generalized inverted exponential distribution.
Due to the importance of the TPXG distribution in lifetime distributions and to our knowledge, this is the first study to consider the RSS design for parameter estimations of the TPXG distribution. Hence, the main focus of this paper is to use RSS design for estimating the TPXG distribution parameters and then use some well-known methods of estimation, including the method of maximum product of spacings, maximum likelihood method, ordinary least square method, method of Cramer and von Mises, weight least square method, and the Anderson–Darling method. Then, the suggested estimators based on the RSS design are compared with their competitors in SRS for the same number of measured observations. A real data set is analyzed to explain the usefulness of the offered estimators. Based on the gained results, the RSS estimators are found to be better than the SRS counterparts in terms of the MSE, bias, and efficiency values for all methods of estimation considered in the study.
The layout of this paper is as follows. The estimation methods of the TPXG distribution parameters are presented in Section 2. A simulation study is conducted to show the superiority of the RSS relative to the SRS estimators in Section 3. In Section 4, the suggested estimators’ usefulness is examined using a real data set fitted to the TPXG distribution. The last section will present the conclusion and remarks.

2. Method of Estimation

Here, based on RSS design, six estimation methods are considered to estimate the λ and β parameters of the TPXG distribution, which are: the maximum likelihood (MLE) method, the maximum product of spacings (MPS) method, ordinary least square (OLS) method, weight least square (WLS) method, Cramer–von Mises (CV) method, and Anderson–Darling (AD) method. In all methods, we denote by { X ( i : k ) j , i = 1 , 2 , , k ;   j = 1 , 2 , , n } the ith order statistics from the ith set of size k of the jth cycle and take them to be the RSS data for X with sample size N = k n .

2.1. MLE Method

Considering an RSS sample of size N = k n , the likelihood function is obtained by
M R S S ( β , λ ) = j = 1 r i = 1 k f ( i : k ) ( x ( i : k ) j , β , λ ) ,
with
f ( i : k ) ( x ( i : k ) j , β , λ ) = k ! ( i 1 ) ! ( k i ) ! [ F ( x ( i : k ) j ) ] i 1 [ 1 F ( x ( i : k ) j ) ] k i f ( x ( i : k ) j ) = k ! ( i 1 ) ! ( k i ) ! 1 e λ x ( i : k ) j λ 2 x ( i : k ) j 2 β 2 + λ x ( i : k ) j β + β + λ β + λ i 1 × e λ x ( i : k ) j λ 2 x ( i : k ) j 2 β 2 + λ x ( i : k ) j β + β + λ β + λ k i × λ 2 β λ x ( i : k ) j 2 + 2 e λ x ( i : k ) j 2 λ + β .
Let the log-likelihood function Ψ R S S = log M R S S ( β , λ ) be
Ψ R S S = j = 1 r i = 1 k log f ( i : k ) ( x ( i : k ) j , λ , β ) = j r i = 1 k log k ! ( i 1 ) ! ( k i ) ! + j = 1 r i = 1 k ( i 1 ) log ( F ( x ( i : k ) j ) ) + j = 1 r i = 1 k ( k i ) log ( 1 F ( x ( i : k ) j ) ) + j = 1 r i = 1 k log f ( x ( i : k ) j ) .
The Ψ R S S λ = 0 and Ψ R S S β = 0 cannot be obtained explicitly and they are not in closed form. Hence, they should be solved numerically to find the MLEs, λ ^ R S S M L E and β ^ R S S M L E of λ and β , respectively.

2.2. Method of MPS

Cheng and Amin [21,22] introduced this method, which depends on maximizing the geometric mean of data spacings. Consider X ( 1 : N ) , X ( 2 : N ) , , X ( N : N ) to be an ordered sample forming a RSS of size N = n k from the TPXG distribution. The uniform spacings are given by
υ i ( β , λ ) = F ( x ( i : N ) | β , λ ) F ( x ( i 1 : N ) | β , λ ) , i = 1 , 2 , , N .
Note that F ( x 0 : N | β , λ ) = 0 and F ( x N + 1 : N | β , λ ) = 1 . It is clear that i = 1 N + 1 υ i ( β , λ ) = 1 .
Let the geometric mean of the spacing be
Q ( β , λ | x ) = i = 1 N + 1 υ i ( β , λ ) 1 N + 1 .
The natural logarithm of (11) is
Υ ( β , λ | x ) = 1 N + 1 i = 1 N + 1 log υ i ( β , λ ) .
The estimators, β ^ R S S M P S and λ ^ R S S M P S , are the values of β and λ , which maximize the geometric mean of spacings. The determination of these estimators can be achieved by determining the solution of the following nonlinear equations:
β Υ ( β , λ ) = 1 N + 1 i = 1 N + 1 1 υ i ( β , λ ) Φ 1 ( x ( i : N ) | β , λ ) Φ 1 ( x ( i 1 : N ) | β , λ ) = 0 ,
λ Υ ( β , λ ) = 1 N + 1 i = 1 N + 1 1 υ i ( β , λ ) Φ 2 ( x ( i : N ) | β , λ ) Φ 2 ( x ( i 1 : N ) | β , λ ) = 0 ,
where
Φ 1 ( x ( i : N ) | β , λ ) = β F ( x ( i : N ) | β , λ ) = λ 2 x ( i : N ) · λ x ( i : N ) + 2 e λ x ( i : N ) 2 β + λ 2 ,
and
Φ 2 ( x ( i : N ) | β , λ ) = λ F ( x ( i : N ) | β , λ ) = x ( i : N ) λ · β x ( i : N ) 2 λ 2 + β 2 x ( i : N ) 2 + β x ( i : N ) + 2 λ + 4 β e x ( i : N ) λ 2 λ + β 2 ,
that can be solved numerically.

2.3. Methods of LS

Well-known results in probability theory indicate that F ( x ( i : N ) ) B e t a ( i , N i + 1 ) , where F is a distribution function, and X ( i : N ) are the ith-order statistic of the sample X 1 , X 2 , , X N . Therefore, E [ F ( x ( i : N ) ) ] = i N + 1 and V a r [ F ( x ( i : N ) ) ] = i ( N i + 1 ) ( N + 1 ) 2 ( N + 2 ) .
Using the expectation and variance, two variants of the least squares methods can be obtained. Swain et al. [23] were the first to use the method of LS for parameter estimations of the beta distribution.

2.3.1. OLS Method

The OLS estimators β ^ O L S R S S and λ ^ O L S R S S of β and λ , respectively, can be found by minimizing the following function, with respect to β and λ :
Ω ( β , λ | x ) = i = 1 N F ( x ( i : N ) | β , λ ) i N + 1 2 = i = 1 N 1 β x ( i : N ) 2 λ 2 2 + β x ( i : N ) λ + λ + β e x ( i : N ) λ λ + β i N + 1 2 ,
Alternatively, we can obtain the estimators by solving simultaneously the nonlinear equations:
i = 1 N F ( x ( i : N ) | β , λ ) i N + 1 Φ 1 ( x ( i : N ) | β , λ ) = 0 ,
i = 1 N F ( x ( i : N ) | β , λ ) i N + 1 Φ 2 ( x ( i : N ) | β , λ ) = 0 ,
where Φ 1 ( x ( i : N ) | β , λ ) and Φ 2 ( x ( i : N ) | β , λ ) are defined as in (12) and (13), respectively.

2.3.2. WLS Method

The WLS estimators of β and λ , say, β ^ W L S R S S and λ ^ W L S R S S , respectively, can be determined by minimizing the following function, with respect to β and λ :
W ( β , λ | x ) = i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) F ( x ( i : N ) | β , λ ) i N + 1 2 = i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) 1 β x ( i : N ) 2 λ 2 2 + β x ( i : N ) λ + λ + β e x ( i : N ) λ λ + β i N + 1 2 ,
Note that these estimators are also the solution to the following nonlinear equations:
i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) F ( x ( i : N ) | β , λ ) i N + 1 Φ 1 ( x ( i : N ) | β , λ ) = 0 ,
i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) F ( x ( i : N ) | β , λ ) i N + 1 Φ 2 ( x ( i : N ) | β , λ ) = 0 ,
where Φ 1 ( x ( i : N ) | β , λ ) and Φ 2 ( x ( i : N ) | β , λ ) are specified as in (12) and (13), respectively.

2.4. Methods of Minimum Distances

Several methods of estimation can be proposed based on the minimization of test statistics between the empirical cumulative distribution and theoretical functions. The Cramer–von Mises and Anderson–Darling methods are considered here. (See D’Agostino and Stephens [24]).

2.4.1. CV Method

The CV β ^ C V R S S and λ ^ C V R S S of β and λ , respectively, can be found by minimizing the following function, with respect to β and λ :
C V ( β , λ ) = 1 12 N + i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) F ( x ( i : N ) ; β , λ ) 2 i 1 2 N 2 = 1 12 N + i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) 1 β x ( i : N ) 2 λ 2 2 + β x ( i : N ) λ + λ + β e x ( i : N ) λ λ + β 2 i 1 2 N 2 ,
Consequently, these estimators are also the solution to the nonlinear equations:
i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) F ( x ( i : N ) | β , λ ) 2 i 1 2 N Φ 1 ( x ( i : N ) | β , λ ) = 0 ,
i = 1 N ( N + 1 ) 2 ( N + 2 ) i ( N i + 1 ) F ( x ( i : N ) | β , λ ) 2 i 1 2 N Φ 2 ( x ( i : N ) | β , λ ) = 0 ,
where Φ 1 ( x ( i : N ) | β , λ ) and Φ 2 ( x ( i : N ) | β , λ ) are given in (12) and (13), respectively.

2.4.2. AD Method

The AD estimates of the TPXG distribution parameters, β and λ , denoted by β ^ A D R S S and λ ^ A D R S S , can be gained by minimizing the following function with respect to β and λ :
A ( β , λ ) = N 1 N i = 1 N ( 2 i 1 ) { log F ( x ( i : N ) | β , λ ) + log F ¯ ( x ( N i + 1 : N ) | β , λ ) } , = N 1 N i = 1 N ( 2 i 1 ) log 1 β x ( i : N ) 2 λ 2 2 + β x ( i : N ) λ + λ + β e x ( i : N ) λ λ + β + log β x ( i : N ) 2 λ 2 2 + β x ( i : N ) λ + λ + β e x ( i : N ) λ λ + β ,
or by simultaneously solving the two equations:
A ( β , λ ) β = i = 1 N ( 2 i 1 ) Φ 1 ( x ( i : N ) | β , λ ) F ( x ( i : N ) | β , λ ) Φ 1 ( x ( N i + 1 : N ) | β , λ ) F ¯ ( x ( N i + 1 : N ) | β , λ ) = 0 ,
and
A ( β , λ ) λ = i = 1 N ( 2 i 1 ) Φ 2 ( x ( i : N ) | β , λ ) F ( x ( i : N ) | β , λ ) Φ 2 ( x ( N i + 1 : N ) | β , λ ) F ¯ ( x ( N i + 1 : N ) | β , λ ) = 0 ,
where Φ 1 ( x ( i : N ) | β , λ ) and Φ 2 ( x ( i : N ) | β , λ ) are specified in (12) and (13), respectively.

3. Simulation and Discussion

In this section, a simulation study is supplemented using R software. Based on different parameters and sample sizes for both designs RSS and SRS, 1000 samples from the TPXG distribution are generated. The simulation is performed assuming that the ranking in RSS design is perfect, i.e., there is no error in ranking. The number of cycles is nominated to be n = 10, 15, and 20, while the set size k is designated as 4, 5, and 6. The SRS size is N = n k , which must have the same size in the RSS design. For the purpose of comparison between SRS and RSS methods, the estimated mean squared error (MSE) and efficiency are deliberated for each estimator as:
M S E ( β ^ ) = 1 1000 i = 1 1000 β i ^ β 2 and M S E ( λ ^ ) = 1 1000 i = 1 1000 λ i ^ λ 2 ,
where β i ^ and λ i ^ denote the estimate of β and λ , respectively, for the ith simulated sample. The efficiency (Eff) values of the RSS estimators with respect to the SRS estimators based on the same sample size are defined by
E f f β ^ R S S , β ^ S R S = M S E β ^ S R S M S E β ^ R S S and E f f ( λ ^ R S S λ ^ S R S ) = M S E λ ^ S R S M S E λ ^ R S S .
For various selections of the parameters, sample sizes and number of cycles, the estimates (ES), estimated MSE, and the Eff values are displayed in Table 1, Table 2, Table 3 and Table 4. The results in Table 1, Table 2, Table 3 and Table 4 indicate that:
  • Most of the efficiency values are larger than 1 for all cases considered in this study, indicating that the RSS estimators perform better than the SRS estimators in all methods based on the same number of measured units;
  • Using the RSS design with an increasing number of cycles, the MSE values decrease. As an example, when k = 5, the MSEs of the MLE estimators of λ are 0.0501, 0.0346, and 0.0270 for n = 10, 15, and 20, respectively;
  • Based on the RSS method, the MSE value decreases as k value increases. For example, when n = 15, the MSEs of estimators of β , using the AD method, are 0.4206, 0.2248, and 0.1693 for k = 4, 5, and 6, respectively;
  • Under the SRS technique, the MSE value decreases when N = n k increases;
  • The SRS estimators of λ perform better than the RSS estimators using the MPS method in some cases, otherwise, the RSS estimators are better;
  • The MLE estimators accomplish superior performance compared to the other estimation methods under both RSS and SRS schemes for all results in the tables.

4. Application to Read Data

The usefulness of the proposed RSS estimators is examined in this section using a well-known real data set, which embodies the waiting times (in minutes) before service of 100 bank customers. These data were studied by Ghitany et al. [25]. The data observations are: 2.6, 2.7, 2.9, 3.1, 3.2, 3.3, 3.5, 3.6, 4.0, 4.1, 4.2, 4.2, 4.3, 4.3, 4.4, 4.4, 4.6, 4.7, 4.7, 4.8, 4.9, 4.9, 5.0, 0.8, 0.8, 1.3, 1.5, 1.8, 1.9, 1.9, 2.1, 5.3, 5.5, 5.7, 5.7, 6.1, 7.1, 7.1, 7.4, 7.6, 7.7, 8.0, 8.2, 8.6, 8.6, 8.6, 8.8, 8.8, 11.0, 11.1, 11.2, 6.2, 6.2, 6.2, 8.9, 8.9, 9.5, 9.6, 9.7, 9.8, 10.7, 10.9, 11.0, 6.3, 6.7, 6.9, 7.1, 7.1, 11.2, 11.5, 11.9, 12.4, 12.5, 12.9, 13.0, 13.1, 13.3, 17.3, 17.3, 18.1, 18.2, 18.4, 18.9, 19.0, 19.9, 20.6, 21.3, 21.4, 21.9, 23.0, 27.0, 31.6, 33.1, 13.6, 13.7, 13.9, 14.1, 15.4, 15.4, 38.5.
The TPXGD distribution is fitted to this data. We considered different criteria in this study, such as the Akaike information criterion (AIC), Bayesian information criterion (BIC), Hannan Quinn Information Criterion (HQIC), Consistent Akaike Information Criterion (CAIC). Details of these criteria can be found in Akaike [26], and Schwarz [27], Hannan and Quinn [28] and Bozdogan [29]. Additionally, Kolmogorov–Smirnov (KS) is obtained for each model.
The formulae for these criteria are: AIC = 2 L + 2 h , CAIC = 2 L + 2 2 h n n h 1 , HQIC = 2 log log ( n ) [ h 2 L ] , BIC = 2 L + h log ( n ) , and KS = sup | F n ( x ) F ( x ) | , F n ( x ) = 1 n i n ł x i x , where h is the number of parameters and n is the sample size and L is the value of the maximum log-likelihood function.
Since the distribution under study has two parameters, for fitting the data, we considered two distributions of two parameters—Darna distribution and Marshall–Olkin Esscher transformed Laplace distribution—and one distribution of one parameter, the inverse length-biased Maxwell distribution. The pdfs of these distributions are mentioned below.
  • Darna distribution with pdf:
    f D D ( x ; θ , α ) = θ 2 α 2 + θ 2 2 α + θ 4 x 2 2 α 3 e θ x α ; x > 0 , α > 0 , θ > 0 .
  • Marshall–Olkin Esscher transformed Laplace distribution (MOETL) with pdf:
    f ( x ; λ , k ) = λ k 1 + k 2 e λ k x , x < 0 e k λ x , x 0
  • Inverse length-biased Maxwell distribution (ILBMD) with pdf:
    f I L B M D ( x ; α ) = 1 2 α 4 x 5 e 1 2 α 2 x 2 , 0 < x < , α > 0 .
The results are reported in Table 5. They show that the TPXG distribution provides a superior fit over other competing continuous models, since it has the smallest values for all measures with smallest values of the Kolmogorov–Smirnov distance; Figure 4 supports this claim.
Total Time on the Test (TTT) plot plays a vital role in selecting the proper model for fitting the underlying data regarding the failure rates. This informs us of the altered forms of the model failure rate. If the plot has a straight line, then the given data have a constant failure rate. The failure rates will be decreased if it is convex and increased if this plot is concave. For the bathtub shape, the TTT plot decreases first and then increases. Whereas, if the TTT plot is concave first and then convex, the failure rates will have an inverted bathtub shape. The TTT and density plots for TPXG distribution for the bank customers’ data are given in Figure 5. The probability–probability (P-P) and quantile–quantile (Q-Q) plots for the TPXG model based on the real data are given in Figure 6. Figure 7 presents the box and Bee Swarm plots for these data.
Based on these data, we take an SRS of size 20, while for the ranked set sampling, a small sample size of k = 5 is considered with number of cycles as n = 4 . Table 6 and Table 7 include the RSS (n = 4 and k = 5) and SRS ( N = 20 ) samples taken from the bank customers data. It is of interest to note here that the SRS and RSS methods are compared based on the same number of measured units. Using the previous methods, we calculate the estimates of β and λ in each design. Here, we assumed that the ranking is perfect. To compare estimators, we considered the previous criteria measures, AIC, BIC, CAIC HQC, and KS. The results are summarized in Table 8.
The findings in Table 8 illustrate that the TPXG parameter estimates, based on the RSS method, are improved compared to their counterparts in SRS in terms of the smallest values of AIC, BIC, CAIC HQC, and KS, using the MLE, MPS, OLS, WLS, CV, and AD.

5. Conclusions

This paper discusses several estimation methods for the TPXG distribution parameters based on RSS and SRS designs. A simulation study is performed to compare the performance of these various estimators, considering the same number of measuring units. A real data set is analyzed to illustrate the usefulness of the suggested estimators. Based on the obtained results, the RSS estimators are better than the SRS estimators in terms of MSE, bias, and efficiency values for all estimation methods considered in the paper. For future work, the topics discussed in this paper can be considered under imperfect ranking (Dell and Clutter [4]) based on RSS or using some of its modifications. For some modifications of the RSS, see the balanced groups RSS by Jemain et al. [30], percentile DRSS by Al-Omari and Jaber [31], and new mixed RSS by Hanandeh et al. [32]. Furthermore, one can use other real data sets for additional investigation of the suggested estimators of the TPXG distribution.

Author Contributions

The authors contributed equally to this work. Formal analysis, S.B.; Validation, A.I.A.-O.; Writing—review & editing, A.I.A.-O. and I.M.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deanship of Scientific Research at King Khalid University through the Research Groups Program under the grant number R.G.P. 2/132/43.

Data Availability Statement

The data used in this study are available within the paper.

Acknowledgments

The authors thank and extend their appreciation to the Deanship of Scientific Research at King Khalid University for funding this work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sen, S.; Maiti, S.S.; Chandra, N. The xgamma distribution: Statistical properties and application. J. Mod. Appl. Stat. Methods 2016, 15, 38. [Google Scholar] [CrossRef]
  2. McIntyre, G.A. A method for unbiased selective sampling using ranked set sampling. Aust. J. Agric. Res. 1952, 3, 385–390. [Google Scholar] [CrossRef]
  3. Takahasi, K.; Wakimoto, K. On unbiased estimates of the population mean based on the sample stratified by means of ordering. Ann. Inst. Stat. Math. 1968, 20, 1–31. [Google Scholar] [CrossRef]
  4. Dell, T.R.; Clutter, J.L. Ranked set sampling theory with order statistics background. Biometrics 1972, 28, 545–555. [Google Scholar] [CrossRef]
  5. Samawi, H.; Abu-Daayeh, H.A.; Ahmed, S. Estimating the population mean using extreme ranked set sampling. Biom. J. 1996, 38, 577–586. [Google Scholar] [CrossRef]
  6. Mutllak, H.A. Median ranked set sampling. J. Appl. Stat. Sci. 1997, 6, 245–255. [Google Scholar]
  7. Al-Saleh, M.F.; Al-Kadiri, M.A. Double-ranked set sampling. Stat. Probab. Lett. 2000, 48, 205–212. [Google Scholar] [CrossRef]
  8. Muttlak, H.A. Modified ranked set sampling methods. Pak. J. Stat.-All Ser. 2003, 19, 315–324. [Google Scholar] [CrossRef]
  9. Al-Nasser, A.D. L-Ranked set sampling: A generalization procedure for robust visual sampling. Commun. Stat. 2007, 6, 33–43. [Google Scholar] [CrossRef]
  10. Haq, A.; Brown, J.; Moltchanova, E.; Al-Omari, A.I. Partial ranked set sampling design. Environmetrics 2013, 24, 201–207. [Google Scholar] [CrossRef]
  11. Zamanzade, E.; Al-Omari, A.I. New ranked set sampling for estimating the population mean and variance. Hacet. J. Math. Stat. 2016, 45, 1891–1905. [Google Scholar] [CrossRef]
  12. Abu-Dayyeh, W.A.; Al-Subh, S.A.; Muttlak, H.A. Logistic parameters estimation using simple random sampling and ranked set sampling data. Appl. Math. Comput. 2004, 150, 543–554. [Google Scholar] [CrossRef]
  13. Al-Omari, A.I.; Benchiha, S.; Almanjahie, I.M. Efficient Estimation of the Generalized Quasi-Lindley Distribution Parameters under Ranked Set Sampling and Applications. Math. Probl. Eng. 2021, 1, 9982397. [Google Scholar] [CrossRef]
  14. Yousef, O.M.; Al-Subh, S.A. Estimation of Gumbel parameters under ranked set sampling. J. Mod. Appl. Stat. Methods 2014, 13, 24. [Google Scholar] [CrossRef]
  15. Akgul, F.G.; Şenoglu, B. Estimation of P (X< Y) using ranked set sampling for the Weibull distribution. Qual. Technol. Quant. Manag. 2017, 14, 296–309. [Google Scholar]
  16. Hussian, M.A. Bayesian and maximum likelihood estimation for Kumaraswamy distribution based on ranked set sampling. Am. J. Math. Stat. 2014, 4, 30–37. [Google Scholar]
  17. Chen, W.; Xie, M.; Wu, M. Parametric estimation for the scale parameter for scale distributions using moving extremes ranked set sampling. Stat. Probab. Lett. 2013, 83, 2060–2066. [Google Scholar] [CrossRef]
  18. Al-Omari, A.I.; Bouza, C.N. Ratio estimators of the population mean with missing values using ranked set sampling. Environmetrics 2015, 26, 67–76. [Google Scholar] [CrossRef]
  19. Al-Omari, A.I. Maximum likelihood estimation in location-scale families using varied L ranked set sampling. RAIRO-Oper. Res. 2021, 55, S2759–S2771. [Google Scholar] [CrossRef]
  20. Hassan, A.S.; Al-Omari, A.I.; Nagy, H.F. The Determination of the order of an autoregression. Iran. J. Sci. Technol. Trans. Sci. 2021, 45, 641–659. [Google Scholar] [CrossRef]
  21. Cheng, R.C.H.; Amin, N.A.K. Maximum product-of-spacings estimation with applications to the log-Normal distribution; Technical Report; Department of Mathematics, University of Wales: Cardiff, UK, 1979. [Google Scholar]
  22. Cheng, R.C.H.; Amin, N.A.K. Estimating parameters in continuous univariate distributions with a shifted origin. J. R. Stat. Soc. Ser. 1983, 45, 394–403. [Google Scholar] [CrossRef]
  23. Swain, J.; Venkatraman, S.; Wilson, J. Least squares estimation of distribution function in Johnson’s translation system. J. Stat. Comput. Simul. 1988, 29, 271–297. [Google Scholar] [CrossRef]
  24. D’Agostino, R.B.; Stephens, M.A. Goodness-of-Fit Techniques; Marcel Dekker: New York, NY, USA, 1986. [Google Scholar]
  25. Ghitany, M.E.; Atieh, B.; Nadarajah, S. Lindley distribution and its application. Math. Comput. Simul. 2008, 78, 493–506. [Google Scholar] [CrossRef]
  26. Akaike, H. A new look at the statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–723. [Google Scholar] [CrossRef]
  27. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
  28. Hannan, E.J.; Quinn, B.G. The Determination of the order of an autoregression. J. R. Stat. Soc. Ser. B 1979, 41, 190–195. [Google Scholar] [CrossRef]
  29. Bozdogan, H. Model selection and Akaike’s information criterion (AIC): The general theory and its analytical extensions. Psychometrika 1987, 52, 345–370. [Google Scholar] [CrossRef]
  30. Jemain, A.A.; Al-Omari, A.I.; Ibrahim, K. Balanced groups ranked set sampling for estimating the population median. J. Appl. Stat. 2008, 17, 39–46. [Google Scholar]
  31. Al-Omari, A.I.; Jaber, K. Percentile double ranked set sampling. J. Math. Stat. 2008, 4, 60–64. [Google Scholar] [CrossRef]
  32. Hanandeh, A.A.; Al-Nasser, A.D.; Al-Omari, A.I. New mixed ranked set sampling variations. Stat 2022, 11, e417. [Google Scholar] [CrossRef]
Figure 1. The pdfs plots of the xgamma distribution for different values of the λ .
Figure 1. The pdfs plots of the xgamma distribution for different values of the λ .
Mathematics 10 03170 g001
Figure 2. The pdfs plots of the TPXG distribution for different values of the parameters.
Figure 2. The pdfs plots of the TPXG distribution for different values of the parameters.
Mathematics 10 03170 g002
Figure 3. Plots of H ( x , β , λ ) of the TPXG distribution for (a) β = 2.1 , λ = 3 , (b) β = 2.2 , λ = 4 , (c) β = 0.3 , λ = 2.5 , (d) β = 10 , λ = 0.3 , (e) β = 9.9 , λ = 5 , (f) β = 0.9 , λ = 16 .
Figure 3. Plots of H ( x , β , λ ) of the TPXG distribution for (a) β = 2.1 , λ = 3 , (b) β = 2.2 , λ = 4 , (c) β = 0.3 , λ = 2.5 , (d) β = 10 , λ = 0.3 , (e) β = 9.9 , λ = 5 , (f) β = 0.9 , λ = 16 .
Mathematics 10 03170 g003
Figure 4. Plots of estimated pdf and CDF for the bank customers data.
Figure 4. Plots of estimated pdf and CDF for the bank customers data.
Mathematics 10 03170 g004
Figure 5. TTT and density plots for TPXG distribution for bank customer data.
Figure 5. TTT and density plots for TPXG distribution for bank customer data.
Mathematics 10 03170 g005
Figure 6. P-P and Q-Q plots for TPXG distribution for bank customers.
Figure 6. P-P and Q-Q plots for TPXG distribution for bank customers.
Mathematics 10 03170 g006
Figure 7. Box and Bee Swarm plots for TPXG distribution for bank customers.
Figure 7. Box and Bee Swarm plots for TPXG distribution for bank customers.
Mathematics 10 03170 g007
Table 1. The ESs, MSEs, and Effs values, using the MLE, MPS, and OLS methods for the TPXG model parameters with β = 1 , λ = 1 .
Table 1. The ESs, MSEs, and Effs values, using the MLE, MPS, and OLS methods for the TPXG model parameters with β = 1 , λ = 1 .
nkSamplingWLSCVAD
EsMSEEffEsMSEEffEsMSEEff
104 β ^ RSS1.28220.89702.83950.89200.61661.73051.10811.01363.0978
SRS1.50182.54711.03481.06701.30573.1398
λ ^ RSS1.00470.03341.48670.88370.06340.91570.94630.04631.1537
SRS1.00680.04970.90340.05800.96520.0534
5 β ^ RSS1.20510.64853.95440.89590.47661.98831.08030.71481.7602
SRS1.42972.56441.02640.94761.19591.2582
λ ^ RSS1.00140.02481.72050.90380.04941.02780.95710.03571.2049
SRS1.00920.04260.91800.05080.96680.0430
6 β ^ RSS1.18110.42842.82850.88140.42841.63331.05260.51381.9028
SRS1.30481.21170.98090.69971.15670.9777
λ ^ RSS1.00960.01911.73110.90560.04790.97600.96500.02921.2111
SRS1.00380.03310.91640.04680.96790.0354
154 β ^ RSS1.19800.52202.24030.93270.50511.38421.06140.60181.7132
SRS1.33151.16931.01310.69921.17281.0311
λ ^ RSS1.00920.02121.54170.92170.04061.09950.95740.03491.0294
SRS1.00550.03270.92080.04470.96600.0359
5 β ^ RSS1.12790.38312.13940.88830.37301.46221.03630.48181.2633
SRS1.22330.81960.97460.54541.09540.6086
λ ^ RSS0.99780.01831.66810.91390.04130.91430.96130.02811.1134
SRS0.99790.03050.92940.03780.96980.0313
6 β ^ RSS1.11020.23672.69700.90080.24361.91601.01010.28571.9564
SRS1.18330.63850.97500.46671.08640.5590
λ ^ RSS1.00570.01172.26640.93590.02841.11220.97100.01851.3982
SRS0.99890.02650.94140.03160.97430.0259
204 β ^ RSS1.15350.38511.94530.92910.36671.41941.04850.44121.5535
SRS1.23230.74920.99260.52051.13590.6854
λ ^ RSS1.00440.01621.95530.92950.03311.08390.96920.02261.2817
SRS1.00220.03170.93780.03590.98110.0290
5 β ^ RSS1.11470.23622.06810.92900.24821.52161.03380.28071.5711
SRS1.15840.48850.96380.37771.07080.4410
λ ^ RSS1.00700.01271.60040.94380.02731.00550.97850.01781.2573
SRS1.00040.02020.94400.02750.97640.0224
6 β ^ RSS1.07090.16262.01860.89820.20071.38661.00670.21401.5474
SRS1.13060.32820.96370.27831.07490.3312
λ ^ RSS1.00140.00901.96980.94070.02420.96700.97700.01451.3105
SRS0.99990.01770.95070.02340.98330.0191
Table 2. The ESs, MSEs, and Effs values, using the WLS, CV, and AD methods for the TPXG model parameters with β = 1 , λ = 1 .
Table 2. The ESs, MSEs, and Effs values, using the WLS, CV, and AD methods for the TPXG model parameters with β = 1 , λ = 1 .
nkSamplingWLSCVAD
EsMSEEffEsMSEEffEsMSEEff
104 β ^ RSS1.14540.96272.39711.50241.84203.26631.18090.92352.1388
SRS1.33232.30771.78926.01651.37461.9752
λ ^ RSS0.96090.04211.17351.02970.04231.26950.97360.03871.1148
SRS0.97440.04941.04720.05370.99220.0432
5 β ^ RSS1.10850.69412.03201.38731.16472.11871.14580.66522.0651
SRS1.25101.41031.56512.46761.28971.3736
λ ^ RSS0.96730.03311.19061.02850.03081.33320.98140.03021.2161
SRS0.97870.03941.03600.04110.99180.0368
6 β ^ RSS1.08500.49152.05351.31000.75832.10901.12090.48932.0556
SRS1.21451.00931.46271.59931.23731.0057
λ ^ RSS0.97620.02651.22981.02700.02531.33540.98810.02391.2891
SRS0.98210.03261.02990.03380.99070.0308
154 β ^ RSS1.11450.64841.61971.34260.99271.64181.14370.62881.6188
SRS1.22991.05011.48191.62991.25351.0180
λ ^ RSS0.97440.03001.08101.02420.02881.09400.98600.02731.0980
SRS0.98030.03251.03070.03150.99010.0299
5 β ^ RSS1.06460.45641.41431.24300.63651.45031.08490.43091.5054
SRS1.13700.64541.32340.92321.15510.6487
λ ^ RSS0.97270.02431.14331.01490.02221.23600.98170.02201.2253
SRS0.98130.02781.02270.02750.98680.0270
6 β ^ RSS1.04340.26252.18301.18620.35232.19111.06840.26292.1731
SRS1.12530.57311.27870.77201.13830.5713
λ ^ RSS0.98370.01481.52181.01790.01421.54580.99170.01391.5096
SRS0.98690.02251.02130.02190.99220.0209
204 β ^ RSS1.08010.43471.67501.24400.60071.71311.10110.43271.5615
SRS1.17840.72811.35681.02911.19020.6757
λ ^ RSS0.98050.01921.39831.01800.01851.47100.98720.01861.3114
SRS0.99220.02681.02990.02720.99880.0244
5 β ^ RSS1.06360.26551.66451.19300.34421.68221.08150.26601.6672
SRS1.10580.44201.24120.57891.11900.4434
λ ^ RSS0.98910.01541.20951.02050.01471.22420.99520.01431.2327
SRS0.98820.01871.01970.01800.99330.0176
6 β ^ RSS1.02880.19281.66841.13720.23531.71421.03990.18491.7156
SRS1.10030.32171.21340.40341.10860.3172
λ ^ RSS0.98570.01221.30901.01370.01031.48230.99090.01051.4428
SRS0.99210.01591.01900.01530.99560.0151
Table 3. The ESs, MSEs, and Effs values, using the MLE, MPS, and OLS methods for the TPXG model parameters with β = 0.5 , λ = 1 .
Table 3. The ESs, MSEs, and Effs values, using the MLE, MPS, and OLS methods for the TPXG model parameters with β = 0.5 , λ = 1 .
nkSamplingWLSCVAD
EsMSEEffEsMSEEffEsMSEEff
104 β ^ RSS0.76810.58632.29260.51270.37421.56080.69420.67922.6051
SRS0.90161.34420.58430.58400.80471.7693
λ ^ RSS1.03600.06201.34530.90250.06931.01320.99910.06371.0801
SRS1.03770.08340.91910.07021.01450.0688
5 β ^ RSS0.69680.39472.46770.47680.25871.69240.66060.63181.5033
SRS0.83740.97400.55570.43770.78780.9498
λ ^ RSS1.01780.05011.31000.90060.06060.95380.99060.05431.1423
SRS1.03070.06570.91800.05781.02330.0621
6 β ^ RSS0.60980.21672.89600.41660.16112.16580.55910.29682.9335
SRS0.76010.62770.53800.34890.73580.8708
λ ^ RSS0.99880.03771.64170.88500.05381.02650.96550.04171.3349
SRS1.01880.06190.92390.05521.01930.0557
154 β ^ RSS0.69830.34942.29470.48950.25501.77920.68710.52041.5954
SRS0.77810.80190.55120.45360.74320.8303
λ ^ RSS1.02500.04711.33930.90960.05860.95041.01090.04991.1248
SRS1.01440.06310.92070.05571.01070.0561
5 β ^ RSS0.60960.20532.14090.42480.15141.88250.56780.25701.8184
SRS0.72030.43950.52750.28510.68020.4673
λ ^ RSS1.00500.03461.53340.89740.04961.02260.97860.03981.1677
SRS1.01990.05300.92880.05071.01090.0465
6 β ^ RSS0.58970.15231.81810.42840.13241.50130.56330.21841.5977
SRS0.63010.27690.46940.19870.65170.3489
λ ^ RSS1.00360.02961.69040.90460.04581.05190.98210.03631.1660
SRS0.99330.05010.91430.04821.01350.0424
204 β ^ RSS0.60630.22371.94660.43840.17591.65290.59710.29241.7466
SRS0.68940.43540.51160.29080.68500.5108
λ ^ RSS0.99970.03811.40770.90020.05210.97070.98880.04451.0531
SRS1.00640.05360.92330.05061.01140.0469
5 β ^ RSS0.56800.12702.16270.42300.11691.65050.54580.17152.1895
SRS0.65400.27480.49520.19300.66690.3756
λ ^ RSS0.99490.02701.66910.90430.04261.01780.97700.03291.2466
SRS1.00240.04500.92590.04341.01130.0411
6 β ^ RSS0.54270.10491.78700.41280.10991.37550.52350.15921.5091
SRS0.57910.18740.45140.15120.57180.2402
λ ^ RSS0.98800.02591.39780.90520.04190.94970.96980.03251.0423
SRS0.98380.03620.91490.03980.98220.0339
Table 4. The ESs, MSEs, and Effs values, using the WLS, CV, and AD methods for the TPXG model parameters with β = 0.5 , λ = 1 .
Table 4. The ESs, MSEs, and Effs values, using the WLS, CV, and AD methods for the TPXG model parameters with β = 0.5 , λ = 1 .
nkSamplingWLSCVAD
EsMSEEffEsMSEEffEsMSEEff
104 β ^ RSS0.68740.67051.55210.95021.33242.34100.70760.64322.2283
SRS0.78571.04071.09353.11910.82931.4332
λ ^ RSS0.99040.06381.04741.07130.08231.08471.00360.06331.0894
SRS1.00970.06681.08270.08931.02550.0689
5 β ^ RSS0.64120.49212.18370.84430.87892.17450.65470.44711.8495
SRS0.80651.07471.04071.91120.80520.8269
λ ^ RSS0.98370.05351.11611.05360.06661.13810.99460.05181.0992
SRS1.02290.05971.08730.07581.03180.0570
6 β ^ RSS0.55390.26842.73500.70730.43862.72280.56920.23632.6108
SRS0.72330.73420.91341.19420.73590.6170
λ ^ RSS0.96270.04171.27611.02320.04801.34640.97460.04051.2523
SRS1.00710.05331.07000.06471.01930.0507
154 β ^ RSS0.67970.51871.55510.85610.83611.51320.67410.42061.8506
SRS0.74310.80670.92181.26520.74880.7784
λ ^ RSS1.00400.05001.09431.06550.06111.07261.00830.04781.0979
SRS1.00650.05471.06220.06561.01300.0525
5 β ^ RSS0.56220.24501.86650.68960.36031.85270.58020.22481.9059
SRS0.68960.45730.83020.66750.69810.4285
λ ^ RSS0.97410.03901.15011.02710.04411.18360.98670.03771.1327
SRS1.01190.04491.06220.05221.02090.0427
6 β ^ RSS0.55250.19101.66900.66460.26751.67290.56120.16931.7132
SRS0.63500.31880.75210.44750.63390.2900
λ ^ RSS0.97690.03531.12071.02630.03801.17430.98530.03341.1265
SRS1.00360.03961.04960.04461.00730.0376
204 β ^ RSS0.59420.27781.86450.72080.41411.79320.59350.25181.8232
SRS0.69600.51800.82870.74250.68740.4590
λ ^ RSS0.98590.04191.07681.03610.04831.08820.98990.04021.0705
SRS1.01290.04511.05900.05261.01440.0431
5 β ^ RSS0.53160.15522.30080.62350.21122.28780.54310.15072.0644
SRS0.66460.35700.77060.48330.65880.3112
λ ^ RSS0.96910.03271.16081.00940.03521.17480.97710.03191.0960
SRS1.00960.03791.05050.04131.01210.0350
6 β ^ RSS0.52150.14061.58230.60400.18011.58280.52750.13341.5362
SRS0.57070.22250.65220.28500.57090.2049
λ ^ RSS0.97030.03001.03731.00880.03091.05760.97530.02950.9944
SRS0.98050.03111.01540.03270.98430.0293
Table 5. The goodness-of-fit tests for data sets.
Table 5. The goodness-of-fit tests for data sets.
DensityAICCAICBICHQICKS
TXPG644.6638644.7875649.8741646.77250.4893
DD662.2815662.4053667.4919664.39030.5924
MOETL662.1237662.2474667.334664.23240.6170
ILBMD978.2727978.3135980.8779979.32710.4961
Table 6. RSS sample taken from the bank customers data for n = 4 and k = 5.
Table 6. RSS sample taken from the bank customers data for n = 4 and k = 5.
1st-Cycle:1.98.97.415.438.5
2nd-Cycle:2.11.98.67.413.3
3rd-Cycle:4.86.38.917.318.2
4th-Cycle:0.88.615.418.911.0
Table 7. SRS sample taken form the data set of size 20.
Table 7. SRS sample taken form the data set of size 20.
21.92.94.37.121.36.24.22.712.517.3
18.46.138.511.131.615.433.17.61.513.9
Table 8. Estimates, AIC, BIC, CAIC HQC, and KS in SRS and RSS design using MLE, MPS, OLS, WLS, CV, and AD.
Table 8. Estimates, AIC, BIC, CAIC HQC, and KS in SRS and RSS design using MLE, MPS, OLS, WLS, CV, and AD.
MethodDesignAICBICAICSHQCKS
MLESRS673.043678.253673.167675.1520.223
RSS648.315653.526648.439650.4240.102
MPSSRS680.050685.261680.174682.1590.228
RSS656.563661.773656.686658.6710.125
OLSSRS678.286683.496678.409680.3940.226
RSS650.269655.479650.393652.3780.098
WLSSRS678.174683.384678.298680.2830.235
RSS649.969655.179650.092652.0770.095
CVSRS672.151677.362672.275674.2600.226
RSS647.083652.293647.207649.1920.096
ADSRS674.156679.366674.280676.2650.228
RSS648.987654.198649.111651.0960.097
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al-Omari, A.I.; Benchiha, S.; Almanjahie, I.M. Efficient Estimation of Two-Parameter Xgamma Distribution Parameters Using Ranked Set Sampling Design. Mathematics 2022, 10, 3170. https://doi.org/10.3390/math10173170

AMA Style

Al-Omari AI, Benchiha S, Almanjahie IM. Efficient Estimation of Two-Parameter Xgamma Distribution Parameters Using Ranked Set Sampling Design. Mathematics. 2022; 10(17):3170. https://doi.org/10.3390/math10173170

Chicago/Turabian Style

Al-Omari, Amer Ibrahim, SidAhmed Benchiha, and Ibrahim M. Almanjahie. 2022. "Efficient Estimation of Two-Parameter Xgamma Distribution Parameters Using Ranked Set Sampling Design" Mathematics 10, no. 17: 3170. https://doi.org/10.3390/math10173170

APA Style

Al-Omari, A. I., Benchiha, S., & Almanjahie, I. M. (2022). Efficient Estimation of Two-Parameter Xgamma Distribution Parameters Using Ranked Set Sampling Design. Mathematics, 10(17), 3170. https://doi.org/10.3390/math10173170

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