Next Article in Journal
Study on the Influence of an Under-Crossing Parallel Double-Line Shield Tunnel on the Existing Tunnel Structure
Next Article in Special Issue
Sensitivity of Survival Analysis Metrics
Previous Article in Journal
Enhancing Decomposition Approach for Solving Multi-Objective Dynamic Non-Linear Programming Problems Involving Fuzziness
Previous Article in Special Issue
Weighted Competing Risks Quantile Regression Models and Variable Selection
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Survival Analysis of the PRC Model from Adaptive Progressively Hybrid Type-II Censoring and Its Engineering Applications

by
Ahmed Elshahhat
1,*,
Osama E. Abo-Kasem
2 and
Heba S. Mohammed
3
1
Faculty of Technology and Development, Zagazig University, Zagazig 44519, Egypt
2
Department of Statistics, Faculty of Commerce, Zagazig University, Zagazig 44519, Egypt
3
Department of Mathematical Sciences, College of Science, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(14), 3124; https://doi.org/10.3390/math11143124
Submission received: 13 June 2023 / Revised: 7 July 2023 / Accepted: 13 July 2023 / Published: 14 July 2023
(This article belongs to the Special Issue Statistical Methods and Models for Survival Data Analysis)

Abstract

:
A new two-parameter statistical model, obtained by compounding the generalized-exponential and exponential distributions, called the PRC lifetime model, is explored in this paper. This model can be easily linked to other well-known six-lifetime models; namely the exponential, log-logistic, Burr, Pareto and generalized Pareto models. Adaptive progressively hybrid Type-II censored strategy, used to increase the efficiency of statistical inferential results and save the total duration of a test, has become widely used in various sectors such as medicine, biology, engineering, etc. Via maximum likelihood and Bayes inferential methodologies, given the presence of such censored data, the challenge of estimating the unknown parameters and some reliability time features, such as reliability and failure rate functions, of the PRC model is examined. The Markov-Chain Monte Carlo sampler, when the model parameters are assumed to have independent gamma density priors, is utilized to produce the Bayes’ infer under the symmetric (squared-error) loss of all unknown subjects. Asymptotic confidence intervals as well as the highest posterior density intervals of the unknown parameters and the unknown reliability indices are also created. An extensive Monte Carlo simulation is implemented to investigate the accuracy of the acquired point and interval estimators. Four various optimality criteria, to select the best progressive censored design, are used. To demonstrate the applicability and feasibility of the proposed model in a real-world scenario, two data sets from the engineering sector; one based on industrial devices and the other on aircraft windshield, are analyzed. Numerical evaluations showed that the PRC model furnishes a superior fit compared to seven other models in the literature, including: alpha-power exponential, log-logistic, Nadarajah–Haghighi, generalized-exponential, Weibull, gamma and exponential lifetime distributions. The findings demonstrate that, in order to obtain the necessary estimators, the Bayes’ paradigm via Metropolis–Hastings sampler is recommended compared to its competitive likelihood approach.

1. Introduction

Compounding of several discrete and lifetime distributions is a common strategy for producing new distributions in lifetime modelling. Via compounding the generalized exponential and exponential distributions, Popović et al. [1] introduced a new two-parameter lifetime model, called the generalized exponential distribution. They also derived various mathematical and reliability properties of this model.
Henceforth, due to the fact that several extended exponential models have been proposed in the literature, we will call this model PRC, which is an abbreviation of its authors’ names. Suppose Y, a random variable of lifetime of a test item(s), follows PRC ( σ ) , where σ = ( μ , γ ) . Then, the probability density function (PDF) f ( · ) and cumulative distribution function (CDF) F ( · ) of Y, are given by
f ( y ; σ ) = μ e μ y γ ( 1 e μ y ) 1 γ 1 log ( 1 e μ y ) 2 , y > 0 , μ , γ > 0 ,
and
F ( y ; σ ) = 1 γ 1 log ( 1 e μ y ) 1 ,
respectively. Usually, a reliability practitioner uses two reliability time indices, namely: reliability function (RF) R ( · ) and hazard rate function (HRF) h ( · ) , from any lifespan model as key aspects for assessing capability.
Thus, besides μ and γ , we investigate both R ( · ) and h ( · ) of the PRC model as unknown time parameters, for t > 0 , are given by
R ( t ; σ ) = 1 1 γ 1 log ( 1 e μ t ) 1 , t > 0 ,
and
h ( t ; σ ) = μ e μ t ( 1 e μ t ) ( log ( 1 e μ t ) ) 1 γ 1 log ( 1 e μ t ) 1 ,
respectively. Figure 1 depicts the PDF and HRF shapes of the PRC ( σ ) distribution based on several options of its parameters. It exhibits that the density (1) function has a unimodal, decreasing, or increasing shape, while the failure rate (4) function has a bathtub-shaped or decreasing shape.
From (1), Popović et al. [1] showed also that the PRC distribution could be investigated as a new extended model of five widely known distributions in the literature, namely: log-logistic, generalized Pareto, Pareto, Burr and exponential models. Taking Y ( μ ) = log 1 e μ Y , the PRC model can be easily linked to:
  • Exponential distribution with scale parameter one if setting X = log 1 γ 1 Y ( μ ) .
  • Log-logistic distribution with shape β ( = 1 ) parameter and scale λ ( = μ ) parameter if setting X = Y ( μ ) .
  • Burr distribution with shape ( β 1 , β 2 ) = ( 2 , 1 ) parameters if setting X = γ 1 Y ( μ ) .
  • Pareto distribution with shape β ( = 1 ) parameter if setting X = 1 γ 1 Y ( μ ) .
  • Generalized Pareto distribution without location parameter, shape β ( = 1 ) parameter and scale λ ( = μ ) parameter if setting X = Y ( μ ) .
In the reliability sector, as a result of time constraints or a lack of money, life-testing investigations typically conclude before all of the items fail. The observations resulting from these circumstances are known as the censored sample.
The adaptive progressive hybrid Type-II (AT2-PH) censored method, proposed by Ng et al. [2], has received a lot of attention from numerous authors since it allows for the construction of quite effective statistical analyses. According to this mechanism, on a test at time zero, n independent experimental items units are considered. The threshold point T > 0 , the progressive censoring design R = ( R 1 , R 2 , , R m ) , and the number of failures m ( < n ) are predetermined. Specifically, the fixed values of some of progressive stages R i for i = 1 , 2 , , m , may be reassigned consequently during the experiment. It should be mentioned here that T is an ideal full test time due to the fact that it allows the testing to pass the preassigned time T. At the moment of the first failure (say Y 1 : m : n ) recorded, R 1 surviving items are randomly selected and extracted from the staying n 1 items. Again, at the moment of the second failure Y 2 : m : n recorded, R 2 items of the staying n R 1 2 items are randomly selected and withdrawn out of the test and so on. If Y m : m : n < T , the test proceeds using the same R i s and stops at Y m : m : n . This situation is similar to the same traditional Type-II progressive (T2-P) censored strategy. Otherwise, if T < Y m : m : n such as Y d : m : n < T < Y d + 1 : m : n , we reset the removal scenario by setting R i = 0 for i = d + 1 , , m 1 , where d is the number of failed items recorded before time T. Then, the test ends at mth failure occur and all staying surviving subjects, i.e., R m = n m i = 1 d R i are removed. However, let y = { ( Y 1 : m : n , R 1 ) , , ( Y d : m : n , R d ) , ( Y d + 1 : m : n , 0 ) , , ( Y m : m : n , R m ) } be an AT2-PH censored sample from a continuous population, with PDF f ( · ) and CDF F ( · ) , hence, the joint PDF for the AT2-PH failure times of size m can be expressed as
L ( σ | y ) = K i = 1 m f ( y i ; σ ) i = 1 d [ R ( y i ; σ ) ] R i [ R ( y m ; σ ) ] R m * ,
where K is a constant. The most significant benefit of this technique is that it allows us to obtain the effective number of failures m while also ensuring that the overall test duration does not deviate too much from the pre-specified period T. Various investigations based on AT2-PH have been conducted; readers can refer for example to Nassar et al. [3], Panahi and Moradi [4], Elshahhat and Nassar [5], Panahi and Asadi [6], Alotaibi et al. [7], Alotaibi et al. [8], Ateya et al. [9] and references cited therein.
The originality of this work stems from the fact that, in the presence of incomplete sampling, it is the first time since the PRC distribution’s establishment that two maximum likelihood and Bayesian methods for its parameters of life have been compared. The impetus for this work stems from the relevance of the suggested censoring mechanism in boosting the efficacy of statistical inference when compared to Type-I (time) and Type-II (failure) censoring mechanisms. However, we are motivated to perform this work for two reasons: (i) the HRF of the PRC distribution exhibits a bathtub-shaped or decreasing shape, which is a preferable occurrence in many practical domains; (ii) the usefulness of AT2-PH censoring plan is that it provides flexibility in terminating experiments at a predetermined time and lowering total test length while retaining the desirable qualities of progressive design in practical reliability investigations. As far as we know, there is no treatment of inferred aspects of the RRC distribution, especially in reliability practice. The goal of this work, using AT2-PH censoring plan, is to address this gap by proving that the PRC lifespan model may be employed as a survival model. Thus, the present study has five objectives, which can be viewed as:
  • The issue of estimating the model parameters ( μ , γ ) as well as the reliability indices ( R ( t ) , h ( t ) ) of the PRC distribution using maximum likelihood and Bayes estimation techniques from adaptive progressively hybrid Type-II censoring is addressed.
  • Via independent gamma conjugate density priors, the Bayes estimates against the squared-error loss (SEL) of μ , γ , R ( t ) and h ( t ) are evaluated through the Monte Carlo Markov-chain (MCMC) techniques.
  • Approximate confidence interval (ACI) and highest posterior density (HPD) interval estimates of μ , γ , R ( t ) and h ( t ) are also obtained.
  • As anticipated, the analytical solutions of α or μ developed by the proposed estimation approaches cannot be represented in closed expressions, so in R programming software, ‘ maxLik ’ (by Henningsen and Toomet [10]) and ‘ coda ’ (by Plummer et al. [11]) packages are recommended to evaluate the offered estimates.
  • Determine the most efficient progressively designed solution that delivers a substantial amount of information on the model parameter(s) of interest, depending on four distinct optimality criteria.
  • Through a series of numerical comparisons, we assess the efficiency of the offered estimate in terms of root-mean-squared error, average relative absolute bias, average interval length and coverage probability simulated values.
  • Two real-world engineering-data-sets-based applications demonstrate the PRC distribution’s capacity to fit various data types and adapt the offered methodologies to actual practical circumstances.
The remaining sections of the paper are: Non-Bayesian and Bayesian estimations are provided in Section 2 and Section 3, respectively. In Section 4, Monte Carlo outcomes are presented. Engineering applications are examined in Section 5. In Section 6, criteria for specifying the best progressive censoring are depicted. Lastly, the study concludes in Section 7.

2. Likelihood Estimation

The likelihood estimation method is used to evaluate the parameters of a probability distribution by maximizing the likelihood (or log-likelihood) function in order to make the observed data most probable for the statistical model. In general, the Expectation-Maximization (EM) approach is used to discover the local maximum likelihood parameters of a statistical model when latent variables are present or data are missing or incomplete; see Zhang et al. [12]. For brevity, the main advantages of the EM algorithm are: (i) it is guaranteed that the likelihood will increase on each iteration; (ii) the E-Step and M-Step are both easy for many problems; and (iii) the M-Step solution is often found in closed form. Unfortunately, the main drawbacks of the EM algorithm are: (i) its convergence is very slow; (ii) it achieves convergence to the local optima only; and (iii) it requires both forward and backward probabilities.
However, the likelihood methodology is considered in this section to provide the point ML and ACI estimates of μ , γ , R ( t ) and h ( t ) . Suppose y , is an AT2-PH censored sample taken from the PRC population whose PDF (1) and CDF (2) are given. Substituting (1) and (2) into (5), where y i is used, for instance, in place of y i : m : n , the likelihood Function (5) becomes
L ( σ | y ) γ 1 μ e μ y ¯ m 1 ξ 1 ( y m ; μ , γ ) R m × i = 1 d 1 ξ 1 ( y i ; μ , γ ) R i i = 1 m ( 1 e μ y i ) 1 ξ 2 ( y i ; μ , γ ) ,
where ξ ( y i ; μ , γ ) = 1 γ 1 log 1 e μ y i for i = 1 , 2 , , m .
As a result, the natural logarithm of (6), symbolized by ( · ) log L ( · ) , is
σ y m log γ 1 μ e μ y ¯ i = 1 m log ( 1 e μ y i ) 2 i = 1 m log ξ ( y i ; μ , γ ) + i = 1 d R i log 1 ξ 1 ( y i ; μ , γ ) + R m log 1 ξ 1 ( y m ; μ , γ ) .
Consequently the MLEs ( μ ^ , γ ^ ) of ( μ , γ ) , from (7), can be derived directly by solution of the following normal formulas
μ = m ( μ 1 y ¯ ) i = 1 m y i e μ y i ( 1 e μ y i ) 1 2 i = 1 m ξ ( y i ; μ , γ ) ξ 1 ( y i ; μ , γ ) + i = 1 d R i ξ ( y i ; μ , γ ) ξ 2 ( y i ; μ , γ ) ( 1 ξ 1 ( y i ; μ , γ ) ) 1 + R m ξ ( y m ; μ , γ ) ξ 2 ( y m ; μ , γ ) ( 1 ξ 1 ( y m ; μ , γ ) ) 1
and
γ = m γ 1 2 i = 1 m ξ ( y i ; μ , γ ) ξ 1 ( y i ; μ , γ ) + i = 1 d R i ξ ( y i ; μ , γ ) ξ 2 ( y i ; μ , γ ) ( 1 ξ 1 ( y i ; μ , γ ) ) 1 + R m ξ ( y m ; μ , γ ) ξ 2 ( y m ; μ , γ ) ( 1 ξ 1 ( y m ; μ , γ ) ) 1 ,
where ξ ( y i ; μ , γ ) = γ 1 y i e μ y i 1 e μ y i 1 and ξ ( y i ; μ , γ ) = γ 2 log 1 e μ y i .
It is obvious, from (8) and (9), that there are no closed-form expressions for the offered MLEs μ ^ or γ ^ . Given an AT2-PH censored sample, we recommend employing the Newton–Raphson technique by installing the ‘ maxLik ’ package in R 4.2.2 software ( R Foundation for Statistical Computing, Vienna, Austria) to evaluate the acquired μ ^ and γ ^ .
Once the values of μ ^ and γ ^ obtained, via replacing ( μ , γ ) by ( μ ^ , γ ^ ) , we obtain the MLEs of R ( t ) (3) and h ( t ) (4), at a distinct t > 0 , respectively, as follows:
R ^ ( t ) = 1 1 γ ^ 1 log ( 1 e μ ^ t ) 1 ,
and
h ^ ( t ) = μ ^ e μ ^ t ( 1 e μ ^ t ) ( log ( 1 e μ ^ t ) ) 1 γ ^ 1 log ( 1 e μ ^ t ) 1 .
To build the 100 ( 1 α ) % ACIs of μ , γ , R ( t ) and h ( t ) , the asymptotic behavior of each MLE is used. Following a large sample theory, the asymptotic distribution of the MLEs ( μ ^ , γ ^ ) is the normal distribution with mean ( σ ) and 2 × 2 variance/covariance I 1 ( σ ) matrix. Following Lawless [13], using the observed Fisher information, we utilize I 1 ( σ ^ ) to estimate I 1 ( σ ) as
I 1 ( σ ^ ) = I 11 I 12 I 21 I 22 ( μ , γ ) = ( μ ^ , γ ^ ) 1 = ν ^ 11 ν ^ 12 ν ^ 21 ν ^ 22 ,
where, in Appendix A, the Fisher elements I i j for i , j = 1 , 2 of (10) are reported.
Hence, the 100 ( 1 α ) % ACIs of μ and γ can be obtained as follows:
μ ^ z α 2 ν ^ 11 and γ ^ z α 2 ν ^ 22 ,
respectively, where z α 2 is the ( α 2 ) th percentile of standard normal distribution.
Additionally, to create the 100 ( 1 α ) % ACIs of R ( t ) and h ( t ) , the variances ν ^ R and ν ^ h of the MLEs R ^ ( t ) and h ^ ( t ) of R ( t ) and h ( t ) , respectively, must be derived first. To achieve this aim, the delta technique, which is a generic strategy for estimating confidence intervals for any unknown parametric function, is considered. Following Greene [14], we estimate ν ^ R and ν ^ h of R ^ ( t ) and h ^ ( t ) as
ν ^ R [ M R I 1 ( σ ) M R ] ( μ ^ , γ ^ ) and ν ^ h [ M h I 1 ( σ ) M h ] ( μ ^ , γ ^ ) ,
respectively, where M R = ( R μ R γ ) and M h = ( h μ h γ ) .
As a result, the respective ACIs of R ( t ) and h ( t ) at a confidence level 100 ( 1 α ) are obtained, respectively, by
R ^ ( t ) z α 2 ν ^ R and h ^ ( t ) z α 2 ν ^ h .

3. Bayes’ Estimation

The Bayes’ paradigm of statistical analysis has received substantial attention in recent decades as an efficient and useful alternative to the conventional ML approach. This section deals with deriving the Bayes estimators of μ , γ , R ( t ) and h ( t ) in addition to their HPD intervals. To achieve this goal, we assume that the PRC parameters μ and γ are stochastically independent and have gamma PDFs G a m m a ( θ 1 , β 1 ) and G a m m a ( θ 2 , β 2 ) , respectively.
Thus, the joint prior PDF (say ω ( · ) ) of μ and γ is
ω ( μ , γ ) μ θ 1 1 γ θ 2 1 e ( β 1 μ + β 2 γ ) , μ , γ > 0 ,
where θ j , β j > 0 , j = 1 , 2 are the hyperparameters.
Loss functions play an important part in the Bayesian framework because they may be used to detect the overestimation and underestimating of a research of interest. We thus consider the most important loss, called the SEL function (say l ( · ) ), which is given by
l τ , τ ˜ = τ ˜ τ 2 ,
where τ ˜ denotes the symmetric Bayes estimates of τ . Using (13), the Bayes estimator τ ˜ is provided directly by the posterior mean of τ .

3.1. Posterior Distribution

Substituting (6) with (12) into the continuous Bayes’ theorem, the posterior PDF (say Ω ( · ) ) of μ and γ can be written as
Ω σ y = G 1 μ m + θ 1 1 e μ β 1 + y ¯ γ θ 2 m 1 e β 2 γ 1 ξ 1 ( y m ; μ , γ ) R m × i = 1 d 1 ξ 1 ( y i ; μ , γ ) R i i = 1 m ( 1 e μ y i ) 1 ξ 2 ( y i ; μ , γ ) ,
where G = 0 0 L ( σ | y ) × ω ( μ , γ ) d μ d γ .
It is obvious, from (14), that the Bayes estimators μ ˜ and γ ˜ of μ and γ , respectively, cannot be formulated in closed form. As a result, the Bayes estimate may be approximated using the MCMC approach. To adopt this technique, the respective full conditional distributions of μ and γ take the following forms
Ω 1 μ γ , y μ m + θ 1 1 e μ β 1 + y ¯ 1 ξ 1 ( y m ; μ , γ ) R m i = 1 d 1 ξ 1 ( y i ; μ , γ ) R i i = 1 m ( 1 e μ y i ) 1 ξ 2 ( y i ; μ , γ ) ,
and
Ω 2 γ μ , y γ θ 2 m 1 e β 2 γ 1 ξ 1 ( y m ; μ , γ ) R m i = 1 d 1 ξ 1 ( y i ; μ , γ ) R i i = 1 m ξ 2 ( y i ; μ , γ ) .
As a result, the full conditional PDFs of μ and γ cannot be reduced to any known density. So, acquiring samples straightforward of μ and γ from (14) and (15), respectively, is unattainable. Plotting a random AT2-PH censored sample taken from the PRC ( 0.5 , 0.2 ) population, when ( n , m , T ) = ( 80 , 40 , 5 ) and R is uniform censoring, Figure 2 indicates that the full conditional PDFs of μ and γ are similar to the normal density. Subsequently, we consider the normal PDF as the proposal model to derive the Bayes estimates as well as to derive the HPD intervals of μ , γ , R ( t ) or h ( t ) ; see Gelman et al. [15] and Lynch [16] for further details.

3.2. Metropolis–Hastings Algorithm

To offer the Bayes point and HPD interval estimates, we propose to employ the following Metropolis–Hastings steps:
Step 1: 
Set the starting points of ( μ , γ ) say ( μ ( 0 ) , γ ( 0 ) ) = ( μ ^ , γ ^ ) .
Step 2: 
Set k = 1 .
Step 3: 
Obtain μ * and γ * from (14) and (15) using the normal distributions N ( μ ^ , ν ^ 11 ) and N ( γ ^ , ν ^ 22 ) , respectively.
Step 4: 
Obtain Ψ i , i = 1 , 2 as
Ψ 1 = Ω 1 μ * γ ( k 1 ) , y Ω 1 μ ( k 1 ) γ ( k 1 ) , y and Ψ 2 = Ω 2 γ * μ ( k ) , y Ω 2 γ ( k 1 ) μ ( k ) , y ,
Step 5: 
Obtain u i , i = 1 , 2 from a uniform U ( 0 , 1 ) distribution.
Step 6: 
If u 1 min { 1 , Ψ 1 } , set μ ( k ) = μ * ; else, set μ ( k ) = μ ( k 1 ) .
Step 7: 
If u 2 min { 1 , Ψ 2 } , set γ ( k ) = γ * ; else, set γ ( k ) = γ ( k 1 ) .
Step 8: 
Obtain R ( k ) ( t ) and h ( k ) ( t ) for t > 0 , respectively, as
R ( k ) ( t ) = 1 1 γ ( k ) 1 log ( 1 e μ ( k ) t ) 1 ,
and
h ( k ) ( t ) = μ ( k ) e μ ( k ) t ( 1 e μ ( k ) t ) ( log ( 1 e μ ( k ) t ) ) 1 γ ( k ) 1 log ( 1 e μ ( k ) t ) 1 .
Step 9: 
Set k = k + 1 .
Step 10: 
Repeat Steps 2–8 B times and ignore the first B draws as ‘burn-in’ to obtain
μ ( B + 1 ) , γ ( B + 1 ) , R ( B + 1 ) ( t ) , h ( B + 1 ) ( t ) , , μ ( B ) , γ ( B ) , R ( B ) ( t ) , h ( B ) ( t ) .
Step 11: 
Get the Bayes estimates of μ , γ , R ( t ) , or h ( t ) (say τ ) as
τ ˜ = 1 B B k = B + 1 M τ ( k ) ,
Step 12: 
Create the 100 ( 1 α ) % HPD interval of μ , γ , R ( t ) or h ( t ) (say τ ) as:
(a)
Sort the MCMC draws of τ ( k ) for k = B + 1 , B + 2 , , B as
τ ( B + 1 ) , τ ( B + 2 ) , , τ ( B ) .
(b)
Obtain the 100 ( 1 α ) % two-sided HPD interval of τ as
( τ k * , τ k * + 1 α B B ) ,
where k * = B + 1 , B + 2 , , M is chosen such that
τ k * + [ 1 α B B ] τ k * = min 1 k α ( B B ) ( τ ( k + [ 1 γ B B ] ) τ ( k ) ) .
where the largest integer less than or equal to e is denoted by [ e ] , Chen and Shao [17].

4. Numerical Comparisons

To assess the validity of the obtained point and interval estimations of the PRC’s parameters and reliability features such as RF R ( t ) and HRF h ( t ) , for t > 0 , presented in the preceding sections, based on various scenarios, this section demonstrates extensive Monte Carlo simulations.

4.1. Simulation Scenarios

This subsection exhibits several censoring designs, in turn, to examine the behavior of acquired estimators, based on the proposed censored data. To do this, based on various choices of T (threshold point), n (total number of items on the test), m (size of censored sample) and R (progressive design), we replicate 1000 AT2-PH samples from different sets of PRC ( μ , γ ) , namely; Set-A: PRC (0.2, 0.5) and Set-B: PRC (1.5, 1). Fixing t = 0.1 , the plausible values of ( R ( t ) , h ( t ) ) are taken as (0.8869, 0.2854) and (0.6634, 1.5825) using Set-A and -B, respectively. Taking T ( = 2 , 5 ) and n ( = 40 , 80 ) , we obtain m as a failure percentage (FP%) for each n such as m n × 100 % (=50, 80)%. Clearly, the test is terminated when FP% achieves the pre-specified 50% (or 80%). Further, three different scenarios of progressive patterns R are used, namely:
Scheme - 1 : R = n m , 0 * ( m 1 ) , Scheme - 2 : R = 0 * 0.5 m 1 , n m , 0 * 0.5 m , and Scheme - 3 : R = 0 * ( m 1 ) , n m ,
where 0 * ( m 1 ) represents 0 is repeated m 1 times.
Specifically, to obtain a AT2-PH censored sample of size m from the PRC model, after assigning T, n and R , do the following procedure:
Step 1: 
Set the specified actual values of PRC ( μ , γ ) .
Step 2: 
Obtain an ordinary T2-P censored sample as:
(a)
Simulate κ independent items (say κ 1 , κ 2 , , κ m ) from uniform U ( 0 , 1 ) distribution.
(b)
Set ϱ i = κ i i + l = m i + 1 m R l 1 , for i = 1 , 2 , , m .
(c)
Set U i = 1 ϱ m ϱ m 1 ϱ m i + 1 for i = 1 , 2 , , m .
(d)
Collect a T2-P censored sample (with size m) from PRC ( μ , γ ) distribution by setting
Y i = μ 1 log 1 e γ ( 1 u i 1 ) , i = 1 , 2 , , m .
Step 3: 
Determine d at predetermined T.
Step 4: 
Discard the remaining sample ( Y d + 2 , , Y m ) when T < Y m .
Step 5: 
Use the truncated f ( y ; μ , γ ) 1 F ( y d + 1 ; μ , γ ) 1 distribution to obtain ( Y d + 2 , , Y m ) order statistics with size n d l = 1 d R l 1 .
To do the required computations of the acquired Bayes point (or interval) estimates, according to Kundu [18], we assign different priors to the values of the hyperparameters ( θ 1 , θ 2 , β 1 , β 2 ) as follows:
  • For Set-A: Pr.1:(1,2.5,5,5) and Pr.2:(2,5,10,10);
  • For Set-B: Pr.1:(7.5,5,5,5) and Pr.2:(15,10,10,10).
Following the Bayes’ MCMC technique described in Section 3, to remove of the affection for the starting choice, the first 2000 (from the total 12,000 MCMC iterations) are removed. One the target 1000 AT2-PH samples collected, by installing ‘ maxLik ’ and ‘ coda ’ packages (by Henningsen and Toomet [10] and Plummer et al. [11], respectively) in R 4.2.2 software, the maximum likelihood and Bayes MCMC estimates (in addition to their 95% ACI and HPD interval estimates) of μ , γ , R ( t ) and h ( t ) are evaluated.
Though MCMC is a powerful technique for evaluating complex models and posterior distributions in statistical programming, it needs thorough verification and adjustment to assure its validity and efficiency. Before proceeding, in Bayes’ history, various plotting tools for visualizing the convergence of MCMC drawn from the posterior density of the parameters of a studied model are suggested. There are two MCMC diagnostics, namely:
(i)
Brooks–Gelman–Rubin (BGR): It measures a chain’s convergence by comparing the variances within and between chains.
(ii)
Trace (with thinning): It shows the sampled values per chain and node throughout iterations.
The thinning procedure entails selecting separate points from the sample at each predetermined phase in order to produce an independent sample; see for example Riabiz et al. [19]. Here, we take each 5th point (for instance) for thinning procedure. For instance, from both given sets A and B, we draw 12,000 samples from the joint posterior PDF of μ and γ using Scheme-1, Pr.1, ( n , m ) = ( 40 , 20 ) and T = 2 (as an example); see Figure 3. It highlights that the MCMC draws displayed in each chain are sufficiently mixed, and the length of the draws removed at the start of each chain (burn) is sufficient to reduce autocorrelation. It also shows there are no differences among the variances within and between chains. As a result, the acquired estimations are reliable.
Practically, the average estimates (Av.Es) for the theoretical point estimators of μ , γ , R ( t ) and h ( t ) (say ρ ) derived from frequentist (or Bayes) approach are given by
ρ ˇ ¯ q = 1 1000 i = 1 1000 ρ ˇ q ( i ) ,
where ρ ˇ ( i ) is the computed estimate of ρ at ith sample, ρ 1 = μ , ρ 2 = γ , ρ 3 = R ( t ) and ρ 4 = h ( t ) .
Comparing the acquired point estimates of ρ is carried out based on their root mean squared-errors (RMSEs) and average relative absolute biases (ARABs) as
RMSE ( ρ ˇ q ) = 1 1000 i = 1 1000 ρ ˇ ( i ) ρ 2 ,
and
ARAB ( ρ ˇ q ) = 1 1000 i = 1 1000 1 ρ q ρ ˇ q ( i ) ρ q ,
respectively.
On the other hand, the acquired 95% ACI and HPD interval estimates obtained of ρ are compared based on their average confidence lengths (ACLs) and coverage percentages (CPs) as
ACL ( 1 α ) % ( ρ ) = 1 1000 i = 1 1000 U ρ ˇ q ( i ) L ρ ˇ q ( i ) ,
and
CP ( 1 α ) % ( ρ ) = 1 1000 i = 1 1000 J L ρ ˇ q ( i ) ; U ρ ˇ q ( i ) ρ ,
respectively, where J ( · ) refers the indicator operator and ( L ( · ) , U ( · ) ) refers the (lower-limit, upper-limit) of the interval estimator.

4.2. Simulation Results

A heatmap (or heat-map) is a popular tool for data visualization. It is described as a visual representation of data in which individual values in a matrix are shown by colors. So, via a heatmap tool, this subsection presents the simulated RMSE, ARAB, ACL and CP values of μ , γ , R ( t ) and h ( t ) in Figure 4, Figure 5, Figure 6 and Figure 7. To better distinguish the methods, in Figure 4, Figure 5, Figure 6 and Figure 7, the proposed estimation methods are defined on the ‘ x - axis ’ line while the proposed censored test (referred by “T-n[FP%]-Scheme”) is defined on the ‘ y - axis ’ line. For each elegant heatmap, the values of the RMSE, ARAB, ACL, or CP are specified using colors from yellow (lowest values) to red (highest values). In each heat-map, for distinguishing, using Pr.1 (as an example), the Bayes estimator using Pr.1 is symbolized by “BE-Pr.1” and the HPD estimator is symbolized by “HPD-Pr.1”. In the Supplementary Materials, the numerical outcomes of μ , γ , R ( t ) and h ( t ) are provided.
From Figure 4, Figure 5, Figure 6 and Figure 7, for the smallest RMSE, ARAB and ACL values as well as the largest values of CP, we draw several comments on the behavior of σ , R ( t ) and h ( t ) as follows:
  • All acquired estimates of μ , γ , R ( t ) and h ( t ) show superior behavior and provide more efficient results.
  • As n (or m) grows, for μ , γ , R ( t ) or h ( t ) , the offered Bayesian and classical estimates perform well. A similar pattern is also noted when the total number of the removal items, i.e., l = 1 m R l , decreases.
  • As T i , i = 1 , 2 increases, it can be noted that
    (i)
    For Set-A: The RMSE, ARAB and ACL values developed from the maximum likelihood (or Bayes’) approach of μ decrease while those of γ , R ( t ) and h ( t ) increase.
    (ii)
    For Set-B: The RMSE, ARAB and ACL values developed from the maximum likelihood (or Bayes’) approach of μ and γ decrease while those of R ( t ) and h ( t ) increase.
    (iii)
    The opposite result of (i)–(ii) is observed when comparing the acquired ACI (or HPD interval) estimators in relation to their CPs.
  • As μ and γ increase, it is clear that
    (i)
    For Set-A (or Set-B): The RMSE and ARAB values created from the maximum likelihood (or Bayes’) estimates of μ , γ , R ( t ) or h ( t ) increase.
    (ii)
    For Set-A (or Set-B): The ACL values created from the asymptotic (or HPD interval) estimates of μ , γ , R ( t ) or h ( t ) increase. The opposite pattern is also recorded if acquired ACI (or HPD interval) estimators are compared in terms of their CPs.
  • Comparing the proposed schemes 1, 2 and 3, it is observed that
    (i)
    For Set-A: The offered estimates of μ , γ and R ( t ) perform satisfactorily via Scheme-1 ‘left-censoring’ whereas those of h ( t ) perform satisfactorily via Scheme-3 ‘right-censoring’ than others.
    (ii)
    For Set-B: The offered estimates of μ and h ( t ) perform satisfactorily via Scheme-3 ‘right-censoring’ whereas those of γ and R ( t ) perform satisfactorily via Scheme-1 ‘left-censoring’ than others.
  • From Set-A (or Set-B), since Pr.2’s variance is smaller than Pr.1’s variance, the Bayes’ evaluations employing Prior-2 produce more reliable estimations than others for all unknown parameters.
  • Due to the fact that the offered Bayes’ estimates involved more priority information compared to the competitive frequentist estimates, the Bayes estimates of μ , γ , R ( t ) or h ( t ) showed better behavior compared to others.
  • As a consequence, in the case of data created using an adaptive Type-II gradually hybrid technique, the Bayes M-H methodology is advised for assessing the PRC parameters or its reliability time features.

5. Engineering Applications

In this section, we will give two real data sets from the engineering sector to analyze the capacity to adapt the research objectives to actual situations and to demonstrate the flexibility and superiority of the suggested lifespan model.

5.1. Industrial Devices Data Analysis

This application provides an analysis of the lifetimes of fifty industrial devices put on life test at time zero, given by Aarset [20] and reanalyzed by Elshahhat and Nassar [5]; see Table 1. The main feature of this type of data is that it presents a bathtub-shaped failure rate.
Just before going any further, we shall fit the PRC model to the complete industrial devices data set along with seven other statistical models in the literature as its competitors, namely: (1) alpha-power exponential ( APE ( μ , γ ) ) by Mahdavi and Kundu [21]; (2) log-logistic ( LoL ( μ , γ ) ) by Johnson et al. [22]; (3) Nadarajah–Haghighi NH ( μ , γ ) by Nadarajah and Haghighi [23]; (4) generalized exponential GE ( μ , γ ) by Gupta and Kundu [24]; (5) Weibull W ( μ , γ ) by Weibull [25]; (6) gamma G ( μ , γ ) by Johnson et al. [22]; and (7) exponential ( E ( μ ) ) by Johnson et al. [22] distributions.
To highlight the best model, besides the Kolmogorov–Smirnov (K–S) distance (along its p-value), several measures of fit are taken into account, namely: negative log-likelihood (NL), Akaike (A), Bayesian (B), consistent Akaike (CA), Hannan–Quinn (HQ), Anderson–Darling (A–D) and Cramér–von Mises (CvM). By installing the ‘ AdequacyModel ’ package, by Marinho et al. [26] in R 4.2.2 software, all suggested criteria are evaluated through the ML approach. So, Table 2 displays the MLEs (along with their standard–errors (St.Es)) of μ and γ , as well as the fitted results of the PRC model and other competing models. It shows that the PRC model fits the industrial devices data adequately and that the PRC model is the best compared to other distributions because it provides the lowest value of all proposed statistics and the highest p-value.
Additionally, based on the full industrial devices data set, the superiority of the PRC model is examined via various graphical plot tools, namely: (a) probability–probability (PP); (b) estimated PDFs; and (c) estimated/empirical RFs; see Figure 8. It supports the same findings provided in Table 2. To show the existence and uniqueness characteristics of the acquired ML estimates of the PRC parameters μ and γ , based on the full industrial devices data set, the contour plot of the natural logarithm of the likelihood function for different options of μ and γ is shown in Figure 8d. It shows that the MLEs μ ^ and γ ^ exist and are unique.
Henceforward, to run any additional calculations based on the industrial devices data, we propose considering these values as suitable starting points. We recommend using the acquired estimates μ ^ 0.0384 and γ ^ 0.2160 as starting points for any future computations based on the data from industrial devices. Moreover, Figure 8e demonstrates the estimated/empirical scaled–TTT transform of the PRC lifetime distribution and indicates that the PRC failure rate model provides a bathtub-shaped failure rate.
Now, to evaluate the offered estimators in the presence of industrial devices data, various AT2-PH samples of size m = 25 , using different choices of T and R , are created and listed in Table 3. For each artificial sample, the maximum likelihood and Bayes’ MCMC estimates (with their St.Es) in addition to the associated 95% two-sided ACI/HPD-interval estimates (with their lengths) of μ , γ , R ( t ) and h ( t ) (at t = 1 ) are obtained; see Table 4.
Since no a priori information about the PRC parameters μ and γ is available, the Bayes’ evaluations are developed based on 50,000 MCMC iterations, the first 10,000 of which are considered as ‘burn-in’. It can be shown that the point and interval estimates obtained using the MCMC method are quite close to those created using the maximum likelihood approach. Remaining 40,000 draws of μ , γ , R ( t ) and h ( t ) are used to evaluate several properties of each unknown parameter including: mean, mode, (1st, 2nd, 3rd) quartiles denoted by Q i , i = 1 , 2 , 3 , standard deviation (St.D) and skewness (Skew.); see Table 5. It appears that the findings of central tendency are relatively near to one other.
Trace and MCMC frequencies (with Gaussian kernel) plots, based on the staying 40,000 MCMC variates of μ , γ , R ( t ) and h ( t ) , are shown in Figure 9. In each MCMC plot, to distinguish, the Bayes estimate (or sample mean) and HPD interval limits are displayed with horizontal solid and dashed lines, respectively. It reveals that the MCMC process converges successfully and that deleting the first 10,000 samples as burn-in is an adequate size to remove the influence of the starting values. It is also evident from the estimates that the acquired estimates of μ , γ , R ( t ) or h ( t ) are fairly symmetric.

5.2. Aircraft Windshield Data Analysis

A big aircraft’s windscreen is a sophisticated piece of equipment made up of multiple layers of material, including a highly tough outer skin with a heated layer just behind it. These objects’ failures are not structural failures and do not cause aircraft damage, but they do require that the windscreen be repaired. This application provides an analysis of the remaining 65 service times of windshields that had not failed at the time of observation. This application provides an analysis of 63 observations of windshields’ service times for a particular windshield model; see Murthy et al. [27]. In Table 6, aircraft windshield service times (AWSTs) are displayed.
As described in Section 5.1, just like our comparison of the PRC model and its competitive (including: APE, LoL, NH, GE, W, G and E) distributions, here we shall calculate the suggested measures of model selection (including: NL, A, B, CA, HQ, A–D, CvM and K–S (p-value)) based on the complete AWST data; see Table 7. It demonstrates that the PRC model gives the lowest values with regard to the NL, A, B, CA, HQ, A–D, CvM and K–S statistics, except for the highest p-value among all fitted lifetime models. As a result, from the AWST data, the PRC model is the best choice. Further, Figure 10 displays the PP, estimated PDFs, estimated/empirical RFs, contour plot of the log-likelihood function and estimated/empirical scaled–TTT transform plots. As we anticipated, Figure 10a–c supports the same facts provided by Table 7. Figure 8d shows that the acquired ML estimates μ ^ 1.2648 and γ ^ 0.0829 existed and are unique. These offered estimates can be easily utilized as starting points for any forthcoming assessments in presence of the complete AWST data. Furthermore, Figure 8e indicates that the AWST data have an increasing failure rate, which is one of the common PRC failure rates.
To obtain the acquired estimates of μ , γ , R ( t ) and h ( t ) , using m = 33 along with different selects of T and R , three artificial AT2-PH samples are collected from the AWST data and provided in Table 8. Using S[i] for i = 1 , 2 , 3 , the frequentist and Bayes estimates (with their St.Es) as well as the 95% ACI and HPD-interval estimates (with their lengths) of μ , γ , R ( t ) and h ( t ) (at t = 0.1 ) are calculated; see Table 9. Recall that the acquired Bayes point and HPD interval estimates, based on 50,000 MCMC iterations (when the first 10,000 iterations are taken as ‘burn-in’), are carried out using non-informative priors. As a result, Table 9 shows that the point (or interval) estimates of μ , γ , R ( t ) or h ( t ) obtained by considered maximum likelihood and Bayes’ estimation methodologies are high similar to each other. This finding is due to the fact that the MCMC calculations are implemented based on improper gamma densities. Again, using S[i] for i = 1 , 2 , 3 created from AWST data, the same properties reported in Table 5 are re-computed and presented in Table 10. The results of the central tendency, reported in Table 10, appear to be quite close to one another and support the same findings listed in Table 5.
Based on S[i], i = 1 , 2 , 3 created from AWST data, Figure 11 displays both trace and MCMC frequencies plots of μ , γ , R ( t ) and h ( t ) based on their staying 40,000 MCMC iterations. It indicates that the simulated MCMC draws of μ , γ , R ( t ) or h ( t ) are mixed sufficiently and behave fairly symmetrically. Also, the facts displayed in Figure 11 support the same facts displayed in Table 10.

6. Optimal Progressive Designs

In the area of reliability, the investigator can choose the best effective censoring technique from a set of all accessible progressive patterns in order to provide the most information on the parameter(s) under research as feasible; see for additional details Ng et al. [28]. Depending on unit capacity, experimental facilities and budgetary restrictions, when n, m and T i , i = 1 , 2 are predetermined, the ideal censoring design R = ( R 1 , R 2 , , R m ) is obtained. However, different criteria have been established in the literature, and several results on the optimal censoring fashions have also been suggested; see for example Pradhan and Kundu [29]; Sen et al. [30]; Ashour et al. [31]; Elshahhat and Abu El Azm [32] and references cited therein. Table 11 lists several metrics to help us determine the best censoring strategy.
From Table 11, it is clear that the proposed optimum Q 1 , Q 2 and Q 4 , criteria aim to minimize the determinant, trace and the variance of the logarithmic-ML estimate of the ζ th quantile of the estimated variance–covariance I 1 ( μ ^ , γ ^ ) matrix, while criterion Q 3 aims to maximize the observed Fisher information I ( μ ^ , γ ^ ) matrix. Obviously, the best progressive pattern must correspond to the highest value of Q 3 and the lowest value of other criteria Q i for i = 1 , 2 , 4 .
Subsequently, from (2), the logarithm of the quantile of the PRC lifetime distribution Q ζ is
log ( Q ^ ζ ) = log ( μ ) + log log 1 e γ ( 1 ζ 1 ) ( μ ^ , γ ^ ) , 0 < ζ < 1 .
Again, performing the delta method, the approximated variance (say v a r ^ ( · ) ) of log ( Q ^ ζ ) is given by
v a r ^ ( log ( Q ^ ζ ) ) = Σ log ( Q ζ ) I 1 ( ϑ ) Σ log ( Q ζ ) ( μ ^ , γ ^ ) ,
where
Σ log ( Q ^ ζ ) = μ log ( Q ζ ) log ( Q ζ ) , γ log ( Q ζ ) ( μ ^ , γ ^ ) .

6.1. Optimum for Industrial Devices

In this part, using the three created samples S[i] for i = 1 , 2 , 3 , which are reported in Table 3, we shall evaluate the optimum criteria in Table 11 in turn to suggest the best progressive censoring plans from the industrial devices data. However, based on the acquired MLEs μ ^ and γ ^ obtained from industrial devices data, the optimum criteria Q i for i = 1 , 2 , 3 , 4 are evaluated; see Table 12.
As a result, Table 12 shows that:
  • According to O i , i = 1 , 3 , 4 ; the progressive design used in S[3] is the optimum censoring than others.
  • According to O 2 ; the progressive design used in S[2] is the optimum censoring than others.

6.2. Optimum for Aircraft Windshield

In this part, based on the generated samples S[i] for i = 1 , 2 , 3 reported from the aircraft windshield data, we shall suggest the optimal progressive censored plan. However, from Table 8 and Table 11, the optimum criteria are evaluated; see Table 13.
Table 13 shows, according to all optimum criteria O i , i = 1 , 2 , 3 , 4 , the progressive design used in sample S[3] is the best compared to others. Additionally, Table 12 and Table 13 support the same recommended censoring mechanisms as set out in Section 4.
Ultimately, we found that the suggested infer techniques, which use both industrial devices and aircraft windshield data sets, give a good illustration of the unknown parameters along with the reliability properties of the PRC lifetime model when an adaptive progressive hybrid Type-II censored sample is created.

7. Concluding Remarks

A new two-parameter lifetime distribution called the PRC lifetime model obtained by compounding the generalized-exponential and exponential distributions, which allows a bathtub-shaped or decreasing failure rate shape, has been explored in the presence of incomplete data collected from an adaptive progressively hybrid Type-II censored mechanism. In this study, using the proposed censoring plan, we took into account both maximum likelihood (non-Bayesian) and Bayesian estimations of the parameter, reliability and hazard functions of the PRC distribution. The asymptotic confidence intervals of each unknown quantity using the asymptotic distribution of the frequentist estimates have been obtained. The Newton–Raphson technique, via ‘ maxLik ’ package, has been utilized in turn to evaluate the point and interval estimates developed by the likelihood method. Squared-error loss as well as independent gamma prior functions have been considered to create the Bayes estimates in addition to their HPD interval estimates. The Metropolis–Hastings approximation technique, via ‘ coda ’ package, has also been utilized to approximate the Bayes estimates and also to construct the associated HPD intervals. Four accuracy metrics in the proposed Monte Carlo simulation experiments, namely: root-mean-squared error, average absolute relative bias, average interval length and coverage probability, have been used to evaluate the effectiveness of the suggested estimation approaches. Using four optimality metrics, the optimum progressive censored designs have been proposed. Two applications using the lifetimes of industrial devices and the remaining service times of windshields have been analyzed to exemplify the offered methodologies. The main results of the applications examined demonstrate that the suggested PRC lifespan model may be employed as an appropriate model when compared to the other seven competing models. As a future study, it may be useful to extend the proposed estimation methodologies to accelerated tests, competing risks, or other censoring designs. We hope that the approaches provided in this paper will be valuable to statisticians, reliability practitioners, clinicians and anyone else who needs to perform this kind of life test.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/math11143124/s1, Table S1: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of μ for Set-A; Table S2: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of γ for Set-A; Table S3: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of R ( t ) for Set-A; Table S4: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of h ( t ) for Set-A; Table S5: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of μ for Set-B; Table S6: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of γ for Set-B; Table S7: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of R ( t ) for Set-B; Table S8: Av.Es (1st column), RMSEs (2nd column) and ARABs (3rd column) of h ( t ) for Set-B; Table S9: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of μ for Set-A; Table S10: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of γ for Set-A; Table S11: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of R ( t ) for Set-A; Table S12: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of h ( t ) for Set-A; Table S13: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of μ for Set-B; Table S14: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of γ for Set-B; Table S15: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of R ( t ) for Set-B; Table S16: The ACLs (1st column) and CPs (2nd column) of 95% ACI/HPD intervals of h ( t ) for Set-B.

Author Contributions

Methodology, A.E., O.E.A.-K. and H.S.M.; Funding acquisition, H.S.M.; Software, A.E.; Supervision O.E.A.-K.; Writing—original draft, A.E., O.E.A.-K. and H.S.M.; Writing—review and editing A.E. and H.S.M. 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 like to express their thanks to the editor and anonymous referees for helpful comments and observations. The authors would like to thank the Princess Nourah bint Abdulrahman University for supporting (PNURSP2023R175) this work.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Differentiating (7) with regard to μ and γ , the Fisher’s items I i j , i , j = 1 , 2 are
I 11 = i = 1 m μ 1 ( ( y i ( ( 2 δ ( y i ; μ , γ ) 1 + μ y i ( ( 1 + 2 δ ( y i ; μ , γ ) 2 2 δ ( y i ; μ , γ ) 1 ) e μ y i ( 1 e μ y i ) 1 + 1 2 δ ( y i ; μ , γ ) 1 ) 1 ) × e μ y i ( 1 e μ y i ) 1 1 ) μ 1 ( 1 + μ y i ( ( 2 δ ( y i ; μ , γ ) 1 1 ) e μ y i ( 1 e μ y i ) 1 1 ) ) ) ) + i = 1 d ( R i y i 2 ( ( η ( y i ; μ , γ ) ( γ { 1 η ( y i ; μ , γ ) 1 } η ( y i ; μ , γ ) 2 ) 2 { 1 η ( y i ; μ , γ ) 1 } ) η ( y i ; μ , γ ) e μ y i × [ γ { 1 η ( y i ; μ , γ ) 1 } ( 1 e μ y i ) η ( y i ; μ , γ ) 2 ] 2 + [ γ { 1 η ( y i ; μ , γ ) 1 } ( 1 e μ y i ) η ( y i ; μ , γ ) 2 ] 1 ) e μ y i ) + R m y m 2 { ( η ( y m ; μ , γ ) ( γ [ 1 η ( y m ; μ , γ ) 1 ] η ( y m ; μ , γ ) 2 ) 2 [ 1 η ( y m ; μ , γ ) 1 ] ) η ( y m ; μ , γ ) e μ y m + [ γ [ 1 η ( y m ; μ , γ ) 1 ] ( 1 e μ y m ) η ( y m ; μ , γ ) 2 ] 1 } e μ y m ,
I 22 = i = 1 m R i γ 1 δ ( y i ; μ , γ ) log ( 1 e μ y i ) ( γ 2 γ 2 δ 2 ( y i ; μ , γ ) [ 1 [ γ 1 δ ( y i ; μ , γ ) ] 1 ] ) 2 { γ 1 δ ( y i ; μ , γ ) × ( 2 γ ( 1 [ γ 1 δ ( y i ; μ , γ ) ] 1 ) + log ( 1 e μ y i ) ( γ 2 δ 2 ( y i ; μ , γ ) ) 1 ) + 2 log ( 1 e μ y i ) ( 1 [ γ 1 δ ( y i ; μ , γ ) ] 1 ) } + γ 1 i = 1 d ( γ 1 ( 1 + 2 log ( 1 e μ y i ) δ 1 ( y i ; μ , γ ) ) + 2 log ( 1 e μ y i ) δ 2 ( y i ; μ , γ ) ) R m γ 1 δ ( y i ; μ , γ ) log ( 1 e μ y m ) ( γ 2 [ 1 [ γ 1 δ ( y m ; μ , γ ) ] 1 ] γ 2 δ 2 ( y m ; μ , γ ) ) 2 ( γ 1 δ ( y i ; μ , γ ) × ( 2 ( γ [ 1 [ γ 1 δ ( y m ; μ , γ ) ] 1 ] ) + log ( 1 e μ y m ) [ γ 2 δ 2 ( y m ; μ , γ ) ] 1 ) + 2 ( [ 1 [ γ 1 δ ( y m ; μ , γ ) ] 1 ] log ( 1 e μ y m ) ) )
and
I 12 = i = 1 m ( ( 2 ( y i e μ y i ( ( 1 e μ y i ) ( γ ( 1 γ 1 log ( 1 e μ y i ) ) ) 2 ) 1 ) ) ) + i = 1 d ( R i y i ( ( ( log ( 1 e μ y i ) ( γ δ ( y i ; μ , γ ) ) 1 1 ) γ δ 1 ( y i ; μ , γ ) + 1 ) δ ( y i ; μ , γ ) + 2 ( ( 1 γ δ 1 ( y i ; μ , γ ) ) γ 1 log ( 1 e μ y i ) ) ) ( 1 e μ y i ) δ ( y i ; μ , γ ) e μ y i ( γ ( 1 γ δ 1 ( y i ; μ , γ ) ) ( 1 e μ y i ) δ 2 ( y i ; μ , γ ) ) 2 ) + R m y m ( 1 e μ y m ) γ 1 δ ( y m ; μ , γ ) e μ y m ( γ ( 1 γ δ 1 ( y i ; μ , γ ) ) ( 1 e μ y m ) γ 2 δ 2 ( y m ; μ , γ ) ) 2 × ( ( ( log ( 1 e μ y m ) ( γ γ 1 δ ( y m ; μ , γ ) ) 1 1 ) γ δ 1 ( y m ; μ , γ ) + 1 ) γ 1 δ ( y m ; μ , γ ) + 2 ( ( 1 γ δ 1 ( y i ; μ , γ ) ) γ 1 log ( 1 e μ y m ) ) ) ,
where η ( y i ; μ , γ ) = γ 1 ( 1 log ( 1 e μ y i ) ) and δ ( y i ; μ , γ ) = γ ( 1 γ 1 log ( 1 e μ y i ) ) .

References

  1. Popović, B.V.; Ristić, M.M.; Cordeiro, G.M. A two-parameter distribution obtained by compounding the generalized exponential and exponential distributions. Mediterr. J. Math. 2016, 13, 2935–2949. [Google Scholar] [CrossRef]
  2. Ng, H.K.T.D.; Chan, P.S. Statistical analysis of exponential lifetimes under an adaptive Type-II progressive censoring scheme. Nav. Res. Logist. 2009, 56, 687–698. [Google Scholar] [CrossRef] [Green Version]
  3. Nassar, M.; Abo-Kasem, O.; Zhang, C.; Dey, S. Analysis of Weibull distribution under adaptive Type-II progressive hybrid censoring scheme. J. Indian Soc. Probab. Stat. 2018, 19, 25–65. [Google Scholar] [CrossRef]
  4. Panahi, H.; Moradi, N. Estimation of the inverted exponentiated Rayleigh distribution based on adaptive Type II progressive hybrid censored sample. J. Comput. Appl. Math. 2020, 364, 112345. [Google Scholar] [CrossRef]
  5. 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]
  6. Panahi, H.; Asadi, S. On adaptive progressive hybrid censored Burr type III distribution: Application to the nano droplet dispersion data. Qual. Technol. Quant. Manag. 2021, 18, 179–201. [Google Scholar] [CrossRef]
  7. Alotaibi, R.; Nassar, M.; Elshahhat, A. Computational Analysis of XLindley Parameters Using Adaptive Type-II Progressive Hybrid Censoring with Applications in Chemical Engineering. Mathematics 2022, 10, 3355. [Google Scholar] [CrossRef]
  8. Alotaibi, R.; Elshahhat, A.; Rezk, H.; Nassar, M. Inferences for Alpha Power Exponential Distribution Using Adaptive Progressively Type-II Hybrid Censored Data with Applications. Symmetry 2022, 14, 651. [Google Scholar] [CrossRef]
  9. Ateya, S.F.; Amein, M.M.; Mohammed, H.S. Prediction under an adaptive progressive type-II censoring scheme for Burr Type-XII distribution. Commun. Stat. Theory Methods 2022, 51, 4029–4041. [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. Plummer, M.; Best, N.; Cowles, K.; Vines, K. ‘coda’: Convergence diagnosis and output analysis for MCMC. R News 2006, 6, 7–11. [Google Scholar]
  12. Zhang, Y.M.; Wang, H.; Bai, Y.; Mao, J.X.; Chang, X.Y.; Wang, L.B. Switching Bayesian dynamic linear model for condition assessment of bridge expansion joints using structural health monitoring data. Mech. Syst. Signal Process. 2021, 160, 107879. [Google Scholar] [CrossRef]
  13. Lawless, J.F. Statistical Models and Methods for Lifetime Data, 2nd ed.; John Wiley and Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  14. Greene, W.H. Econometric Analysis, 4th ed.; Prentice-Hall: New York, NY, USA, 2000. [Google Scholar]
  15. Gelman, A.; Carlin, J.B.; Stern, H.S.; Rubin, D.B. Bayesian Data Analysis, 2nd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2004. [Google Scholar]
  16. Lynch, S.M. Introduction to Applied Bayesian Statistics and Estimation for Social Scientists; Springer: New York, NY, USA, 2007. [Google Scholar]
  17. Chen, M.H.; Shao, Q.M. Monte Carlo estimation of Bayesian credible and HPD intervals. J. Comput. Graph. Stat. 1999, 8, 69–92. [Google Scholar]
  18. 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]
  19. Riabiz, M.; Chen, W.Y.; Cockayne, J.; Swietach, P.; Niederer, S.A.; Mackey, L.; Oates, C.J. Optimal thinning of MCMC output. J. R. Stat. Soc. Ser. Stat. Methodol. 2022, 84, 1059–1081. [Google Scholar] [CrossRef]
  20. Aarset, M.V. How to identify a bathtub hazard rate. IEEE Trans. Reliab. 1987, 36, 106–108. [Google Scholar] [CrossRef]
  21. Mahdavi, A.; Kundu, D. A new method for generating distributions with an application to exponential distribution. Commun. Stat. Theory Methods 2017, 46, 6543–6557. [Google Scholar] [CrossRef]
  22. Johnson, N.; Kotz, S.; Balakrishnan, N. Continuous Univariate Distributions, 2nd ed.; Wiley: New York, NY, USA, 1994. [Google Scholar]
  23. Nadarajah, S.; Haghighi, F. An extension of the exponential distribution. Statistics 2011, 45, 543–558. [Google Scholar] [CrossRef]
  24. Gupta, R.D.; Kundu, D. Generalized exponential distribution: Different method of estimations. J. Stat. Comput. Simul. 2001, 69, 315–337. [Google Scholar] [CrossRef]
  25. Weibull, W. A statistical distribution function of wide applicability. J. Appl. Mech. 1951, 18, 293–297. [Google Scholar] [CrossRef]
  26. Marinho, P.R.D.; Silva, R.B.; Bourguignon, M.; Cordeiro, G.M.; Nadarajah, S. AdequacyModel: An R package for probability distributions and general purpose optimization. PLoS ONE 2019, 14, e0221487. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Murthy, D.N.P.; Xie, M.; Jiang, R. Weibull Models; Wiley Series in Probability and Statistics; Wiley: Hoboken, NJ, USA, 2004. [Google Scholar]
  28. 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]
  29. 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]
  30. Sen, T.; Tripathi, Y.M.; Bhattacharya, R. Statistical inference and optimum life testing plans under Type-II hybrid censoring scheme. Ann. Data Sci. 2018, 5, 679–708. [Google Scholar] [CrossRef]
  31. Ashour, S.K.; El-Sheikh, A.A.; Elshahhat, A. Inferences and optimal censoring schemes for progressively first-failure censored Nadarajah-Haghighi distribution. Sankhyā Indian J. Stat. 2022, 84, 885–923. [Google Scholar] [CrossRef]
  32. 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]
Figure 1. The PDF (left) and HRF (right) shapes of the PRC model.
Figure 1. The PDF (left) and HRF (right) shapes of the PRC model.
Mathematics 11 03124 g001
Figure 2. Plots of the conditional PDFs of μ (left) and γ (right).
Figure 2. Plots of the conditional PDFs of μ (left) and γ (right).
Mathematics 11 03124 g002
Figure 3. Trace (with Gaussian kernel) and BGR plots of μ and γ in Monte Carlo simulation.
Figure 3. Trace (with Gaussian kernel) and BGR plots of μ and γ in Monte Carlo simulation.
Mathematics 11 03124 g003
Figure 4. Heatmap for the inferred results of μ .
Figure 4. Heatmap for the inferred results of μ .
Mathematics 11 03124 g004
Figure 5. Heatmap for the inferred results of γ .
Figure 5. Heatmap for the inferred results of γ .
Mathematics 11 03124 g005
Figure 6. Heatmap for the inferred results of R ( t ) .
Figure 6. Heatmap for the inferred results of R ( t ) .
Mathematics 11 03124 g006
Figure 7. Heatmap for the inferred results of h ( t ) .
Figure 7. Heatmap for the inferred results of h ( t ) .
Mathematics 11 03124 g007
Figure 8. The PP (a), Fitted PDF (b) Fitted RF (c), Contour of log-likelihood (d,e) scaled–TTT transform from industrial devices data.
Figure 8. The PP (a), Fitted PDF (b) Fitted RF (c), Contour of log-likelihood (d,e) scaled–TTT transform from industrial devices data.
Mathematics 11 03124 g008
Figure 9. Density (left) and Trace (right) plots of μ , γ , R ( t ) and h ( t ) from industrial devices data.
Figure 9. Density (left) and Trace (right) plots of μ , γ , R ( t ) and h ( t ) from industrial devices data.
Mathematics 11 03124 g009
Figure 10. The PP (a), Fitted PDF (b) Fitted RF (c), Contour of log-likelihood (d), (e) scaled–TTT transform from AWST data.
Figure 10. The PP (a), Fitted PDF (b) Fitted RF (c), Contour of log-likelihood (d), (e) scaled–TTT transform from AWST data.
Mathematics 11 03124 g010aMathematics 11 03124 g010b
Figure 11. Density (left) and Trace (right) plots of μ , γ , R ( t ) and h ( t ) from AWST data.
Figure 11. Density (left) and Trace (right) plots of μ , γ , R ( t ) and h ( t ) from AWST data.
Mathematics 11 03124 g011
Table 1. Failure times of 50 industrial devices.
Table 1. Failure times of 50 industrial devices.
0.10.211111236
7111218181818182132
36404545475055606363
67676767727579828283
84848485858585858686
Table 2. Results of fit of the PRC and its competitive models from industrial devices data.
Table 2. Results of fit of the PRC and its competitive models from industrial devices data.
ModelMLE(St.E)NLABCAK–S (p-Value)
HQA–DCvM
μ γ
PRC0.0384 (0.0059)0.2160 (0.0881)234.7804473.561477.385473.8160.1484 (0.2209)
475.0172.133980.32793
APE6.5479 (17.940)0.0307 (0.0117)241.1075486.215490.039486.4700.1587 (0.1609)
487.6712.568960.41510
LoL1.0880 (0.1334)29.937 (6.7031)251.0748506.150509.974506.4050.2408 (0.0061)
507.6064.033820.69540
NH2.2030 (0.5870)0.0074 (0.0024)238.3326480.665484.489480.9210.1738 (0.0976)
482.1222.555770.41148
GE0.7804 (0.1351)0.0187 (0.0036)239.9734483.947487.771484.2020.2045 (0.0305)
485.4032.941970.48371
W0.9487 (0.1195)44.793 (6.9216)240.9797485.959489.783486.2150.1937 (0.0469)
487.4163.000240.49469
G0.7991 (0.1375)0.0175 (0.0041)240.1683484.337488.161484.5920.2024 (0.0333)
485.7932.962490.48758
E0.0219 (0.0031)-241.068484.135486.047484.2190.1915 (0.0512)
484.8642.954540.48605
Table 3. Various AT2-PH samples from industrial devices data.
Table 3. Various AT2-PH samples from industrial devices data.
Sample T ( d ) R R m Censored Data
S[1]0.5 (2) ( 5 5 , 0 20 ) 150.1, 0.2, 1, 1, 1, 1, 2, 3, 6, 7, 11, 12, 18, 18, 18, 21, 32, 36, 40, 45, 45, 47, 50, 55, 60
S[2]15 (13) ( 0 10 , 5 5 , 0 10 ) 100.1, 0.2, 1, 1, 1, 1, 1, 2, 3, 6, 7, 11, 12, 18, 18, 18, 21, 32, 36, 40, 47, 50, 55, 60, 63
S[3]60 (24) ( 0 20 , 5 5 ) 50.1, 0.2, 1, 1, 1, 1, 1, 2, 3, 6, 7, 11, 12, 18, 18, 18, 18, 18, 21, 32, 36, 40, 45, 50, 67
Table 4. Estimates of μ , γ , R ( t ) and h ( t ) from industrial devices data.
Table 4. Estimates of μ , γ , R ( t ) and h ( t ) from industrial devices data.
SamplePar.MLEMCMCACIHPD
Est.St.EEst.St.ELowerUpperLengthLowerUpperLength
S[1] μ 0.02370.00730.00940.03800.02860.02330.00340.01670.03000.0133
γ 0.45770.20200.06180.85360.79180.45750.00500.44740.46710.0197
R ( 1 ) 0.89130.03700.81870.96380.14510.89190.00390.88430.89940.0150
h ( 1 ) 0.02860.00840.01220.04510.03290.02830.00200.02430.03220.0079
S[2] μ 0.02750.00790.01210.04290.03080.02710.00360.02020.03430.0141
γ 0.40590.16800.07660.73520.65860.40550.00500.39560.41520.0196
R ( 1 ) 0.89890.03270.83490.96290.12800.89940.00350.89280.90640.0136
h ( 1 ) 0.02760.00780.01240.04290.03040.02740.00190.02380.03110.0073
S[3] μ 0.02080.00720.00670.03490.02820.02040.00330.01410.02680.0128
γ 0.45730.19570.07370.84090.76720.45700.00500.44710.46660.0195
R ( 1 ) 0.89460.03420.82770.96160.13390.89530.00400.88770.90320.0155
h ( 1 ) 0.02690.00730.01250.04120.02870.02650.00210.02270.03060.0079
Table 5. Characteristics of μ , γ , R ( t ) and h ( t ) from industrial devices data.
Table 5. Characteristics of μ , γ , R ( t ) and h ( t ) from industrial devices data.
SamplePar.MeanMode Q 1 Q 2 Q 3 St.DSkew.
S[1] μ 0.023350.022960.021020.023300.025620.003400.09963
γ 0.457520.450710.454160.457520.460890.00500 0.01970
R ( 1 ) 0.891900.890210.889260.891830.894430.003840.12186
h ( 1 ) 0.028340.027870.026980.028340.029720.00202 0.02828
S[2] μ 0.027080.025140.024610.027020.029540.003610.07485
γ 0.405540.400200.402110.405510.408940.005010.02838
R ( 1 ) 0.899440.902300.897070.899390.901770.003470.09654
h ( 1 ) 0.027380.026110.026120.027390.028650.00187 0.02191
S[3] μ 0.020420.018200.018150.020340.022650.003280.10599
γ 0.456990.455230.453580.456960.460390.005010.02319
R ( 1 ) 0.895330.893550.892600.895250.897970.003980.14052
h ( 1 ) 0.026550.025510.025170.026550.027950.00203 0.04074
Table 6. Aircraft windshield service times.
Table 6. Aircraft windshield service times.
0.0460.1400.1500.2480.2800.3130.3890.4870.6220.900
0.9520.9961.0031.0101.0851.0921.1521.1831.2441.249
1.2621.4361.4921.5801.7191.7941.9151.9201.9631.978
2.0532.0652.1172.1372.1412.1632.1832.2402.3412.435
2.4642.5432.5922.6002.6702.7172.8192.8202.8782.950
3.0033.1023.3043.4833.5003.6223.6653.6954.0154.628
4.8064.8815.140
Table 7. Results of fit of the PRC and its competitive models from AWST data.
Table 7. Results of fit of the PRC and its competitive models from AWST data.
ModelMLE(St.E)NLABCAK–S (p-Value)
HQA–DCvM
μ γ
PRC1.2649 (0.1551)0.0830 (0.0329)99.6012203.202207.489203.4020.0652 (0.9356)
204.8880.227010.03013
APE21.339 (14.195)0.8551 (0.0972)100.430204.860209.146205.0600.1167 (0.3315)
206.5460.602200.09895
LoL2.1352 (0.2312)1.8015 (0.1809)108.742221.485225.771221.6850.1199 (0.3012)
223.1701.978840.32711
NH7.3061 (2.9223)0.0447 (0.0192)100.886205.773210.059205.9730.1483 (0.1127)
207.4580.443580.07303
GE1.8978 (0.3402)0.6921 (0.0942)103.547211.093215.380211.2930.1438 (0.1340)
212.7791.231530.20347
W1.6327 (0.1686)2.3124 (0.1866)100.318204.636208.922204.8360.1078 (0.4270)
206.3220.629800.10395
G1.9085 (0.3148)0.9152 (0.1725)102.833209.665213.951209.8650.1387 (0.1613)
211.3511.111740.18364
E0.4796 (0.0604)-109.299220.597222.740220.6630.2078 (0.0073)
221.4401.126430.18614
Table 8. Various AT2-PH samples from AWST data.
Table 8. Various AT2-PH samples from AWST data.
Sample T ( d ) R R m Censored Data
S[1]0.1 (1) ( 10 3 , 0 30 ) 200.046, 0.140, 0.150, 0.248, 0.280, 0.313, 0.389, 0.487, 0.622, 0.900, 0.952,
0.996, 1.003, 1.010, 1.092, 1.152, 1.183, 1.244, 1.262, 1.436, 1.492, 1.580,
1.719, 1.915, 1.920, 1.978, 2.053, 2.065, 2.117, 2.137, 2.141, 2.163, 2.240
S[2]1.2 (17) ( 0 15 , 10 3 , 0 15 ) 100.046, 0.140, 0.150, 0.248, 0.280, 0.313, 0.389, 0.487, 0.622, 0.900, 0.952,
0.996, 1.003, 1.010, 1.085, 1.092, 1.183, 1.244, 1.249, 1.262, 1.436, 1.580,
1.719, 1.794, 1.978, 2.053, 2.141, 2.163, 2.341, 2.435, 2.543, 2.592, 2.600
S[3]2.2 (33) ( 0 30 , 10 3 ) 00.046, 0.140, 0.150, 0.248, 0.280, 0.313, 0.389, 0.487, 0.622, 0.900, 0.952,
0.996, 1.003, 1.010, 1.085, 1.092, 1.152, 1.183, 1.244, 1.249, 1.262, 1.436,
1.492, 1.580, 1.719, 1.794, 1.915, 1.920, 1.963, 1.978, 2.053, 2.065, 2.117
Table 9. Estimates of μ , γ , R ( t ) and h ( t ) from AWST data.
Table 9. Estimates of μ , γ , R ( t ) and h ( t ) from AWST data.
SamplePar.MLEMCMCACIHPD
Est.St.EEst.St.ELowerUpperLengthLowerUpperLength
S[1] μ 1.28650.24110.81391.75900.94501.28620.00501.27631.29580.0195
γ 0.09120.04540.00220.18010.17800.09070.00490.08120.10050.0192
R ( 0.1 ) 0.95870.01690.92550.99190.06640.95890.00210.95460.96300.0084
h ( 0.1 ) 0.18320.06490.05600.31040.25440.18220.00950.16380.20090.0371
S[2] μ 1.33420.23740.86901.79950.93051.33400.00501.32411.34350.0195
γ 0.08050.03810.00570.15520.14950.07990.00490.07050.08960.0191
R ( 0.1 ) 0.96280.01460.93410.99140.05730.96300.00220.95870.96720.0085
h ( 0.1 ) 0.16730.05700.05560.27900.22340.16620.00980.14770.18600.0383
S[3] μ 1.24500.24970.75561.73450.97891.24480.00501.23491.25430.0195
γ 0.08560.04340.00060.17060.17000.08510.00490.07570.09480.0191
R ( 0.1 ) 0.96160.01590.93050.99280.06230.96180.00210.95760.96590.0082
h ( 0.1 ) 0.16800.05930.05180.28430.23250.16700.00930.14930.18540.0361
Table 10. Characteristics of μ , γ , R ( t ) and h ( t ) from AWST data.
Table 10. Characteristics of μ , γ , R ( t ) and h ( t ) from AWST data.
SamplePar.MeanMode Q 1 Q 2 Q 3 St.DSkew.
S[1] μ 1.286181.288101.282791.286171.289520.004960.01673
γ 0.090650.082970.087300.090650.093970.004910.02521
R ( 0.1 ) 0.958900.962220.957450.958890.960350.00213 0.01195
h ( 0.1 ) 0.182150.167520.175720.182160.188580.009470.01234
S[2] μ 1.333961.321701.330571.333951.337300.004960.01546
γ 0.079920.069190.076600.079920.083210.004870.02718
R ( 0.1 ) 0.963010.964750.961540.963010.964490.00217 0.01377
h ( 0.1 ) 0.166210.143930.159580.166220.172820.009760.01415
S[3] μ 1.244771.232511.241371.244761.248110.004960.01716
γ 0.085100.074350.081770.085090.088390.004870.02597
R ( 0.1 ) 0.961850.963560.960430.961850.963280.00210 0.01292
h ( 0.1 ) 0.167000.145880.160750.167000.173220.009210.01331
Table 11. Criteria for optimum censoring.
Table 11. Criteria for optimum censoring.
CriterionGoal
O 1 Minimize det ( I 1 ( μ ^ , γ ^ ) )
O 2 Minimize trace ( I 1 ( μ ^ , γ ^ ) )
O 3 Maximize trace ( I ( μ ^ , γ ^ ) )
O 4 Minimize v a r ^ ( log ( Q ^ ζ ) )
Table 12. Optimum progressive plans from industrial devices data.
Table 12. Optimum progressive plans from industrial devices data.
Sample O 1 O 2 O 3 O 4
ζ 0.30.60.9
S[1]8.72 × 10 7 0.0408546,849.3850.0820126.346750.602
S[2]7.71 × 10 7 0.0282936,712.7138.7364102.774558.193
S[3]7.49 × 10 7 0.0383651,235.3354.2833186.1831316.36
Table 13. Optimum progressive plans from AWST data.
Table 13. Optimum progressive plans from AWST data.
Sample O 1 O 2 O 3 O 4
ζ 0.30.60.9
S[1]2.89 × 10 5 0.060182082.760.032950.041280.14746
S[2]2.46 × 10 5 0.064252613.760.030970.045830.18699
S[3]2.17 × 10 5 0.057802658.200.029400.041190.13977
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.; Abo-Kasem, O.E.; Mohammed, H.S. Survival Analysis of the PRC Model from Adaptive Progressively Hybrid Type-II Censoring and Its Engineering Applications. Mathematics 2023, 11, 3124. https://doi.org/10.3390/math11143124

AMA Style

Elshahhat A, Abo-Kasem OE, Mohammed HS. Survival Analysis of the PRC Model from Adaptive Progressively Hybrid Type-II Censoring and Its Engineering Applications. Mathematics. 2023; 11(14):3124. https://doi.org/10.3390/math11143124

Chicago/Turabian Style

Elshahhat, Ahmed, Osama E. Abo-Kasem, and Heba S. Mohammed. 2023. "Survival Analysis of the PRC Model from Adaptive Progressively Hybrid Type-II Censoring and Its Engineering Applications" Mathematics 11, no. 14: 3124. https://doi.org/10.3390/math11143124

APA Style

Elshahhat, A., Abo-Kasem, O. E., & Mohammed, H. S. (2023). Survival Analysis of the PRC Model from Adaptive Progressively Hybrid Type-II Censoring and Its Engineering Applications. Mathematics, 11(14), 3124. https://doi.org/10.3390/math11143124

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