Next Article in Journal
Phase Conventions in Hadron Physics from the Perspective of the Quark Model
Previous Article in Journal
Separable Reversible Data Hiding for Encrypted 3D Meshes Based on Self-Organized Blocking and Most Significant Bit Prediction
Previous Article in Special Issue
A Generalized Iterated Tikhonov Method in the Fourier Domain for Determining the Unknown Source of the Time-Fractional Diffusion Equation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Analysis of the Comprehensive Lifetime Performance Index for Exponentiated Fréchet Lifetime Distribution Products with Multi-Components

Department of Statistics, Tamkang University, Tamsui, Taipei 251301, Taiwan, China
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(8), 1060; https://doi.org/10.3390/sym16081060
Submission received: 18 July 2024 / Revised: 14 August 2024 / Accepted: 14 August 2024 / Published: 16 August 2024
(This article belongs to the Special Issue Computational Mathematics and Its Applications in Numerical Analysis)

Abstract

:
The lifetime performance index is commonly used in the manufacturing industry to evaluate the performance of the capabilities of the production process. For products with multiple components, the comprehensive lifetime performance index, which is a monotonically increasing function of the overall process yield, is used to relate to each individual lifetime performance index. For products where the lifetime of the ith component follows an exponentiated Fréchet lifetime distribution, we examine the maximum likelihood estimators for both the comprehensive and individual lifetime performance indices based on the progressive type I interval-censored samples, deriving their asymptotic distributions. By specifying the target level for the comprehensive lifetime performance index, we can set the desired level for individual indices. A testing procedure, using the maximum likelihood estimator as the test statistic, was developed to determine if the comprehensive lifetime performance index meets the target. Given that the lifetime distribution is asymmetric, this study pertains to asymmetrical probability distributions and their applications across diverse fields. We illustrate the power analysis of this testing procedure with figures and summarize key findings. Finally, we demonstrate the application of this testing algorithm with a practical example involving two components to verify if the overall production process achieves the assigned target level.

1. Introduction

Process capability analysis plays a crucial role in statistical quality control applications. It helps producers determine if their products meet the specified quality standards and requirements. The process capability index (PCI) is a tool used to measure and quantify this capability. Numerous widely recognized PCIs have been proposed in the literature to enhance product quality, including C p ,   C p k ,   C p m ,   C p m k ,   a n d   C L (Juran [1], Kane [2], Hsiang and Taguchi [3], Pearn et al. [4], and Montgomery [5]). Given that product lifetime is a “the-larger-the-better” quality characteristic, the PCI for unilateral tolerance, C L , is utilized to assess its performance, also known as the lifetime performance index. This investigation focuses on the evaluation of the comprehensive lifetime performance of products with multi-components by using the lifetime performance index C L . Using this index, Tong et al. [6] employed the uniformly minimum variance unbiased estimator (UMVUE) in the manufacturing of products with a single component. They used this estimator as the test statistic in their hypothesis testing procedure for exponential product lifetimes, relying on the complete sample data set. In real-world situations, researchers often face constraints that prevent them from consistently observing the lifetimes of all tested items. These constraints can include time limitations, budgetary and material restrictions, errors by typists or recorders, and mechanical or experimental difficulties. As a result, incomplete data, such as progressive censoring data, may be collected. For insights on handling such data, refer to Balakrishnan [7], Aggarwala [8], Balakrishnan and Aggarwala [9], Wu et al. [10], Sanjel and Balakrishnan [11], and Lee et al. [12]. To facilitate data collection, the progressive type I interval censoring scheme was adopted in numerous quality control applications. Wu and Lin [13] developed a testing algorithmic procedure to assess the lifetime performance index of products with the Weibull distribution. Wu et al. [14] proposed a computational algorithm for the evaluation of the lifetime performance index of products with Burr XII lifetime distribution under progressive type I interval censoring. Wu and Chang [15] assessed the lifetime performance index of products with lifetime following the exponentiated Fréchet lifetime distribution based on the progressive type I interval-censored sample. Previous studies have focused on the production of products containing a single component. For products with multiple components produced in multiple production lines, Wu and Chiang [16] assessed the overall lifetime performance index of Weibull products in multiple production lines. Wu et al. [17] investigated the comprehensive lifetime performance index of Burr XII products in multiple production lines. This investigation is an extension of the case study by Wu and Chang [15], where they extend beyond products with one component produced in a single production line to products with multiple components produced in multiple production lines. This extended research proposes a testing procedure for the comprehensive lifetime performance index for products with multiple components. Results show that their lifetimes follow the exponentiated Fréchet distribution based on the progressive type I interval-censored sample. Since the lifetime distribution was asymmetric, that study explored the asymmetrical probability distributions and their practical applications across various fields.
In this study, the comprehensive lifetime performance index proposed by Wu et al. [17] for products with multiple components manufactured across multiple production lines is utilized to evaluate whether the comprehensive lifetime performance of multiple components following the exponentiated Fréchet distributions meets a predetermined target level. The maximum likelihood estimator (MLE) for the lifetime performance index is used as the test statistic for the testing procedure, using a progressive type I interval-censored sample. Regarding the maximum likelihood estimation, Rytgaard and Van der Laan [18] investigated the target maximum likelihood estimation for causal inference in survival and competing risk analysis. Luca et al. [19] proposed an evolutionary optimization method for maximum likelihood and approximate maximum likelihood estimation of discrete latent variable models. For recent literature on reliability statistics, Goldengorin et al. [20] provide a comprehensive exploration of reliability in engineering and statistical contexts. Zhuang et al. [21] investigated the life prediction for a two-phase degradation model based on a reparameterized inverse Gaussian process. Xu et al. [22] proposed the Bayesian reliability assessment of a permanent magnet brake for a small sample size.

2. The Lifetime Performance Index and the Maximum Likelihood Estimation

In this section, we utilize the comprehensive lifetime performance index for the lifetime model of exponentiated Fréchet distribution. Section 2.1 clarifies the monotonic relationship among the overall process yield, the comprehensive lifetime performance index, and the individual lifetime performance indices. In Section 2.2, the maximum likelihood estimator and the asymptotic distribution are derived for all lifetime performance indices.

2.1. The Relationship between the Comprehensive Lifetime Performance Index and Individual Lifetime Performance Index

We consider products having d components produced in d production lines. For the ith component, the product lifetime U i follows an exponentiated Fréchet (EF) lifetime distribution with the scale parameter θ i and the shape parameter δ i , as proposed by Pedro et al. [23]. The probability density function (p.d.f.), the cumulative distribution function (c.d.f.), and the hazard function (h.f.) are defined as follows:
f U i u = θ i δ i u δ i 1 e u δ i ( 1 e u δ i ) θ i 1 , u > 0 , δ i > 0 , θ i > 0 ,
F U i u = 1 ( 1 e u δ i ) θ i , u > 0 , δ i > 0 , θ i > 0
and
r U i u = f U i u 1 F U i u = θ i δ i u δ i 1 e u δ i ( 1 e u δ i ) 1 ,
where θ i is the scale parameter, δ i is the shape parameter, and i = 1, …, d.
The p.d.f. for   δ i = 2, 4, 6 under θ i = 1 is depicted in Figure 1a and for θ i = 1, 3, 5 under   δ i = 1 is depicted in Figure 1b. The h.f. for δ i = 2, 4, 6 under θ i = 1 is depicted in Figure 2a and for θ i = 1, 3, 5 under   δ i = 1 is depicted in Figure 2b. From Figure 1a, b, we observe that the shape of the p.d.f. changes when the shape parameter δ i changes. However, when parameter θ i changes, it only affects the degree of data dispersion. Additionally, Figure 2a, b show that as the scale parameter θ i increases, the hazard rate appears to increase, especially for larger values of the shape parameter   δ i , where the increase is more pronounced. For the estimation of different censoring schemes, refer to Pedro et al. [23] and Tabassum et al. [24].
Let Y i = ln 1 e U i δ i , then we have U i = ln 1 e Y i 1 δ i with Jacobian of J = δ i u δ i 1 e u δ i 1 . We can find the probability density function (p.d.f.) of Y i as
f Y i y = f U i ln 1 e Y i 1 δ i J = θ i e θ i y ,   y > 0 , θ i > 0 .
The related cumulative distribution function (c.d.f.) and the hazard function (h.f.) of   Y i are obtained as follows:
F Y i y = 1 e θ i y , y > 0 , θ i > 0 ,
r Y i y = f Y i y 1 F Y i y = θ i .
The lifetime of components of products is a “the-larger-the-better” quality characteristic. A longer lifetime indicates better product quality, and products with superior quality have higher market appeal. To assess this characteristic, the lifetime performance index proposed by Montgomery [5] is employed. After assessing consumer demands, manufacturers set a lower specification limit L U i for the product’s lifetime.
According to Montgomery [5], the lifetime performance index for the production of the ith component is defined as follows:
C L i = μ i L U i   σ i ,
where μ i denotes the mean lifetime of the ith component, σ i denotes the standard deviation of the lifetime of the ith component, and L U i is the specified lower specification limit for the lifetime of the ith component.
The mean and standard deviation of the new lifetime variable is calculated using Y i are μ i = E Y i = 1 θ i , σ i = V a r ( Y i ) = 1 θ i . The lower specification limit becomes L i = ln 1 e L U i δ i . The lifetime performance index for the lifetime of the ith component becomes
C L i = μ i L i σ i = 1 θ i L i 1 θ i = 1 θ i L i .
Equation (8), shows that when   1 θ i > L i , C L i > 0 ; conversely, when 1 θ i < L i , C L i < 0 . It can also be noted that as the hazard rate θ i decreases, the process capability index C L i increases. This demonstrates that the process capability index   C L i can be used to evaluate lifetime performance.
A product in the ith production line is classified as a conforming product only if its lifetime exceeds L U i . The proportion of conforming products in the ith production process is defined as the process yield, which is obtained as follows:
P r i = P U i L U i = P Y i L i = e θ i L i = e C L i 1 , < C L i < 1 .
Equation (9) shows that there is a strict increasing relationship between the process yield P r i   and the lifetime performance index C L i . When the process yield P r i is higher, the lifetime performance index C L i   increases accordingly. This once again demonstrates that the lifetime performance index   C L i can be used to evaluate the product’s lifetime performance.
Assuming the lifetimes of d components of products produced in d independent production lines, the overall process yield P r   is defined as
P r = P Y i L i ,   i = 1,2 , , d = e i = 1 d θ i L i = e i = 1 d C L i d , < C L i < 1 .
We employ the comprehensive lifetime performance index CT defined by Wu et al. [17] in the following equation:
P r = e x p C T 1 , < C T < 1 .
As shown in Equation (11), the comprehensive lifetime performance index C T still exhibits a strictly increasing relationship with the overall process yield P r . Furthermore, based on Equations (10) and (11), we can find the relationship between the comprehensive lifetime performance index C T and the individual lifetime performance index C L i   as follows:
C T = i = 1 d C L i d + 1 = i = 1 d C L i d 1 , < C T < 1 .
A reasonable assumption in this context would be that the capabilities of the production lines are equally important for quality engineers across all components of products. The rational setup of equal individual life performance indices is defined as C L 1 = C L 2 = = C L d = C L . Substituting this back into Equation (12), we obtain the relationship between C T   and C L as
C L = C T + d 1 d , < C T < 1 .
The comprehensive lifetime performance index C T , the related overall process yield P r , and the corresponding values of C L are displayed in Table 1 for d = 2 , 3 , 4 , 5 , 6 , 7 . As Table 1 shows, if we require P r   to surpass 0.9753, then C T must exceed 0.975. Regarding the target value of C T = 0.975, we can find the corresponding values for the individual lifetime performance index C L   to exceed 0.9875, 0.9917, 0.9938, 0.9950, 0.9958 ,   a n d   0.9964   f o r   d = 2 , 3 , 4 , 5 , 6 , 7 .

2.2. The Maximum Likelihood Estimators for the Life Performance Indices

The progressive type I interval censoring is addressed as follows: Let n products be subjected to a life test with m inspection intervals. This censoring scheme, used in survival analysis and reliability engineering, addresses scenarios where the exact failure times are not observed; instead, the number of failure units Xi is recorded at the ith inspection time point ti, i = 1, …, m. Here, (t1, …, tm) are the predetermined inspection time points, and tm is the ending time of this experiment. When the number of failure units Xi is observed at the ith inspection time point ti, we remove Ri products from the unexamined units under the removal probability pi, i = 1, …, m. When the experiment is completed, we collect the progressive type I interval-censored sample (X1, …, Xm) under the censoring scheme of (R1, …, Rm). We will utilize this sample set to make inferences on the comprehensive lifetime performance index and each individual lifetime performance index. Assuming the products have d components produced in d independent production lines, we observe X i 1 , , X i m as the progressive type I interval-censored sample for the lifetimes of the ith component of products at the time points of t 1 , , t m under the progressive censoring scheme of   R i 1 , , R i m . The likelihood function of this censored sample is
L θ i j = 1 m F U i t j F U i t j 1 X i j 1 F U i t j R i j            j = 1 m 1 e t j 1 δ i θ i 1 e t j δ i θ i X i j 1 e t j δ i θ i R i j   j = 1 m 1 1 e t j δ i 1 e t j 1 δ i θ i X i j 1 e t j δ i θ i X i j 1 e t j δ i θ i R i j .
Firstly, taking the logarithm of the likelihood function yields the log-likelihood function as follows:
ln L θ i = j = 1 m X i j ln 1 1 e t j δ i 1 e t j 1 δ i θ i θ i j = 1 m R i j ln 1 e t j δ i + X i j ln ( 1 e t j 1 δ i ) .
Next, differentiating Equation (15) with respect to the parameter θ i , setting the derivative equal to zero, yields the log-likelihood equation as follows:
  d d θ i l n L θ i = j = 1 m X i j ln ( 1 e t j 1 δ i ) ln ( 1 e t j δ i ) 1 1 e t j δ i 1 e t j 1 δ i θ i 1 1 e t j δ i 1 e t j 1 δ i θ i   j = 1 m R i j ln 1 e t j δ i + X i j ln ( 1 e t j 1 δ i ) 0 .
The maximum likelihood estimator for the parameter θ i , denoted by θ ^ i , can be obtained by solving Equation (16). However, since Equation (16) does not have a closed-form solution, numerical methods are required to solve it. The maximum likelihood estimator k ^ i has an asymptotic normal distribution. The definition of the Fisher’s information number is as follows:
I θ i = E d 2 l n L θ i d θ i 2 .
Taking the second derivative of log-likelihood function yields
d 2 d θ i 2 l n L θ i = j = 1 m X i j ln 1 e t j 1 δ i ln 1 e t j δ i 2 1 e t j δ i 1 e t j 1 δ i θ i 1 1 e t j δ i 1 e t j 1 δ i θ i 2 .
We assume that the number of failures X i j under the progressive type I interval-censored scenario follows a binomial distribution, which is denoted by the following distribution:
X i j X i , j 1 , , X i 1 , R i , j 1 , , R i 1 ~ B i n o m i a l n l = 1 j 1 X i l l = 1 j 1 R i l , q i j ,
where q i j = F t j F t j 1 1 F t j 1 = 1 1 e t j δ i 1 e t j 1 δ i θ i , j = 1 , , m , i = 1 , , d .
Likewise, we assume that the number of progressive censorings at the ith observation time point follows a binomial distribution, which is denoted by the following distribution:
R i j R i , j 1 , , R i 1 , X i , j 1 , , X i 1 ~ B i n o m i a l n l = 1 j 1 X i l l = 1 j 1 R i l , p j ,
where p j represents the progressive removal probabilities.
Through Equations (18) and (19), the expected values are obtained as follows:
E X i j = E E X i j X i , j 1 , , X i 1 , R i , j 1 , , R i 1 n q i j l = 1 j 1 ( 1 q i l ) ( 1 p i l ) , j = 1 , , m .
E R i j = E E R i j R i , j 1 , , R i 1 , X i , j 1 , , X i 1 = n p i j 1 q i j l = 1 j 1 ( 1 q i l ) ( 1 p i l ) ,   j = 1 , , m .
Based on the above results, the Fisher’s information number can be derived as follows:
I θ i = E d 2 l n L ( θ i ) d θ i 2 = j = 1 m n q i j ln 1 e t j 1 δ i ln 1 e t j δ i 2 1 e t j δ i 1 e t j 1 δ i θ i 1 1 e t j δ i 1 e t j 1 δ i θ i 2 l = 1 j 1 1 p i l 1 q i l .
Using the information obtained from Equation (22), we can deduce that the asymptotic normal distribution of θ ^ i as θ ^ i d n N ( θ i , I 1 θ i ) .
For the convenience of data collection, this paper considers the special case where time intervals are all equal. Specifically, the time length in the jth interval is set as t j t j 1 = t , j = 1 , , m , where t is the equal time interval. Simplifying Equation (16) using this assumption, we get
  d d θ i l n L θ i = j = 1 m X i j ln ( 1 e ( j 1 t ) δ i ) ln ( 1 e ( j t ) δ i ) 1 1 e ( j t ) δ i 1 e ( j 1 t ) δ i θ i 1 1 e ( j t ) δ i 1 e ( j 1 t ) δ i θ i j = 1 m R i j ln 1 e j t δ i + X i j ln 1 e j 1 t δ i 0 .
By solving Equation (23), we obtain the maximum likelihood estimator for the parameter θ i , denoted as θ ^ i .
According to the invariance property of the maximum likelihood estimator, we can obtain the maximum likelihood estimator for the lifetime performance index C L i as follows:
C ^ L i = 1 θ ^ i L i .
By using the delta method, we can obtain
C ^ L i d n N C L i , L i 2 V θ ^ i .
According to the invariance property of maximum likelihood estimators, we can obtain the maximum likelihood estimator for the comprehensive lifetime performance index C T as follows:
  C ^ T = i = 1 d C ^ L i d 1 , < C T < 1 .
It is evident that the asymptotic variance of   C ^ T is i = 1 d V C ^ L i , and the estimator of the asymptotic variance is i = 1 d V ^ C ^ L i . Apparently, the asymptotic distribution for C ^ T is
C ^ T   d n   N C T , i = 1 d V ^ C ^ L i .

3. The Testing Algorithm Procedure for the Comprehensive Lifetime Performance Index

In Section 3.1, the testing procedure to determine whether the comprehensive lifetime performance index meets the target level and the related power analysis are presented. Section 3.2 includes a numerical example to showcase the testing procedure.

3.1. The Testing Procedure for the Comprehensive Lifetime Performance Index and Power Analysis

Firstly, let us assume that the manufacturer sets a target level c 0   for the comprehensive lifetime performance of products with d components. Based on this target level c 0 , the appropriate target level c 0 * for the individual lifetime performance index C L can be obtained from Table 1. We can then establish a null and alternative hypothesis as follows:
H0: 
C L i c 0 * , for some i (indicating the process is not capable).
vs.
H1: 
C L i > c 0 * ,   i = 1 , , d (indicating the process is capable).
This hypothesis test is known as the intersection–union test (IUT) in Casella and Berger [25].
As proposed by Wu et al. [13], the ith test of this hypothesis test is formulated as follows:
H 0 i : C L i c 0 *     a n d   H 1 i : C L i > c 0 *
We use the maximum likelihood estimator of the individual lifetime performance index C L i as C ^ L i = 1 θ ^ i L i to as the test statistic. Let α = α 1 d be the significance level for the ith test, where the critical value is determined as C L i 0 = 1 Z 1 α w θ i 0 + θ i 0 L i , and where w θ i 0 = I 1 θ i 0 as per Wu and Chang [11]. The rejection region for the ith test is given by R i = C ^ L i C ^ L i > C L i 0 . For the composite test with H 0 : C L i c 0 * , for some i and H 1 : C L i > c 0 * ,   i = 1 , , d , the overall rejection region is determined as R = i = 1 d R i = i = 1 d C ^ L i     C ^ L i > C L i 0 such that the size of this test is still α.
The proposed testing algorithmic procedure about C T can be constructed as Algorithm 1:
Algorithm 1. The testing procedure for C T .
Step 1: For two components with specified lower limits L1 and L2, gather the progressive type I interval-censored samples   X i 1 , , X i m at predetermined inspection times t 1 , , t m , using censoring schemes of R i 1 , , R i m from the exponentiated Fréchet distribution.
Step 2: For a defined process yield Pr and the component count d, locate the appropriate values for CT (denoted as c0) and CL (denoted as c 0 * ) from Table 1 to meet each component’s required level. Following this, we outline the testing hypothesis H 0 : C L i c 0 * for some i vs. H 1 : C L i > c 0 * , i = 1 , , d .
Step 3: Quantify the test statistic C ^ L i .
Step 4: For the specified significance level α, compute the critical value C L i 0 = 1 L i θ i 0 + Z 1 α w θ i 0 , where θ i 0 = 1 c 0 * L i , α1/d, w θ i 0 = I 1 θ i 0 and I i 1 θ i 0 is defined in Equation (22).
Step 5: Assess C ^ L i in relation to C L i 0 . If C ^ L i   e x c e e d s   C L i 0 for all i, it can be concluded that the comprehensive lifetime performance index for products with d components produced in d production lines has surpassed the target level c0.
Furthermore, the test power function h c 1 of the testing procedure at the point of C T = d C L d 1 = c 1 > c 0 or C L = c 1 + d 1 d = c 1 * in the alternative hypothesisis obtained as follows:
h c 1 = P C ^ L i > C L i 0 c 1 * = 1 θ i 1 L i   ,   i = 1 , , d = P 1 θ ^ i L i > 1 L i θ i 0 + Z 1 α w θ i 0 θ i 0 = 1 c 0 * L i , θ i 1 = 1 c 1 * L i , i = 1 , , d = i = 1 d P θ ^ i < θ i 0 + Z 1 α w θ i 0 θ i 0 = 1 c 0 * L i , θ i 1 = 1 c 1 * L i = i = 1 d P θ ^ i θ i 1 w θ i 1 < θ i 0 + Z 1 α w θ i 0 θ i 1 w θ i 1 θ i 0 = 1 c 0 * L i , θ i 1 = 1 c 1 * L i = i = 1 d Φ θ i 0 θ i 1 + Z 1 α w θ i 0 w θ i 1 θ i 0 = 1 c 0 * L i , θ i 1 = 1 c 1 * L i ,
where Φ ( . ) is the c.d.f. for the standard normal distribution, θ i 0 = 1 c 0 * L i ,   θ i 1 = 1 c 1 * L i , i = 1 , , d . For the case of L1 = … = Ld = L, this yields the following equation for the power function
h c 1 = Φ θ 0 θ 1 + Z 1 α w θ 0 w θ 1 θ 0 = 1 c 0 * L i , θ 1 = 1 c 1 * L i d
To test H 0 : C T 0.8 , the powers h(c1) are computed using Equation (28) and they are presented in Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8 and Table A9 for d = 2, 3, 4 with α = 0.050, 0.075, 0.100, respectively, for c1 = 0.8, 0.825, 0.85, 0.875, 0.9, 0.925, 0.95, 0.975; m = 2, 4, 6, n = 25, 50, 75,100; and p = 0.050, 0.075, 0.100 with L = 0.05, T = 0.5. The power values depicted in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 serve to exemplify various standard cases. Our observations are as follows: (1) In Figure 3, decreasing p enhances the power when d = 1, n = 25, m = 4, and α = 0.05. This behavior is consistent across different combinations of d, n, m, and α. (2) In Figure 4, as n increases, the power rises with d = 1, m = 4, p = 0.05, and α = 0.05 held constant. Similar patterns emerge across various combinations of d, n, p, and α. (3) In Figure 5, decreasing m enhances the power for fixed values of d = 1, α = 0.05, n = 25 and p = 0.05. (4) In Figure 6, the power increases with α under d = 1, n = 25, m = 4, and p = 0.05. This trend is observed across various combinations of d, n, m, and p. (5) In Figure 7, the power improves as d increases with n = 25, m = 4, p = 0.05, and α = 0.05 held constant. This trend remains consistent across different combinations of n, m, p, and α. The improvement in power is getting not significant when d is greater than 3. (6) Across Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7, the power increases with the value of c1 for any combination of d, n, m, p, and α.

3.2. Practical Example

We consider the case of products with two components (d = 2) in this section. The lifetime data for the first component consists of the failure times (number of cycles) of n = 36 electrical appliances in Lawless [26]. The data of 36 failure times U1j, j = 1, …, 36, is listed as follows: 0.001, 0.0035, 0.0049, 0.017, 0.0329, 0.0381, 0.0707, 0.0958, 0.1062, 0.1167, 0.1594, 0.1925, 0.199, 0.2223, 0.3270, 0.2400, 0.2451, 0.2471, 0.2551, 0.2565, 0.2568, 0.2702, 0.2761, 0.2831, 0.3034, 0.3034, 0.3059, 0.3113, 0.3214, 0.3478, 0.3504, 0.4329, 0.6367, 0.6976, 0.7846, 1.3403. The p-value of the G test based on the Gini statistic (see Gail and Gastwirth [27]) is a function of the values of δ1. The Gini statistic is obtained by G n = j = 1 n 1 j n j Y 1 ( j + 1 ) Y 1 ( j ) n 1 j = 1 n n j + 1 Y 1 ( j ) Y 1 ( j 1 ) , where Y 1 ( j ) = ln 1 e U 1 j δ 1 , j = 1, …, 36. Figure 8 depicts the p-values vs. the values of δ1. From this figure, it reveals that the p-value peaks at 0.961 when δ1 = 0.6, so δ1 is concluded to be 0.6.
The data for the lifetimes of the second component consists of n = 50 failure times in Aarset [28]. The data of U2j, j = 1, …, 50, is listed as follows: 0.001, 0.002, 0.010, 0.010, 0.010, 0.010, 0.020, 0.030, 0.060, 0.070, 0.110, 0.120, 0.180, 0.180, 0.180, 0.180, 0.180, 0.210, 0.320, 0.360, 0.400, 0.450, 0.450, 0.470, 0.500, 0.550, 0.600, 0.630, 0.630, 0.670, 0.670, 0.670, 0.670, 0.670, 0.720, 0.750, 0.790, 0.820, 0.820, 0.830, 0.840, 0.840, 0.840, 0.850, 0.850, 0.850, 0.850, 0.850, 0.860, 0.860. Similarly, Figure 9 shows the p-values vs. the values of δ2. Figure 9 reveals that the p-value peaks at 0.9881 when δ2 = 0.87, so δ2 is concluded to be 0.87.
High p-values indicate that two datasets fit the exponentiated Fréchet distribution very well. We use these two datasets to illustrate how to implement the testing procedure proposed in Section 3.1 with an end time of T = 2 and the number of inspections m = 8 under the pre-specified removal percentages of (p1, p2, p3, p4, p5, p6, p7, p8) = (0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 0.05, 1.00, respectively).
The testing procedure is implemented in the following order according to the Algorithm 1 given in Section 3.1:
Step 1: For two components with specified lower limits L1 = L2 = L = 0.05, observe the progressive type I interval-censored sample ( X 11 , X 12 , X 13 , X 14 , X 15 , X 16 , X 17 , X 18 ) = (18, 13, 2, 0, 0, 1, 0, 0, respectively) and ( X 21 , X 22 , X 23 , X 24 ,   X 25 , X 26 , X 27 , X 28 ) = ( 19 ,   6 ,   9 ,   12 ,   0 ,   0 ,   0 ,   0 , respectively) for each component at the preset times (t1, …, t8) = ( 0.25 ,   0.5 ,   0.75 ,   1.0 ,   1.25 ,   1.5 ,   1.75 ,   2.0 , respectively) with censoring schemes of ( R 11 , R 12 , R 13 , R 14 , R 15 , R 16 , R 17 , R 18 ) = ( 1 ,   0 ,   1 ,   0 ,   0 ,   0 ,   0 ,   0 , respectively and ( R 21 , R 22 , R 23 , R 24 , R 25 , R 26 , R 27 , R 28 ) = ( 1 ,   1 ,   1 ,   0 ,   0 ,   0 ,   0 ,   1 , respectively) from the exponentiated Fréchet distribution.
Step 2: For a defined conforming rate Pr = 0.9048 and the component count d = 2, the appropriate values for c0 = 0.9 and c 0 * = 0.95 are found from Table 1. The value of c 0 * = 0.95 is the required level for each individual lifetime performance index. We then outline the testing hypotheses H 0 : C L i 0.95 for some i vs. H 1 : C L i > 0.95 , i = 1 , 2 .
Step 3: Obtain the maximum likelihood estimators θ ^ 1 = 7.376792 and θ ^ 2 = 5.062244 for two components. Compute the values of the test statistics C ^ L 1 = 0.9823 and C ^ L 2 = 0.9969 , respectively.
Step 4: For the level of significance α = 0.05, we have α = 0.05 1 2 = 0.2236 . Compute θ 10 = 1 c 0 * L 1 = 1 0.95 0.0024 = 20.8333 and θ 20 = 1 c 0 * L 2 = 1 0.95 0.0006 = 83.3333 . Then we can calculate the critical values of C L 1 0 = 0.9578 and C L 2 0 = 0.9578 .
Step 5: Since C ^ L 1 = 0.9823 > C L 1 0 = 0.9578 and C ^ L 2 = 0.9969 > C L 2 0 = 0.9578 , we can infer that both individual lifetime performance indices reach the desired target values, so that the comprehensive lifetime performance index for products with two components surpasses the target level c0 = 0.9.

4. Conclusions

Assessing lifetime performance indices for products, particularly those with exponentiated Fréchet lifetime distribution, is essential across various manufacturing sectors. For products composed of multiple components produced on separate production lines, each with lifetimes following an exponentiated Fréchet distribution, we employed a comprehensive lifetime performance index, which is a monotonically increasing function of both the overall process yield and the individual lifetime performance indices. Using the progressive type I interval-censored samples, we developed a hypothesis testing procedure to evaluate whether the comprehensive index meets the target level by testing multiple hypotheses for each individual lifetime performance index. We analyzed the effects of the number of inspection intervals, sample size, removal probability, significance level, and the number of components on test power and demonstrated the proposed testing process with a numerical example involving two-component products for a specified target value and a given level of significance. In the last few years, developments related to the detailed study of families of classical probability distributions and their modifications have been proposed in relation to the saturation characteristic regarding the Hausdorff distance. We will study this topic for exponentiated Fréchet distributions in the future. On the other hand, we will also work on the optimal experimental design related to the proposed testing procedure for the comprehensive lifetime performance index to obtain the given test power or to find the minimum total experimental cost under a specific cost structure. The current lifetime performance index has some limitations in addressing asymmetrical distributions. Future research should consider the new lifetime performance index by replacing the mean with the median.

Author Contributions

Conceptualization, S.-F.W.; methodology, S.-F.W.; software, S.-F.W. and H.-C.Y.; validation, H.-C.Y.; formal analysis, S.-F.W.; investigation, S.-F.W. and H.-C.Y.; resources, S.-F.W.; data curation, S.-F.W. and H.-C.Y.; writing—original draft preparation, S.-F.W. and H.-C.Y.; writing—review and editing, S.-F.W.; visualization, H.-C.Y.; supervision, S.-F.W.; project administration, S.-F.W.; funding acquisition, S.-F.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science and Technology Council, Taiwan, NSTC 113-2118-M-032-002.

Data Availability Statement

Data is available in a publicly accessible repository. The data presented in this study is openly available in Lawless [26] and Aarset [28].

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

  θ i the scale parameter for the exponentiated Fréchet lifetime distribution;
δ i the shape parameter for the exponentiated Fréchet lifetime distribution;
C L i the lifetime performance index for the ith component;
P r i the process yield for the ith component;
C T the comprehensive lifetime performance index;
P r the overall process yield.

Appendix A

Table A1. The power h(c1) at d = 2, α = 0.01.
Table A1. The power h(c1) at d = 2, α = 0.01.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.01000.01840.03550.07190.15130.32170.64000.9670
0.0750.01000.01820.03480.06980.14590.30960.62130.9609
0.10.01000.01800.03410.06780.14060.29760.60210.9537
500.050.01000.02710.07360.19170.44050.78650.98401.0000
0.0750.01000.02670.07180.18590.42790.77290.98121.0000
0.10.01000.02630.07000.18010.41510.75870.97801.0000
750.050.01000.03580.11930.33470.69280.95590.99971.0000
0.0750.01000.03530.11610.32530.67910.95070.99961.0000
0.10.01000.03470.11300.31580.66510.94490.99951.0000
1000.050.01000.04490.17050.47670.85020.99271.00001.0000
0.0750.01000.04410.16590.46490.84000.99141.00001.0000
0.10.01000.04330.16140.45300.82910.98981.00001.0000
4250.050.01000.01910.03790.07750.16210.33760.65210.9672
0.0750.01000.01860.03580.07160.14710.30490.60220.9490
0.10.01000.01810.03390.06610.13310.27390.55080.9236
500.050.01000.02820.07800.20200.45410.79220.98351.0000
0.0750.01000.02720.07290.18580.41980.75540.97501.0000
0.10.01000.02620.06810.17060.38640.71560.96321.0000
750.050.01000.03740.12560.34740.70200.95620.99961.0000
0.0750.01000.03570.11680.32150.66510.94120.99931.0000
0.10.01000.03420.10860.29670.62680.92250.99861.0000
1000.050.01000.04680.17870.48970.85480.99251.00001.0000
0.0750.01000.04460.16590.45750.82670.98841.00001.0000
0.10.01000.04250.15380.42580.79540.98271.00001.0000
6250.050.01000.01900.03710.07470.15390.31720.61760.9538
0.0750.01000.01810.03400.06580.13160.26830.53710.9132
0.10.01000.01740.03120.05800.11220.22460.45690.8517
500.050.01000.02780.07550.19270.43190.76600.97691.0000
0.0750.01000.02620.06780.16840.37910.70340.95791.0000
0.10.01000.02470.06090.14680.32980.63560.92830.9997
750.050.01000.03660.12080.33140.67670.94500.99931.0000
0.0750.01000.03410.10760.29200.61640.91520.99811.0000
0.10.01000.03190.09580.25590.55470.87570.99521.0000
1000.050.01000.04570.17140.46900.83490.98941.00001.0000
0.0750.01000.04230.15210.41890.78550.98000.99991.0000
0.10.01000.03920.13490.37130.72990.96470.99981.0000
Table A2. The power h(c1) at d = 2, α = 0.05.
Table A2. The power h(c1) at d = 2, α = 0.05.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.05000.08870.15980.28590.49030.75550.95870.9999
0.0750.05000.08810.15760.28080.48120.74520.95440.9998
0.10.05000.08740.15540.27570.47210.73440.94970.9998
500.050.05000.11720.25840.50210.79310.96970.99971.0000
0.0750.05000.11610.25430.49370.78420.96660.99961.0000
0.10.05000.11490.25020.48520.77490.96320.99951.0000
750.050.05000.14290.35010.66740.92480.99691.00001.0000
0.0750.05000.14130.34430.65820.91950.99641.00001.0000
0.10.05000.13970.33860.64880.91370.99581.00001.0000
1000.050.05000.16730.43430.78410.97420.99971.00001.0000
0.0750.05000.16530.42730.77560.97160.99971.00001.0000
0.10.05000.16320.42020.76680.96870.99961.00001.0000
4250.050.05000.09050.16430.29300.49690.75660.95660.9998
0.0750.05000.08870.15820.27900.47250.72850.94390.9997
0.10.05000.08690.15230.26540.44830.69900.92860.9993
500.050.05000.11970.26460.50930.79520.96850.99961.0000
0.0750.05000.11650.25330.48650.77070.95940.99931.0000
0.10.05000.11340.24230.46400.74480.94820.99881.0000
750.050.05000.14610.35740.67330.92490.99661.00001.0000
0.0750.05000.14160.34170.64810.90980.99501.00001.0000
0.10.05000.13720.32630.62250.89270.99261.00001.0000
1000.050.05000.17100.44220.78830.97390.99971.00001.0000
0.0750.05000.16530.42310.76500.96630.99941.00001.0000
0.10.05000.15970.40430.74060.95710.99901.00001.0000
6250.050.05000.08970.16100.28440.48010.73530.94610.9997
0.0750.05000.08680.15160.26290.44210.68890.92150.9991
0.10.05000.08410.14280.24290.40560.64050.88990.9976
500.050.05000.11800.25780.49400.77700.96110.99931.0000
0.0750.05000.11310.24040.45830.73610.94320.99841.0000
0.10.05000.10840.22410.42400.69290.91980.99651.0000
750.050.05000.14360.34760.65580.91340.99521.00001.0000
0.0750.05000.13670.32320.61530.88620.99131.00001.0000
0.10.05000.13020.30030.57470.85430.98490.99991.0000
1000.050.05000.16780.43000.77180.96800.99951.00001.0000
0.0750.05000.15890.40020.73320.95330.99881.00001.0000
0.10.05000.15060.37170.69260.93410.99741.00001.0000
Table A3. The power h(c1) at d = 2, α = 0.1.
Table A3. The power h(c1) at d = 2, α = 0.1.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.10000.16970.28550.46280.69180.90110.99271.0000
0.0750.10000.16870.28250.45680.68390.89550.99171.0000
0.10.10000.16770.27940.45080.67580.88950.99061.0000
500.050.10000.21250.41090.67890.90650.99261.00001.0000
0.0750.10000.21080.40600.67150.90140.99171.00001.0000
0.10.10000.20910.40110.66390.89590.99071.00001.0000
750.050.10000.24900.51320.81110.97280.99951.00001.0000
0.0750.10000.24680.50710.80440.97050.99941.00001.0000
0.10.10000.24450.50090.79740.96800.99931.00001.0000
1000.050.10000.28210.59810.89000.99231.00001.00001.0000
0.0750.10000.27930.59130.88470.99131.00001.00001.0000
0.10.10000.27660.58440.87900.99030.99991.00001.0000
4250.050.10000.17180.28970.46710.69260.89820.99161.0000
0.0750.10000.16900.28130.45100.67120.88230.98841.0000
0.10.10000.16610.27310.43510.64930.86480.98421.0000
500.050.10000.21540.41600.68190.90550.99191.00001.0000
0.0750.10000.21070.40260.66170.89100.98900.99991.0000
0.10.10000.20610.38940.64110.87510.98520.99991.0000
750.050.10000.25250.51870.81290.97200.99941.00001.0000
0.0750.10000.24630.50200.79430.96520.99901.00001.0000
0.10.10000.24020.48540.77480.95710.99851.00001.0000
1000.050.10000.28610.60360.89100.99191.00001.00001.0000
0.0750.10000.27850.58510.87610.98910.99991.00001.0000
0.10.10000.27100.56650.85990.98550.99991.00001.0000
6250.050.10000.17030.28460.45600.67620.88500.98871.0000
0.0750.10000.16580.27150.43080.64180.85710.98181.0000
0.10.10000.16150.25920.40670.60720.82600.97190.9999
500.050.10000.21270.40730.66730.89390.98940.99991.0000
0.0750.10000.20540.38630.63480.86870.98320.99981.0000
0.10.10000.19850.36640.60220.84030.97450.99961.0000
750.050.10000.24880.50750.79920.96650.99911.00001.0000
0.0750.10000.23910.48130.76830.95360.99821.00001.0000
0.10.10000.23000.45590.73590.93730.99651.00001.0000
1000.050.10000.28150.59100.87990.98960.99991.00001.0000
0.0750.10000.26970.56160.85420.98380.99981.00001.0000
0.10.10000.25840.53260.82590.97580.99951.00001.0000
Table A4. The power h(c1) at d = 3, α = 0.01.
Table A4. The power h(c1) at d = 3, α = 0.01.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.01000.02090.04490.09830.21390.44150.77910.9916
0.0750.01000.02070.04410.09570.20740.42890.76560.9899
0.10.01000.02050.04320.09310.20090.41620.75140.9880
500.050.01000.03050.08880.23390.51190.83850.98971.0000
0.0750.01000.03010.08670.22750.49980.82810.98811.0000
0.10.01000.02960.08470.22110.48740.81710.98621.0000
750.050.01000.04000.13890.38010.73540.96350.99971.0000
0.0750.01000.03940.13540.37060.72360.95950.99961.0000
0.10.01000.03870.13200.36100.71130.95510.99951.0000
1000.050.01000.04970.19310.51620.86700.99261.00001.0000
0.0750.01000.04890.18820.50490.85830.99141.00001.0000
0.10.01000.04800.18330.49340.84910.99011.00001.0000
4250.050.01000.02110.04530.09830.21120.43110.76250.9888
0.0750.01000.02050.04300.09140.19410.39800.72450.9825
0.10.01000.01990.04080.08490.17810.36580.68380.9733
500.050.01000.03060.08880.23120.50200.82580.98701.0000
0.0750.01000.02950.08330.21440.46970.79600.98111.0000
0.10.01000.02850.07810.19850.43770.76340.97311.0000
750.050.01000.04020.13830.37460.72370.95790.99951.0000
0.0750.01000.03850.12910.34930.69100.94520.99911.0000
0.10.01000.03680.12050.32490.65690.92960.99841.0000
1000.050.01000.04990.19170.50840.85750.99081.00001.0000
0.0750.01000.04750.17880.47810.83240.98671.00001.0000
0.10.01000.04530.16650.44810.80460.98110.99991.0000
6250.050.01000.02070.04360.09290.19700.40200.72740.9826
0.0750.01000.01980.04010.08270.17200.35180.66280.9669
0.10.01000.01900.03690.07370.14970.30520.59430.9414
500.050.01000.02980.08450.21730.47370.79830.98131.0000
0.0750.01000.02810.07640.19230.42350.74620.96761.0000
0.10.01000.02660.06900.16980.37570.68900.94710.9998
750.050.01000.03890.13090.35310.69440.94590.99911.0000
0.0750.01000.03630.11740.31480.64040.92030.99781.0000
0.10.01000.03400.10520.27930.58460.88720.99511.0000
1000.050.01000.04810.18120.48230.83470.98691.00001.0000
0.0750.01000.04460.16190.43510.79020.97740.99991.0000
0.10.01000.04140.14460.38990.74040.96290.99961.0000
Table A5. The power h(c1) at d = 3, α = 0.05.
Table A5. The power h(c1) at d = 3, α = 0.05.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.05000.09590.18170.33120.55750.81490.97531.0000
0.0750.05000.09520.17930.32590.54900.80680.97280.9999
0.10.05000.09450.17690.32050.54040.79840.97010.9999
500.050.05000.12530.28250.53950.81940.97400.99971.0000
0.0750.05000.12410.27830.53150.81180.97150.99961.0000
0.10.05000.12280.27400.52350.80380.96890.99951.0000
750.050.05000.15150.37290.69010.92940.99651.00001.0000
0.0750.05000.14980.36720.68170.92470.99601.00001.0000
0.10.05000.14820.36150.67300.91980.99541.00001.0000
1000.050.05000.17610.45400.79420.97280.99951.00001.0000
0.0750.05000.17400.44720.78650.97030.99941.00001.0000
0.10.05000.17190.44030.77850.96770.99931.00001.0000
4250.050.05000.09590.18080.32700.54740.80200.97020.9999
0.0750.05000.09400.17440.31300.52480.77940.96220.9999
0.10.05000.09210.16830.29930.50220.75550.95240.9997
500.050.05000.12510.28020.53190.80900.96960.99951.0000
0.0750.05000.12190.26890.51060.78770.96180.99921.0000
0.10.05000.11870.25800.48930.76500.95250.99871.0000
750.050.05000.15120.36940.68140.92270.99551.00001.0000
0.0750.05000.14670.35410.65840.90890.99371.00001.0000
0.10.05000.14230.33930.63480.89350.99121.00001.0000
1000.050.05000.17560.44940.78600.96910.99931.00001.0000
0.0750.05000.16990.43120.76450.96150.99901.00001.0000
0.10.05000.16440.41330.74200.95250.99841.00001.0000
6250.050.05000.09440.17550.31470.52650.78000.96200.9998
0.0750.05000.09150.16600.29330.49100.74190.94580.9996
0.10.05000.08870.15700.27320.45650.70160.92500.9990
500.050.05000.12250.27070.51280.78890.96180.99921.0000
0.0750.05000.11760.25360.47940.75290.94650.99821.0000
0.10.05000.11290.23750.44710.71480.92710.99661.0000
750.050.05000.14750.35630.66060.90960.99371.00001.0000
0.0750.05000.14070.33300.62340.88470.98940.99991.0000
0.10.05000.13420.31090.58600.85600.98300.99991.0000
1000.050.05000.17100.43380.76650.96190.99901.00001.0000
0.0750.05000.16230.40550.73070.94710.99791.00001.0000
0.10.05000.15420.37850.69340.92850.99611.00001.0000
Table A6. The power h(c1) at d = 3, α = 0.1.
Table A6. The power h(c1) at d = 3, α = 0.1.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.10000.17930.31060.50410.73540.92350.99501.0000
0.0750.10000.17830.30750.49840.72860.91930.99441.0000
0.10.10000.17720.30440.49260.72160.91490.99371.0000
500.050.10000.22220.43280.70150.91440.99261.00001.0000
0.0750.10000.22050.42800.69470.91000.99181.00001.0000
0.10.10000.21870.42320.68780.90540.99090.99991.0000
750.050.10000.25840.53000.81920.97170.99921.00001.0000
0.0750.10000.25610.52410.81310.96960.99911.00001.0000
0.10.10000.25380.51820.80670.96720.99901.00001.0000
1000.050.10000.29090.60940.89000.99050.99991.00001.0000
0.0750.10000.28810.60290.88500.98950.99991.00001.0000
0.10.10000.28530.59640.87980.98850.99991.00001.0000
4250.050.10000.17890.30810.49720.72450.91500.99351.0000
0.0750.10000.17600.29980.48200.70580.90260.99131.0000
0.10.10000.17320.29170.46690.68660.88900.98851.0000
500.050.10000.22140.42880.69310.90710.99080.99991.0000
0.0750.10000.21680.41590.67460.89430.98810.99991.0000
0.10.10000.21230.40330.65580.88030.98460.99981.0000
750.050.10000.25730.52490.81140.96810.99901.00001.0000
0.0750.10000.25120.50910.79430.96140.99851.00001.0000
0.10.10000.24530.49350.77640.95360.99781.00001.0000
1000.050.10000.28960.60370.88360.98880.99991.00001.0000
0.0750.10000.28220.58630.86940.98560.99981.00001.0000
0.10.10000.27490.56880.85420.98170.99971.00001.0000
6250.050.10000.17650.30080.48310.70620.90220.99111.0000
0.0750.10000.17210.28810.45940.67600.88030.98641.0000
0.10.10000.16780.27600.43650.64530.85570.97971.0000
500.050.10000.21760.41740.67580.89450.98800.99991.0000
0.0750.10000.21050.39750.64610.87210.98220.99981.0000
0.10.10000.20370.37840.61630.84710.97430.99951.0000
750.050.10000.25220.51090.79540.96150.99841.00001.0000
0.0750.10000.24290.48630.76700.94880.99721.00001.0000
0.10.10000.23400.46240.73720.93330.99531.00001.0000
1000.050.10000.28340.58830.87040.98570.99981.00001.0000
0.0750.10000.27200.56070.84600.97920.99961.00001.0000
0.10.10000.26110.53360.81950.97060.99911.00001.0000
Table A7. The power h(c1) at d = 4, α = 0.01.
Table A7. The power h(c1) at d = 4, α = 0.01.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.01000.02240.05090.11540.25300.50770.83610.9961
0.0750.01000.02220.05000.11250.24610.49570.82590.9953
0.10.01000.02200.04900.10960.23930.48360.81510.9945
500.050.01000.03230.09680.25490.54270.85600.99101.0000
0.0750.01000.03190.09470.24840.53130.84720.98971.0000
0.10.01000.03140.09250.24190.51970.83780.98821.0000
750.050.01000.04200.14770.39740.74680.96370.99961.0000
0.0750.01000.04130.14420.38820.73600.96010.99951.0000
0.10.01000.04070.14070.37880.72480.95610.99941.0000
1000.050.01000.05180.20150.52640.86670.99131.00001.0000
0.0750.01000.05100.19660.51560.85860.99011.00001.0000
0.10.01000.05010.19170.50470.85010.98871.00001.0000
4250.050.01000.02230.05010.11200.24300.48660.81460.9941
0.0750.01000.02170.04760.10470.22540.45510.78450.9910
0.10.01000.02110.04530.09780.20870.42410.75180.9865
500.050.01000.03200.09450.24600.52370.83830.98791.0000
0.0750.01000.03080.08900.22940.49330.81220.98291.0000
0.10.01000.02970.08370.21350.46310.78370.97641.0000
750.050.01000.04140.14360.38390.72770.95600.99931.0000
0.0750.01000.03970.13460.35980.69760.94410.99891.0000
0.10.01000.03810.12600.33630.66600.92980.99811.0000
1000.050.01000.05100.19560.51010.85190.98861.00001.0000
0.0750.01000.04870.18300.48150.82810.98420.99991.0000
0.10.01000.04650.17100.45310.80200.97840.99991.0000
6250.050.01000.02180.04790.10520.22610.45520.78350.9907
0.0750.01000.02090.04430.09450.20020.40680.73100.9827
0.10.01000.02000.04090.08490.17680.36070.67380.9696
500.050.01000.03100.08950.23020.49380.81160.98261.0000
0.0750.01000.02930.08120.20560.44650.76580.97131.0000
0.10.01000.02770.07380.18310.40090.71520.95490.9999
750.050.01000.03990.13530.36080.69770.94360.99881.0000
0.0750.01000.03740.12190.32410.64770.91990.99741.0000
0.10.01000.03510.10990.28990.59600.88980.99481.0000
1000.050.01000.04900.18390.48250.82800.98400.99991.0000
0.0750.01000.04550.16530.43790.78610.97400.99981.0000
0.10.01000.04240.14830.39520.73960.95970.99941.0000
Table A8. The power h(c1) at d = 4, α = 0.05.
Table A8. The power h(c1) at d = 4, α = 0.05.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.05000.09980.19380.35550.59060.84000.98081.0000
0.0750.05000.09910.19140.35030.58280.83330.97901.0000
0.10.05000.09840.18900.34500.57480.82620.97701.0000
500.050.05000.12900.29290.55350.82660.97410.99961.0000
0.0750.05000.12780.28880.54600.81970.97190.99951.0000
0.10.05000.12660.28460.53830.81260.96960.99941.0000
750.050.05000.15480.38010.69380.92710.99571.00001.0000
0.0750.05000.15320.37460.68580.92270.99511.00001.0000
0.10.05000.15150.36900.67770.91800.99451.00001.0000
1000.050.05000.17890.45740.79100.96910.99931.00001.0000
0.0750.05000.17680.45080.78370.96660.99911.00001.0000
0.10.05000.17470.44420.77610.96400.99901.00001.0000
4250.050.05000.09900.19030.34640.57450.82400.97580.9999
0.0750.05000.09710.18410.33270.55350.80470.96970.9999
0.10.05000.09530.17800.31940.53250.78420.96240.9998
500.050.05000.12760.28680.54010.81220.96870.99941.0000
0.0750.05000.12440.27590.52000.79260.96150.99901.0000
0.10.05000.12130.26540.50000.77200.95310.99851.0000
750.050.05000.15280.37180.67950.91770.99431.00001.0000
0.0750.05000.14850.35730.65790.90450.99221.00001.0000
0.10.05000.14420.34310.63580.88990.98960.99991.0000
1000.050.05000.17640.44750.77770.96380.99891.00001.0000
0.0750.05000.17090.43020.75730.95590.99841.00001.0000
0.10.05000.16550.41320.73600.94670.99761.00001.0000
6250.050.05000.09730.18430.33280.55280.80320.96900.9999
0.0750.05000.09440.17480.31200.51990.77050.95680.9998
0.10.05000.09160.16590.29230.48750.73570.94130.9994
500.050.05000.12460.27630.52000.79190.96090.99901.0000
0.0750.05000.11980.25980.48870.75900.94690.99811.0000
0.10.05000.11520.24430.45810.72420.92960.99651.0000
750.050.05000.14880.35780.65780.90400.99201.00001.0000
0.0750.05000.14210.33560.62300.88030.98750.99991.0000
0.10.05000.13590.31450.58790.85330.98100.99981.0000
1000.050.05000.17130.43080.75720.95560.99831.00001.0000
0.0750.05000.16290.40410.72330.94050.99701.00001.0000
0.10.05000.15510.37850.68810.92200.99481.00001.0000
Table A9. The power h(c1) at d = 4, α = 0.1.
Table A9. The power h(c1) at d = 4, α = 0.1.
c1
mnp0.80.8250.850.8750.90.9250.950.975
2250.050.10000.18440.32360.52430.75470.93200.99571.0000
0.0750.10000.18330.32050.51890.74860.92850.99521.0000
0.10.10000.18220.31740.51340.74230.92490.99471.0000
500.050.10000.22620.44090.70760.91460.99180.99991.0000
0.0750.10000.22450.43620.70130.91060.99110.99991.0000
0.10.10000.22280.43150.69480.90630.99020.99991.0000
750.050.10000.26130.53290.81690.96870.99891.00001.0000
0.0750.10000.25910.52730.81110.96660.99871.00001.0000
0.10.10000.25680.52170.80510.96430.99861.00001.0000
1000.050.10000.29260.60800.88400.98820.99981.00001.0000
0.0750.10000.29000.60180.87920.98710.99981.00001.0000
0.10.10000.28720.59560.87410.98590.99981.00001.0000
4250.050.10000.18280.31810.51310.74020.92250.99421.0000
0.0750.10000.18000.31010.49880.72340.91200.99251.0000
0.10.10000.17720.30220.48450.70610.90050.99031.0000
500.050.10000.22380.43300.69490.90510.98970.99991.0000
0.0750.10000.21940.42070.67770.89320.98690.99991.0000
0.10.10000.21500.40870.66020.88030.98360.99981.0000
750.050.10000.25820.52350.80540.96370.99851.00001.0000
0.0750.10000.25240.50860.78920.95690.99781.00001.0000
0.10.10000.24670.49380.77230.94920.99701.00001.0000
1000.050.10000.28900.59770.87440.98560.99981.00001.0000
0.0750.10000.28190.58120.86060.98210.99961.00001.0000
0.10.10000.27490.56470.84590.97780.99941.00001.0000
6250.050.10000.18010.31010.49820.72210.91080.99221.0000
0.0750.10000.17570.29780.47590.69480.89210.98851.0000
0.10.10000.17160.28610.45430.66710.87140.98341.0000
500.050.10000.21960.42080.67720.89240.98660.99981.0000
0.0750.10000.21270.40190.64950.87160.98090.99971.0000
0.10.10000.20620.38380.62180.84870.97360.99941.0000
750.050.10000.25270.50880.78880.95650.99781.00001.0000
0.0750.10000.24370.48550.76190.94380.99631.00001.0000
0.10.10000.23520.46290.73400.92880.99411.00001.0000
1000.050.10000.28230.58150.86030.98190.99961.00001.0000
0.0750.10000.27140.55540.83670.97480.99921.00001.0000
0.10.10000.26100.52980.81130.96570.99861.00001.0000

References

  1. Juran, J.M. Juran’s Quality Control Handbook, 3rd ed.; McGraw-Hill: New York, NY, USA, 1974. [Google Scholar]
  2. Kane, V.E. Process capability indices. J. Qual. Technol. 1986, 18, 41–52. [Google Scholar] [CrossRef]
  3. Hsiang, T.C.; Taguchi, G. A tutorial on quality control and assurance—The Taguchi methods. In Joint Meetings of the American Statistical Association; InASA Annual Meeting LA: Las Vegas, NV, USA, 1985; p. 188. [Google Scholar]
  4. Pearn, W.L.; Kotz, S.; Johnson, N.L. Distributional and inferential properties of process capability indices. J. Qual. Technol. 1992, 24, 216–231. [Google Scholar] [CrossRef]
  5. Montgomery, D.C. Introduction to Statistical Quality Control; John Wiley and Sons Inc.: New York, NY, USA, 1985. [Google Scholar]
  6. Tong, L.I.; Chen, K.S.; Chen, H.T. Statistical testing for assessing the performance of lifetime index of electronic components with exponential distribution. Int. J. Qual. Reliab. Manag. 2002, 19, 812–824. [Google Scholar] [CrossRef]
  7. Balakrishnan, N. Progressive censoring methodology: An appraisal. Test 2007, 16, 211–259. [Google Scholar] [CrossRef]
  8. Aggarwala, R. Progressive interval censoring: Some mathematical results with applications to inference. Commun. Stat. Theory Methods 2001, 30, 1921–1935. [Google Scholar] [CrossRef]
  9. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods and Applications; Birkhäuser: Boston, MA, USA, 2000. [Google Scholar]
  10. Wu, J.W.; Lee, W.C.; Lin, L.S.; Hong, M.L. Bayesian test of lifetime performance index for exponential products based on the progressively type II right censored sample. J. Quant. Manag. 2011, 8, 57–77. [Google Scholar]
  11. Sanjel, D.; Balakrishnan, N. A Laguerre polynomial approximation for a goodness-of-fit test for exponential distribution based on progressively censored data. J. Stat. Comput. Simul. 2008, 78, 503–513. [Google Scholar] [CrossRef]
  12. Lee, W.C.; Wu, J.W.; Hong, C.W. Assessing the lifetime performance index of products with the exponential distribution under progressively type II right censored samples. J. Comput. Appl. Math. 2009, 231, 648–656. [Google Scholar] [CrossRef]
  13. Wu, S.F.; Lin, M.J. Computational testing algorithmic procedure of assessment for lifetime performance index of products with weibull distribution under progressive type I interval censoring. J. Comput. Appl. Math. 2017, 311, 364–374. [Google Scholar] [CrossRef]
  14. Wu, S.F.; Chen, Z.C.; Chang, W.J.; Chang, C.W.; Lin, C. A hypothesis testing procedure for the evaluation on the lifetime performance index of products with Burr XII distribution under progressive type I interval censoring. Commun. Stat. Simul. Comput. 2018, 47, 2670–2683. [Google Scholar] [CrossRef]
  15. Wu, S.F.; Chang, W.T. The evaluation on the process capability index CL for exponentiated Frech’et lifetime product under progressive type I interval censoring. Symmetry 2021, 13, 1032. [Google Scholar] [CrossRef]
  16. Wu, S.F.; Chiang, K.Y. Assessment of the overall lifetime performance index of Weibull products in multiple production lines. Mathematics 2024, 12, 514. [Google Scholar] [CrossRef]
  17. Wu, S.F.; Kuo, P.H.; Deng, W.S. The computational testing procedure for the comprehensive lifetime performance index of Burr XII products in multiple production lines. Mathematics 2024, 12, 584. [Google Scholar] [CrossRef]
  18. Rytgaard, H.C.W.; Van der Laan, M. Targeted maximum likelihood estimation for causal inference in survival and competing risks analysis. Lifetime Data Anal. 2024, 30, 4–33. [Google Scholar] [CrossRef]
  19. Luca, B.; Fulvia, P.; Francesco, B. Maximum likelihood estimation for discrete latent variable models via evolutionary algorithms. Stat. Comput. 2024, 34, 62. [Google Scholar]
  20. Goldengorin, B.; Golovnya, N.; Kudinov, O.; Volkova, E. Reliability: Measurement, Theory, and Practice. Qual. Reliab. Eng. Int. 2020, 36, 505–506. [Google Scholar]
  21. Zhuang, L.; Xu, A.; Wang, Y.; Tang, Y. Remaining useful life prediction for two-phase degradation model based on reparameterized inverse Gaussian process. Eur. J. Oper. Res. 2024. [Google Scholar] [CrossRef]
  22. Xu, A.; Wang, B.; Zhu, D.; Pang, J.; Lian, X. Bayesian reliability assessment of permanent magnet brake under small sample size. IEEE Trans. Reliab. 2024. [Google Scholar] [CrossRef]
  23. Pedro, L.R.; Francisco, L.; Eduardo, R.; Sanku, D. The Fréchet distribution: Estimation and application—An overview. J. Stat. Manag. Syst. 2020, 23, 549–578. [Google Scholar]
  24. Tabassum, N.S.; Muhammad, A.; Navid, F. On Estimation of the Fr´echet Distribution with Known Shape under Different Sample Schemes. Int. J. Intell. Technol. Appl. Stat. 2014, 7, 353–368. [Google Scholar]
  25. Casella, G.; Berger, R.L. Statistical Inference, 2nd ed.; Duxbury Press: Pacific Grove, CA, USA, 2002. [Google Scholar]
  26. Lawless, J.F. Statistical Models and Methods for Lifetime Data, 2nd ed.; John Wiley and Sons: New York, NY, USA, 2003. [Google Scholar]
  27. Gail, M.H.; Gastwirth, J.L. A scale-free goodness of fit test for the exponential distribution based on the Gini Statistic. J. R. Stat. Soc. B 1978, 40, 350–357. [Google Scholar] [CrossRef]
  28. Aarset, M.V. How to identify a bathtub hazard rate. IEEE Trans. Reliab. 1987, 36, 106–110. [Google Scholar] [CrossRef]
Figure 1. (a) The p.d.f. for   δ i = 2, 4, 6 under θ i = 1. (b) The p.d.f. for θ i = 1, 3, 5 under   δ i = 1.
Figure 1. (a) The p.d.f. for   δ i = 2, 4, 6 under θ i = 1. (b) The p.d.f. for θ i = 1, 3, 5 under   δ i = 1.
Symmetry 16 01060 g001
Figure 2. (a) The h.f. for   δ i = 2, 4, 6 under θ i = 1. (b) The h.f. for θ i = 1, 3, 5 under   δ i = 1.
Figure 2. (a) The h.f. for   δ i = 2, 4, 6 under θ i = 1. (b) The h.f. for θ i = 1, 3, 5 under   δ i = 1.
Symmetry 16 01060 g002
Figure 3. The power curve for d = 1, α = 0.05, n = 25, m = 4, and p = 0.050, 0.075, 0.100.
Figure 3. The power curve for d = 1, α = 0.05, n = 25, m = 4, and p = 0.050, 0.075, 0.100.
Symmetry 16 01060 g003
Figure 4. The power curve for d = 1, α = 0.05, m = 4, p = 0.05, and n = 25, 50, 75, 100.
Figure 4. The power curve for d = 1, α = 0.05, m = 4, p = 0.05, and n = 25, 50, 75, 100.
Symmetry 16 01060 g004
Figure 5. The power curve for d = 1, α = 0.05, n = 25, p = 0.05, and m = 2, 4, 6, 8.
Figure 5. The power curve for d = 1, α = 0.05, n = 25, p = 0.05, and m = 2, 4, 6, 8.
Symmetry 16 01060 g005
Figure 6. The power curve for d = 1, n = 25, m = 4, p = 0.05, and α = 0.01, 0.05, 0.10.
Figure 6. The power curve for d = 1, n = 25, m = 4, p = 0.05, and α = 0.01, 0.05, 0.10.
Symmetry 16 01060 g006
Figure 7. The power curve for n = 25, m = 4, p = 0.05, α = 0.05, and d = 2, 3, 4, 5.
Figure 7. The power curve for n = 25, m = 4, p = 0.05, α = 0.05, and d = 2, 3, 4, 5.
Symmetry 16 01060 g007
Figure 8. The p-values vs. the values of δ1.
Figure 8. The p-values vs. the values of δ1.
Symmetry 16 01060 g008
Figure 9. The p-values vs. the values of δ2.
Figure 9. The p-values vs. the values of δ2.
Symmetry 16 01060 g009
Table 1. The corresponding values of C L for a given value of P r   a n d   C T .
Table 1. The corresponding values of C L for a given value of P r   a n d   C T .
Pr
0.6873 0.7047 0.7225 0.7408 0.7596 0.7788 0.7985 0.8187
CT
d0.6250 0.6500 0.6750 0.7000 0.7250 0.7500 0.7750 0.8000
20.8125 0.8250 0.8375 0.8500 0.8625 0.8750 0.8875 0.9000
30.8750 0.8833 0.8917 0.9000 0.9083 0.9167 0.9250 0.9333
40.9063 0.9125 0.9188 0.9250 0.9313 0.9375 0.9438 0.9500
50.9250 0.9300 0.9350 0.9400 0.9450 0.9500 0.9550 0.9600
60.9375 0.9417 0.9458 0.9500 0.9542 0.9583 0.9625 0.9667
70.9464 0.9500 0.9536 0.9571 0.9607 0.9643 0.9679 0.9714
Pr
0.8395 0.8607 0.8825 0.9048 0.9277 0.9512 0.9753 1.0000
CT
d0.8250 0.8500 0.8750 0.9000 0.9250 0.9500 0.9750 1.0000
20.9125 0.9250 0.9375 0.9500 0.9625 0.9750 0.9875 1.0000
30.9417 0.9500 0.9583 0.9667 0.9750 0.9833 0.9917 1.0000
40.9563 0.9625 0.9688 0.9750 0.9813 0.9875 0.9938 1.0000
50.9650 0.9700 0.9750 0.9800 0.9850 0.9900 0.9950 1.0000
60.9708 0.9750 0.9792 0.9833 0.9875 0.9917 0.9958 1.0000
70.9750 0.9786 0.9821 0.9857 0.9893 0.9929 0.9964 1.0000
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

Wu, S.-F.; Yeh, H.-C. Computational Analysis of the Comprehensive Lifetime Performance Index for Exponentiated Fréchet Lifetime Distribution Products with Multi-Components. Symmetry 2024, 16, 1060. https://doi.org/10.3390/sym16081060

AMA Style

Wu S-F, Yeh H-C. Computational Analysis of the Comprehensive Lifetime Performance Index for Exponentiated Fréchet Lifetime Distribution Products with Multi-Components. Symmetry. 2024; 16(8):1060. https://doi.org/10.3390/sym16081060

Chicago/Turabian Style

Wu, Shu-Fei, and Hsueh-Chien Yeh. 2024. "Computational Analysis of the Comprehensive Lifetime Performance Index for Exponentiated Fréchet Lifetime Distribution Products with Multi-Components" Symmetry 16, no. 8: 1060. https://doi.org/10.3390/sym16081060

APA Style

Wu, S. -F., & Yeh, H. -C. (2024). Computational Analysis of the Comprehensive Lifetime Performance Index for Exponentiated Fréchet Lifetime Distribution Products with Multi-Components. Symmetry, 16(8), 1060. https://doi.org/10.3390/sym16081060

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