Next Article in Journal
New Formulas Involving Fibonacci and Certain Orthogonal Polynomials
Next Article in Special Issue
Daily Semiparametric GARCH Model Estimation Using Intraday High-Frequency Data
Previous Article in Journal
A Fast Circle Detection Algorithm Based on Circular Arc Feature Screening
Previous Article in Special Issue
A Gamma Process with Three Sources of Variability
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Survival Analysis and Applications of Weighted NH Parameters Using Progressively Censored Data

by
Ahmed Elshahhat
1,* and
Heba S. Mohammed
2
1
Faculty of Technology and Development, Zagazig University, Zagazig 44519, Egypt
2
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2023, 15(3), 735; https://doi.org/10.3390/sym15030735
Submission received: 29 January 2023 / Revised: 11 March 2023 / Accepted: 13 March 2023 / Published: 15 March 2023
(This article belongs to the Special Issue Mathematical Models and Methods in Various Sciences)

Abstract

:
A new weighted Nadarajah–Haghighi (WNH) distribution, as an alternative competitor model to gamma, standard half-logistic, generalized-exponential, Weibull, and other distributions, is considered. This paper explores both maximum likelihood and Bayesian estimation approaches for estimating the parameters, reliability, and hazard rate functions of the WNH distribution when the sample type is Type-II progressive censored order statistics. In the classical interval setup, both asymptotic and bootstrap intervals of each unknown parameter are constructed. Using independent gamma priors and symmetric squared-error loss, the Bayes estimators cannot be obtained theoretically. Thus, two approximation techniques, namely: Lindley and Markov-Chain Monte Carlo (MCMC) methods, are used. From MCMC variates, the Bayes credible and highest posterior density intervals of all unknown parameters are also created. Extensive Monte Carlo simulations are implemented to compare the performance of the proposed methodologies. Numerical evaluations showed that the estimates developed by the MCMC sampler performed better than the Lindley estimates, and both behaved significantly better than the frequentist estimates. To choose the optimal censoring scheme, several optimality criteria are considered. Three engineering applications, including vehicle fatalities, electronic devices, and electronic components data sets, are provided. These applications demonstrated how the proposed methodologies could be applied in real practice and showed that the proposed model provides a satisfactory fit compared to three new weighted models, namely: weighted exponential, weighted Gompertz, and new weighted Lindley distributions.

1. Introduction

The Nadarajah–Haghighi (NH) distribution, proposed by Nadarajah and Haghighi [1], has a density shape that can be growing or unimodal as well as a hazard rate that can be increasing, decreasing, or constant. Recently, as a new weighted version of the NH distribution, according to the idea of created weighted distributions, Khan et al. [2] proposed the weighted Nadarajah–Haghighi (WNH) distribution as a competitive model to other lifetime models, such as gamma, exponentiated half-logistic, generalized-exponential, Weibull, etc. They also stated that the WNH distribution is suitable for modeling data from different areas; for example, reliability, survival analysis, forest, ecological, etc. Suppose that X is a lifetime random variable of the test unit(s) following the two-parameter WNH ( δ , μ ) distribution. Then, its probability density function (PDF) f ( · ) and cumulative distribution function (CDF) F ( · ) are given, respectively, as
f ( x ; δ , μ ) = 2 δ μ 1 + μ x δ 1 e 1 1 + μ x δ ( 1 + e 1 1 + μ x δ ) 2 , x > 0 ,
and
F ( x ; δ , μ ) = 1 e 1 1 + μ x δ 1 + e 1 1 + μ x δ , x > 0 ,
where δ > 0 and μ > 0 are the shape and scale parameters, respectively. From (1), Khan et al. [2] showed that the WNH distribution could be considered an extended model of three common distributions, namely: NH, standard half-logistic and zero truncated weighted Weibull models, such as:
  • Nadarajah–Haghighi distribution by multiplying (1) by 1 2 1 + e 1 1 + μ x δ 2 .
  • Standard half-logistic distribution by setting Y = ( 1 + μ x ) δ 1 .
  • Truncated weighted Weibull by setting Y = X + μ 1 .
Further, the reliability function (RF) R ( · ) and hazard function (HF) h ( · ) at distinct time t are
R ( t ; δ , μ ) = 2 ( 1 + e 1 + μ t δ 1 ) 1 , t > 0 ,
and
h ( t ; δ , μ ) = δ μ 1 + μ t δ 1 1 + e 1 1 + μ t δ , t > 0 ,
respectively. From (2), where X WNH ( δ , μ ) , two lifetime models can be obtained by making a useful transformation, namely (i) standard half-logistic distribution when Y = ( 1 + μ X ) δ 1 and (ii) zero truncated weighted Weibull distribution when Y = X + μ 1 . Figure 1 depicts the density and hazard rate shapes of the WNH ( δ , μ ) distribution for several choices of its parameters. It reveals that the PDF (1) is positively skewed while the HF (4) has increasing, decreasing, and constant shapes.
In the past decade, researchers have favored dealing with progressive Type-II censoring (PCTII) over other censoring mechanisms because it (i) reduces the total testing time, (ii) reduces the cost of failed units, (iii) terminates the experiment rapidly, and (iv) allows the researcher to make use of the units removed during the experiment for use in other tests. Briefly, the PCTII can be described as: Suppose that n independent and identical units are subjected to a life-testing experiment at time zero, m is a fixed number of failures and the censoring scheme R = ( R 1 , R 2 , , R m ) is also pre-specified. At the time that the first failure is observed (say X ( 1 ) ), R 1 survival units are randomly removed from the experiment. Following the second observed failure (say X ( 2 ) ), R 2 of survival units are randomly removed from n R 1 1 units, and so on. This mechanism continues until the m t h failure occurs and the experiment is stopped. At the same time, all remaining R m = n m i = 1 m 1 R i units are removed. Let x = { x ( 1 ) < x ( 2 ) < < x ( m ) } be a PCTII sample of size m with censoring R from a continuous population.
Then, the likelihood function of x , where Θ is the parameter vector, can be expressed as
L ( Θ | x ) = K i = 1 m f ( x ( i ) ; Θ ) [ 1 F ( x ( i ) ; Θ ) ] R i ,
where K = n ( n R 1 1 ) ( n i = 1 m ( R i + 1 ) ) .
Various studies on the NH distribution using censored data have been made in the literature, readers may refer to Mohie El-Din et al. [3], Ashour et al. [4], Wu and Gui [5], Elshahhat et al. [6], Elshahhat and Abu El Azm [7], Dey et al. [8], Almarashi et al. [9], among others.
Although the NH distribution has received a lot of attention from several authors to the best of our knowledge, the problem of estimating the WNH parameters and/or the RF and HF via PCTII sampling has yet to be investigated. Therefore, this article offers an analysis of PCTII lifetime data when each test unit follows the WNH distribution. The maximum likelihood estimators (MLEs) with their Fisher information members are evaluated through the standard Newton–Raphson algorithm. The associated asymptotic normality of the MLE and log-MLE is used to construct the asymptotic confidence intervals (ACIs). Further, for all unknown quantities, bootstrap-p and bootstrap-t intervals are also obtained. Using the Lindley and Metropolis–Hastings (M-H) methods, the Bayes estimates against the squared-error loss (SEL) function can be easily approximated. Using Markov chain Monte Carlo (MCMC) steps, the Bayes credible interval (BCI) and the highest posterior density (HPD) interval for each unknown parameter are obtained. To compare two (or more) different progressive censoring plans, using their Fisher information, different optimality criteria are taken into account. To assess the performance of the proposed estimation methodologies, an extensive Monte Carlo simulation is performed. Lastly, three applications based on real-life engineering data are presented to show the superiority and flexibility of the WNH model over three weighted distributions (as competitors), namely: weighted exponential, weighted Gompertz, and weighted Lindley distributions. Additionally, the proposed applications aim to demonstrate the applicability of the acquired estimators in a real-life scenario.
The rest of the article is classified as: Section 2 provides the likelihood inference. Two kinds of bootstrapping are presented in Section 3. Bayes approximation techniques are discussed in Section 4. Simulation results are provided in Section 5. In Section 6, a brief explanation of the optimum progressive censoring is presented. Three real applications are presented in Section 7. In Section 8, some conclusions are offered.

2. Likelihood Inference

This section deals with finding the MLEs of the model parameters ( δ , μ ) and the reliability characteristics ( R ( t ) , h ( t ) ) of the WNH distribution based on PCTII data. Two-sided ( 1 τ ) 100 % ACIs of the same unknown parameters are also obtained here.

2.1. Maximum Likelihood Estimators

Suppose X ( i ) , i = 1 , 2 , , m , ( 1 m < n ) is a PCTII sample acquired from the WNH ( δ , μ ) having a fixed R . Substituting (1) and (2) into (5), where x i is used in the place of X ( i ) , we receive
L δ , μ x δ μ m e i = 1 m R i + 1 ψ i δ , μ m i = 1 1 + μ x i δ 1 1 + e ψ i δ , μ R i + 2 ,
where ψ i δ , μ = 1 1 + μ x i δ , i = 1 , 2 , . . . , m .
The log-likelihood function ( L ( · ) log L ( · ) ) of (6) can be expressed (up to proportional) as
L δ , μ x m log δ μ + i = 1 m R i + 1 ψ i δ , μ + δ 1 i = 1 m log 1 + μ x i i = 1 m R i + 2 log 1 + e ψ i δ , μ .
Differentiating (7) partially with respect to δ and μ , two likelihood equations are obtained as
L δ = m δ + i = 1 m R i + 1 ψ i δ δ , μ + i = 1 m log 1 + μ x i i = 1 m R i + 2 ψ i δ δ , μ e ψ i δ , μ 1 + e ψ i δ , μ 1 ,
and
L μ = m μ + i = 1 m R i + 1 ψ i μ δ , μ + δ 1 i = 1 m x i 1 + μ x i 1 i = 1 m R i + 2 ψ i μ δ , μ e ψ i δ , μ 1 + e ψ i δ , μ 1 ,
where ψ i η ( · ) is the first-partial derivative with respect to η , such as
ψ i δ δ , μ = 1 + μ x i δ log 1 + μ x i and ψ i μ δ , μ = δ x i 1 + μ x i δ 1 for i = 1 , 2 , . . . , m .
It is observed, from (8) and (9), that the MLEs δ ^ and μ ^ are derived in a system of nonlinear equations. Thus, a suitable iterative procedure, such as the Newton–Raphson iterative method, is used. For this reason, the ‘ maxLik ’ package suggested by Henningsen and Toomet [10] is considered. Next, by replacing δ and μ with their MLEs δ ^ and μ ^ , respectively, the MLEs R ^ ( t ) and h ^ ( t ) can be easily derived.
Remark 1. 
The results of Khan et al. [2] in the case of the WNH distribution based on complete sampling can be easily obtained as a special case from (7) by setting n = m and R i = 0 , i = 1 , 2 , , m .

2.2. Asymptotic Interval Estimators

To create the 100 ( 1 τ ) % ACIs of δ , μ , R ( t ) , or h ( t ) , the variance-covariance (VCov) matrix of δ ^ and μ ^ must be obtained first as
I 1 ( δ ^ , μ ^ ) L 11 L 12 L 21 L 22 ( δ ^ , μ ^ ) 1 = σ ^ 11 σ ^ 12 σ ^ 21 σ ^ 22 ,
where, from (7), the elements L i j , i , j = 1 , 2 , are obtained and presented in Appendix A.
Under certain regularity conditions, the asymptotic normality of Θ ^ = ( δ ^ , μ ^ ) T is approximately distributed as bivariate normal Θ ^ N ( Θ , I 1 ( Θ ^ ) ) . Thus, following the normal approximation (NA) of the MLEs δ ^ and μ ^ , two bounds 100 ( 1 τ ) % ACI via NA (ACI-NA) of δ and μ can be obtained, respectively, by
δ ^ z τ / 2 σ ^ 11 and μ ^ z τ / 2 σ ^ 22 ,
where z τ / 2 is the upper ( τ / 2 )th probability point of the standard normal distribution.
Further, to obtain the ACIs of R ( t ) and h ( t ) , the delta method is used to approximate their variances, see Greene [11]. However, from the Fisher information (10), the approximated variances of R ^ ( t ) and h ^ ( t ) are given by
σ ^ 2 R ^ ( t ) C R I 1 ( δ , μ ) C R T δ ^ , μ ^ and σ ^ 2 h ^ ( t ) C h I 1 ( δ , μ ) C h T δ ^ , μ ^ ,
respectively, where C R = R δ R μ and C h = h δ h μ , such as
R δ = 2 1 + μ t δ log 1 + μ t e 1 + μ t δ 1 1 + e 1 + μ t δ 1 2 ,
R μ = 2 δ t 1 + μ t δ 1 e 1 + μ t δ 1 1 + e 1 + μ t δ 1 2 ,
h δ = μ 1 + μ t δ 1 1 + e ψ t ; δ , μ 1 + δ log 1 + μ t 1 1 + μ t δ e ψ t ; δ , μ 1 + e ψ t ; δ , μ 1 ,
and
h μ = δ 1 + μ t δ 1 1 + e ψ t ; δ , μ 1 + μ t δ 1 1 + μ t 1 1 + μ t δ e ψ t ; δ , μ 1 + e ψ t ; δ , μ 1 ,
where ψ t ; δ , μ = 1 1 + μ t δ .
Thus, at the confidence percent 100 ( 1 τ ) % , the ACI-NA of R ( t ) and h ( t ) are obtained, respectively, as
R ^ ( t ) z τ / 2 σ ^ 2 R ^ ( t ) and h ^ ( t ) z τ / 2 σ ^ 2 h ^ ( t )
where z τ / 2 is the percentile of the standard normal distribution with right-tail probability ( τ / 2 ) t h .
In practice, the main problem with the ACI-NA method is that it may give a negative lower bound for a lifetime parameter. For handling this disadvantage, Meeker and Escobar [12] developed the normal log-transformed (NL) approximation for the MLE. They also stated that the ACI has the highest CP based on the NL compared to the other NA. Then, the 100 ( 1 τ ) % ACI via NL (ACI-NL) of δ is given by
δ ^ exp z τ / 2 σ ^ 11 δ ^ ,
where, in a similar fashion, the other 100 ( 1 τ ) % ACI-NL of μ , R ( t ) or h ( t ) can be easily computed.

3. Bootstrapping Estimators

If the sample size m is small, due to the law of large numbers, both ACIs developed by NA and NL methods cannot perform well. In this section, two bootstrapping techniques, namely: bootstrap-p and bootstrap-t algorithms, are obtained for δ , μ , R ( t ) , and h ( t ) . To obtain the 100 ( 1 τ ) % bootstrap-p (or bootstrap-t) of δ , μ , R ( t ) or h ( t ) , we first need to receive a bootstrap sample, as shown in Algorithm 1.
Algorithm 1 Bootstrap Sampling:
Step 1: 
Compute the MLEs of δ and μ , using PCTII sample X 1 , X 2 , , X m .
Step 2: 
Generate a bootstrap sample (say X 1 * , X 2 * , , X m * ) from WNH ( δ ^ , μ ^ ) .
Step 3: 
Calculate the bootstrap samples δ * , μ * , R * ( t ) and h * ( t ) based on X 1 * , X 2 * , , X m * .
Step 4: 
Redo Steps 2-3 D times.
Step 5: 
Set the bootstrap samples in ascending order as
δ ( 1 ) * < < δ ( D ) * ; μ ( 1 ) * < < μ ( D ) * ; R ( 1 ) * ( t ) < < R ( D ) * ( t ) ; h ( 1 ) * ( t ) < < h ( D ) * ( t ) .
Then, using the bootstrap-p and bootstrap-t methods, the 100 ( 1 τ ) % bootstrap intervals of δ , μ , R ( t ) and h ( t ) can be constructed as
(i)
Bootstrap-p (Boot-p) Method: The 100 ( 1 τ ) % boot-p intervals of δ , μ , R ( t ) and h ( t ) (say ϕ ) are given, respectively, by
ϕ ^ * D τ 2 , ϕ ^ * D 1 τ 2 .
(ii)
Bootstrap-t (Boot-t) Method: The 100 ( 1 τ ) % boot-t intervals of δ , μ , R ( t ) and h ( t ) are given by
1.
Define the T-statistic for δ , μ , R ( t ) and h ( t ) (say ϕ ) as
T j ϕ = ( v a r ^ ( ϕ ^ ) ) 1 / 2 ( ϕ ^ j * ϕ ^ ) , j = 1 , 2 , , D ,
where ϕ ^ is the MLE of ϕ and v a r ^ ( ϕ ^ ) is the observed estimated variance of ϕ .
2.
Arrange the T-statistics as T ( 1 ) ϕ , T ( 2 ) ϕ , , T ( D ) ϕ .
3.
The 100 ( 1 τ ) % bootstrap-t interval of ϕ is given by
ϕ ^ + T D τ 2 ϕ v a r ^ ( ϕ ^ ) , ϕ ^ + T D 1 τ 2 ϕ v a r ^ ( ϕ ^ ) .

4. Bayes Estimators

The Bayes’ paradigm provides the possibility of incorporating prior information about the unknown parameter(s) of interest. To establish this method, the WNH parameters δ and μ are assumed to be random variables with some extra information included in their prior distributions. One of the most flexible priors is called the gamma conjugate density prior; see Kundu [13]. Gamma density is fairly straightforward and flexible enough to cover a large variety of the experimenter’s prior beliefs. We, therefore, assumed that δ and μ , having independent gamma PDFs as δ G a m m a ( a 1 , b 1 ) and μ G a m m a ( a 2 , b 2 ) , where a i , b i > 0 , i = 1 , 2 , are chosen to reflect prior knowledge about the WNH ( δ , μ ) parameters. Thus, the joint gamma PDF of δ and μ becomes
π δ , μ δ a 1 1 μ a 2 1 e ( b 1 δ + b 2 μ ) , δ , μ > 0 ,
where a i and b i , i = 1 , 2 , are assumed to be known. On the other hand, the most widely used loss is called the SEL function, ( · ) , and is defined as
( η , η ˜ ) = ( η ˜ η ) 2 ,
where η ˜ denotes the Bayes estimate of η .
From (12), the target Bayes estimate is given by the posterior expectation of η . The joint posterior PDF, by combining (11) with (6), of δ and μ is
Ω δ , μ x = C 1 δ m + a 1 1 μ m + a 2 1 e ( b 1 δ + b 2 μ ) e i = 1 m R i + 1 ψ i δ , μ × m i = 1 1 + μ x i δ 1 1 + e ψ i δ , μ R i + 2 ,
where C = 0 0 π δ , μ × L δ , μ x ̲ d δ d μ . However, using (12), the Bayes estimator for any function of δ and μ , say ρ ( δ , μ ) is given by
ρ ˜ δ , μ = E ρ δ , μ x = C 1 0 0 ρ ( δ , μ ) × π δ , μ × L δ , μ x ̲ d δ d μ .
It is clear, based on (14), that the Bayes estimate of ρ ( δ , μ ) is derived in an expression of the ratio of double integrals.
Then, the Bayes solution for each unknown parameter is not available. We, therefore, propose to use both the Lindley and MCMC methods to approximate the Bayes estimates.

4.1. Lindley’s Approximation

In this subsection, following Lindley [14], the closed expressions for the Bayes estimates ρ ( δ , μ ) are obtained. However, setting ( δ , μ ) = ( β 1 , β 2 ) for simplicity, the approximated posterior expectation (say ρ ˜ L ( β 1 , β 2 ) ) of ρ ( β 1 , β 2 ) is given by
ρ ˜ L ( β 1 , β 2 ) = ρ ( β ^ 1 , β ^ 2 ) + 0.5 ( A + Υ 30 B 12 + Υ 03 B 21 + Υ 21 C 12 + Υ 12 C 21 ) + γ 1 A 12 + γ 2 A 21 ,
where all terms in (15) are evaluated at their α ^ and θ ^ , A = i = 1 2 j = 1 2 υ i j σ i j , A i j = υ i σ i i + υ j σ j i , B i j = ( υ i σ i i + υ j σ i j ) σ i i , C i j = 3 υ i σ i i σ i j + υ j ( σ i i σ j j + 2 σ i j 2 ) , υ i = ρ ( β 1 , β 2 ) / β i , υ i j = 2 ρ ( β 1 , β 2 ) / β i β j , Υ i j = i + j L ( β 1 , β 2 ) / β 1 i β 2 j , i , j = 0 , 1 , 2 , 3 , i + j = 3 , γ = log π ( β 1 , β 2 ) , γ i = γ / β i , and σ i j in the ( i , j ) t h element of the VCov matrix (10).
From (15), the Lindley’s estimators δ ˜ L and μ ˜ L of δ and μ , respectively, as
δ ˜ L = δ ^ + γ ^ 1 σ ^ 11 + γ ^ 2 σ ^ 12 + 0.5 ( σ ^ 11 ( Υ ^ 30 σ ^ 11 + 2 Υ ^ 21 σ ^ 12 + Υ ^ 12 σ ^ 22 ) + σ ^ 21 ( Υ ^ 03 σ ^ 22 + 2 Υ ^ 12 σ ^ 21 + Υ ^ 21 σ ^ 11 ) ) ,
and
μ ˜ L = μ ^ + γ ^ 1 σ ^ 21 + γ ^ 2 σ ^ 22 + 0.5 ( σ ^ 12 ( Υ ^ 30 σ ^ 11 + 2 Υ ^ 21 σ ^ 12 + Υ ^ 12 σ ^ 22 ) + σ ^ 22 ( Υ ^ 03 σ ^ 22 + 2 Υ ^ 12 σ ^ 21 + Υ ^ 21 σ ^ 11 ) ) ,
where γ 1 = 1 δ ( a 1 1 ) b 1 , γ 2 = 1 μ ( a 2 1 ) b 2 , and the items Υ i j , i , j = 0 , 1 , 2 , 3 (for i + j = 3 ) are obtained and reported in Appendix B.
Similarly, using (15), the Lindley’s estimators of the reliability characteristic functions R ˜ L ( · ) and h ˜ L ( · ) of R ( · ) and h ( · ) are given (for t > 0 ), respectively, by
R ˜ L ( t ; δ , μ ) = R ( t ; δ ^ , μ ^ ) + γ ^ 1 ( υ ^ 1 σ ^ 11 + υ ^ 2 σ ^ 21 ) + γ ^ 2 ( υ ^ 2 σ ^ 22 + υ ^ 1 σ ^ 12 ) + 0.5 ( ( υ ^ 11 σ ^ 11 + 2 υ ^ 12 σ ^ 12 + υ ^ 22 σ ^ 22 ) + Υ ^ 30 σ ^ 11 ( υ ^ 1 σ ^ 11 + υ ^ 2 σ ^ 12 ) + Υ ^ 03 σ ^ 22 ( υ ^ 2 σ ^ 22 + υ ^ 1 σ ^ 21 ) + Υ ^ 12 ( 3 υ ^ 2 σ ^ 22 σ ^ 12 + υ ^ 1 ( σ ^ 11 σ ^ 22 + 2 σ ^ 12 2 ) ) + Υ ^ 21 ( 3 υ ^ 1 σ ^ 11 σ ^ 12 + υ ^ 2 ( σ ^ 11 σ ^ 22 + 2 σ ^ 12 2 ) ) ) ,
and
h ˜ L ( t ; δ , μ ) = h ( t ; δ ^ , μ ^ ) + γ ^ 1 ( υ ^ 1 σ ^ 11 + υ ^ 2 σ ^ 21 ) + γ ^ 2 ( υ ^ 2 σ ^ 22 + υ ^ 1 σ ^ 12 ) + 0.5 ( ( υ ^ 11 σ ^ 11 + 2 υ ^ 12 σ ^ 12 + υ ^ 22 σ ^ 22 ) + Υ ^ 30 σ ^ 11 ( υ ^ 1 σ ^ 11 + υ ^ 2 σ ^ 12 ) + Υ ^ 03 σ ^ 22 ( υ ^ 2 σ ^ 22 + υ ^ 1 σ ^ 21 ) + Υ ^ 12 ( 3 υ ^ 2 σ ^ 22 σ ^ 12 + υ ^ 1 ( σ ^ 11 σ ^ 22 + 2 σ ^ 12 2 ) ) + Υ ^ 21 ( 3 υ ^ 1 σ ^ 11 σ ^ 12 + υ ^ 2 ( σ ^ 11 σ ^ 22 + 2 σ ^ 12 2 ) ) ) ,
where υ i , i = 1 , 2 is given in Section 2.2 and υ i j , i = 1 , 2 of R ( t ) , and h ( t ) are obtained and presented in Appendix C.
Although the Lindley estimators of δ , μ , R ( t ) , and h ( t ) are expressed explicitly, the main problem of this approach is its disability to obtain an interval estimate. Consequently, to approximate the Bayes estimates and to construct their BCI and HPD interval estimates, the M-H algorithm is considered.

4.2. M-H Sampler

Since the joint posterior PDF cannot be reduced by analytical steps to any familiar form, the M-H algorithm (which is one of the effective MCMC techniques) is used, see Gelman et al. [15].
First, from (13), the conditional distribution of δ and μ must be obtained, respectively, as
Ω 1 δ μ , x δ m + a 1 1 e b 1 δ e i = 1 m R i + 1 ψ i δ , μ i = 1 m 1 + μ x i δ 1 + e ψ i δ , μ R i + 2 ,
and
Ω 2 μ δ , x μ m + a 2 1 e b 2 μ e i = 1 m R i + 1 ψ i δ , μ i = 1 m 1 + μ x i δ 1 1 + e ψ i δ , μ R i + 2 .
As we anticipated, the conditional distributions (16) and (17) of δ and μ , respectively, cannot be reduced to any standard distribution. Figure 2 shows that the distributions (16) and (17) behave like a normal density. Therefore, to obtain the Bayes MCMC estimates along with their BCI and HPD intervals, we perform the generation process described in Algorithm 2.
Algorithm 2 Markov Chain Monte Carlo Sampling:
Step 1: 
Start with initial guesses δ ( 0 ) = δ ^ and μ ( 0 ) = μ ^ .
Step 2: 
Set j = 1
Step 3: 
Generate δ * and μ * from N ( δ ^ , σ ^ 11 ) and N ( μ ^ , σ ^ 22 ) , respectively.
Step 4: 
Obtain Δ δ = Ω 1 δ * μ ( j 1 ) , x Ω 1 δ ( j 1 ) μ ( j 1 ) , x and Δ μ = Ω 2 μ * δ ( j ) , x Ω 2 μ ( j 1 ) δ ( j ) , x .
Step 5: 
Generate sample variates u 1 and u 2 from the uniform U ( 0 , 1 ) distribution and set
δ ( j ) = δ * , if u 1 min 1 , Δ δ , δ ( j 1 ) , otherwise . and μ ( j ) = μ * , if u 2 min 1 , Δ μ , μ ( j 1 ) , otherwise .
Step 6: 
Obtain R ( j ) ( t ) and h ( j ) ( t ) , for t > 0 by replacing δ and μ with their δ ( j ) and μ ( j ) , respectively.
Step 7: 
Put j = j + 1 .
Step 8: 
Repeat Steps 3-7 B times and obtain δ ( j ) and μ ( j ) for j = 1 , 2 , , B .
Step 9: 
Compute the Bayes estimates of δ , μ , R ( t ) or h ( t ) (say ρ ˜ M H ) under the SEL (12) as
ρ ˜ M H = j = B 0 + 1 B ρ ( j ) B B 0 ,
where B 0 is burn-in.
Step 10: 
Compute the 100 ( 1 τ ) % BCI of ρ by ordering ρ ( j ) , j = B 0 + 1 , B 0 + 2 , , B as ρ ( B 0 + 1 ) , ρ ( B 0 + 2 ) , , ρ ( B ) . Thus the 100 ( 1 τ ) % BCI of ρ is obtained as
ρ ( ( B B 0 ) τ 2 ) , ρ ( ( B B 0 ) ( 1 τ 2 ) ) .
Step 11: 
Compute the 100 ( 1 τ ) % HPD interval of ρ as
ρ ( j * ) , ρ ( j * + B B 0 ( 1 τ ) ) ,
where j * = B 0 + 1 , B 0 + 2 , , B is selected so that
ρ j * + [ B B 0 1 τ ] ρ j * = min 1 j τ ( B B 0 ) ρ ( j + [ ( B B 0 ) 1 τ ] ) ρ ( j ) ) .

5. Numerical Comparisons

To compare the efficiency of the proposed estimators of δ , μ , R ( t ) , and h ( t ) described in the preceding sections, based on 2000 PCTII samples generated from WNH ( 0.3 , 0.1 ) via Algorithm 3, an extensive Monte Carlo simulation is conducted. At distinct time t = 5 , the corresponding actual values of R ( t ) and h ( t ) are taken as 0.935 and 0.012, respectively. Various combinations of the total sample size n and effective sample size m are also considered, such as n ( = 50 , 100 , 150 ) , and m is taken as 50% and 80% for each n. Three different progressive designs, where 0 * ( m 1 ) means that 0 is repeated ( m 1 ) times, are used, namely:
S 1 : R = ( n m , 0 * ( m 1 ) ) , S 2 : R = 0 * m 2 1 , n m , 0 * m 2 , i f m i s e v e n ; R = 0 * ( m 1 ) 2 , n m , 0 * ( m 1 ) 2 , i f m i s o d d , S 3 : R = ( n m ) 2 , 0 * ( m 2 ) , ( n m ) 2 , i f m i s e v e n ; R = ( m 1 ) 2 , 0 * ( m 2 ) , ( m + 1 ) 2 , i f m i s o d d .
To draw a PCTII sample from the WNH distribution for pre-specified values of n, m, and R , we perform the following steps:
Algorithm 3 Progressive Type-II Censored Sampling:
Step 1: 
Create ϵ independent observations with size m as ϵ 1 , ϵ 2 , , ϵ m from uniform U ( 0 , 1 ) distribution.
Step 2: 
For specific n, m, and R , set ν i = ϵ i i + j = m i + 1 m R j 1 , i = 1 , 2 , , m .
Step 3: 
Set u i = 1 ν m ν m 1 ν m i + 1 for i = 1 , 2 , , m .
Step 4: 
Set x i = μ 1 1 log 1 u i 1 + u i 1 δ 1 , i = 1 , 2 , , m , is a PCTII sample from WNH ( δ , μ ) .
Using each simulated sample, the point estimators, such as maximum likelihood, Lindley, and MCMC, are compared using their mean absolute bias (MAB), and root mean squared error (RMSE) simulated values. Further, the 95% interval estimators, called ACI-NA, ACI-NL, Boot-p, Boot-t, BCI, and HPD intervals, are examined via their simulated average confidence length (ACL) and coverage probability (CP) values. Following Algorithm 1, the bootstrapping procedure of δ , μ , R ( t ) , and h ( t ) is repeated 10,000 times.
To compute the Bayes estimates and BCI/HPD intervals of δ , μ , R ( t ) , or h ( t ) , we repeat Algorithm 2 12,000 times and ignore the first 2000 times as burn-in. For both WNH parameters δ and μ , two different informative sets of the hyperparameters ( a i , b i ) , i = 1 , 2 , are used, namely: Prior-1 (P1): ( a 1 , a 2 , b 1 , b 2 ) = ( 1.5 , 0.5 , 5 , 5 ) and Prior-2 (P2): ( a 1 , a 2 , b 1 , b 2 ) = ( 3 , 1 , 10 , 10 ) . These prior sets are specified based on two criteria called prior mean and prior variance in such a way that the prior mean becomes the plausible value of the corresponding unknown parameter. All numerical computations are carried out via R 4.1.2 software utilizing two recommended packages, namely: ‘ coda ’ and ‘ maxLik ’ packages, as suggested by Plummer et al. [16] and Henningsen and Toomet [10], respectively.
A heatmap is a graphical representation of numerical data used to show relationships between two variables of interest so that individual data points contained in a data set are represented via different colors. Therefore, by R software, the heatmap data tool is used to display the simulated results, including: MAB, RMSE, ACL, and CP values of δ , μ , R ( t ) , and h ( t ) , see Figure 3, Figure 4, Figure 5 and Figure 6, respectively. Each heatmap displays the proposed estimation methods and the specified test settings on the “ x-axis ” and “ y-axis ” lines, respectively. On the other hand, all simulation results of δ , μ , R ( t ) , and h ( t ) are also available in Table 1, Table 2, Table 3 and Table 4.
From Figure 3, Figure 4, Figure 5 and Figure 6, we can draw the following observations:
  • Overall, the Bayes MCMC estimates of all unknown parameters perform better than the MLEs in terms of the smallest values of MABs, RMSEs, and ACLs, as well as the highest CP values.
  • As n (or m) increases, the MABs, RMSEs, and ACLs tend to decrease (except CPs increase), which infers that the proposed estimates are asymptotically unbiased and consistent. This result is also reached when the total number of progressive items, i = 1 m R i , decreases.
  • Comparing the proposed point estimation methods on the basis of the smallest MABs and RMSEs, it is clear that the performance of MCMC estimates of all unknown parameters is highly recognizable compared to those obtained based on the Lindley approach, and both exhibit reasonably good performance compared to the likelihood method. It is an anticipated result because the calculated MCMC estimates included additional prior information.
  • Comparing the proposed interval estimation methods in terms of the smallest ACLs and highest CPs, it is observed that (i) ACI-NA performs better than ACI-NL, (ii) Boot-t performs better than Boot-p, and (iii) the HPD interval performs better than BCI. This conclusion holds for all unknown parameters. It is also clear, in most cases, that the CPs of classical (or Bayes) interval estimates are mostly close to (or below) the specified nominal level.
  • In particular, based on the ACL and CP criteria, it is also noted that the HPD interval of all unknown parameters provides the best bounds compared to its competitive intervals.
  • Since the variance in P2 is smaller than the variance in P1, consequently, the Bayes estimates (or BCI/HPD intervals) perform better based on P2 than the others.
  • Comparing the proposed censoring schemes, both point and interval estimation methods of δ provide good behavior under S 1 (where the survival n m units are removed in the first stage) compared with other competing schemes. Moreover, both point and interval estimation methods of μ , R ( t ) , and h ( t ) provide good behavior under S 2 (where the survival n m units are removed in the middle stage) compared to others.
  • Lastly, to evaluate the unknown parameter(s) of the WNH lifetime model in the presence of data collected from Type-II progressive censoring, the Bayes MCMC methodology is recommended.

6. Optimum PCTII Fashions

The issue of determining the optimum progressive censoring (removal pattern design R ) from a group of all possible plans is a great challenge for any practitioner. In this context, independently, Balakrishnan and Aggarwala [17] and Ng et al. [18] first studied this issue through various setups. The main idea of optimal progressive censoring is, for fixed n and m, that the experimenter chooses the scheme R = ( R 1 , R 2 , , R m ) in the sense it provides a significant amount of information of the unknown parameter(s) under study. According to the experimenter’s knowledge about the availability of test items, experimental sources, and test cost, the best censoring plan is selected.
In the literature, for different lifetime models, various works have investigated the problem of optimum censoring; see, for example, Pradhan and Kundu [19], Elshahhat and Nassar [20], Elshahhat and Abu El Azm [7], among others. In Table 5, to select the optimal PCTII plan, several criteria are listed.
From Table 5; regarding criterion 1, the experimenter aims to increase the main diagonal values of the observed Fisher information. Further, regarding criteria 2 and 3, the experimenter aims to reduce the trace and determinant values of the estimated VCov matrix. Under criterion 4, the experimenter aims to minimize the estimated variance in the logarithm of the MLE of ϱ th quantile, log ( T ^ ϱ ) for 0 < ϱ < 1 .
The optimized censoring should coincide with the largest value of criterion 1 and the smallest value of the other criteria. According to criterion 4, using the delta method, the variance of log ( T ^ ϱ ) is given by
v a r ^ ( log ( T ^ ϱ ) ) = [ log ( T ^ ϱ ) ] T I 1 ( δ ^ , μ ^ ) [ log ( T ^ ϱ ) ] ,
where [ log ( T ^ ϱ ) ] T = δ log ( T ^ ϱ ) μ log ( T ^ ϱ ) ( δ = δ ^ , μ = μ ^ ) .

7. Three Engineering Applications

To show the applicability of the proposed estimation methodologies to real phenomena, we will analyze three different actual data sets. The first data set (say Data A), reported by Mann [21] and re-analyzed by Alotaibi et al. [22], consists of the number of vehicle fatalities for 39 counties in South Carolina in 2012. The second data set (say Data B), given by Wang [23] and discussed earlier by Elshahhat and Abu El Azm [7], represents the times to failure of eighteen electronic devices. The last data set (say Data C), taken from Lawless [24] and Mohammed et al. [25], represents the failure times (in minutes) of fifteen electronic components in an accelerated life test. These data sets are reported in Table 6.
To highlight the flexibility and adaptability of the proposed model, based on different criteria, the WNH distribution is compared with three weighted distributions as competitors, namely: weighted exponential (WE) by Gupta and Kundu [26], weighted Gompertz (WG) by Bakouch et al. [27] and new weighted Lindley (WL) by Asgharzadeh et al. [28] distributions. The respective PDFs of the compared WE, WG, and WL distributions (for x > 0 and δ , μ > 0 ) are:
  • f W E ( x ) = δ + 1 δ μ e μ x ( 1 e δ μ x )
  • f W G ( x ) = μ δ 2 1 + μ δ ( e μ x + μ 1 ) e μ x δ ( e μ x 1 )
  • f W L ( x ) = μ 2 ( δ + 1 ) 2 ( 1 + x ) e μ x δ 2 + 2 δ + δ μ ( δ + 1 ) ( 1 e δ μ x )
For each data set, the parameter estimates of δ and μ with their standard errors (SEs) of each competing model are evaluated by the ML method. Moreover, the estimated values of the negative log-likelihood (NL), Akaike information (AI), consistent Akaike information (CAI), Bayesian information (BI), Hannan–Quinn information (HQI) and the Kolmogorov–Smirnov statistic with its P-value are computed and reported in Table 7. Figure 7 represents the probability-probability (PP) plots of the WNH, WE, WG and WL distributions. It indicates that the plotted points of the fitted WNH distribution based on each data set are quite close to the empirical distribution line. It also confirms the same findings in Table 7. Moreover, Figure 8 displays the estimated PDFs with their histograms as well as the estimated RFs for the compared distributions under the three data sets. It is clear that the fitted WNH density captures the empirical histograms well. To sum up, we concluded that the WNH distribution has a favorable superiority over other weighted distributions because it has the smallest values of the fitting criteria based on the given data sets. For more illustration, for each real data set, the empirical/estimated TTT transform plot of the WNH distribution is provided in Figure 9. It exhibits that the given real data sets have a constant hazard rate function.
Next, using various choices of m and R , three different PCTII samples from each complete real data set are generated and listed in Table 8. For each PCTII sample in Table 8, the theoretical results of δ , μ , R ( t ) , and h ( t ) are evaluated numerically via R software. The point estimates (including: MLEs, Lindley, and MCMC) with their SEs as well as the 95% interval estimates (including: ACI-NA, ACI-NL, Boot-p, Boot-t, BCI, and HPD) of δ , μ , R ( t ) , and h ( t ) (at t = 10 ) are presented in Table 9 and Table 10, respectively. To examine the uniqueness and existence of the MLEs of δ and μ for Ai, Bi, and Ci for i = 1 , 2 , 3 samples, the contour plots of the log-likelihood function (7) are shown in Figure 10. It shows that all acquired estimates of δ ^ and μ ^ exist and are unique. For obtaining the Bayes MCMC estimates along with their BCI/HPD intervals based on noninformative gamma prior, the M-H algorithm proposed in Section 4.2 is repeated 30,000 times, and the first 5000 draws are thrown out. From Table 9 and Table 10, in terms of minimum standard errors and interval lengths, it can be seen that the MCMC estimates are more precise than the classical estimates. Moreover, the two interval bounds obtained by asymptotic, bootstrap, and credible intervals are close to each other.
In Figure 11, based on samples A1, B1, and C1 (as an example), the simulated estimates with their histograms for the marginal posterior distribution functions of δ , μ , R ( t ) , and h ( t ) are plotted with the Gaussian line. It is evident that the MCMC estimates behave appropriately with the normal distribution, where (:) represents the mean sample line. In addition, trace plots using A1, B1, and C1 of the same parameters are displayed in Figure 11 where the sample mean, and 95 % BCI/HPD intervals are shown with solid (—), dotted (…), and dashed (- - -) horizontal lines, respectively. It shows that the MCMC iterations based on each generated PCTII sample are mixed well.
Using the generated PCTII samples in Table 8, the optimal progressive censoring plans, with respect to different criteria reported in Table 5, are obtained and reported in Table 11. It is noted, from Table 11, that
  • For Data A; the censoring schemes used in samples A1 (with respect to criteria 1 and 3) and A3 (with respect to criteria 2 and 4) are the optimum plans.
  • For Data B; the censoring schemes used in samples B1 (with respect to criterion 3), B2 (with respect to criterion 1), and B3 (with respect to criteria 2 and 4) are the optimum plans.
  • For Data C; the censoring schemes used in samples C1 (with respect to criteria 2 and 3), C2 (with respect to criterion 1) and C3 (with respect to criterion 4) are the optimum plans.
Table 11 also supports the same results of the proposed censoring plans (S i , i = 1 , 2 , 3 , ) discussed in the simulation section. Briefly, we concluded that the WNH lifetime model provides a good and adequate explanation of the proposed censoring plan for all given real data sets.

8. Concluding Remarks

This article derives both point and interval estimators of the parameters, reliability, and hazard functions of the new, two-parameter weighted Nadarajah–Haghighi distribution based on Type-II progressive censored sample via the maximum likelihood and Bayesian estimation methods. Further, two parametric bootstrap confidence intervals of the same unknown parameters have been derived. Since the Bayes estimators cannot be written in closed forms, Lindley and MCMC approximations have been investigated. Numerical comparisons have been performed to judge the behavior of the proposed estimators, and they showed that the Bayes estimates using the Metropolis–Hastings algorithm behave well compared to Lindley’s procedure, and both behave much better than the frequentist estimates. Three real data sets taken from several areas, namely: insurance, reliability, and accelerated life testing, have been analyzed to show the usefulness of the proposed model in real practice. All theoretical results have been evaluated utilizing two useful packages in R software, namely: ‘ maxLik ’ and ‘ coda ’ packages. As a result, the Bayes MCMC paradigm is recommended to estimate the model parameters and/or the reliability characteristic of the weighted Nadarajah–Haghighi model in the presence of a sample obtained from Type-II progressive censoring. As a future work, one can extend the results and methodologies discussed here to accelerate life tests or other censoring plans. We also hope that the proposed methodologies will be beneficial to researchers and reliability practitioners.

Author Contributions

Methodology, A.E. and H.S.M.; Funding acquisition, H.S.M.; Software, A.E.; Supervision A.E.; Writing—original draft, H.S.M.; Writing—review and editing A.E. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R175), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article.

Acknowledgments

The authors would also like to express their gratitude to the Editor-in-Chief and anonymous referees for their constructive comments and suggestions, which significantly improved the paper. Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R175), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Fisher elements L i j , i , j = 1 , 2 from (7) with respect to δ and μ are
L 11 = m δ 2 + i = 1 m R i + 1 ψ i δ δ , μ i = 1 m R i + 2 κ ( x i ; δ , μ ) ψ i δ δ , μ ψ i δ δ , μ 2 κ ( x i ; δ , μ ) 1 ,
L 22 = m μ 2 + i = 1 m R i + 1 ψ i μ δ , μ δ 1 i = 1 m x i 2 1 + μ x i 2 i = 1 m R i + 2 κ ( x i ; δ , μ ) ψ i μ δ , μ ψ i μ δ , μ 2 κ ( x i ; δ , μ ) 1 ,
and
L 12 = i = 1 m R i + 1 ψ i δ μ δ , μ + i = 1 m x i 1 + μ x i 1 i = 1 m R i + 2 κ ( x i ; δ , μ ) ψ i δ μ δ , μ ψ i δ δ , μ ψ i μ δ , μ κ ( x i ; δ , μ ) 1 ,
where ψ i η is the second-partial derivative with respect η to for i = 1 , 2 , , m as
ψ i δ δ , μ = 1 + μ x i δ log 1 + μ x i 2 , ψ i μ δ , μ = δ δ 1 x i 2 1 + μ x i δ 2 ,
ψ i δ μ δ , μ = x i 1 + μ x i δ 1 1 + δ log 1 + μ x i and κ ( x i ; δ , μ ) = e ψ ( x i ; δ , μ ) ( 1 + e ψ ( x i ; δ , μ ) ) 1 .

Appendix B

Lindley’s elements Υ i j , i , j = 0 , 1 , 2 , 3 , i + j = 3 from (7) with respect to δ and μ are
Υ 12 = i = 1 m ( R i + 1 ) x i 2 ( 1 + μ x i ) δ 2 [ 1 2 δ δ ( δ 1 ) log ( μ x i ) ] i = 1 m x i 2 ( 1 + μ x i 2 ) + i = 1 m ( R i + 2 ) x i 2 κ ( x i ; δ , μ ) { δ ( 1 + μ x i ) δ κ ( x i ; δ , μ ) × [ { ( δ 1 ) ( 1 + μ x i ) δ 2 δ ( 1 + μ x i ) 2 ( δ 1 ) } log ( μ x i ) + ( 1 + μ x i ) δ 1 [ ( 1 + μ x i ) 1 + δ log ( μ x i ) ( 1 + μ x i ) δ 1 κ ( x i ; δ , μ ) ] ] ( 1 + μ x i ) δ 2 + δ ( 1 + μ x i ) 1 [ ( 1 + μ x i ) δ 1 ( 1 + μ x i ) 2 δ 1 ] + δ log ( μ x i ) [ ( δ 1 ) ( 1 + μ x i ) δ 2 ( 2 δ 1 ) ( 1 + μ x i ) 2 ( δ 1 ) ] + δ ( 1 + μ x i ) δ 1 [ ( 1 + μ x i ) 1 + δ log ( μ x i ) ( 1 + μ x i ) δ 1 κ ( x i ; δ , μ ) ] δ ( 1 + μ x i ) δ 1 ( 1 κ ( x i ; δ , μ ) ) { ( 1 + μ x i ) δ [ ( 1 + μ x i ) 1 + δ log ( μ x i ) ( 1 + μ x i ) δ 1 κ ( x i ; δ , μ ) ] + δ log ( μ x i ) [ ( 1 + μ x i ) δ 1 ( 1 + μ x i ) 2 δ 1 ] } } ,
Υ 21 = i m ( R i + 1 ) x i log ( μ x i ) ( 1 + μ x i ) δ 1 ( 2 + δ log ( μ x i ) ) + i m ( R i + 2 ) x i log ( μ x i ) κ ( x i ; δ , μ ) { × ( 1 + μ x i ) 2 δ κ ( x i ; δ , μ ) [ 1 + δ log ( μ x i ) { 1 ( 1 + μ x i ) δ ( 1 κ ( x i ; δ , μ ) ) } ] + ( 1 + μ x i ) δ 1 [ 2 + δ log ( μ x i ) ( 1 + μ x i ) δ κ ( x i ; δ , μ ) ] + δ log ( μ x i ) ( 1 + μ x i ) δ 1 [ 1 2 ( 1 + μ x i ) δ ] ( 1 κ ( x i ; δ , μ ) ) ( 1 + μ x i ) δ { ( 1 + μ x i ) δ 1 [ 1 + δ log ( μ x i ) ( 1 + μ x i ) δ κ ( x i ; δ , μ ) ] + δ log ( μ x i ) ( 1 + μ x i ) δ 1 [ 1 ( 1 + μ x i ) δ ] ( 1 + μ x i ) 2 δ 1 } } ,
Υ 30 = 2 m δ 3 + i = 1 m R i + 1 ψ i δ + i = 1 m R i + 2 κ ( x i ; δ , μ ) × { [ ( 1 + μ x i ) 2 δ ( κ ( x i ; δ , μ ) + ( 1 + μ x i ) δ 1 ) ] [ ( 1 + μ x i ) δ ( κ ( x i ; δ , μ ) 1 ) ] + [ ( 1 + μ x i ) 3 δ ( κ ( x i ; δ , μ ) + 2 ( 1 + μ x i ) δ 1 ) ] κ ( x i ; δ , μ ) + ( 1 + μ x i ) δ 2 ( 1 + μ x i ) 2 δ } ,
Υ 03 = 2 m μ 3 + i = 1 m R i + 1 ψ i μ + δ i = 1 m x i 3 R i + 2 κ ( x i ; δ , μ ) × { ( δ 1 ) ( δ 2 ) ( 1 + μ x i ) δ 3 + δ { [ ( 1 + μ x i ) δ 1 ( κ ( x i ; δ , μ ) 1 ) ] × [ ( 1 + μ x i ) δ 2 ( δ 1 ) + δ ( 1 + μ x i ) 2 ( δ 1 ) ( κ ( x i ; δ , μ ) 1 ) ] + [ 2 ( δ 1 ) ( 1 + μ x i ) 2 ( δ 1 ) 1 + δ ( 1 + μ x i ) 2 ( δ 1 ) + δ 1 ( κ ( x i ; δ , μ ) 1 ) ] × κ ( x i ; δ , μ ) 2 ( δ 1 ) ( 1 + μ x i ) 2 ( δ 1 ) 1 ) } } ,
where ψ i η is the third-partial derivative with respect η to for i = 1 , 2 , , m as
ψ i δ δ , μ = 1 + μ x i δ log 1 + μ x i 3 and ψ i μ δ , μ = δ ( δ 1 ) ( δ 2 ) x i 3 1 + μ x i δ 3 .

Appendix C

The second-order partial derivatives of R ( t ) from (3) with respect to δ and μ are
υ 11 = 2 ( log ( μ t ) ) 2 ( 1 + μ t ) 2 δ e ( 1 + μ t ) δ 1 ( 1 + e ( 1 + μ t ) δ 1 ) 2 [ 1 + ( 1 + μ t ) δ 2 ϑ ( t ; δ , μ ) ] ,
υ 22 = 2 δ t 2 e ( 1 + μ t ) δ 1 ( 1 + e ( 1 + μ t ) δ 1 ) 2 × { ( δ 1 ) ( 1 + μ t ) δ 2 + δ ( ( 1 + μ t ) 2 ( δ 1 ) 2 ( ( ( 1 + μ t ) 2 ( δ 1 ) ϑ ( t ; δ , μ ) ) ) ) } ,
and
υ 12 = 2 t e ( 1 + μ t ) δ 1 ( 1 + e ( 1 + μ t ) δ 1 ) 2 [ ( 1 + μ t ) δ ( ( 1 + μ t ) 1 2 ( δ log ( μ t ) ( 1 + μ t ) δ 1 ϑ ( t ; δ , μ ) ) ) + δ log ( μ t ) ( ( 1 + μ t ) 2 δ 1 + ( 1 + μ t ) δ 1 ) ] .
where ϑ ( t ; δ , μ ) = e ( 1 + μ t ) δ 1 ( 1 + e ( 1 + μ t ) δ 1 ) 1 .
Similarly, the second-order partial derivatives of h ( t ) from (4) with respect to δ and μ are
υ 11 = μ log ( t μ ) ( 1 + e ψ ( t ; δ , μ ) ) 1 { 2 ( 1 + μ t ) δ 1 + κ ( t ; δ , μ ) [ ( 1 + μ t ) 2 δ 1 + ( 1 + μ t ) 3 δ 1 + δ log ( t μ ) ( 1 + μ t ) 2 δ 1 { ( 1 + μ t ) δ κ ( t ; δ , μ ) + 1 } ] + δ log ( t μ ) { ( 1 + μ t ) δ 1 + ( 1 + μ t ) 3 δ 1 κ ( t ; δ , μ ) [ κ ( t ; δ , μ ) + 2 ( 1 + μ t ) δ 1 ] } } ,
υ 22 = δ { ( ( 1 + μ t ) δ 1 + t μ [ ( δ 1 ) ( 1 + μ t ) δ 2 + t δ 2 ( 1 + μ t ) 2 ( δ 1 ) κ ( t ; δ , μ ) ] ) ( 1 + μ t ) δ 1 + ( 1 + μ t ) 2 ( δ 1 ) } + μ δ 2 t 3 ( 1 + e ψ ( t ; δ , μ ) ) 1 { ( δ 1 ) ( δ 2 ) ( 1 + μ t ) δ 3 + δ κ ( t ; δ , μ ) ( 2 ( δ 1 ) ( 1 + μ t ) 2 ( δ 1 ) 1 + δ [ ( 1 + μ t ) 2 ( δ 1 ) + δ 1 κ ( t ; δ , μ ) ( 1 + μ t ) 2 ( δ 1 ) + δ 1 ] ) } + 2 t δ ( δ 1 ) ( 1 + μ t ) δ 2 ,
and
υ 12 = ( ( 1 + μ t ) δ 1 + δ log ( μ t ) ( ( 1 + μ t ) 2 δ 1 κ ( t ; δ , μ ) + ( 1 + μ t ) δ 1 ) ) + μ t ( 1 + e 1 ( 1 + μ t ) δ ) { 1 + δ ( 1 + e 1 ( 1 + μ t ) δ ) [ log ( μ t ) { ( ( 1 + μ t ) 2 ( δ 1 ) ( 2 δ 1 ) δ ( 1 + μ t ) δ 1 κ ( t ; δ , μ ) ( ( 1 + μ t ) 2 δ 1 ( 1 + μ t ) 2 δ 1 κ ( t ; δ , μ ) ) ) + ( δ 1 ) ( 1 + μ t ) δ 2 } + ( 1 + μ t ) δ 1 ( ( 1 + μ t ) δ κ ( t ; δ , μ ) + 1 ) ( 1 + μ t ) 1 + ( ( 1 + μ t ) δ 1 + δ log ( μ t ) ( ( 1 + μ t ) 2 δ 1 κ ( t ; δ , μ ) + ( 1 + μ t ) δ 1 ) ) ( 1 + μ t ) δ 1 κ ( t ; δ , μ ) ] } ,
where κ ( t ; δ , μ ) = e ψ ( t ; δ , μ ) ( 1 + e ψ ( t ; δ , μ ) ) 1 .

References

  1. Nadarajah, S.; Haghighi, F. An extension of the exponential distribution. Statistics 2011, 45, 543–558. [Google Scholar] [CrossRef]
  2. Khan, M.N.; Saboor, A.; Cordeiro, G.M.; Nazir, M.; Pescim, R.R. A weighted Nadarajah and Haghighi distribution. Upb Sci. Bull. Ser. Appl. Math. Phys. 2018, 80, 133–140. [Google Scholar]
  3. Mohie El-Din, M.M.; Abu-Youssef, S.E.; Ali, N.S.; Abd El-Raheem, A.M. Classical and Bayesian inference on progressive-stress accelerated life testing for the extension of the exponential distribution under progressive type-II censoring. Qual. Reliab. Eng. Int. 2017, 33, 2483–2496. [Google Scholar] [CrossRef]
  4. Ashour, S.K.; El-Sheikh, A.A.; Elshahhat, A. Inferences and optimal censoring schemes for progressively first-failure censored Nadarajah-Haghighi distribution. Sankhya A 2020, 84, 885–923. [Google Scholar] [CrossRef]
  5. Wu, M.; Gui, W. Estimation and prediction for Nadarajah-Haghighi distribution under progressive type-II censoring. Symmetry 2021, 13, 999. [Google Scholar] [CrossRef]
  6. Elshahhat, A.; Alotaibi, R.; Nassar, M. Inferences for Nadarajah–Haghighi Parameters via Type-II Adaptive Progressive Hybrid Censoring with Applications. Mathematics 2022, 10, 3775. [Google Scholar] [CrossRef]
  7. Elshahhat, A.; Abu El Azm, W.S. Statistical reliability analysis of electronic devices using generalized progressively hybrid censoring plan. Qual. Reliab. Eng. Int. 2022, 38, 1112–1130. [Google Scholar] [CrossRef]
  8. Dey, S.; Wang, L.; Nassar, M. Inference on Nadarajah–Haghighi distribution with constant stress partially accelerated life tests under progressive type-II censoring. J. Appl. Stat. 2022, 49, 2891–2912. [Google Scholar] [CrossRef]
  9. Almarashi, A.M.; Algarni, A.; Okasha, H.; Nassar, M. On reliability estimation of Nadarajah–Haghighi distribution under adaptive type-I progressive hybrid censoring scheme. Qual. Reliab. Eng. Int. 2022, 38, 817–833. [Google Scholar] [CrossRef]
  10. Henningsen, A.; Toomet, O. ‘maxLik’: A package for maximum likelihood estimation in R. Comput. Stat. 2011, 26, 443–458. [Google Scholar] [CrossRef]
  11. Greene, W.H. Econometric Analysis, 7th ed.; Pearson Prentice-Hall: Upper Saddle River, NJ, USA, 2012. [Google Scholar]
  12. Meeker, W.Q.; Escobar, L.A. Statistical Methods for Reliability Data; John Wiley and Sons: New York, NY, USA, 1998. [Google Scholar]
  13. Kundu, D. Bayesian inference and life testing plan for the Weibull distribution in presence of progressive censoring. Technometrics 2008, 50, 144–154. [Google Scholar] [CrossRef]
  14. Lindley, D.V. Approximate Bayesian methods. Trab Estad. 1980, 31, 223–245. [Google Scholar] [CrossRef]
  15. Gelman, A.; Carlin, J.B.; Stern, H.S.; Dunson, D.B.; Vehtari, A.; Rubin, D.B. Bayesian Data Analysis, 2nd ed.Chapman and Hall/CRC: Boca Raton, FL, USA, 2004. [Google Scholar]
  16. Plummer, M.; Best, N.; Cowles, K.; Vines, K. ‘coda’: Convergence diagnosis and output analysis for MCMC. R News 2006, 6, 7–11. [Google Scholar]
  17. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods, and Applications; Springer Science and Business Media, Birkhäuser: Boston, MA, USA, 2000. [Google Scholar]
  18. Ng, H.K.T.; Chan, P.S.; Balakrishnan, N. Optimal progressive censoring plans for the Weibull distribution. Technometrics 2004, 46, 470–481. [Google Scholar] [CrossRef]
  19. Pradhan, B.; Kundu, D. Inference and optimal censoring schemes for progressively censored Birnbaum-Saunders distribution. J. Stat. Plan. Inference 2013, 143, 1098–1108. [Google Scholar] [CrossRef]
  20. Elshahhat, A.; Nassar, M. Bayesian survival analysis for adaptive Type-II progressive hybrid censored Hjorth data. Comput. Stat. 2021, 36, 1965–1990. [Google Scholar] [CrossRef]
  21. Mann, S.P. Introductory Statistics; Wiley Publisher: Hoboken, NJ, USA, 2016. [Google Scholar]
  22. Alotaibi, R.; Rezk, H.; Elshahhat, A. Computational Analysis for Fréchet Parameters of Life from Generalized Type-II Progressive Hybrid Censored Data with Applications in Physics and Engineering. Symmetry 2023, 15, 348. [Google Scholar] [CrossRef]
  23. Wang, F.K. A new model with bathtub-shaped failure rate using an additive Burr XII distribution. Reliab. Eng. Syst. Saf. 2000, 70, 305–312. [Google Scholar] [CrossRef]
  24. Lawless, J.F. Statistical Models and Methods for Lifetime Data, 2nd ed.Wiley Publisher: Hoboken, NJ, USA, 2003. [Google Scholar]
  25. Mohammed, H.S.; Nassar, M.; Alotaibi, R.; Elshahhat, A. Analysis of Adaptive Progressive Type-II Hybrid Censored Dagum Data with Applications. Symmetry 2022, 14, 2146. [Google Scholar] [CrossRef]
  26. Gupta, R.D.; Kundu, D. A new class of weighted exponential distributions. Statistics 2009, 43, 621–634. [Google Scholar] [CrossRef]
  27. Bakouch, H.S.; El-Bar, A.; Ahmed, M.T. A new weighted Gompertz distribution with applications to reliability data. Appl. Math. 2017, 62, 269–296. [Google Scholar] [CrossRef]
  28. Asgharzadeh, A.; Bakouch, H.S.; Nadarajah, S.; Sharafi, F. A new weighted Lindley distribution with application. Braz. J. Probab. Stat. 2016, 30, 1–27. [Google Scholar] [CrossRef]
Figure 1. Density (left) and hazard rate (right) shapes of the WNH distribution.
Figure 1. Density (left) and hazard rate (right) shapes of the WNH distribution.
Symmetry 15 00735 g001
Figure 2. Conditional PDFs of δ and μ .
Figure 2. Conditional PDFs of δ and μ .
Symmetry 15 00735 g002
Figure 3. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for δ .
Figure 3. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for δ .
Symmetry 15 00735 g003
Figure 4. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for μ .
Figure 4. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for μ .
Symmetry 15 00735 g004
Figure 5. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for R ( t ) .
Figure 5. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for R ( t ) .
Symmetry 15 00735 g005
Figure 6. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for h ( t ) .
Figure 6. Heatmaps for the simulation results of MAB (a), RMSE (b), ACL (c), and CP (d) for h ( t ) .
Symmetry 15 00735 g006
Figure 7. The PP plots of the WNH distribution and its competitive models under three real data sets.
Figure 7. The PP plots of the WNH distribution and its competitive models under three real data sets.
Symmetry 15 00735 g007
Figure 8. The estimated PDFs (left) and estimated RFs (right) of the competing models under the three real data sets.
Figure 8. The estimated PDFs (left) and estimated RFs (right) of the competing models under the three real data sets.
Symmetry 15 00735 g008
Figure 9. TTT-Transform plots of the WNH distribution under the three real data sets.
Figure 9. TTT-Transform plots of the WNH distribution under the three real data sets.
Symmetry 15 00735 g009
Figure 10. Contour plots of the log-likelihood function of δ and μ from the three data sets.
Figure 10. Contour plots of the log-likelihood function of δ and μ from the three data sets.
Symmetry 15 00735 g010
Figure 11. Histograms (top) and Trace (bottom) plots of δ , μ , R ( t ) , and h ( t ) from the real data sets.
Figure 11. Histograms (top) and Trace (bottom) plots of δ , μ , R ( t ) , and h ( t ) from the real data sets.
Symmetry 15 00735 g011
Table 1. The point estimation results of δ and μ .
Table 1. The point estimation results of δ and μ .
Par. ( n , m ) SchemeMLELindleyMCMC
Prior → 1212
MABRMSEMABRMSEMABRMSEMABRMSEMABRMSE
δ (50,25)S10.05700.08320.04160.05710.03090.04890.03330.05610.02550.0417
S20.06400.10250.04910.08940.04010.07450.03570.05830.02610.0442
S30.11070.59800.07660.09920.06300.08410.04560.06440.03950.0630
(50,40)S10.04230.05750.03550.04810.02840.04090.02630.04320.01920.0380
S20.04260.05920.03600.04920.02900.04110.02650.04730.02130.0387
S30.04780.06870.03990.05590.03030.05000.03330.05270.02930.0412
(100,50)S10.03670.04930.03230.04330.02740.03650.02720.03880.02090.0354
S20.03810.05250.03420.04670.03030.04060.03070.04090.02160.0397
S30.05180.07400.04600.06310.03410.04370.03990.04550.02580.0427
(100,80)S10.02810.03640.02590.03350.02350.03040.02280.03210.01890.0300
S20.02760.03570.02550.03840.02340.03010.02860.03600.02070.0309
S30.03200.04190.02950.04180.02640.03450.03170.04140.02120.0327
(150,75)S10.02760.03620.02530.03690.02290.03000.02180.03370.01290.0302
S20.02990.03920.02790.04030.02590.03370.02260.03640.01370.0314
S30.04050.05430.03780.05070.03220.04210.02910.04230.01490.0366
(150,120)S10.02210.02820.02090.02670.01970.02520.02000.02600.01180.0242
S20.02240.02920.02130.02770.02010.02620.02010.02680.01330.0244
S30.02400.03130.02280.02980.02130.02770.02080.02850.01370.0267
μ (50,25)S10.06280.09500.05860.09650.05310.09310.04100.06050.03840.0542
S20.05110.07410.04770.07870.03930.06180.02550.05710.03280.0513
S30.07160.17100.23770.18400.13350.14710.05610.07080.05550.0641
(50,40)S10.04700.06800.04400.07690.03540.05950.03620.05450.02510.0484
S20.04490.06550.04330.06860.03520.05970.02530.05200.02040.0434
S30.05070.07370.04740.08050.03720.06070.04030.06270.03080.0542
(100,50)S10.04340.06180.04150.05250.03050.03870.03470.05070.02870.0336
S20.03500.04680.03540.04730.02900.03530.02460.04450.02830.0317
S30.04290.05910.03960.05540.03320.04720.03910.05130.03020.0405
(100,80)S10.03400.04710.03370.04590.02930.03750.03340.04300.02260.0326
S20.03160.04280.03200.04330.02850.03490.02610.04020.02180.0305
S30.03540.04750.03510.04680.03020.03820.03400.04430.02830.0375
(150,75)S10.03360.04580.03370.04560.02910.03730.02570.03290.02070.0266
S20.02790.03670.02820.03750.02580.03350.02150.03240.01960.0249
S30.03760.05190.03610.04870.02890.03750.03390.03310.02680.0286
(150,120)S10.02750.03600.02780.03690.02570.03340.02340.02950.01290.0232
S20.02540.03310.02570.03410.02400.03130.02050.02360.01260.0215
S30.02760.03680.02770.03740.02540.03340.02440.03180.01390.0261
Table 2. The point estimation results of R ( t ) and h ( t ) .
Table 2. The point estimation results of R ( t ) and h ( t ) .
Par. ( n , m ) SchemeMLELindleyMCMC
Prior → 1212
MABRMSEMABRMSEMABRMSEMABRMSEMABRMSE
R ( t ) (50,25)S10.02350.03010.01640.02550.01520.02240.01590.02370.01880.0217
S20.01880.02400.01480.01890.01220.01580.01320.01640.01060.0122
S30.02190.02870.01520.02060.01470.02010.01490.01960.01450.0167
(50,40)S10.01920.02450.01560.01900.01200.01520.01530.01780.01170.0139
S20.01730.02220.01470.01830.01180.01510.01290.01600.01020.0120
S30.01860.02360.01500.01880.01190.01490.01460.01630.01060.0126
(100,50)S10.01710.02180.01480.01840.01210.01480.01330.01620.01070.0121
S20.01330.01690.01220.01540.01070.01210.00870.01110.00750.0109
S30.01500.01900.01240.01550.01080.01360.00900.01220.00900.0117
(100,80)S10.01350.01740.01240.01590.01110.01380.01220.01330.01020.0128
S20.01260.01610.01170.01490.00980.01340.00820.01070.00710.0088
S30.01360.01700.01230.01540.01080.01330.01210.01130.00870.0106
(150,75)S10.01350.01730.01240.01580.01100.01370.01130.01530.00350.0047
S20.01070.01350.01010.01280.00940.01180.00770.00970.00260.0046
S30.01320.01680.01170.01470.01010.01250.00870.01050.00400.0047
(150,120)S10.01120.01420.01070.01350.01000.01250.01050.01180.00350.0045
S20.01010.01270.00970.01220.00910.01140.00510.00670.00210.0037
S30.01090.01370.01030.01300.00950.01200.00830.01020.00260.0043
h ( t ) (50,25)S10.00380.00470.00260.00420.00240.00350.00270.00310.00270.0031
S20.00310.00380.00240.00300.00200.00250.00230.00300.00190.0018
S30.00340.00420.00250.00320.00220.00310.00250.00310.00210.0025
(50,40)S10.00300.00380.00240.00300.00190.00240.00260.00300.00190.0023
S20.00280.00350.00230.00290.00190.00240.00230.00260.00180.0015
S30.00310.00390.00240.00300.00190.00240.00240.00200.00200.0021
(100,50)S10.00280.00350.00230.00290.00200.00240.00220.00270.00180.0025
S20.00220.00270.00190.00240.00160.00200.00170.00220.00150.0016
S30.00240.00300.00190.00240.00180.00220.00190.00230.00160.0023
(100,80)S10.00220.00280.00200.00250.00180.00220.00200.00240.00160.0021
S20.00210.00260.00190.00240.00170.00210.00150.00200.00140.0015
S30.00220.00270.00190.00240.00170.00210.00170.00220.00150.0016
(150,75)S10.00220.00280.00200.00250.00180.00220.00170.00210.00130.0017
S20.00170.00220.00160.00200.00150.00190.00140.00180.00100.0013
S30.00210.00260.00180.00230.00160.00200.00160.00190.00120.0015
(150,120)S10.00190.00230.00170.00220.00160.00200.00150.00200.00080.0012
S20.00160.00210.00150.00190.00150.00180.00130.00150.00060.0012
S30.00180.00220.00160.00210.00150.00190.00150.00160.00050.0011
Table 3. The interval estimation results of δ and μ .
Table 3. The interval estimation results of δ and μ .
( n , m ) SchemeACI-NABoot-pBCI-P1BCI-P2ACI-NABoot-pBCI-P1BCI-P2
ACI-NLBoot-tHPD-P1HPD-P2ACI-NLBoot-tHPD-P1HPD-P2
Par. δ μ
ACLCPACLCPACLCPACLCPACLCPACLCPACLCPACLCP
(50,25)S10.2690.9671.2420.9360.1650.9640.1190.9780.3380.8481.4110.9150.1730.9190.1620.919
0.2780.9490.4940.9390.1550.9870.1080.9800.4810.9470.4080.9200.1580.9540.1460.968
S20.2920.9721.5830.9080.1750.9520.1350.9750.2600.8450.8060.9190.1700.9200.1370.919
0.3040.9400.6240.9610.1640.9820.1250.9740.3330.9280.2810.9140.1600.9600.1310.962
S30.7560.9651.8540.8990.1880.9540.1460.9670.3760.8331.4320.9130.2510.9100.1630.916
0.7860.9780.9560.9840.1770.9900.1330.9890.5660.9560.4070.9160.2240.9460.1470.956
(50,40)S10.2030.9750.6160.9280.1490.9640.1110.9820.2430.8550.9160.9180.1690.9250.1290.925
0.2060.9530.2630.9700.1360.9880.1040.9830.3030.9440.3000.9270.1580.9650.1170.973
S20.2070.9730.6490.9200.1600.9610.1170.9780.2340.8730.7700.9200.1450.9350.1140.940
0.2110.9550.2810.9720.1490.9890.1100.9810.2860.9460.2690.9240.1320.9750.1040.979
S30.2300.9700.8250.9400.1720.9670.1680.9810.2650.8600.9560.9150.1780.9260.1350.929
0.2360.9530.3470.9710.1610.9890.1520.9830.3350.9400.3090.9190.1680.9690.1220.973
(100,50)S10.1740.9570.4430.9260.0960.9530.0910.9670.2180.8780.6770.9290.1540.9420.1360.946
0.1760.9470.2000.9660.0860.9850.0860.9790.2580.9380.2450.9310.1400.9750.1320.975
S20.1810.9610.4820.9310.1630.9580.1140.9720.1760.8940.5190.9300.1330.9510.1100.959
0.1830.9430.2190.9630.1470.9830.1070.9740.1980.9450.2020.9380.1240.9820.1040.987
S30.2580.9661.0340.9740.3170.9820.1650.9870.2270.8910.7570.9270.1660.9480.1570.956
0.2660.9690.4200.9780.3000.9870.1510.9960.2750.9570.2660.9300.1520.9740.1460.976
(100,80)S10.1340.9620.3160.9400.0860.9630.0840.9750.1690.8980.4760.9320.1330.9540.1160.963
0.1360.9560.1480.9720.0830.9880.0800.9850.1880.9490.1880.9390.1260.9840.1040.985
S20.1360.9630.3190.9480.1320.9680.1130.9780.1590.9060.4130.9400.1230.9630.1030.972
0.1370.9600.1500.9760.1180.9920.1070.9880.1750.9410.1700.9440.1160.9800.1010.981
S30.1510.9540.3820.9510.2070.9670.1960.9720.1740.8970.5120.9300.1540.9530.1210.962
0.1520.9440.1780.9640.1900.9840.1800.9760.1960.9450.2020.9360.1400.9810.1110.979
(150,75)S10.1370.9600.3310.9410.1290.9630.0930.9740.1740.9070.4770.9340.1320.9600.0970.970
0.1390.9550.1550.9710.1160.9870.0820.9850.1950.9530.1890.9420.1240.9790.0850.982
S20.1430.9610.3520.9360.1310.9610.0980.9730.1390.9110.3510.9410.1300.9660.0980.976
0.1450.9510.1660.9670.1190.9840.0860.9820.1510.9420.1490.9450.1220.9840.0850.985
S30.1990.9650.5860.9640.1330.9740.0990.9830.1860.8820.5480.9330.1350.9470.1020.951
0.2020.9590.2620.9750.1220.9910.0870.9880.2120.9430.2110.9370.1250.9800.0890.980
(150,120)S10.1070.9580.2400.9490.1040.9660.0930.9740.1360.9130.3430.9390.1280.9660.0950.977
0.1080.9540.1160.9710.0980.9890.0810.9840.1460.9470.1470.9530.1150.9810.0830.986
S20.1100.9520.2430.9410.1100.9590.0960.9670.1280.9240.3070.9450.1160.9750.0940.987
0.1100.9500.1170.9680.0990.9870.0840.9810.1360.9480.1330.9590.1110.9840.0830.988
S30.1200.9610.2650.9400.1340.9630.0990.9740.1410.9180.3610.9360.1300.9670.0950.980
0.1210.9510.1270.9680.1200.9840.0860.9820.1520.9490.1520.9380.1190.9840.0830.984
Table 4. The interval estimation results of R ( t ) and h ( t ) .
Table 4. The interval estimation results of R ( t ) and h ( t ) .
( n , m ) SchemeACI-NABoot-pBCI-P1BCI-P2ACI-NABoot-pBCI-P1BCI-P2
ACI-NLBoot-tHPD-P1HPD-P2ACI-NLBoot-tHPD-P1HPD-P2
Par. δ μ
ACLCPACLCPACLCPACLCPACLCPACLCPACLCPACLCP
(50,25)S10.1150.8870.2430.8600.0610.9000.0450.9090.0180.8910.0330.8360.0110.9110.0080.923
0.1150.8850.0980.8630.0610.9020.0450.9100.0200.9360.0150.8480.0110.9210.0080.930
S20.0900.8930.1860.8750.0500.9100.0360.9190.0140.8960.0290.8640.0110.9280.0070.938
0.0900.8910.0810.8860.0470.9150.0340.9220.0150.9190.0130.8780.0100.9320.0060.953
S30.1060.8910.2110.8710.0580.9070.0420.9170.0160.8920.0290.8380.0100.9120.0070.926
0.1060.8920.0860.8670.0550.9090.0410.9150.0180.9210.0130.8410.0100.9290.0070.937
(50,40)S10.0930.8950.1950.9000.0520.9240.0410.9340.0150.9020.0270.8550.0100.9270.0080.939
0.0930.8940.0840.8950.0510.9210.0400.9300.0160.9350.0120.8640.0100.9440.0080.952
S20.0860.9070.1760.9290.0350.9450.0310.9520.0140.9150.0260.8620.0070.9370.0060.950
0.0860.9050.0780.9340.0340.9470.0310.9560.0150.9380.0120.8730.0060.9500.0060.956
S30.0910.9030.1890.9180.0450.9300.0350.9410.0150.9040.0270.8590.0070.9300.0070.943
0.0910.9030.0810.9030.0440.9380.0340.9440.0160.9250.0120.8630.0070.9410.0070.952
(100,50)S10.0820.9060.1740.9110.0560.9360.0450.9450.0130.9070.0250.8570.0100.9310.0080.943
0.0820.9080.0760.9160.0550.9390.0430.9490.0140.9320.0120.8860.0100.9550.0080.961
S20.0650.9260.1350.9270.0440.9540.0270.9640.0100.9270.0200.8990.0090.9630.0070.976
0.0650.9250.0620.9330.0420.9570.0270.9660.0110.9430.0090.9030.0090.9690.0060.980
S30.0760.9210.1530.9250.0500.9480.0390.9550.0120.9240.0210.8650.0100.9440.0080.956
0.0760.9180.0690.9290.0490.9510.0370.9630.0120.9420.0100.8770.0090.9600.0080.968
(100,80)S10.0660.9210.1380.9200.0400.9480.0370.9590.0110.9210.0210.9040.0080.9620.0050.955
0.0660.9190.0640.9280.0390.9510.0360.9600.0110.9450.0100.9070.0080.9650.0050.976
S20.0610.9210.1250.9420.0320.9590.0190.9690.0100.9240.0190.9040.0060.9640.0040.974
0.0610.9220.0580.9360.0310.9570.0190.9660.0100.9410.0090.9110.0050.9700.0030.982
S30.0650.9230.1390.9370.0400.9530.0220.9610.0100.9270.0210.9000.0080.9630.0050.965
0.0650.9210.0640.9310.0390.9560.0220.9650.0110.9430.0100.8960.0070.9680.0050.980
(150,75)S10.0680.9260.1430.9340.0340.9580.0140.9670.0110.9310.0210.8940.0070.9620.0040.961
0.0680.9250.0640.9400.0340.9600.0130.9700.0110.9470.0100.8960.0070.9690.0030.966
S20.0520.9260.1080.9510.0270.9670.0120.9760.0090.9310.0170.9210.0050.9730.0030.980
0.0520.9260.0510.9550.0270.9690.0110.9780.0090.9470.0080.9280.0050.9790.0030.984
S30.0630.9140.1300.9420.0280.9560.0130.9640.0100.9130.0190.9010.0060.9570.0030.972
0.0630.9160.0600.9490.0280.9600.0120.9700.0100.9350.0090.9060.0060.9710.0030.982
(150,120)S10.0540.9250.1100.9410.0340.9610.0140.9710.0090.9280.0170.9110.0070.9680.0030.982
0.0540.9250.0520.9490.0310.9650.0130.9750.0090.9420.0080.9180.0070.9770.0030.984
S20.0500.9420.1030.9600.0230.9770.0110.9820.0080.9430.0160.9290.0050.9830.0030.989
0.0500.9410.0490.9660.0230.9800.0100.9840.0080.9490.0080.9350.0040.9870.0030.990
S30.0540.9360.1110.9510.0250.9720.0110.9820.0090.9390.0170.9160.0050.9790.0030.985
0.0540.9360.0530.9570.0250.9760.0100.9830.0090.9450.0080.9210.0050.9840.0030.988
Table 5. Optimum criteria of PCTII fashions.
Table 5. Optimum criteria of PCTII fashions.
CriterionGoalFormula
1 Maximize Trace { I ( δ ^ , μ ^ ) } Maximize { ( L 11 + L 22 ) }
2 Minimize Trace { I 1 ( δ ^ , μ ^ ) } Minimize { σ ^ 11 + σ ^ 22 }
3 Minimize Det { I 1 ( δ ^ , μ ^ ) } Minimize { σ ^ 11 σ ^ 22 σ ^ 12 σ ^ 21 }
4 Minimize v a r ^ { log ( T ^ ϱ ) } Minimize v a r ^ log μ 1 1 log 1 ϱ 1 + ϱ 1 δ 1
Table 6. The failure times of three real-life data sets.
Table 6. The failure times of three real-life data sets.
DataTimes
A12344566899991012
121313131415161617172020222326
273133484850515268
B5112131467598122145165196224245293321
330350420
C1.45.16.310.812.118.519.722.223.030.637.346.353.959.866.2
Table 7. Parameter estimates (SEs) and summary fit for the three data sets.
Table 7. Parameter estimates (SEs) and summary fit for the three data sets.
DataModel δ ^ μ ^ NLAICAIBIHQIKS(p-Value)
AWNH0.9853(0.3123)0.0734(0.0415)154.062312.125312.458315.452313.3180.0967(0.859)
WE0.9438(0.5143)0.0978(0.0820)155.563315.126315.459318.453316.3190.1351(0.475)
WG1.9629(0.8424)0.0379(0.0441)183.312370.623370.957373.950371.8170.1534(0.318)
WL1.0323(0.7126)0.1168(0.0886)245.940495.880496.213499.207497.0730.1534(0.318)
BWNH1.2847(0.5563)0.0054(0.0035)110.007224.015224.815225.796224.2600.1157(0.947)
WE0.9656(0.4975)0.0080(0.0044)112.999229.997230.797231.778230.2430.1843(0.516)
WG1.5689(0.9971)0.0055(0.0037)126.963255.926256.726257.707256.1710.2807(0.095)
WL1.5683(0.9768)0.0141(0.1129)160.595325.191325.991326.971325.4360.2015(0.404)
CWNH1.6947(1.6346)0.0224(0.0309)63.9266131.853132.853133.269131.3830.1284(0.939)
WE1.7791(1.4654)0.0514(0.0694)64.6491133.298134.298134.714133.2830.1230(0.956)
WG1.5548(1.5448)0.0288(0.0536)68.6265141.253142.253142.669141.2380.2000(0.522)
WL1.8443(1.5011)0.0969(0.0972)106.630217.261218.261218.677217.2460.1996(0.524)
Table 8. Artificial PCTII samples from the real data sets.
Table 8. Artificial PCTII samples from the real data sets.
Data ( n , m ) SampleSchemePCTII Data
A(39,13)A1 ( 26 , 0 * 12 ) 1, 2, 4, 5, 6, 9, 10, 12, 12, 13, 16, 22, 27
A2 ( 0 * 6 , 26 , 0 * 6 ) 1, 2, 3, 4, 4, 5, 6, 6, 9, 12, 15, 26, 33
A3 ( 0 * 12 , 26 ) 1, 2, 3, 4, 4, 5, 6, 6, 8, 9, 9, 9, 9
B(18,10)B1 ( 8 , 0 * 9 ) 5, 11, 21, 31, 46, 98, 122, 165, 224, 293
B2 ( 0 * 4 , 4 , 4 , 0 * 4 ) 5, 11, 21, 31, 46, 75, 98, 122, 165, 245
B3 ( 0 * 9 , 8 ) 5, 11, 21, 31, 46, 75, 98, 122, 145, 165
C(15,8)C1 ( 7 , 0 * 7 ) 1.4, 5.1, 10.8, 12.1, 19.7, 23.0, 37.3, 46.3
C2 ( 0 * 4 , 7 , 0 * 3 ) 1.4, 5.1, 6.3, 10.8, 12.1, 19.7, 22.2, 37.3
C3 ( 0 * 7 , 7 ) 1.4, 5.1, 6.3, 10.8, 12.1, 18.5, 19.7, 22.2
Table 9. Point estimates (SEs) of δ , μ , R ( t ) , and h ( t ) from the three data sets.
Table 9. Point estimates (SEs) of δ , μ , R ( t ) , and h ( t ) from the three data sets.
DataSampleParameterMLELindleyMCMC
AA1 δ 0.5897(0.3913)0.63260.5887(0.0098)
μ 0.0960(0.1079)0.10520.0944(0.0095)
R ( 5 ) 0.8706(0.0480)0.88420.8730(0.0119)
h ( 5 ) 0.0272(0.0087)0.02380.0267(0.0027)
A2 δ 0.5248(0.2854)0.57210.5237(0.0099)
μ 0.1278(0.1221)0.15790.1264(0.0099)
R ( 5 ) 0.8531(0.0483)0.87930.8549(0.0188)
h ( 5 ) 0.0304(0.0089)0.02440.0300(0.0036)
A3 δ 0.3515(0.1614)0.40140.3503(0.0099)
μ 0.2031(0.1882)0.27760.2017(0.0100)
R ( 5 ) 0.8612(0.0467)0.87400.8625(0.0106)
h ( 5 ) 0.0258(0.0072)0.02530.0255(0.0022)
BB1 δ 0.7141(0.4439)0.85260.7138(0.0050)
μ 0.0240(0.0308)0.03270.0234(0.0042)
R ( 5 ) 0.9579(0.0276)0.96990.9589(0.0886)
h ( 5 ) 0.0086(0.0056)0.00720.0084(0.0189)
B2 δ 1.6471(1.9045)1.91701.6471(0.0010)
μ 0.0064(0.0103)0.00410.0063(0.0008)
R ( 5 ) 0.9735(0.0128)0.97660.9738(0.1032)
h ( 5 ) 0.0055(0.0027)0.00600.0054(0.0218)
B3 δ 0.4845(0.3724)0.45160.4845(0.0010)
μ 0.0256(0.0391)0.02740.0256(0.0010)
R ( 5 ) 0.9699(0.0221)0.97040.9700(0.0994)
h ( 5 ) 0.0060(0.0042)0.00700.0060(0.0212)
CC1 δ 0.6951(0.6325)0.75480.6950(0.0010)
μ 0.0614(0.0973)0.06850.0613(0.0010)
R ( 5 ) 0.8981(0.0599)0.89540.8982(0.0276)
h ( 5 ) 0.0217(0.0116)0.02290.0216(0.0056)
C2 δ 1.1384(1.3816)1.14711.1384(0.0010)
μ 0.0421(0.0741)0.04260.0421(0.0010)
R ( 5 ) 0.8791(0.0599)0.87200.8792(0.0090)
h ( 5 ) 0.0276(0.0130)0.02940.0276(0.0008)
C3 δ 1.2189(3.1401)1.22811.2189(0.0010)
μ 0.0376(0.1245)0.03300.0376(0.0010)
R ( 5 ) 0.8836(0.0693)0.88240.8837(0.0031)
h ( 5 ) 0.0266(0.0136)0.02590.0266(0.0008)
Table 10. Interval estimates of δ , μ , R ( t ) , and h ( t ) from the three real data sets.
Table 10. Interval estimates of δ , μ , R ( t ) , and h ( t ) from the three real data sets.
DataSamplePar.ACI-NABoot-pBCI
ACI-NLBoot-tHPD
LowerUpperLengthLowerUpperLengthLowerUpperLength
AA1 δ 0.00001.35671.35670.25611.09810.84200.56970.60790.0382
0.16062.16502.00430.25551.09670.84120.56900.60710.0381
μ 0.00000.30760.30760.04170.17880.13710.07600.11270.0366
0.01060.86940.85880.04160.17860.13700.07580.11230.0366
R ( 5 ) 0.77660.96470.18820.81450.92150.10700.85030.89610.0458
0.78150.97000.18850.81270.92030.10760.85030.89600.0458
h ( 5 ) 0.01010.04430.03420.02180.05070.02890.02160.03190.0103
0.01450.05100.03650.02150.05090.02940.02170.03190.0102
A2 δ 0.00001.08411.08410.22790.97710.74930.50440.54260.0382
0.18071.52371.34290.22740.97580.74850.50400.54210.0381
μ 0.00000.36710.36710.05550.23800.18250.10730.14560.0383
0.01970.83100.81130.05580.23740.18160.10740.14570.0383
R ( 5 ) 0.75840.94780.18940.78140.88860.10720.83450.87550.0409
0.76340.95320.18980.78050.88740.10690.83390.87460.0407
h ( 5 ) 0.01300.04780.03480.02320.05160.02840.02550.03450.0090
0.01720.05390.03670.02290.05100.02810.02570.03460.0090
A3 δ 0.03520.66790.63270.15260.65460.50190.33150.36960.0381
0.14300.86450.72160.15230.65370.50140.33170.36960.0379
μ 0.00000.57190.57190.08820.37810.34310.18230.22120.0388
0.03301.24841.21540.08800.37760.28960.18310.22180.0387
R ( 5 ) 0.76960.95280.18320.84390.92340.07950.84890.87570.0268
0.77430.95790.18360.84320.92200.07880.84860.87540.0267
h ( 5 ) 0.01160.04000.02830.01710.04260.02550.02300.02830.0054
0.01490.04470.02980.01780.04250.02470.02270.02810.0053
BB1δ0.00001.58401.58400.31001.32961.01950.70410.72360.0195
0.21122.41472.20350.30941.32781.01850.70420.72370.0195
μ 0.00000.08440.08440.01040.04470.03430.01560.03160.0160
0.00190.29750.29550.01080.04470.03390.01560.03150.0159
R ( 5 ) 0.90380.99780.09400.86750.92150.05400.94470.97240.0277
0.90530.99350.08820.86560.91930.05370.94490.97240.0276
h ( 5 ) 0.00120.01950.01830.00180.01960.01780.00550.01140.0059
0.00250.03050.02800.00170.01940.01770.00560.01140.0058
B2 δ 0.00005.38005.38001.01524.12253.10731.64521.64910.0040
0.170815.88415.7131.01304.11713.10421.64511.64910.0039
μ 0.00000.02650.02650.00280.01990.01710.00470.00820.0035
0.00030.14950.14930.00260.01940.01680.00470.00800.0033
R ( 5 ) 0.94850.99850.05000.95440.99890.04450.96660.98050.0139
0.94880.99880.05000.95320.99720.04400.96680.98050.0137
h ( 5 ) 0.00030.01070.01040.00130.01060.00930.00390.00690.0030
0.00210.01420.01200.00140.01040.00900.00400.00690.0029
B3 δ 0.00541.21431.20890.21040.90210.69170.48230.48640.0041
0.10742.18532.07790.20990.90090.69100.48240.48640.0040
μ 0.00000.10230.10230.00780.08260.07480.02370.02760.0041
0.00130.50990.50860.00760.08210.07450.02360.02750.0039
R ( 5 ) 0.92660.99340.06680.94390.99230.04840.96780.97220.0044
0.92750.99170.06420.94120.98890.04770.96790.97230.0044
h ( 5 ) 0.00020.01420.01400.00170.01260.01090.00540.00640.0010
0.00150.02350.02200.00160.01240.01080.00550.00640.0009
CC1 δ 0.00001.93461.93460.30181.29420.99240.69310.69710.0040
0.11684.13574.01890.30111.29250.99140.69310.69700.0039
μ 0.00000.25210.25210.00260.19290.19030.05920.06330.0041
0.00271.37341.37070.00240.19220.18980.05940.06330.0039
R ( 5 ) 0.78080.98410.20330.81750.95150.13400.89500.90140.0063
0.78810.98970.20160.81670.95030.13360.89510.90130.0062
h ( 5 ) 0.00100.04430.04330.00240.02760.02520.02090.02230.0014
0.00760.06160.05400.00230.02710.02480.02090.02230.0014
C2 δ 0.00003.84643.84640.45153.12252.67101.13631.14030.0040
0.105512.28512.1790.45053.11842.66791.13651.14040.0039
μ 0.00000.18740.18740.01830.12990.11160.03980.04400.0042
0.00131.32601.32460.01810.12920.11110.04010.04400.0039
R ( 5 ) 0.76170.99660.23490.83440.97860.14420.87350.88490.0114
0.76920.98480.21560.83300.97660.14360.87370.88490.0112
h ( 5 ) 0.00200.05310.05110.00130.02740.02610.02610.02900.0029
0.01090.06960.05870.00130.02740.02600.02610.02900.0029
C3 δ 0.00007.37347.37340.02574.02483.99911.21681.22080.0041
0.0078190.04190.030.02573.99403.96831.21691.22090.0040
μ 0.00000.28160.28160.01780.17860.16080.03550.03960.0041
0.000524.56824.5680.01770.17800.16030.03570.03960.0039
R ( 5 ) 0.74780.98970.24190.78190.96280.18090.87740.88970.0124
0.75770.98670.22900.78250.96080.17830.87750.88970.0122
h ( 5 ) 0.00160.05320.05160.02170.06460.04290.02480.02820.0033
0.00980.07220.06250.02150.06350.04200.02500.02820.0031
Table 11. Optimum PCTII plans from the real data sets.
Table 11. Optimum PCTII plans from the real data sets.
DataSampleCriterion
1234
ϱ 0.30.60.9
AA11403.9330.1647720.00011716.70746148.93994784.368
A2789.15720.0963500.00012227.44207805.426662029.04
A3436.62850.0614550.00014113.49689143.38244387.135
BB117524.800.1979760.000008405.23051260.4217284.561
B2477786.50.6273940.000012907.68686271.593321362.3
B311544.300.1401860.000011337.3565742.78142647.492
CC12303.0250.4094670.00017853.18096198.70514162.369
C26946.6003.8757240.00057323.53281161.66423427.930
C36278.9671.9144030.00027622.3980360.40023632.9782
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

Elshahhat, A.; Mohammed, H.S. Survival Analysis and Applications of Weighted NH Parameters Using Progressively Censored Data. Symmetry 2023, 15, 735. https://doi.org/10.3390/sym15030735

AMA Style

Elshahhat A, Mohammed HS. Survival Analysis and Applications of Weighted NH Parameters Using Progressively Censored Data. Symmetry. 2023; 15(3):735. https://doi.org/10.3390/sym15030735

Chicago/Turabian Style

Elshahhat, Ahmed, and Heba S. Mohammed. 2023. "Survival Analysis and Applications of Weighted NH Parameters Using Progressively Censored Data" Symmetry 15, no. 3: 735. https://doi.org/10.3390/sym15030735

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