Next Article in Journal
Local Refinement and Adaptive Strategy for a System of Free Boundary Power Options with High Order Compact Differencing
Next Article in Special Issue
An Enhanced Spatial Capture Model for Population Analysis Using Unidentified Counts through Camera Encounters
Previous Article in Journal
New Applications of Faber Polynomials and q-Fractional Calculus for a New Subclass of m-Fold Symmetric bi-Close-to-Convex Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Analysis of Type-II Generalized Progressively Hybrid Alpha-PIE Censored Data and Applications in Electronic Tubes and Vinyl Chloride

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.
Axioms 2023, 12(6), 601; https://doi.org/10.3390/axioms12060601
Submission received: 14 May 2023 / Revised: 11 June 2023 / Accepted: 15 June 2023 / Published: 16 June 2023
(This article belongs to the Special Issue Stochastic and Statistical Analysis in Natural Sciences)

Abstract

:
A new Type-II generalized progressively hybrid censoring strategy, in which the experiment is ensured to stop at a specified time, is explored when the lifetime model of the test subjects follows a two-parameter alpha-power inverted exponential (Alpha-PIE) distribution. Alpha-PIE’s parameters and reliability indices, such as reliability and hazard rate functions, are estimated via maximum likelihood and Bayes estimation methodologies in the presence of the proposed censored data. The estimated confidence intervals of the unknown quantities are created using the normal approximation of the acquired classical estimators. The Bayesian estimators are also produced using independent gamma density priors under symmetrical (squared-error) loss. The Bayes’ estimators and their associated highest posterior density intervals cannot be calculated theoretically since the joint likelihood function is derived in a complicated form, but they can potentially be assessed using Monte Carlo Markov-chain algorithms. We next go through four optimality criteria for identifying the best progressive design. The effectiveness of the suggested estimation procedures is assessed using Monte Carlo comparisons, and certain recommendations are offered. Ultimately, two different applications, one focused on the failure times of electronic tubes and the other on vinyl chloride, are analyzed to illustrate the effectiveness of the proposed techniques that may be employed in real-world scenarios.

1. Introduction

Inverted exponential distribution has been popularly used in several sectors such as engineering, biology, medicine, and others. It has an inverted bathtub failure rate, which describes the relative failure rate, which initially increases, comes to a head after some time, and then decreases over time. A new simple two-parameter alpha power-inverted exponential distribution (denoted as Alpha-PIE( α , μ )) has been proposed by Ünal et al. [1] as an extension of the conventional inverted exponential distribution. They also stated, by making use of two lifetime data sets, that the Alpha-PIE model offers superior fits to some other lifetime models, namely, Lindley, inverted exponential, generalized inverted exponential, and inverted Rayleigh distributions. However, a lifetime random variable X is said to have the Alpha-PIE ( α , μ )  distribution, where  ϑ = ( α , μ ) T  is a vector of the model parameters, if its probability density function (PDF)  f ( · ) ; cumulative distribution function, say  F ( · )  (CDF); reliability function, say  R ( · )  (RF); and hazard rate function, say  h ( · )  (HRF), are provided by
f ( x ; ϑ ) = μ log ( α ) α 1 x 2 exp μ x α exp μ x , x > 0 , α 1 ,
F ( x ; ϑ ) = α exp μ x 1 α 1 ,
R ( x ; ϑ ) = α α 1 1 α exp μ x 1 , x > 0 ,
and
h ( x ; ϑ ) = μ log ( α ) x 2 exp μ x α exp μ x 1 1 α exp μ x 1 ,
respectively, where  α > 0  and  μ > 0  are the shape and scale parameters, respectively. Obviously, when we put  α 1 , the alpha-PIE distribution is reduced to the inverted exponential distribution by Keller et al. [2]. When we fix  μ = 1  and consider different values of  α , in Figure 1, several shapes for the PDF and HRF of the Alpha-PIE distribution are demonstrated. This shows that the density of the alpha-PIE distribution is unimodal, while its HRF is a monotonically increasing function. Recently, in the presence of complete data, Amjad et al. [3] derived various Bayes estimators of the alpha-PIE parameters and stated that the alpha-PIE model behaved better than the inverted exponential and generalized inverted exponential models.
In a reliability context, the failure time data of the experimental objects are frequently not totally accessible, so minimizing both the expense and duration of the experiment is critical for any researcher. Progressive Type-II censoring (PTIIC) was proposed to allow the experimenter to remove items at times other than the stop time. This technique may be advantageous when some of the live test objects that are removed early on can be used again for future tests or when a balance between observing even some extreme lifetimes and reducing testing duration is desired, one may refer to Balakrishnan and Cramer [4] for additional details. Recently, if PCS-II is applied, because of the lengthy lifespan of many items, particularly electronics, the total experimental time might be quite long. This is the main drawback of PTIIC sampling. Therefore, to ensure that the experiment is completed at a specified optimal time without losing the ability to collect an effective number of observable failures, a generalized Type-II progressively hybrid censored (GTIIPHC) strategy was introduced by Lee et al. [5]. This strategy starts with putting n independent items into a test at time zero, specifying two threshold points of failures,  T i , i = 1 , 2  and the desired number m, where ( 1 < m n ), and designing the progressive censoring  R = ( R 1 , R 2 , , R m ) . When the experimenter records the first failure (say  X 1 : m : n ),  R 1  items are randomly chosen from  n 1  and removed from the test; next, at the second failure (say  X 2 : m : n ),  R 2  items are randomly chosen from  n R 1 2  and out of the test, and so on. At  T * = max min X m : m : n , T 2 , T 1 , where  0 < T 1 < T 2 , the test is stopped and all remaining survival subjects are removed. Note that ( d 1 , d 2 ) represents the size of recorded failures up to ( T 1 , T 2 ).
If  X m : n < T 1 , the test is stopped at  T 1  (Case 1) without any additional removals. If  T 1 < X m : m : n < T 2 , similar to the PTIIC, the test is stopped at  X m : m : n  (Case-2); otherwise, the test is stopped at  T 2  (Case-3). Thus, the partitioner will collect one of the samples:
{ X , R } = { ( X 1 : m : n , R 1 ) , , ( X m 1 : m : n , R m 1 ) , ( X m : m : n , 0 ) , , ( X d 1 : n , 0 ) } ; Case- 1 , { ( X 1 : m : n , R 1 ) , , ( X d 1 : n , R d 1 ) , , ( X m 1 : m : n , R m 1 ) , ( X m : m : n , R m ) } ; Case- 2 , { ( X 1 : m : n , R 1 ) , , ( X d 1 : n , R d 1 ) , , ( X d 2 1 : n , R d 2 1 ) , ( X d 2 : n , R d 2 ) } ; Case- 3 .
Let  { X , R }  be the censoring members of GTIIPHC in a distribution with a PDF and CDF. Therefore, the joint likelihood function, say  L q ( · ) , of the proposed censoring mechanism can be expressed as
L q ( ϑ | X ) = K q q ( T τ ; ϑ ) i = 1 D q f ( x i : m : n ; ϑ ) 1 F ( x i : m : n ; ϑ ) R i , q = 1 , 2 , 3 ,
where  q ( · )  is a composite term of reliability functions and  τ = 1 , 2 .
It is critical to note that the GTIIPHC modifies the Type-II progressive hybrid censoring (PHC), proposed by Childs et al. [6], by guaranteeing that the test is completed at the specified time  T 2 . Thus,  T 2  is the greatest duration that the examiner is willing to allow the test to continue. Table 1 lists the GTIIPHC notations. Furthermore, from (5), six sampling plans can be introduced and provided in Table 2. A diagrammatic demonstration of GTIIPHC sampling is depicted in Figure 2.
In the context of GTIIPHC data, several researchers have carried out important research on the statistical estimation of unknown parameter(s) and/or reliability time functions in various lifetime models; for example, Ashour and Elshahhat [11] studied both frequentist and Bayes estimators of the Weibull parameters; Ateya and Mohammed [12] studied the prediction issue of the Burr-XII failure times; Seo [13] discussed the Bayesian inference of Weibull’s model; Cho and Lee [14] analyzed the competing risks from exponential data; Nagy et al. [15] pointed out different estimates of the Burr-XII parameters; Wang et al. [16] derived various estimators of the Kumaraswamy parameters; Elshahhat et al. [17] addressed the Nadarajah–Haghighi parameters; later, Alotaibi et al. [18] estimated the Fréchet Parameters.
Although there are many studies that give a mathematical treatment to the proposed distribution, they do not shed light on the application aspects of the alpha-PIE distribution, especially in reliable practice. To our understanding, no study has been performed to evaluate the parameters or reliability features of the alpha-PIE lifetime model when an incomplete dataset is available. Thus, to resolve this issue, our objectives in this work are fourfold:
  • Derive the maximum likelihood estimators (MLEs) in addition to their two-sided approximate confidence intervals (ACIs), using observed Fisher’s information, of the alpha-PIE parameters  α  and  μ  or any associated function such as  R ( t )  and  h ( t ) .
  • Derive the Bayes’ estimators in addition to their two-sided highest posterior density (HPD) intervals, under independent gamma priors assumption, of  α μ R ( t ) , and  h ( t )  using the squared error loss (SEL) function.
  • To select the best progressive censoring patterns among various competing strategies, several criteria of optimality are proposed.
  • Via extensive Monte Carlo simulations, on the basis of four accuracy criteria, namely, (i) root mean squared-errors, (ii) mean relative absolute biases, (iii) average confidence lengths, and (iv) coverage percentages, the performance of the acquired estimators is examined. Additionally, two real-word applications from the engineering and chemistry sectors, to evaluate how the offered approaches operate in practice and to choose the best censoring strategy, are examined.
The organization of the article is as follows: Section 2 and Section 3 provide the point and interval inferences using frequentist and Bayes approaches. Simulation results are obtained and discussed in Section 4. Section 5 examines two applications of actual data sets. Optimum criteria of progressive patterns are presented in Section 6. Lastly, Section 7 lists the study’s conclusions.

2. Classical Inference

This section investigates the maximum likelihood and approximate asymptotic interval estimators of  α μ R ( t )  and  h ( t ) .

2.1. Maximum Likelihood Estimators

Suppose  x = { x i : m : n , R i }  for  i = 1 ,   ,   d 2  is a GTIIPHC sample of size  d 2  collected from the alpha-PIE population with PDF (1) and CDF (2). Substitute (1) and (2) into (5). Utilizing  x i  in place of  x i : m : n , the likelihood function of GTIIPHC (5), ignoring the constant term, can be expressed as
L q ( ϑ | x ) ξ q T τ ; ϑ α α 1 n μ log ( α ) α D q e i = 1 D q η x i , R i ; ϑ , q = 1 , 2 , 3 ,
where  ξ 1 T 1 ; ϑ = α α 1 1 α exp ( μ T 1 1 ) 1 R d 1 + 1 * ξ 3 T 2 ; ϑ = α α 1 1 α exp ( μ T 2 1 ) 1 R d 2 + 1 * ξ 2 T τ ; ϑ = 1  and  η x i , R i ; ϑ = μ x i 1 exp ( μ x i 1 ) log α R i log 1 α exp ( μ x i 1 ) 1 .
Consequently, the natural logarithm of (6) becomes
log L q ( ϑ | x ) ζ q T τ ; ϑ + n log α α 1 + D q log μ log ( α ) α i = 1 D q η x i , R i ; ϑ ,
where  ζ 1 T 1 ; ϑ = R d 1 + 1 * log α α 1 1 α exp ( μ T 1 1 ) 1 ,   ζ 2 T τ ; ϑ = 0 , and  ζ 3 T 2 ; ϑ = R d 2 + 1 * log α α 1 1 α exp ( μ T 2 1 ) 1 .
To acquire the MLEs  α ^  and  μ ^  of  α  and  μ , respectively, by differentiating (7) with regard to  α  and  μ , we obtain the following normal equations:
α log L q ϑ x = ζ q T τ ; ϑ + 1 α D q log 1 ( α ) 1 n α 1 i = 1 D q η x i , R i ; ϑ ,
and
μ log L q ϑ x = ζ q T τ ; ϑ + D q μ 1 i = 1 D q η x i , R i ; ϑ ,
where  ζ 2 T τ ; ϑ = ζ 2 T τ ; ϑ = 0 ,   ζ q T τ ; ϑ = R d τ + 1 * ( α 1 ) R T τ ; α , μ { α exp ( μ T τ 1 ) 1 e μ T τ 1 ( α 1 ) 1 [ α exp ( μ T τ 1 ) 1 1 ] } , and  ζ q T τ ; ϑ = R d τ + 1 * ( α 1 ) T τ R T τ ; α , μ α exp ( μ T τ 1 ) log ( α ) e μ T τ 1 .
It is more clear, due to the complex expressions in (8) and (9), that the MLEs  α ^  and  μ ^  cannot be derived in closed form. Applying the Newton–Raphson procedure via the ’ maxLik ’ package, proposed by Henningsen and Toomet [19], the offered estimators can be directly evaluated. Once  α ^  and  μ ^  are obtained, replacing  α  and  μ  in (3) and (4) by their MLEs  α ^  and  μ ^ , the MLEs  R ^ ( t )  and  h ^ ( t )  of  R ( t )  and  h ( t ) , for a mission time  t > 0 , are
R ^ ( t ) = α ^ α ^ 1 1 α ^ exp μ ^ t 1 ,
and
h ^ ( t ) = μ ^ log ( α ^ ) t 2 exp μ ^ t α ^ exp μ ^ t 1 1 α ^ exp μ ^ t 1 ,
respectively.
In addition to obtaining a point estimate for an unknown parameter, it is also beneficial to obtain a limit of values that may include the real parameter with a given degree of confidence, this methodology is called interval estimation. Therefore, in the next subsection, we shall derive the ACIs of all unknown quantities of our interest.

2.2. Asymptotic Interval Estimators

To construct the  ( 1 ϵ ) 100 %  ACIs for  α μ R ( t )  or  h ( t ) , the asymptotic  2 × 2  variance-covariance (inverted Fisher information) matrix must be first obtained. Following mild regularity conditions, we can say that  ϑ ^  has an approximately normal distribution, i.e.,  ϑ N ( ϑ , I 1 ( ϑ ) ) . By replacing  α  and  μ  by  α ^  and  μ ^ , respectively, one can estimate  I 1 ( ϑ )  by  I 1 ( ϑ ^ )  (see Lawless [20]), as
I 1 ( ϑ ^ ) F 11 F 12 F 21 F 22 1 = v ^ 11 v ^ 12 v ^ 21 v ^ 22 ,
where  ϑ ^ = ( α ^ , μ ^ ) T  and  F i j , i , j = 1 , 2  are presented in Appendix A. Therefore, the two-sided  100 ( 1 ϵ ) 100 %  ACIs for  α  and  μ  are given by
α ^ z ϵ / 2 v ^ 11 and μ ^ z ϵ / 2 v ^ 22 ,
respectively, where  z ϵ / 2  refers to the upper  ϵ / 2  percentage level for the distribution of the standard normal.
On the other hand, to build the two-sided  100 ( 1 ϵ ) 100 %  ACIs of RF  R ( t )  and HRF  h ( t ) , following Greene [21], the delta method is considered to derive the estimated variances  v ^ R  and  v ^ h  of  R ^ ( t )  and  h ^ ( t ) , respectively, as
v ^ R = C R T I 1 ( ϑ ) C R ( α ^ , μ ^ ) and v ^ h = C h T I 1 ( ϑ ) C h ( α ^ , μ ^ ) ,
where  C R T = [ R ( t ) α R ( t ) μ ]  and  C h T = [ h ( t ) α h ( t ) μ ] .
Therefore, the  ( 1 ϵ ) 100 %  ACIs of  R ( t )  and  h ( t )  are given by
R ^ ( t ) z ϵ / 2 v ^ R and h ^ ( t ) z ϵ / 2 v ^ h ,
respectively.

3. Bayes Inference

In this section, the Bayes estimators of  α μ R ( t ) , and  h ( t ) , along with their HPD intervals, based on the SEL, are obtained.

3.1. Prior Functions

In practice, the Bayesian technique might leverage additional prior information, such as historical data or expertise in the statistical inferential process, to obtain more accurate estimates for tests with small sample sizes or when censored data are available. Regarding Bayes’ reliability analysis, several recent works have addressed this issue; see, for example, Chen and Ye [22], Wang et al. [23], Luo et al. [24], and Luo and Xu [25]. To acquire the Bayes point estimator of  α μ R ( t )  or  h ( t ) , we suppose that  α  and  μ  are independent and distributed with gamma (G) density priors, i.e.,  α G ( a 1 , b 1 )  and  μ G ( a 2 , b 2 ) , respectively. We employ gamma priors, which are regarded as having greater adaptability than others, and adjust them to the parameters’ support. Furthermore, the independent gamma priors are plain and unambiguous, potentially avoiding many complex inferential concerns. Then, the joint prior PDF of  α  and  μ , denoted by  ρ ( · ) , is
ρ ( α , μ ) α a 1 1 μ a 2 1 exp ( α b 1 + μ b 2 ) ,
where  a i > 0  and  b i > 0  for  i = 1 , 2 ,  must be known. In the next subsections, we derive the Bayes point and HPD interval estimators of  α μ R ( t ) , or  h ( t ) .

3.2. Bayes Estimators

From (5) and (11), the joint posterior PDF (say  Q ( · ) ) of  α  and  μ  becomes
Q q ϑ x α n + a 1 D q 1 μ D q a 2 1 α 1 n log ( α ) D q ξ q T τ ; ϑ e b 1 α + b 2 μ + i = 1 D q η x i , R i ; ϑ , q = 1 , 2 , 3 ,
where its normalizing term, say  Q q * , is given by
Q q * = 0 0 α n + a 1 D q 1 μ D q a 2 1 α 1 n log ( α ) D q ξ q T τ ; ϑ e b 1 α + b 2 μ + i = 1 D q η x i , R i ; ϑ , d α d μ .
Under the SEL, the Bayes estimate (say  π ˜ ( · ) ) of any parametric function of  α  and  μ  (say  π ( · ) ) is known as the posterior mean of (12); for example,
π ˜ ( α , μ ) = 0 0 π ( α , μ ) Q q ϑ x d α d μ .
It is clear, from (13), that the Bayes estimate of  α  or  μ  cannot be extracted analytically because the posterior expectation involves double integrals in both the numerator and the denominator, and Bayes expressions of  α  and  μ  thus cannot be easily extracted. As a result, the Bayes Markov chain Monte Carlo (MCMC) process is used to create samples from which the offered Bayes estimates and HPD intervals can be obtained.
To start the MCMC technique, from (12), the full conditionals PDFs of  α  and  μ  are
Q q α α x , μ α n + a 1 D q 1 α 1 n log ( α ) D q ξ q T τ ; ϑ e b 1 α + i = 1 D q η x i , R i ; ϑ ,
and
Q q μ μ x , α μ D q a 2 1 ξ q T τ ; ϑ e b 2 μ + i = 1 D q η x i , R i ; ϑ ,
respectively.
Since Eqns. (14) and (15) of  α  and  μ , respectively, cannot be expressed analytically with any statistical model. Thus, the Metropolis-Hastings (M-H) algorithm discussed by Gelman et al. [26] and Lynch [27] is considered for this purpose. The following sampling process of M-H is adopted:
Step 1:
Set initial values  α ( 0 ) = α ^  and  μ ( 0 ) = μ ^ .
Step 2:
Set  j = 1
Step 3:
Obtain  α *  and  μ *  from  N ( α ^ , v ^ 11 )  and  N ( μ ^ , v ^ 22 ) , respectively.
Step 4:
Obtain  ϕ α = π ρ α ( α * μ [ j 1 ] , x ) π ρ α ( α [ j 1 ] μ [ j 1 ] , x )  and  ϕ μ = π ρ μ μ * α [ j ] , x π ρ μ μ [ j 1 ] α [ j ] , x .
Step 5:
Obtain  u 1  and  u 2  from uniform  U ( 0 , 1 )  distribution.
Step 6:
If  u α min 1 , ϕ α  and  u μ min 1 , ϕ μ  set  α [ j ] = α *  and  μ [ j ] = μ * ; else set  α [ j ] = α [ j 1 ]  and  μ [ j ] = μ ( j 1 ) , respectively.
Step 7:
Put  j = j + 1 .
Step 8:
Redo Steps 3–7 for  B  times and obtain  α [ j ]  and  μ [ j ]  for  j = 1 ,   2 ,   , B .
Step 9:
Obtain the RF (3) and HRF (4) using  ( α [ j ] , μ [ j ] ) , j = 1 ,   2 ,   , B , at  t > 0 , respectively as
R [ j ] ( t ) = α [ j ] α [ j ] 1 1 α [ j ] exp μ [ j ] t 1 1 ,
and
h [ j ] ( t ) = μ [ j ] log ( α [ j ] ) t 2 exp μ [ j ] t 1 α [ j ] exp μ [ j ] t 1 1 1 α [ j ] exp μ [ j ] t 1 1 .
Step 10:
Obtain the Bayes estimate  π ˜ ( · )  of  π ( · ) , after eliminating the first  B *  samples as burn-in, as
π ˜ ( α , μ ) = 1 B B * j = B * + 1 B π [ j ] ( α , μ ) .
Step 11:
Obtain the HPD interval of  π ( · )  via ordering  π [ j ] ( α , μ )  for  j = B * + 1 , , B . Following Chen and Shao [28], the  ( 1 ϵ ) 100 %  HPD interval of  π ( α , μ )  is given by
π j * , π j * + 1 ϵ B B * ,
where  j * = B * + 1 , , B  is specified such that
δ j * + [ 1 ϵ B B * ] π j * = min 1 j ϵ ( B B * ) ( π ( j + [ 1 ϵ B B * ] ) π [ j ] ) ,
where the highest integer less (or equal) than to y is denoted by  [ y ] .

4. Monte Carlo Simulations

This section deals with comparing the behavior of the acquired frequentist and Bayes estimators of  α μ R ( t )  and  h ( t )  obtained in the proceeding sections via extensive Monte Carlo simulations. This goal is developed based on large 1000 GPHC-T-II samples generated from two different groups of the alpha-PIE parameters, namely Set-1, Alpha-PIE ( 1.2 , 0.2 ) , and Set-2, Alpha-PIE ( 1.5 , 0.5 ) . For distinct time  t = 0.1 , from Sets 1 and 2, the plausible values of ( R ( t ) h ( t ) ) are (0.87509, 2.89008) and (0.99453, 0.27545), respectively. Several choices of  T i , i = 1 , 2  (thresholds), n (total test items), m (effective test items), and  R  (progressive pattern) are also used, namely  ( T 1 , T 2 ) = ( 0.1 , 0.3 )  and (0.4,0.8),  n = ( 50 , 80 ) . For each n, the value of m is determined as failure percentages (FPs) such as  m n = (40,80)%. To assess the removal mechanism in GPHC-T-II, for each group  ( n , m ) , three patterns of  ( R 1 , R 2 , , R m )  are considered:
Scheme- 1 : R 1 = n m , R i = 0 f o r i 1 , Scheme- 2 : R m 2 = n m , R i = 0 f o r i m 2 , Scheme- 3 : R m = n m , R i = 0 f o r i m .
As soon as 1000 GPHC-T-II data are collected, all suggested point (or interval) estimators created by maximum likelihood (or Bayes) inferential approaches of  α μ R ( t )  and  h ( t )  are evaluated via  R  4.2.2 programming software. Therefore, we recommend installing two statistical packages in  R , namely (i) ’ maxLik ’ package (by Henningsen and Toomet [19]) to calculate the maximum likelihood estimates along with their 95% ACI estimates, (ii) ’ coda ’ package (by Plummer et al. [29]) to evaluate the Bayes MCMC estimates for the same unknown parameters, as well as corresponding HPD interval estimations. Taking the classical estimates of  α  and  μ  as starting points, to compute the acquired Bayes point/interval estimates of  α μ R ( t )  or  h ( t ) , we replicated the MCMC sampler 12,000 times and ignored the first 2000 times as burn-in. To demonstrate the performance of the Bayes findings on different gamma priors, following Kundu [30], two sets of the hyperparameters  ( a i , b i ) , i = 1 , 2  are utilized, called
(a)
For Set-1:
  • Prior-1: ( a 1 , a 2 ) = ( 6 , 1 ) and b i = 5 for i = 1 , 2 ;
  • Prior-2: ( a 1 , a 2 ) = ( 12 , 2 ) and b i = 10 for i = 1 , 2 .
(b)
For Set-2:
  • Prior-1: ( a 1 , a 2 ) = ( 7.5 , 2.5 ) and b i = 5 for i = 1 , 2 ;
  • Prior-2: ( a 1 , a 2 ) = ( 15 , 5 ) and b i = 10 for i = 1 , 2 .
To examine the convergence status of the collected iterations of  α μ R ( t ) , and  h ( t )  developed from the Bayes MCMC algorithm, for Set-1 along with Prior-1,  ( T 1 , T 2 , n , m ) = ( 0.1 , 0.3 , 50 , 20 )  and Scheme-1 as an example, the trace and autocorrelation plots are provided in Figure 3. Obviously, it is evident that the Markovian graphs of  α μ R ( t )  and  h ( t )  are satisfactorily mixed, and thus, the results of the acquired points (or interval) become more significant.
However, for each test, the average point estimates (Av.Es) of  α  are given by
α ˇ ¯ = 1 1000 i = 1 1000 α ˇ ( i ) ,
where  α ˇ ( i )  denotes the calculated estimate of  α  at the ith sample.
The acquired point estimates of  α  are compared based on their root mean squared-errors (RMSEs) and mean relative absolute biases (MRABs) as
RMSE ( α ˇ ) = 1 1000 i = 1 1000 α ˇ ( i ) α 2 ,
and
MRAB ( α ˇ ) = 1 1000 i = 1 1000 1 α α ˇ ( i ) α ,
respectively.
Further, the acquired interval estimates of  α  are compared with regard to their average confidence lengths (ACLs) and coverage percentages (CPs) as
ACL ( 1 ϵ ) % ( α ) = 1 1000 i = 1 1000 U α ˇ ( i ) L α ˇ ( i ) ,
and
CP ( 1 ϵ ) % ( α ) = 1 1000 i = 1 1000 1 L α ˇ ( i ) ; U α ˇ ( i ) α ,
respectively, where  1 ( · )  is the indicator function,  ( L ( · ) , U ( · ) )  denote the (lower,upper) bounds of  ( 1 ϵ ) %  ACI (or HPD) interval of  α . In a similar way, the Av.E, RMSE, BRAB, ACL, and CP values of  μ R ( t ) , or  h ( t )  can be easily computed.
A heat map is a tool for graphically representing numerical data. Via  R  4.2.2 software with the heat-map programming tool, the calculated criteria (including RMSEs, MRABs, ACLs, and CPs) of  α μ R ( t ) , and  h ( t )  are represented in Figure 4, Figure 5, Figure 6 and Figure 7, respectively. As supplementary materials, all numerical values of  α μ R ( t ) , or  h ( t )  are reported. To distinguish, for each plot in Figure 4, Figure 5, Figure 6 and Figure 7, the proposed approaches are displayed on the ‘x-axis’ line, whereas the censoring settings are displayed on the ‘y-axis’. In addition, as notations, the Bayes estimates (for Prior-I (say P1) as an example) are referred to as BE-P1, while the HPD interval is referred to as HPD-P1.
In terms of the smallest RMSE, MRAB, and ACL values as well as the largest CP values, from Figure 4, Figure 5, Figure 6 and Figure 7, the following comments are made:
  • All acquired point and interval estimates of  α μ R ( t )  or  h ( t )  have good behavior; this is a general note.
  • As n(or FP%) increases, all results of all unknown parameters of life perform satisfactorily. A similar point is also true when the spacing between n and m is reduced.
  • As  T i , i = 1 , 2  grow, for both Sets 1 and 2, the RMSEs, MRABs, and ACLs of  α μ R ( t )  and  h ( t )  narrowed down, while their CPs increase.
  • As anticipated, due to the gamma information, the Bayes point (or HPD interval) estimates of  α μ R ( t ) , or  h ( t )  behave better compared to the others.
  • All Bayesian computations performed based on Prior-2 provide more accurate results than those obtained based on Prior-1. This finding is due to the fact that the associated variance of Prior-2 is less than the associated variance of Prior-1.
  • Comparing the suggested schemes 1, 2, and 3, for both Sets 1 and 2, it is seen that the point (or interval) estimates of  δ  have good results when all survival items  n m  are removed at the first stage (i.e., Scheme-1) and of  μ R ( t ) , and  h ( t )  at the last stage (i.e., Scheme-3).
  • In summary, it is advised to use MCMC samples to estimate the model parameters and reliability features of the alpha-PIE lifetime model when Type-II generalized progressively hybrid censored data are available.

5. Real Applications

To show the applicability of the suggested estimating methodologies and to illustrate the ability to adapt the study objectives to real-world situations, this section demonstrates two different applications from the engineering and chemistry fields. In each application, we prove that the offered model furnishes a good fit compared to five other competitive models in literature.

5.1. Electronic Tubes

An electron tube is a device that conducts electricity via electrons through a vacuum or a gas within a sealed glass or metal container and has a variety of popular applications, such as radio and television. For further applications of electron tubes, one may refer to Rosebury [31]. This application provides an analysis of a data set, taken from the engineering area, representing twenty lifetime values (in hundred hours) of electronic tubes; see Table 3. This data set has been first reported by Dixit and Nooghabi [32] and later discussed by Ibrahim et al. [33].
First, before proceeding, to highlight the utility of the offered model as competitors based on the complete data on electronic tubes, the alpha-PIE distribution is compared to five other inverted distributions (for  x > 0 α > 0  is a shape parameter and  μ > 0  is a scale parameter), namely: inverted exponential (IE( μ )) by Keller et al. [2]; inverted Lindley (IL( μ )) by Sharma et al. [34]; inverted Weibull (IW( α , μ )) by Keller et al. [35]; inverted gamma (IG( α , μ )) by Glen [36]; and inverted Nadarajah–Haghighi (INH( α , μ )) by Tahir et al. [37]. The comparison of fits of the alpha-PIE (or simply, APIE) and other competitive models is made using several statistics, called negative log-likelihood (N-L), Akaike (A), Bayesian (B), consistent Akaike (C-A), Hannan-Quinn (H-Q), and Kolmogorov–Smirnov (K-S) statistic with its p-value. These measures are evaluated through the maximum likelihood estimates of  α  and  μ . Using Table 3, the MLEs (along with associated standard-errors (St.Ers)) of  α  and  μ , as well as the fitted values of N-L, A, B, C-A, H-Q, and K-S(p-value) are obtained and reported in Table 4. This indicates that the Alpha-PIE distribution has the lowest values of all given goodness-of-fit statistics, except that it has the highest p-value. This conclusion implies that the Alpha-PIE lifetime model provides a best fit than other distributions.
Graphically, in Figure 8, the probability–probability (PP) plots of alpha-PIE and their competitive distributions are displayed. Figure 8 supports the same numerical findings listed in Table 4. Moreover, in Figure 9, three different plots are also considered, called (i) fitted densities with relative histograms of given data; (ii) estimated/empirical reliability functions of alpha-PIE, IE, IL, IW, IG, and INH distributions; and (iii) estimated/empirical total time on test (TTT) transform plots.
Each sub-plot in Figure 9 shows that the Alpha-PIE distribution provides a good fit compared to other competing models and indicates that the complete electronic tube data have an upside-down bathtub-shaped failure rate. In Figure 9, except for the alpha-PIE model, the fitted objects of other competitive models can be seen far away from the empirical ones. This observation is due to the fact that the estimated p-values of all compared models were lower than 0.5.
Additionally, to show the existence and uniqueness of  α ^  and  μ ^ , the contour of the log-likelihood function under complete electronic tube data with respect to various choices of  α  and  μ  are also plotted and displayed in Figure 9, which provides evidence that the MLEs of  α ^ 21.826  and  μ ^ 0.4149  exist and are unique; in addition, they are considered starting points for forthcoming numerical evaluations.
To illustrate our acquired estimates of  α μ R ( t ) , and  h ( t )  from the total data points of electronic tubes, various artificial GPHC-T-II samples of size  m = 10  are obtained based on different choices of times  T i , i = 1 , 2  and censoring plans  R i , i = 1 , 2 , , m , namely,  Sch [ 1 ] : ( 1 10 ) Sch [ 2 ] : ( 2 4 , 0 5 , 2 ) , and  Sch [ 3 ] : ( 2 , 0 5 , 2 4 ) ; see Table 5. Here, for brevity, the censoring plan  ( 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 )  is referred to as  ( 1 10 ) . For each simulated sample in Table 5, Table 6 represents the point estimates (along their St.Ers) and the interval estimates (along their interval widths (IWs)) of  α μ R ( t ) , and  h ( t )  at time  t = 1 . Since we lack the prior information on  α  and  μ  from the electronic tubes data set, the improper gamma priors (where  a i = b i = 0 , i = 1 , 2 ) are used. Then, the Bayes estimates, along with their HPD intervals, are developed by repeating the MCMC technique 50,000 times and eliminating the first 10,000 times as burn-in. It can be seen, from Table 6, that the approximate Bayes estimates of  α μ R ( t )  and  h ( t )  created from the MCMC sampler behave better than others in terms of the lowest St.Er and IW values.
To emphasize the convergence of MCMC sequences of  α μ R ( t )  and  h ( t ) , from each sample  S 1  generated by  Sch [ i ] , i = 1 , 2 , 3  (as examples), Figure 10 displays the density and trace diagrams from 40,000 MCMC draws of  α μ R ( t )  and  h ( t ) . To distinguish THEM, for each plot in Figure 10, the Bayes estimate (posterior mean) and two 95% HPD interval bounds are represented by solid and dashed lines, respectively. It is clear, from Figure 10, that the MCMC strategy converges well and that the recommended size of the burn-in sample is adequate to ignore the impact of the starting points. The figure also demonstrates that the drawn estimates of  α  and  μ  are fairly symmetrical. Figure 10 also shows that the drawn estimates of  R ( t )  and  h ( t )  are negatively and positively skewed, respectively.

5.2. Vinyl Chloride

Vinyl chloride is an organochloride, flammable gas with a sickly sweet odor, which is colorless at ambient temperature and is a recognized human carcinogen that burns readily. In 1835, it was first discovered by Justus von Liebig and Henri Victor Regnault. It is not produced naturally and must be produced industrially for its commercial uses, for example, pipes, packaging materials, and coatings for wire and cable; see Ware [38]. In this application, from clean up-gradient monitoring wells, we shall consider a data set consisting of thirty-four data points of vinyl chloride; see Table 7. This data set was given by Bhaumik et al. [39] and later, Elshahhat and Elemary [40] as well as Alotaibi et al. [41] also discussed it.
To monitor the adaptability and flexibility of the alpha-PIE distribution, from the complete vinyl chloride data, we compare its fits with several other competitive models reported in Section 5.1. For this purpose, the MLEs (with their St.Ers) of  α  and  μ , N-L, A, B, C-A, H-Q, and K-S(p-value) of alpha-PIE and its competing distributions are calculated and provided in Table 8. Furthermore, the corresponding PP plot of the alpha-PIE distribution and its five competitive models are shown in Figure 11. In Figure 12, the fitted density and reliability function of the alpha-PIE distribution and all studied models, as well as the TTT transform plot, are displayed. Figure 11 and Figure 12 indicate confirmed the same results presented in Table 8. Moreover, Figure 12 shows that a decreasing failure rate is appropriate for modeling vinyl chloride data. In addition, Figure 12 shows that the MLEs  α ^ 21.364  and  μ ^ 0.2524  developed from the vinyl chloride data exist and are unique. Therefore, in the forthcoming classical (or Bayesian) evaluations, these values are taken as initial guesses.
To assess the efficacy of the suggested inferential techniques in the chemical area, from the complete vinyl chloride data, several artificial GPHC-T-II samples (when  m = 17  and different choices of  T i , i = 1 , 2  and  R i , i = 1 , 2 , , m  namely:  Sch [ 1 ] : ( 1 17 ) Sch [ 2 ] : ( 2 8 , 0 8 , 1 )  and  Sch [ 3 ] : ( 1 , 0 8 , 2 8 ) ) are created and listed in Table 9. From each GPHC-T-II sample reported Table 9, the maximum likelihood and Bayes MCMC estimates (along with their St.Ers) and the asymptotic and HPD interval estimates (along with their IWs) of  α μ R ( t )  or  h ( t )  (at  t = 0.5 ) are computed and provided in Table 10.
By ignoring the first 10,000 iterations (burn-in) from the full MCMC 50,000 iterations, the Bayes’ estimates are approximated when the hyperparameter values are selected to be zero. Because we lacked additional historical information from the given data, which led to little to no difference between the proposed classical and Bayesian estimates, Table 10 stated that the offered estimates of  α μ R ( t )  or  h ( t )  behave similarly and seem to be close to one another.
From each sample  S 1  obtained by  Sch [ i ] , i = 1 , 2 , 3  (as examples), to display the convergence of 40,000 MCMC draws of  α μ R ( t )  and  h ( t ) , the density and trace pots of all unknown quantities are plotted and shown in Figure 13. This figure supports the same numerical findings reported in Table 10 and indicates that all MCMC iterations of  α μ R ( t )  or  h ( t )  converged satisfactorily. It also shows that the simulated MCMC estimates of  α  and  μ  are almost symmetrical, while those of  R ( t )  and  h ( t )  are negatively and positively skewed, respectively.
As a result, the reported discussions of real data sets representing 20 electron tubes and 34 vinyl chloride tubes, show that the derived estimators perform well in practical situations and proposed lifetime distribution is a suitable choice for modeling these data in the presence of Type-II generalized progressive hybrid censored data.

6. Optimal Progressive Designs

In regard to the context of reliability, the experimenter can decide to select the most effective censoring strategy from a group of all available progressive designs in order to offer as much information as possible about the parameter(s) under investigation. Independently, this issue was first addressed by Balakrishnan and Aggarwala [42] and Ng et al. [43]. The ideal censoring fashion  R = ( R 1 , R 2 , , R m )  such as  j = 1 m R j = n m , can be proposed, and the plausible choices of n, m and  T i , i = 1 , 2  are predetermined in advance depending on unit capacity, experimentation facilities, and budgetary constraints. In the literature, several criteria, however, have been introduced, and numerous results on the optimum censoring systems have been investigated; for examples, see Pradhan and Kundu [44]; Sen et al. [45]; Elshahhat and Rastogi [46]; Ashour et al. [47]; and Elshahhat and Abu El Azm [48]. In Table 11, a list of commonly used metrics to help us choose the best censoring strategy is presented.
It should be noted that the objective of the studied criteria  Q i , i = 2 , 3 , 4  is to minimize the trace, determinant, and variance of the logarithmic MLE of the  ϱ th quantile of the estimated variance–covariance  I 1 ( · )  matrix, while the objective of criterion  Q 1  is to maximize the observed Fisher information  F i i , i = 1 , 2  elements.
Obviously, the best censoring pattern must correspond to the largest value of  Q 1  and the smallest value of  Q i  for  i = 2 , 3 , 4 . Specifically, from (2), the logarithmic of the alpha-PIE lifetime distribution  T ϱ  is given by
log ( T ^ ϱ ) = log μ log { log 1 ( α ) log [ ϱ ( α 1 ) + 1 ] } 1 ( α ^ , μ ^ ) , 0 < ϱ < 1 .
Again, utilizing the delta method, the approximated variance of  log ( T ^ ϱ )  is given by
v ^ ( log ( T ^ ϱ ) ) = Σ log ( T ϱ ) T I 1 ( ϑ ) Σ log ( T ϱ ) ( α ^ , μ ^ ) ,
where
Σ log ( T ^ ϱ ) T = α log ( T ϱ ) log ( T ϱ ) , μ log ( T ϱ ) ( α ^ , μ ^ ) .

6.1. Optimum from Electronic Tubes

This subsection aims to determine the best progressive censoring based on the generated samples, which are created from the electronic tube data, reported in Table 5. However, from Table 5, Table 6 and Table 11, the optimum criteria are evaluated; see Table 12.
Table 12 supports the same recommended censoring schemes considered in Section 4 and shows that
  • According to  O i , i = 1 , 2 , 3 , the design of Sch[2] (in Sample  S 1 ) and the design of Sch[3] (in Samples  S 2  and  S 3 ) are the optimum censoring plans compared to others.
  • According to  O 4 , the design of Sch[1] (in Sample  S 1 ), the design of Sch[3] (in Sample  S 2 ) and the design of Sch[2] (in Sample  S 3 ) are the optimum censoring plans compared to others.

6.2. Optimum from Vinyl Chloride

In this subsection, from the vinyl chloride data, we shall propose an optimal progressive censored plan based on the generated samples reported in Table 9. However, from Table 9, Table 10 and Table 11, the optimum criteria are evaluated and presented in Table 13.
Table 13 shows that
  • According to  O 1 , the design of Sch[1] (in Sample  S 1 ) and the design of Sch[2] (in Samples  S i , i = 2 , 3 ) are the optimum censoring plans compared to others.
  • According to  O i , i = 2 , 3 , the design of Sch[1] (in Samples  S i , i = 1 , 3 ) are the optimum censoring plans compared to others.
  • According to  O i , i = 2 , 3 , the designs Sch[3] and Sch[2] (in Sample  S 2 ), respectively, are the optimum censoring plans compared to others.
  • According to  O 4 , the design of Sch[3] (in Samples  S i , i = 1 , 2 ) and the design of Sch[2] (in Samples  S 3 ) are the optimum censoring plans compared to others.
Finally, it is clear that the optimum progressive censoring plans suggested in this section support the same findings reported in Section 4.

7. Concluding Remarks

This work considers the problem of statistical inference of new alpha power-inverted exponential parameters of life using generalized Type-II progressively hybrid censored data. Classically, the maximum likelihood and asymptotic confidence interval estimates of the model parameters and any related time function have been derived using Newton–Raphson optimization process via the “ maxLik ” language. The joint posterior density has been derived in a nonlinear form because the formula of the likelihood function has a complex expression. Independent gamma priors, to derive the Bayes and HPD interval estimates, are considered. Metropolis–Hastings sampler via the “ coda ” language is also recommended to simulate the MCMC samples of the same unknown parameters. Extensive simulation experiments, based on various options of n, m R , and  T i , i = 1 , 2 , have been conducted to judge the behavior of the offered estimates. These studies show that the MCMC methodology performs quite well compared to the maximum likelihood approach. An optimum progressive censoring has also been presented using several optimality criterion measurements. Two scenarios based on real-world data sets from the engineering and chemical sectors are examined to highlight the superiority of the proposed model and how the provided estimates can be performed in practice. These applications show that the derived estimators perform well in practical situations and that the proposed lifetime distribution is a good choice compared to several models in the literature, namely, the inverted exponential, inverted Lindley, inverted Weibull, inverted gamma, and inverted Nadarajah–Haghighi distributions. We also believe that the findings and methodology presented here will be useful to reliability technicians, statisticians, and/or other scientists.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/axioms12060601/s1, Table S1: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  α  from Set-1; Table S2: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  α  from Set-2; Table S3: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  μ  from Set-1; Table S4: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  μ  from Set-2; Table S5: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  R ( t )  from Set-1; Table S6: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  R ( t )  from Set-2; Table S7: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  h ( t )  from Set-1; Table S8: The Av.Es (1 s t  column), RMSEs (2 n d  column) and MRABs (3 r d  column) of  h ( t )  from Set-2; Table S9: The ACLs (1 s t  column) and CPs (2 n d  column) of 95% ACI/HPD intervals of  α ; Table S10: The ACLs (1 s t  column) and CPs (2 n d  column) of 95% ACI/HPD intervals of  μ ; Table S11: The ACLs (1 s t  column) and CPs (2 n d  column) of 95% ACI/HPD intervals of  R ( t ) ; Table S12: The ACLs (1 s t  column) and CPs (2 n d  column) of 95% ACI/HPD intervals of  h ( t ) .

Author Contributions

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

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Differentiating (7) with regard to  α  and  μ , Fisher’s items  F i j , i , j = 1 , 2  are
F 11 = i = 1 D q ( R i [ ( ( α 1 ) ( 1 α 1 ( α exp ( ( μ x i 1 ) ) 1 ) ) ) 1 ( α e μ x i 1 2 e μ x i 1 [ e μ x i 1 1 ] ( α 1 ) { ( α exp ( ( μ x i 1 ) ) 1 ) e μ x i 1 α 1 ( α exp ( ( μ x i 1 ) ) 1 ) } ) ( ( α 1 ) ( 1 α 1 ( α exp ( ( μ x i 1 ) ) 1 ) ) ) 2 ( 1 α exp ( ( μ x i 1 ) ) e μ x i 1 { ( α exp ( ( μ x i 1 ) ) 1 ) e μ x i 1 α 1 ( α exp ( ( μ x i 1 ) ) 1 ) } ) ] + i = 1 D q ( ( ( α exp ( ( μ x i 1 ) ) 1 α 1 α ) 1 log ( α ) + α exp ( ( μ x i 1 ) ) 2 + α exp ( ( μ x i 1 ) ) 2 [ e μ x i 1 1 ] × log ( α ) ) e μ x i 1 α exp ( ( μ x i 1 ) ) ( ( α 1 ) α 1 α ) 1 log ( α ) + α 2 ) ) [ α exp ( ( μ x i 1 ) ) log ( α ) ] 1 ( ( α exp ( ( μ x i 1 ) ) 1 ) e μ x i 1 × log ( α ) + α exp ( ( μ x i 1 ) ) α 1 α ) 1 log ( α ) ) ( α exp ( ( μ x i 1 ) ) 1 + ( α exp ( ( μ x i 1 ) ) 1 ) e μ x i 1 log ( α ) ) [ α exp ( ( μ x i 1 ) ) log ( α ) ] 2 ) R d τ + 1 * ( ( α ( exp ( ( μ T τ 1 ) ) 2 ) e μ T τ 1 ( e μ T τ 1 1 ) ( α 1 ) { α exp ( ( μ T τ 1 ) ) 1 e μ T τ 1 α 1 ( α exp ( ( μ T τ 1 ) ) 1 ) } ) × [ ( α 1 ) ( 1 α 1 ( α exp ( ( μ T τ 1 ) ) 1 ) ) ] 1 ( 1 α exp ( ( μ T τ 1 ) ) 1 e μ T τ 1 { α exp ( ( μ T τ 1 ) ) 1 e μ T τ 1 α 1 ( α exp ( ( μ T τ 1 ) ) 1 ) } [ ( α 1 ) ( 1 α 1 ( α exp ( ( μ T τ 1 ) ) 1 ) ) ] 2 ) ,
F 22 = i = 1 D q [ ( ( x i μ α exp ( μ x i 1 ) ) 1 ( α exp ( μ x i 1 ) + α exp ( μ x i 1 ) ( 1 μ x i 1 ) μ x i 1 e μ x i 1 log ( α ) α exp ( μ x i 1 ) ( 1 + α exp ( μ x i 1 ) ) ) ×   log ( α ) e μ x i 1 + α exp ( μ x i 1 ) ) + ( α exp ( μ x i 1 ) μ α exp ( μ x i 1 ) e μ x i 1 x i 1 log ( α ) ) ( α exp ( μ x i 1 ) ( 1 μ x i 1 ) μ α exp ( μ x i 1 ) × x i 1 e μ x i 1 log ( α ) ) ( μ α exp ( μ x i 1 ) ) 2 ] i = 1 D q [ R i ( α exp ( μ x i 1 ) ( 1 + α exp ( μ x i 1 ) e μ x i 1 log ( α ) ) [ x i 2 ( 1 ( α 1 ) 1 × ( α exp ( μ x i 1 ) 1 ) ) ( α 1 ) ] 1 + α ( 2 e μ x i 1 ) e μ x i 1 log ( α ) ( x i ( 1 ( α 1 ) 1 ( α exp ( μ x i 1 ) 1 ) ) ( α 1 ) ) 2 ) e μ x i 1 log ( α ) ] R d τ + 1 * e μ T τ 1 log ( α ) [ α exp ( ( μ T τ 1 ) ) ( 1 + α exp ( ( μ T τ 1 ) ) exp ( ( μ T τ 1 ) ) log ( α ) ) ( T τ 2 ( 1 ( α 1 ) 1 ( α exp ( ( μ T τ 1 ) ) 1 ) ) × ( α 1 ) ) 1 + α 2 e μ T τ 1 e μ T τ 1 log ( α ) ( T τ ( 1 ( α 1 ) 1 ( α exp ( ( μ T τ 1 ) ) 1 ) ) ( α 1 ) ) 2 ] ,
and
F 12 = i = 1 D q [ R i ( ( ( ( α 1 ) 1 e μ x i 1 ) α exp ( μ x i 1 ) log ( α ) α exp ( μ x i 1 ) 1 ) [ ( α 1 ) ( 1 ( α 1 ) 1 ( α exp ( μ x i 1 ) 1 ) ) ] 1 α exp ( μ x i 1 ) [ α exp ( μ x i 1 ) 1 e μ x i 1 ( α 1 ) 1 ( α exp ( μ x i 1 ) 1 ) ] log ( α ) [ ( α 1 ) ( 1 ( α 1 ) 1 ( α exp ( μ x i 1 ) 1 ) ) ] 2 ) ×   x i 1 e μ x i 1 ] + i = 1 D q ( ( α exp ( μ x i 1 ) ( α exp ( μ x i 1 ) 1 e μ x i 1 log ( α ) + α exp ( μ x i 1 ) [ α 1 ( α 1 ) 1 log ( α ) ] ) ( log ( α ) ) 2 ×   ( α exp ( μ x i 1 ) log ( α ) ) 2 ( α exp ( μ x i 1 ) 1 + α exp ( μ x i 1 ) 1 e μ x i 1 log ( α ) + α exp ( μ x i 1 ) [ α 1 ( α 1 ) 1 log ( α ) ] ) ×   α exp ( μ x i 1 ) ) x i 1 e μ x i 1 ) R d τ + 1 * ( ( ( ( α 1 ) 1 α exp ( μ T τ 1 ) α exp ( μ T τ 1 ) 1 e μ T τ 1 ) log ( α ) α exp ( μ T τ 1 ) 1 ) ×   [ ( α 1 ) ( 1 ( α 1 ) 1 ( α exp ( μ T τ 1 ) 1 ) ) ] 1 α exp ( μ T τ 1 ) [ α exp ( μ T τ 1 ) 1 e μ T τ 1 ( α 1 ) 1 ( α exp ( μ T τ 1 ) 1 ) ] × log ( α ) [ ( α 1 ) ( 1 ( α 1 ) 1 ( α exp ( μ T τ 1 ) 1 ) ) ] 2 ) T τ 1 e μ T τ 1 .

References

  1. Ünal, C.; Cakmakyapan, S.; Özel, G. Alpha power inverted exponential distribution: Properties and application. Gazi Univ. J. Sci. 2018, 31, 954–965. [Google Scholar]
  2. Keller, A.Z.; Kamath, A.R.R.; Perera, U.D. Reliability analysis of CNC machine tools. Reliab. Eng. 1982, 3, 449–473. [Google Scholar] [CrossRef]
  3. Amjad, N.; Taufique, M.; Saud, N. Bayes estimation of the alpha power inverted exponential parameters under various approximation techniques. Pak. J. Stat. 2022, 38, 99–112. [Google Scholar]
  4. Balakrishnan, N.; Cramer, E. The Art of Progressive Censoring; Springer: Berlin/Heidelberg, Germany; Birkhäuser: New York, NY, USA, 2014. [Google Scholar]
  5. Lee, K.; Sun, H.; Cho, Y. Exact likelihood inference of the exponential parameter under generalized Type II progressive hybrid censoring. J. Korean Stat. Soc. 2016, 45, 123–136. [Google Scholar] [CrossRef]
  6. Childs, A.; Chandrasekar, B.; Balakrishnan, N. Exact Likelihood Inference for an Exponential Parameter under Progressive Hybrid Censoring Schemes. In Statistical Models and Methods for Biomedical and Technical Systems; Vonta, F., Nikulin, M., Limnios, N., Huber-Carol, C., Eds.; Birkhäuser: Boston, MA, USA, 2008; pp. 319–330. [Google Scholar]
  7. Kundu, D.; Joarder, A. Analysis of Type-II progressively hybrid censored data. Comput. Stat. Data Anal. 2006, 50, 2509–2528. [Google Scholar] [CrossRef]
  8. Epstein, B. Truncated life tests in the exponential case. Ann. Math. Stat. 1954, 25, 555–564. [Google Scholar] [CrossRef]
  9. Childs, A.; Chandrasekar, B.; Balakrishnan, N.; Kundu, D. Exact likelihood inference based on Type-I and Type-II hybrid censored samples from the exponential distribution. Ann. Inst. Stat. Math. 2003, 55, 319–330. [Google Scholar] [CrossRef]
  10. Epstein, B.; Sobel, M. Life testing. J. Am. Stat. Assoc. 1953, 48, 486–502. [Google Scholar] [CrossRef]
  11. Ashour, S.; Elshahhat, A. Bayesian and non-Bayesian estimation for Weibull parameters based on generalized Type-II progressive hybrid censoring scheme. Pak. J. Stat. Oper. Res. 2016, 12, 213–226. [Google Scholar]
  12. Ateya, S.; Mohammed, H. Prediction under Burr-XII distribution based on generalized Type-II progressive hybrid censoring scheme. J. Egypt. Math. Soc. 2018, 26, 491–508. [Google Scholar]
  13. Seo, J.I. Objective Bayesian analysis for the Weibull distribution with partial information under the generalized Type-II progressive hybrid censoring scheme. Commun. Stat.-Simul. Comput. 2020, 51, 5157–5173. [Google Scholar] [CrossRef]
  14. Cho, S.; Lee, K. Exact Likelihood Inference for a Competing Risks Model with Generalized Type II Progressive Hybrid Censored Exponential Data. Symmetry 2021, 13, 887. [Google Scholar] [CrossRef]
  15. Nagy, M.; Bakr, M.E.; Alrasheedi, A.F. Analysis with applications of the generalized Type-II progressive hybrid censoring sample from Burr Type-XII model. Math. Probl. Eng. 2022, 2022, 1241303. [Google Scholar] [CrossRef]
  16. Wang, L.; Zhou, Y.; Lio, Y.; Tripathi, Y.M. Inference for Kumaraswamy Distribution under Generalized Progressive Hybrid Censoring. Symmetry 2022, 14, 403. [Google Scholar] [CrossRef]
  17. Elshahhat, A.; Mohammed, H.S.; Abo-Kasem, O.E. Reliability Inferences of the Inverted NH Parameters via Generalized Type-II Progressive Hybrid Censoring with Applications. Symmetry 2022, 14, 2379. [Google Scholar] [CrossRef]
  18. 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]
  19. Henningsen, A.; Toomet, O. maxLik: A package for maximum likelihood estimation in R. Comput. Stat. 2011, 26, 443–458. [Google Scholar] [CrossRef]
  20. Lawless, J.F. Statistical Models and Methods for Lifetime Data, 2nd ed.; John Wiley and Sons: Hoboken, NJ, USA, 2003. [Google Scholar]
  21. Greene, W.H. Econometric Analysis, 4th ed.; Prentice-Hall: Scottsville, NY, USA, 2000. [Google Scholar]
  22. Chen, P.; Ye, Z.S. Estimation of field reliability based on aggregate lifetime data. Technometrics 2017, 59, 115–125. [Google Scholar] [CrossRef]
  23. Wang, P.; Tang, Y.; Bae, S.J.; He, Y. Bayesian analysis of two-phase degradation data based on change-point Wiener process. Reliab. Eng. Syst. Saf. 2018, 170, 244–256. [Google Scholar] [CrossRef]
  24. Luo, C.; Shen, L.; Xu, A. Modelling and estimation of system reliability under dynamic operating environments and lifetime ordering constraints. Reliab. Eng. Syst. Saf. 2022, 218, 108136. [Google Scholar] [CrossRef]
  25. Luo, C.; Xu, L. Online-to-offline on the railway: Optimization of on-demand meal ordering on high-speed railway. Transp. Res. Part C Emerg. Technol. 2023, 152, 104143. [Google Scholar] [CrossRef]
  26. 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]
  27. Lynch, S.M. Introduction to Applied Bayesian Statistics and Estimation for Social Scientists; Springer: New York, NY, USA, 2007. [Google Scholar]
  28. 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]
  29. Plummer, M.; Best, N.; Cowles, K.; Vines, K. CODA: Convergence diagnosis and output analysis for MCMC. R News 2006, 6, 7–11. [Google Scholar]
  30. 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]
  31. Rosebury, F. Handbook of Electron Tube and Vacuum Techniques; American Institute of Physics Melville: New York, NY, USA, 1993. [Google Scholar]
  32. Dixit, U.J.; Nooghabi, M.J. Estimation of parameters of gamma distribution in the presence of outliers in right censored samples. Aligar J. Stat. 2011, 31, 17–29. [Google Scholar]
  33. Ibrahim, H.A.; Mahmoud, M.R.; Khalil, F.A.; El-Kelany, G.A. TL-moments for Type-I Censored Data with an Application to the Weibull Distribution. Math. Comput. Appl. 2018, 23, 47. [Google Scholar] [CrossRef] [Green Version]
  34. Sharma, V.K.; Singh, S.K.; Singh, U.; Agiwal, V. The inverse Lindley distribution: A stress-strength reliability model with application to head and neck cancer data. J. Ind. Prod. Eng. 2015, 32, 162–173. [Google Scholar] [CrossRef]
  35. Keller, A.Z.; Goblin, M.T.; Farnworth, N.R. Reliability analysis of commercial vehicle engines. Reliab. Eng. 1985, 10, 15–25. [Google Scholar] [CrossRef]
  36. Glen, A. On the inverse gamma as a survival distribution. J. Qual. Technol. 2011, 43, 158–166. [Google Scholar] [CrossRef]
  37. Tahir, M.H.; Cordeiro, G.M.; Ali, S.; Dey, S.; Manzoor, A. The inverted Nadarajah-Haghighi distribution: Estimation methods and applications. J. Stat. Comput. Simul. 2018, 88, 2775–2798. [Google Scholar] [CrossRef]
  38. Ware, G.W. Vinyl Chloride. In Reviews of Environmental Contamination and Toxicology; Ware, G.W., Ed.; Springer: New York, NY, USA, 1988; Volume 107. [Google Scholar]
  39. Bhaumik, D.K.; Kapur, K.; Gibbons, R.D. Testing parameters of a gamma distribution for small samples. Technometrics 2009, 51, 326–334. [Google Scholar] [CrossRef]
  40. Elshahhat, A.; Elemary, B.R. Analysis for Xgamma parameters of life under Type-II adaptive progressively hybrid censoring with applications in engineering and chemistry. Symmetry 2021, 13, 2112. [Google Scholar] [CrossRef]
  41. 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]
  42. Balakrishnan, N.; Aggarwala, R. Progressive Censoring Theory, Methods and Applications; Birkhäuser: Boston, MA, USA, 2000. [Google Scholar]
  43. 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]
  44. 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]
  45. 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]
  46. Elshahhat, A.; Rastogi, M.K. Estimation of parameters of life for an inverted Nadarajah–Haghighi distribution from Type-II progressively censored samples. J. Indian Soc. Probab. Stat. 2021, 22, 113–154. [Google Scholar] [CrossRef]
  47. Ashour, S.K.; El-Sheikh, A.A.; Elshahhat, A. Inferences and optimal censoring schemes for progressively first-failure censored Nadarajah-Haghighi distribution. Sankhyā A Indian J. Stat. 2022, 84, 885–923. [Google Scholar] [CrossRef]
  48. 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. Plots for the PDF (a) and HRF (b) of the alpha-PIE distribution when  μ = 1 .
Figure 1. Plots for the PDF (a) and HRF (b) of the alpha-PIE distribution when  μ = 1 .
Axioms 12 00601 g001
Figure 2. Diagram of GTIIPHC sampling.
Figure 2. Diagram of GTIIPHC sampling.
Axioms 12 00601 g002
Figure 3. Autocorrelation (top) and trace (bottom) plots for MCMC draws of  α μ R ( t )  and  h ( t ) .
Figure 3. Autocorrelation (top) and trace (bottom) plots for MCMC draws of  α μ R ( t )  and  h ( t ) .
Axioms 12 00601 g003
Figure 4. Heat-map for the Monte Carlo outputs of  α . (a) Set 1; (b) Set 2.
Figure 4. Heat-map for the Monte Carlo outputs of  α . (a) Set 1; (b) Set 2.
Axioms 12 00601 g004
Figure 5. Heat-map for the Monte Carlo outputs of  μ . (a) Set 1; (b) Set 2.
Figure 5. Heat-map for the Monte Carlo outputs of  μ . (a) Set 1; (b) Set 2.
Axioms 12 00601 g005
Figure 6. Heat-map for the Monte Carlo outputs of  R ( t ) . (a) Set 1; (b) Set 2.
Figure 6. Heat-map for the Monte Carlo outputs of  R ( t ) . (a) Set 1; (b) Set 2.
Axioms 12 00601 g006
Figure 7. Heat-map for the Monte Carlo outputs of  h ( t ) . (a) Set 1; (b) Set 2.
Figure 7. Heat-map for the Monte Carlo outputs of  h ( t ) . (a) Set 1; (b) Set 2.
Axioms 12 00601 g007
Figure 8. The PP plots of the alpha-PIE and its competitive models from electronic tube data.
Figure 8. The PP plots of the alpha-PIE and its competitive models from electronic tube data.
Axioms 12 00601 g008
Figure 9. Fitted PDFs (a); fitted RFs (b); scaled TTT (c); contour (d); plots from electronic tubes data.
Figure 9. Fitted PDFs (a); fitted RFs (b); scaled TTT (c); contour (d); plots from electronic tubes data.
Axioms 12 00601 g009
Figure 10. Density (left) and trace (right) plots of  α μ R ( t )  and  h ( t )  from electronic tubes data.
Figure 10. Density (left) and trace (right) plots of  α μ R ( t )  and  h ( t )  from electronic tubes data.
Axioms 12 00601 g010
Figure 11. The PP plots of the alpha-PIE and its competitive distributions from vinyl chloride data.
Figure 11. The PP plots of the alpha-PIE and its competitive distributions from vinyl chloride data.
Axioms 12 00601 g011
Figure 12. Fitted PDFs (a); fitted RFs (b); scaled-TTT (c); contour (d) plots from vinyl chloride data.
Figure 12. Fitted PDFs (a); fitted RFs (b); scaled-TTT (c); contour (d) plots from vinyl chloride data.
Axioms 12 00601 g012
Figure 13. Density (left) and Trace (right) plots of  α μ R ( t )  and  h ( t )  from vinyl chloride data.
Figure 13. Density (left) and Trace (right) plots of  α μ R ( t )  and  h ( t )  from vinyl chloride data.
Axioms 12 00601 g013
Table 1. The GTIIPHC notations.
Table 1. The GTIIPHC notations.
q K q D q q ( T τ ; ϑ ) R d τ + 1 *
1 Π i = 1 d 1 k = i m R k + 1 d 1 1 F ( T 1 ) R d 1 + 1 * n d 1 k = 1 m 1 R k
2 Π i = 1 m k = i m R k + 1 m10
3 Π i = 1 d 2 k = i m R k + 1 d 2 1 F ( T 2 ) R d 2 + 1 * n d 2 k = 1 d 2 R k
Table 2. Special sampling plans from GTIIPHC.
Table 2. Special sampling plans from GTIIPHC.
PlanAuthor(s)Setting
Type-I PHCKundu and Joarder [7] T 1 0
Type-II PHCChilds et al. [6] T 2
Type-I HybridEpstein [8] T 1 0 R i = 0 , i = 1 ,   2 ,   ,   m 1 ,  and  R m = n m
Type-II HybridChilds et al. [9] T 2 R i = 0 , i = 1 ,   2 ,   ,   m 1 ,  and  R m = n m
Type-I censoringEpstein and Sobel [10] T 1 = 0 m = n R i = 0 , i = 1 ,   2 ,   ,   m 1 ,  and  R m = n m
Type-II censoringEpstein and Sobel [10] T 1 = 0 T 2 R i = 0 , i = 1 ,   2 ,   ,   m 1 ,  and  R m = n m
Table 3. Lifetimes of 20 electronic tubes.
Table 3. Lifetimes of 20 electronic tubes.
0.14150.34840.39940.41740.59371.10451.73231.83482.34672.4651
2.61552.74253.13563.22593.41773.55513.56813.72879.28179.3208
Table 4. Outputs of the fitting alpha-PIE and its competitive models from electronic tube data.
Table 4. Outputs of the fitting alpha-PIE and its competitive models from electronic tube data.
ModelMLE(St.Er)N-LABC-AH-QK-S (p-Value)
α μ
APIE21.826 (55.135)0.4149 (0.2867)43.379890.759692.751191.465591.14840.1962 (0.375)
IE-0.9052 (0.2024)45.028292.056493.052292.278792.25080.2930 (0.051)
IL-1.2990 (0.2240)46.411894.823695.819495.045995.01800.3265 (0.0212)
IW0.8531 (0.1361)1.0019 (0.2383)44.484792.969494.960993.675393.35810.2342 (0.1897)
IG0.8692 (0.2390)0.7869 (0.2871)44.893293.786395.777894.492294.17510.2677 (0.0935)
INH0.5622 (0.1660)2.7258 (1.7357)43.516791.033593.024991.739391.42220.2033 (0.3340)
Table 5. Artificial GPHC-T-II samples from electronic tube data.
Table 5. Artificial GPHC-T-II samples from electronic tube data.
SchemeSample T 1 ( d 1 ) T 2 ( d 2 ) Generated Data R * T * *
Sch[1] S 1 3.8 (10)4.5 (10)0.1415, 0.3994, 0.4174, 0.5937, 1.1045, 1.8348, 2.4651, 2.7425, 3.2259, 3.417713.8
S 2 2.5 (6)3.8 (10)0.1415, 0.4174, 0.5937, 1.7323, 2.3467, 2.4651, 2.6155, 3.2259, 3.5551, 3.568103.5681
S 3 1.5 (4)3.3 (9)0.1415, 0.4174, 0.5937, 1.1045, 1.7323, 1.8348, 2.6155, 2.7425, 3.225923.3
Sch[2] S 1 2.8 (10)3.5 (10)0.1415, 0.4174, 0.5937, 1.1045, 1.7323, 1.8348, 2.3467, 2.4651, 2.6155, 2.742522.8
S 2 1.8 (5)3.6 (10)0.1415, 0.3484, 0.5937, 1.1045, 1.7323, 2.3467, 2.4651, 2.6155, 3.1356, 3.555103.5551
S 3 1.4 (4)3.6 (9)0.1415, 0.3994, 0.4174, 1.1045, 1.8348, 2.3467, 2.4651, 3.4177, 3.555133.6
Sch[3] S 1 2.8 (10)3.2 (10)0.1415, 0.3994, 0.4174, 0.5937, 1.1045, 1.7323, 1.8348, 2.3467, 2.6155, 2.742522.8
S 2 1.2 (6)3.4 (10)0.1415, 0.3484, 0.3994, 0.4174, 0.5937, 1.1045, 1.7323, 1.8348, 2.6155, 3.225903.2259
S 3 1.9 (6)3.5 (9)0.1415, 0.4174, 0.5937, 1.1045, 1.7323, 1.8348, 2.3467, 2.7425, 3.417733.5
Table 6. Bayesian and classical estimates of  α μ R ( t )  and  h ( t )  from electronic tube data.
Table 6. Bayesian and classical estimates of  α μ R ( t )  and  h ( t )  from electronic tube data.
SchemeSamplePar.MLEMCMCACIHPD
Est.St.ErEst.St.ErLowerUpperIWLowerUpperIW
Sch[1] S 1 α 28.05419.39527.9530.14170.000066.06766.06727.76028.1490.3884
μ 0.43910.15200.37420.09710.14110.73700.59600.23490.51540.2805
R ( 0.5 ) 0.88920.05420.85000.05890.78310.99540.21230.75990.92470.1649
h ( 0.5 ) 0.40430.15780.50980.15660.09490.71370.61870.30980.75220.4424
S 2 α 55.49617.19655.3950.143121.79289.20067.40855.19355.5880.3945
μ 0.50420.15090.43460.10380.20830.80000.59170.28470.58460.2999
R ( 0.5 ) 0.93890.03440.91380.03830.87160.99230.12070.85890.96170.1029
h ( 0.5 ) 0.25000.11790.32950.11970.01900.48100.46200.17700.50370.3267
S 3 α 48.97617.34948.8750.141914.97382.97968.00548.68249.0730.3905
μ 0.46760.14250.39950.10070.18840.74680.55840.26120.54770.2865
R ( 0.5 ) 0.92480.04050.89420.04600.84540.99470.14930.82530.95210.1267
h ( 0.5 ) 0.29650.13240.38880.13670.03710.55590.51880.21380.59460.3808
Sch[2] S 1 α 45.39817.54345.2970.143011.01579.78268.76745.09545.4880.3926
μ 0.38080.12210.32540.08770.14160.62000.47840.20050.46100.2605
R ( 0.5 ) 0.88870.05570.84830.06300.77960.99790.21840.75600.93280.1768
h ( 0.5 ) 0.40840.16570.51780.16890.08370.73320.64940.29370.77860.4849
S 2 α 36.00817.98335.9090.13990.761071.25470.49335.71236.1050.3924
μ 0.41450.13780.35190.09500.14450.68450.54000.21730.49020.2729
R ( 0.5 ) 0.89200.05510.85160.06130.78410.99890.21480.75880.92690.1681
h ( 0.5 ) 0.39690.16330.50690.16490.07680.71710.64020.28900.74620.4572
S 3 α 30.43718.63830.3330.14470.000066.96666.96630.12530.5230.3977
μ 0.43780.15320.37140.09890.13770.73800.60040.22710.51760.2905
R ( 0.5 ) 0.89300.05550.85310.06000.78420.99190.20770.75960.92740.1678
h ( 0.5 ) 0.39330.16380.50160.16080.07230.71440.64210.29550.74800.4525
Sch[3] S 1 α 34.01318.49333.9130.13960.000070.25970.25933.72834.1090.3815
μ 0.48380.14960.41850.09870.19060.77700.58650.27250.56210.2896
R ( 0.5 ) 0.91460.04230.88450.04590.83160.99750.16590.81520.94300.1278
h ( 0.5 ) 0.32810.13210.41510.13120.06910.58710.51790.24170.61260.3709
S 2 α 17.02711.35316.9260.14260.000039.27939.27916.716117.1080.3920
μ 0.44270.15320.37870.09690.14250.74290.60040.24420.52590.2817
R ( 0.5 ) 0.86150.06160.81840.06490.74060.98230.24170.72350.90720.1838
h ( 0.5 ) 0.48310.16680.59070.16080.15620.80990.65370.37560.83650.4609
S 3 α 49.67712.06149.5780.138326.03873.31547.27749.39249.7670.3754
μ 0.50610.14990.43490.10530.21240.79990.58760.29040.59130.3009
R ( 0.5 ) 0.93560.03560.90920.04000.86590.99670.13080.84990.95830.1084
h ( 0.5 ) 0.26080.12010.34350.12330.02540.49630.47090.18960.52700.3374
Table 7. Data point (in milligrams/liter) of vinyl chloride.
Table 7. Data point (in milligrams/liter) of vinyl chloride.
0.10.10.20.20.40.40.40.50.50.50.60.60.80.90.91.01.1
1.21.21.31.82.02.02.32.42.52.72.93.24.05.15.36.88.0
Table 8. Outputs of the fitting alpha-PIE and its competitive models from vinyl chloride data.
Table 8. Outputs of the fitting alpha-PIE and its competitive models from vinyl chloride data.
ModelMLE(St.Er)N-LABC-AH-QK-S (p-Value)
α μ
APIE21.364 (52.796)0.2524 (0.1725)57.2457118.4915121.5442118.8786119.53250.0952 (0.918)
IE-0.5725 (0.0982)59.1930120.3860121.9124120.5110120.90660.1470 (0.454)
IL-0.8774 (0.1127)61.8136125.6272127.1535125.7522126.14770.1908 (0.168)
IW0.8805 (0.1093)0.6539 (0.1347)58.6266121.2532124.3059121.6403122.29420.1134 (0.774)
IG0.9002 (0.1904)0.5154 (0.1434)59.0659122.1319125.1846122.5190123.17290.1310 (0.604)
INH0.6089 (0.1518)1.4368 (0.7440)57.5539119.1079122.1606119.4950120.14900.1004 (0.883)
Table 9. Artificial GPHC-T-II samples from vinyl chloride data.
Table 9. Artificial GPHC-T-II samples from vinyl chloride data.
SchemeSample T 1 ( d 1 ) T 2 ( d 2 ) Generated Data R * T * *
Sch[1] S 1 4.2 (17)4.8 (17)0.1, 0.2, 0.4, 0.4, 0.5, 0.6, 0.8, 0.9, 0.9, 1.0, 1.1, 1.2, 1.3, 1.8, 2.3, 2.4, 4.014.2
S 2 1.5 (12)5.5 (17)0.1, 0.2, 0.4, 0.4, 0.4, 0.5, 0.6, 0.6, 0.9, 1.0, 1.1, 1.2, 2.0, 2.0, 2.5, 2.7, 5.105.1
S 3 1.1 (12)1.3 (15)0.1, 0.2, 0.4, 0.5, 0.5, 0.5, 0.6, 0.6, 0.8, 0.9, 0.9, 1.0, 1.1, 1.2, 1.241.3
Sch[2] S 1 2.2 (11)6.9 (17)0.1, 0.2, 0.4, 0.5, 0.5, 0.6, 0.8, 0.9, 1.2, 1.3, 2.0, 2.3, 2.7, 4.0, 5.1, 5.3, 6.816.9
S 2 1.7 (13)5.4 (17)0.1, 0.2, 0.5, 0.5, 0.6, 0.6, 0.8, 0.9, 0.9, 1.0, 1.1, 1.2, 1.2, 2.3, 2.4, 4.0, 5.305.3
S 3 1.1 (12)1.3 (14)0.1, 0.2, 0.4, 0.4, 0.4, 0.5, 0.6, 0.6, 0.8, 0.9, 0.9, 1.0, 1.1, 1.241.3
Sch[3] S 1 2.6 (17)2.8 (17)0.1, 0.2, 0.2, 0.4, 0.4, 0.4, 0.5, 0.5, 0.5, 0.6,0.9, 0.9, 1.2, 1.2, 1.3, 2.0, 2.522.6
S 2 1.4 (13)4.4 (17)0.1, 0.2, 0.2, 0.4, 0.4, 0.4, 0.5, 0.5, 0.5, 0.6, 0.9, 1.1, 1.3, 1.8, 2.4, 2.9, 4.004.0
S 3 0.8 (10)1.4 (16)0.1, 0.2, 0.2, 0.4, 0.4, 0.4, 0.5, 0.5, 0.5, 0.6, 0.9, 1.0, 1.1, 1.2, 1.2, 1.331.4
Table 10. Bayesian and classical estimates of  α μ R ( t )  and  h ( t )  from vinyl chloride data.
Table 10. Bayesian and classical estimates of  α μ R ( t )  and  h ( t )  from vinyl chloride data.
SchemeSamplePar.MLEMCMCACIHPD
Est.St.ErEst.St.ErLowerUpperIWLowerUpperIW
Sch[1] S 1 α 60.3968.48960.3490.083243.75777.03433.27860.21860.4850.2670
μ 0.28620.06520.26410.04960.15840.41400.25560.17750.35150.1741
R ( 0.5 ) 0.84660.05080.82160.04910.74710.94610.19900.73740.90030.1629
h ( 0.5 ) 0.53260.13950.59560.12590.25920.80600.54670.39650.82260.4262
S 2 α 41.36812.57041.3200.084416.73066.00649.27541.18841.4590.2718
μ 0.30180.07290.27770.05210.15890.44470.28590.18880.36730.1785
R ( 0.5 ) 0.83510.05320.80890.05020.73090.93930.20840.72380.88580.1620
h ( 0.5 ) 0.55820.14150.62240.12530.28090.83540.55450.41740.83010.4127
S 3 α 183.766.8550183.740.0554170.33197.2026.871183.64183.830.1951
μ 0.18970.04210.17980.03210.10720.27210.16490.12320.24130.1181
R ( 0.5 ) 0.81160.05810.79070.05210.69760.92550.22790.69790.87760.1797
h ( 0.5 ) 0.64670.15130.69500.12740.35010.94330.59320.47170.92380.4521
Sch[2] S 1 α 50.61117.84750.5620.085615.63285.58969.95750.42750.7010.2738
μ 0.30780.07800.28100.05480.15490.46070.30580.18940.37410.1847
R ( 0.5 ) 0.85220.05230.82500.05010.74970.95470.20500.74120.90150.1602
h ( 0.5 ) 0.51460.14440.58370.12910.23150.79770.56620.37620.79730.4211
S 2 α 104.5911.9356104.550.070581.194127.9846.787104.43104.670.2364
μ 0.23520.05470.21950.04080.12800.34230.21430.14820.29580.1476
R ( 0.5 ) 0.83330.05570.80960.05110.72410.94240.21840.71530.88960.1743
h ( 0.5 ) 0.57850.15060.63620.12890.28340.87360.59020.42760.87570.4481
S 3 α 671.1611.8650670.21663.50647.91694.4246.510664.45675.7811.3371
μ 0.12770.02860.18150.05840.07170.18370.11200.14100.22730.0863
R ( 0.5 ) 0.77060.06650.85970.09380.64020.90100.26080.80090.91010.1092
h ( 0.5 ) 0.77180.16040.53340.25280.45731.08620.62890.38510.70140.3163
Sch[3] S 1 α 21.7389.60921.6890.08642.905040.57237.66721.55421.8280.2732
μ 0.31150.07880.28760.05150.15710.46590.30880.20050.37940.1789
R ( 0.5 ) 0.79680.05720.77000.05170.68480.90880.22410.68250.85350.1710
h ( 0.5 ) 0.64930.14000.71050.12000.37490.92370.54880.50810.91300.4049
S 2 α 20.3789.718220.3300.08491.330939.42538.09520.19620.4700.2736
μ 0.35290.09090.32480.05710.17490.53100.35620.22620.42030.1940
R ( 0.5 ) 0.82300.05350.79720.04790.71820.92790.20970.71470.86890.1542
h ( 0.5 ) 0.58340.13610.64500.11540.31660.85020.53360.46380.84030.3765
S 3 α 24.2849.353924.2340.08635.951042.61836.66724.10324.3740.2706
μ 0.29220.07200.26990.04950.15110.43320.28210.18680.35970.1728
R ( 0.5 ) 0.78870.05840.76140.05380.67430.90310.22880.66850.84740.1788
h ( 0.5 ) 0.66980.14180.73160.12390.39190.94780.55590.52680.94500.4182
Table 11. Optimality criteria for the best censoring pattern.
Table 11. Optimality criteria for the best censoring pattern.
CriterionAim
O 1 Maximize trace( I ( ϑ ^ ) )
O 2 Minimize trace( I 1 ( ϑ ^ ) )
O 3 Minimize det( I 1 ( ϑ ^ ) )
O 4 Minimize  v ^ ( log ( T ^ ϱ ) )
Table 12. Optimum progressive plans from electronic tubes data.
Table 12. Optimum progressive plans from electronic tubes data.
SampleScheme O 1 O 2 O 3 O 4
ϱ 0.30.60.9
S 1 Sch[1]57.813376.196.50690.10300.694518.150
Sch[2]73.399307.774.19310.10510.701918.223
Sch[3]54.285342.026.30050.12210.820021.367
S 2 Sch[1]46.118295.746.41260.18601.233831.917
Sch[2]58.574323.425.34210.11280.756819.711
Sch[3]60.542128.922.20100.07560.516013.586
S 3 Sch[1]52.612301.005.71180.15271.016226.337
Sch[2]46.011347.396.60290.11830.796920.809
Sch[3]52.697145.483.16190.17591.171430.362
Table 13. Optimum progressive plans from vinyl chloride data.
Table 13. Optimum progressive plans from vinyl chloride data.
SampleScheme O 1 O 2 O 3 O 4
ϱ 0.30.60.9
S 1 Sch[1]240.8272.0740.29930.03710.24586.3534
Sch[2]185.93318.511.71300.04390.29237.5784
Sch[3]211.5092.3400.43660.02410.16424.3147
S 2 Sch[2]211.13158.020.74850.03480.23286.0573
Sch[2]338.57142.460.42080.03430.22405.7400
Sch[3]162.5894.4510.58100.03030.20615.4156
S 3 Sch[3]566.1546.9930.08300.02620.16864.2879
Sch[2]1224.9140.780.11490.01960.12273.0752
Sch[3]240.0587.5000.36450.02270.15404.0421
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. Statistical Analysis of Type-II Generalized Progressively Hybrid Alpha-PIE Censored Data and Applications in Electronic Tubes and Vinyl Chloride. Axioms 2023, 12, 601. https://doi.org/10.3390/axioms12060601

AMA Style

Elshahhat A, Abo-Kasem OE, Mohammed HS. Statistical Analysis of Type-II Generalized Progressively Hybrid Alpha-PIE Censored Data and Applications in Electronic Tubes and Vinyl Chloride. Axioms. 2023; 12(6):601. https://doi.org/10.3390/axioms12060601

Chicago/Turabian Style

Elshahhat, Ahmed, Osama E. Abo-Kasem, and Heba S. Mohammed. 2023. "Statistical Analysis of Type-II Generalized Progressively Hybrid Alpha-PIE Censored Data and Applications in Electronic Tubes and Vinyl Chloride" Axioms 12, no. 6: 601. https://doi.org/10.3390/axioms12060601

APA Style

Elshahhat, A., Abo-Kasem, O. E., & Mohammed, H. S. (2023). Statistical Analysis of Type-II Generalized Progressively Hybrid Alpha-PIE Censored Data and Applications in Electronic Tubes and Vinyl Chloride. Axioms, 12(6), 601. https://doi.org/10.3390/axioms12060601

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