Next Article in Journal
Comprehensive Neural Cryptanalysis on Block Ciphers Using Different Encryption Methods
Next Article in Special Issue
A Network Reliability Analysis Method for Complex Real-Time Systems: Case Studies in Railway and Maritime Systems
Previous Article in Journal
A Range Query Scheme for Spatial Data with Shuffled Differential Privacy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fuzzy Evaluation Model for Lifetime Performance Using Type-I Censoring Data

1
Department of Finance, Chaoyang University of Technology, Taichung 413310, Taiwan
2
Department of Industrial Engineering and Management, National Chin-Yi University of Technology, Taichung 411030, Taiwan
3
Department of Business Administration, Chaoyang University of Technology, Taichung 413310, Taiwan
4
Department of Business Administration, Asia University, Taichung 413305, Taiwan
*
Authors to whom correspondence should be addressed.
Mathematics 2024, 12(13), 1935; https://doi.org/10.3390/math12131935
Submission received: 16 April 2024 / Revised: 13 June 2024 / Accepted: 19 June 2024 / Published: 21 June 2024

Abstract

:
As global warming becomes increasingly serious, humans start to consider how to coexist with the natural environment. People become more and more aware of environmental protection and sustainable development. Therefore, in the pursuit of economic growth, it has become a consensus that enterprises should be responsible for the social and ecological environment. Regarding the manufacturing of electronic devices, as long as both component production quality and assembly quality are ensured, consumers can be provided with high-quality, safe, and efficient products. In light of this trend, enhancing product availability and reliability can reduce costs and carbon emissions resulting from repairing or replacing components, thus becoming a vital factor for corporate and environmental sustainability. Accordingly, enterprises enhance their economic benefits as well as have the effects of energy conservation and waste reduction by extending products’ service lifetime and increasing their added value. According to several studies, it takes a long time to retrieve electronic products’ lifetime data. Moreover, acquiring complete samples is often challenging. Consequently, when analyzing real cases, samples are usually collected using censoring techniques. The type-I right censoring data is suitable for industrial processes. Thus, this study utilized type-I right censoring sample data to estimate the lifetime performance index. It usually takes a large amount of time to access lifetime data for electronic products and it is often impossible to obtain complete samples since the size of the sample is usually small. Hence, to avoid misjudgment caused by sampling errors, this study followed suggestions from existing research and applied fuzzy tests built on confidence intervals to establish a fuzzy evaluation model for the lifetime performance index. This model helps relevant electronic industries not only evaluate the lifetime of their electronic components but also instantly seize opportunities for improvement.

1. Introduction

According to a number of studies, Taiwan’s electronics industry owns a complete ecological chain in the world supply chain system involving information and communications technology, driven by a high clustering effect. As a result, it has played a pivotal role in the world electronics industry [1,2,3,4]. In the production of electronic products, it is essential to ensure both production quality and assembly quality of all components. This assures consumers of receiving high-quality, safe, and efficient products, thereby satisfying their demands [5,6]. Additionally, plenty of studies have highlighted the increasing seriousness of global warming, urging humans to contemplate how to coexist with the natural environment. Meanwhile, awareness of environmental protection and sustainable development is gradually growing as well [7,8]. Consequently, corporate social responsibility (CSR) has emerged as the business philosophy that global enterprises promote. In the pursuit of economic growth, enterprises have reached a consensus that they must take greater responsibility for the social and ecological environment [9,10]. Under such a trend, boosting product availability and reliability can cut down on costs and carbon emissions caused by repairs or replacement of components, which is one of the key factors in the sustainable operation of enterprises and the environment. In this era of artificial intelligence, constant changes in technology have made manufacturing techniques sophisticated and complicated. Therefore, enterprises strive to extend the service lifetime of products as well as to increase the added value of products to increase their overall competitiveness and economic benefits.
As noted by some studies, the lifetime T of an electronic component has an exponential distribution with mean θ   1 [11,12,13]. Thus, the probability density function (p.d.f.) is f T ( t ) = θ   e θ t for t > 0 and the cumulative function (c.d.f.) of T is F T ( t ) =   1 e θ t for t > 0. Building upon this assumption, some studies have proposed a relative lifetime performance index to offer the industry a more intuitive and user-friendly tool in practice [14,15,16]. According to Chiou [17] and Chen and Yu [16], the relative electronic component lifetime is denoted with X = T/L, where L represents the period of warranty. Thus, the relative electronic component lifetime X possesses a Gamma distribution with two parameters, 1 and λ L , expressed as X~ G ( 1 ,   λ L ) . The second parameter λ L = L θ denotes the relative lifetime performance index for an electronic component where parameter θ refers to the electronic component’s mean time between failures (MTBF). Therefore, the relative electronic component lifetime X has an exponential distribution with the mean λ L . Its probability density function (p.d.f.) is f X x = λ L 1 × exp x / λ L for x > 0 and the cumulative distribution function (c.d.f.) is F X x = 1 exp x / λ L , for x > 0. According to Chen et al. [18], the relative failure rate is denoted by λ L 1 , and the reliability function of the relative lifetime X is R X ( x ) = exp x / λ L for x > 0.
Based on Chen et al. [18], product reliability ( p r ) represents the probability showing that the relative lifetime X is greater than 1, such that p r = e λ L 1 , indicating its direct mathematical relationship with the relative lifetime performance index λ L . Obviously, when the value of θ   1 (= μ T ) is larger than that of L, then the relative lifetime performance index is larger. Additionally, p r increases as the value of the indicator λ L rises. Therefore, the relative lifetime performance index λ L can reveal the lifetime performance of the electronic component in a reasonable manner.
If the relative lifetime performance index λ L is employed as an evaluation tool to establish an evaluation model for the electronic component lifetime, it will help the electronics industry improve the lifetime performance of electronic components through evaluation. Since this index has unknown parameters, it must be evaluated using sample data. However, Chiou [17] stated that it is time-consuming to obtain electronic products’ lifetime data. In addition, often complete samples cannot be obtained. For example, not all observation values can be acquired due to constraints such as time, cost, or human negligence. Accordingly, when analyzing an actual case, it is common to encounter difficulty in obtaining a complete sample. Thus, samples are usually collected using censoring techniques. Censoring falls into three categories: type-I censoring, type-II censoring, and random censoring [19]. In many cases, data of type-II and type-I censoring are frequently applied to engineering [19]. This paper explores the study by Miller [19], emphasizing the right censoring type. This involves dividing the obtained data into two parts: the first part represents the individual’s failure time X (focusing on the relative failure time), while the second part represents the censoring time C. If the failure time X of the experiment itself is less than or equivalent to the censoring time C, then the data for X is the data for uncensored observation. If the failure time X of the experiment itself is greater than the censoring time C, then X refers to the data for censored observation, symbolized by X + . Concerning the censored observation data ( X + ), the corresponding censoring time C is used as the imputation value of the individual’s failure time X, and this imputation value is called censored data. When the censoring time C is the order statistic X(m), where m ≤ n, the censoring type is called type-II right censoring. As the censoring time C is the fixed censoring time xc, the censoring type is called type-I right censoring. When the censoring time C is a random variable (independent failure time X), the censoring type is called random right censoring. According to several studies, of the above three censoring types, random right censoring arises in medical applications with clinical trials or animal studies. Conversely, type-I right censoring data is suitable for industrial processes [19,20,21,22]. Therefore, in this paper, we utilize type-I right censoring sample data to estimate the lifetime performance index. As mentioned above, it is laborious to retrieve electronic products’ lifetime data.
If the sample size is large enough, the confidence interval of index λ L is used to statistically test the life performance of electronic products. Since the length of the confidence interval will shorten as the sample size increases, the risk of misjudgment can be reduced. If the sample size is not large enough, the length of the confidence interval will be too long, which will make the statistical test error too large and increase the risk of misjudgment. Moreover, obtaining a complete sample is often challenging, particularly in the case of a small sample size [23]. In order to solve the risk of misjudgment caused by too-long confidence intervals caused by small samples, this study follows some suggestions from previous studies and adopts fuzzy tests built on confidence intervals to establish a fuzzy evaluation model for the lifetime performance index [16,24]. Fuzzy testing based on confidence intervals is not traditional fuzzy testing. Its advantage is that it can incorporate past accumulated data and expert experience. Therefore, the risk of misjudgment caused by too-long confidence intervals can be reduced. Obviously, since obtaining electronic product lifetime data is very time-consuming, in the case of a small sample, this model can assist relevant electronics industries with evaluating the lifetime of their electronic components as well as seizing opportunities for improvement in real time.
The remaining sections are organized in the following order. In Section 2, this study estimates the lifetime performance index and derives its 100 ( 1 α )% confidence intervals using type-I right censoring sample data. Subsequently, in Section 3, these 100 ( 1 α )% confidence intervals are employed to establish a fuzzy evaluation model for the lifetime performance index, and then testing rules are established. In Section 4, this study proposes an empirical case to clarify the fuzzy evaluation model so as to benefit its application and promotion in the industry. Finally, Section 5 gives the conclusions.

2. Confidence Intervals of the Relative Lifetime Performance Index

As mentioned above, the relative lifetime X of electronic components is distributed as an exponential distribution with mean λ L , where λ L denotes the relative lifetime performance index for an electronic component as follows:
λ L =   μ T L = θ   1 L = 1 L θ ,
When reliability is adopted to analyze an empirical case, data collection is incomplete for external factors or human factors. Consequently, the results obtained after analysis are often questionable. Considering labor costs and time, applying censoring types to data collection is one of the methods that can improve this shortcoming. The type-I progressive censoring scheme has become relatively widespread for analyzing lifetime data of highly reliable products [25,26,27]. In this paper, we took type-I right censoring (i.e., given censoring time) into consideration. According to Equation (2), as the type of censoring data for fixed C = x c is called type-I right censoring, then Y i can be rewritten as follows:
Y i = X i   , X i     x c   ( uncensored   data ) x c   ,   X i >   x c   ( censored   data ) ,
where i = 1, 2, …, n. Let δ i = I   ( X i   x c ) be an indicator function, expressed as follows:
δ i = I   ( X i   x c )   =   1   ,   X i   x c 0   ,   X i >   x c   .
Next, let mu be the number of uncensored data, defined as follows:
m u = i = 1 n δ i   =   i = 1 n I   ( X i   x c )    .
In fact, mu follows a Binomial distribution, denoted by mu~B(n, p), where p = P(Xxc). According to the probability density function and the cumulative distribution function of X, we can find the distribution of estimate λ ^ L using type-I right censoring. Let the likelihood function of pair ( y i , δ i ) be as follows:
L ( y i , δ i ) = f X ( y i ) δ i R X ( y i ) 1 δ i ,   δ i   =   0.1 ,
where y1, y2, …, yn are the type-I censored observation data for Y1, Y2, …, Yn. Let λ L 1 = θ . Then, the likelihood function of type-I right censoring is expressed as follows:
L ( y 1 ,   y 2 , , y n ; δ 1 ,   δ 2 , , δ n ) = i =   1 n L ( y i ; δ i ) = x i    U f X ( y i ) x i   C * R X ( y i )   = θ m u e U θ * x i e C * θ *   x c = θ m u e   U θ *   x i + C * θ *   x c   = θ * m u e   i = 1 n θ *   t i ,
where U represents the uncensored data and C* represents the censored data. mu represents the number of uncensored data as well as a random variable. Using Equation (6), the first- and second-order derivatives can be found below.
d d θ * ln L ( θ *   |   y 1 ,   y 2 , , y n ) = d d θ * ln ( θ * ) m u   e i = 1 n θ *   y i   =   m u θ * i = 1 n y i ;
d 2 d θ * 2 ln L ( θ *   |   y 1 ,   y 2 , , y n ) =     m u θ * 2 .
Setting Equation (7) to 0, we can find the estimator of θ as follows:
θ ^ * =   m u i = 1 n Y i .
Using Equation (8), the Fisher information of θ ^ * is defined as follows:
I ( θ * )   = E d 2 d θ * 2 ln L ( θ *   |   Y 1 ,   Y 2 , , Y n ) =   m u θ * 2 .
Based on the large-sample theorem, the asymptotic variance of the estimate θ ^ * can be expressed as:
Var ( θ ^ * )   =   I ( θ * ) 1   =   θ * 2 m u .
Then, the asymptotic distribution of estimate θ ^ * is presented as follows:
m u ( θ ^ * θ * )   N 0 ,   θ * 2 .
When xc → ∞ and mun, the estimate θ ^ * is defined as
θ ^ *   n i = 1 n X i = 1 X ¯ .
Given that the censoring time is xc and (Y1, Y2, …, Yn) is a random sample, the number of uncensored data m u = i = 1 n δ i in Equation (4) is unknown, and it follows a Binomial distribution (n, p) with p = P(X ≤ xc) = 1 e   θ * x c . Next, we can get the expected value E(mu) = np = n 1 e   θ * x c .
As noted above, this study needs to estimate the lifetime performance index λ L using the sample data since its real value is unknown. Let the electronic component lifetime T have an exponential distribution with mean θ 1 . Since the relative lifetime X = T/L, X is distributed as an exponential distribution with mean ( θ * ) 1 . Assume that (X1, X2, …, Xn) is a random sample set of X, while the type-I right censoring data (Y1, Y2, …, Yn) is a random sample set of Y (given censoring time xc).
In Equation (9), using the property of the invariance for the maximum likelihood estimator (MLE), we can find the MLE of index λ L as follows:
λ ^ L = 1 θ ^ * = i = 1 n Y i / m u .
Let the lifetime performance index λ L be the function h( θ ), h( θ ) = θ 1   . Using Equation (12) and the delta method (Casella and Berger [28]), we can get the asymptotic normal distribution for the estimator of λ ^ L = h( θ ^ * ) as follows:
m u ( h ( θ ^ * ) h ( θ * ) )   =   m u ( λ ^ L λ L )   N 0 ,   1 θ * 2 ,
where the asymptotic variance of λ ^ L (=h( θ ^ * )) is 1 m u   θ * 2 = λ L 2 / m u . Thus, we have the following equation:
Z   =   h ( θ ^ * ) h ( θ * ) 1 m u   θ * = λ ^ L λ L λ L m u     N 0 ,   1 .
Given
lim n E λ ^ L λ L 2 = lim n V a r λ ^ L = lim m u λ L 2 m u   =   0 ,
the estimator of λ ^ L is unbiased and consistent with λ L , displayed as follows:
E λ ^ L =   λ L .
According to the Lehmann–Scheffé theorem, λ ^ L represents the uniformly minimum variance unbiased estimator (UMVUE) of λ L . Using Equation (16), we can get the asymptotic standard normal distribution. The random variable Z, as noted above, is displayed in the following equation:
p   z α / 2 Z z α / 2   =   p   z α / 2 λ ^ L λ L λ L / m u z α / 2 = p   z α / 2 m u λ ^ L λ L 1 z α / 2 = 1 α .
Thus, the (1 − α) × 100% confidence interval of λ L is L ( λ ^ L ) , U ( λ ^ L ) , where the functions L ( λ ^ L ) and U ( λ ^ L ) are written as follows:
L ( λ ^ L ) = λ ^ L 1 + z α / 2 / m u  
and
U ( λ ^ L ) = λ ^ L 1 z α / 2 / m u   ,
where z α / 2 represents the upper α / 2 quantile of the standard normal distribution. As mentioned above, the fuzzy test proposed in this article is not a traditional fuzzy test. It is based on the confidence interval. Under the same sample size, the shorter the length of the confidence interval, the smaller the error of the interval estimate, and the uniformly minimum variance unbiased estimator (UMVUE) is the estimator with the shortest confidence interval among the unbiased estimators, which can improve the accuracy of the test. Let the confidence interval length of λ L be
l λ ^ L = λ ^ L 1 1 z α / 2 / m u 1 1 + z α / 2 / m u .
Then the expected value of l λ ^ L is written as follows:
E l λ ^ L = λ L 1 1 z α / 2 / m u 1 1 + z α / 2 / m u .
The electronic component lifetime (T) has an exponential distribution with mean θ 1 , while the relative lifetime X (=T/L) has an exponential distribution with mean θ * 1 . Let P   k be the kth percentile of the relative lifetime X, k = 1, 2, …,99, such that P(X P   k ) = k%. Then, the kth percentile ( P   k ) is expressed as follows:
P   k =   ln ( 100 k ) % θ * .
Given the confidence interval level (1 − α) × 100% = 95%, sample size n = 100, relative lifetime performance index λ L = 1(1)5 and censoring time xc = P   k with k = 10(10)100, the number of uncensored data mu is unknown, and then mu is replaced with E(mu). When the censoring time is denoted as xc = P   k , then the value of E(mu) is received as follows:
E ( m u ) = n p = n ×   ( X   x c ) = n ×   P ( Y   P   k )   =   n   ×   ( k ) % .
When the number of samples, n = 100, is fixed, based on Equation (22), we can get the distribution for E l λ ^ L , as shown in Figure 1. As λ L is fixed, the larger the censoring time xc (= P   k ) is, the smaller the average confidence interval length of E l λ ^ L is, which means the estimate of the index λ L is better. That is, when the number of uncensored data mu is higher, the estimated indicator λ L is better.
When the k value of quantile P   k is fixed, both the sample number n and the value of E(mu) become larger. When other conditions remain unchanged, the sample number n gets larger, whereas the value of E l λ ^ L gets smaller, indicating the estimate of the index λ L is better, as shown in Figure 2.
Based on Equations (23) and (25), the expected value E l λ ^ L of the confidence interval length for the relative lifetime performance index λ L is related to the number of uncensored observations mu. When mu is larger, then E l λ ^ L is smaller, that is, the estimate of index λ L is better. This result can be verified and illustrated in Figure 1 and Figure 2. These findings, demonstrated by interval estimates in statistics, carry significant implications. Consequently, these factors can be dominated by research to boost interval estimates’ accuracy.

3. Fuzzy Testing Method Built on Confidence Intervals

According to [16], when the sample size is large enough, due to the short length of the confidence interval, the statistical test for the lifetime performance of electronic products can be directly performed using the confidence interval of index λ L . Since it is very time-consuming to obtain the lifetime data of electronic products, the sample size is usually not large enough and the confidence interval is wide. In order to reduce the risk of misjudgment caused by small samples, this article will use the confidence interval of index λ L derived from Section 3 to construct the confidence interval-based fuzzy test. The observed value of this UMVUE is 50% of the quantile of this confidence interval. As mentioned above, the length of its confidence interval is the shortest estimator among unbiased estimators, which can not only improve the accuracy of statistical tests but also improve fuzzy tests based on confidence intervals.
Then, this section presents a fuzzy testing method—null hypothesis H 0 : λ L    e versus alternative hypothesis H 1 : λ L   <  e—to evaluate if the lifetime performance index λ L has reached the desired level. Before performing this method, we examine the testing rules of statistical hypotheses as follows:
(1)
Given that λ ^ L     d c , do not reject H 0 and conclude that λ L    e.
(2)
Given that λ ^ L   <   d c , reject H 0 and conclude that λ L e.
In the rules, the critical value d 0 is obtained using Equation (16). Then, the significant level α is received and expressed as
α   = m i n   P λ ^ L < d c λ L H 0 =   P Z < d c λ L λ L / m u λ L = e .
Next, we obtain the critical value d c as follows:
d c   =   e z α e m u   .
Using Equations (2), (9), and (14), the observed value of estimate λ ^ L is calculated using the following equation:
λ ^ L 0 = 1 θ ^ *   =   i = 1 n y i /   m u .
Based on Buckley’s approach [29], Equations (20) and (21), the α - cuts of the triangular fuzzy number λ ^ L 0 is denoted by the following equation [16,17]:
λ ˜ L 0 α = λ L 01 α , λ L 02 α ,   for   0.01 α 1 λ L 01 0.01 , λ L 02 0.01 ,   for   0 α 0.01 ,
where
λ L 01 α = λ ^ L 0 1 + z α / 2 / m u  
and
λ L 02 α = λ ^ L 0 1 z α / 2 / m u   .
Clearly, the value of λ L 01 α is not equivalent to that of λ L 02 α given that α < 1. As α = 1 and − z α / 2 = z 1 α / 2 , we obtain the result λ L 01 1 λ L 02 1 = λ ^ L 0 , where λ ^ L 0 =   λ ^ L 0 / 1 + z 0.5 / m u . Therefore, let λ L 0 = λ ^ L 0 , such that the α - cuts of the triangular fuzzy number λ ˜ L 0 * is defined as
λ ˜ L 0 * α = λ L 01 * α , λ L 02 * α ,   for   0.01 α 1 λ L 01 * 0.01 , λ L 02 * 0.01 ,   for   0 α 0.01 ,
where
λ L 01 * α = 1 1 + z α / 2 / m u   λ ^ L 0  
and
λ L 02 * α = 1 1 z α / 2 / m u λ ^ L 0   .
Obviously, the value of λ L 01 α is equal to that of λ L 02 α with α = 1, and the newly transformed triangular fuzzy number is represented as λ ˜ L 0 * * = Δ λ L 0 * * , λ M 0 * * , λ R 0 * * , where λ M 0 * * = λ ^ L 0 ,
λ L 0 * * = 1 1 + z 0.005 / m u   λ ^ L 0  
and
λ R 0 * * = 1 1 z 0.005 / m u   λ ^ L 0   .
Furthermore, the membership function of the fuzzy number λ ˜ L 0 * * is denoted by
g ( x ) = 0   ,   i f   x < λ L 0 * * 2 1 F Z m u λ ^ L 0 x 1 , i f   λ L 0 * * x <   λ ^ L 0 1   , i f   x = λ ^ L 0 2 1 F Z m u 1 λ ^ L 0 x , i f   λ ^ L 0 < x λ R 0 * * 0   ,   i f   λ R 0 * * < x ,
where F Z is a cumulative distribution function of the standard normal distribution. Similar to the fuzzy number λ ˜ L 0 * , the α - cuts of triangular-shaped fuzzy critical value number d ˜ c is expressed as
d ˜ c α = d c 1 α , d c 2 α ,   f o r   0.01 α 1 d c 1 0.01 , d c 2 0.01 ,   f o r   0 α 0.01 ,
where
d c 1 α = 1 1 + z α / 2 / m u   d c  
and
d c 2 α = 1 1 z α / 2 / m u   d c   .
Clearly, the value of d c 1 α is equivalent to that of d c 2 α given that α = 1  Then, a newly transformed triangular fuzzy number is denoted by d ˜ c   =   Δ   d c L , d c M , d c R , where d c M   =   d c ,
d c L = 1 1 + z 0.005 / m u   d c  
and
d c R = 1 1 z 0.005 / m u   d c   .
Furthermore, the membership function of the fuzzy number d ˜ c is represented as
g c ( x ) = 0 ,   i f   x < d c L 2 1 F Z m u d c x 1 , i f   d c L x <   d c 1 , i f   x = d c 2 1 F Z m u 1 d c x ,   i f   d c < x d c R 0 ,   i f   d c R < x ,
where F Z is a cumulative distribution function of the standard normal distribution. Subsequently, the curves of g ( x ) and g c ( x ) are depicted below:
Based on Buckley [29] and Chiou [17], we let A R be the area in the curved graph of g c ( x ) and be defined as
A R = x , α d c 1 α x d c 2 α , 0 α 1 .
Buckley [29], Chen and Yu [16], and Chiou [17] have stated that it is complicated to utilize integration to find the area of set A R . A Riemann sum is performed in the study to find the area of the block. The procedures are described as follows:
(1) We divide A R with n = 100 into 100 horizontal blocks. (2) Each block is calculated using a quasi-trapezoid area. (3) We find the sum of the areas for these 100 horizontal blocks. Suppose that j = [100 ×   α ], and j = 0, 1, 2, …, 100 for 0   α   1. [100 ×   α ] represents the largest integer, less than or equivalent to 100 ×   α . Similarly, α   =  j × 0.01, j = 0, 1, 2, …, 100. Set A R is cut into 100 trapezoidal blocks by these 101 horizontal lines. Therefore, the area of the jth block is denoted by
A R   j = x , α d c 1 0.01 × j x d c 2 α 0.01 × j , 0.01 × j 1 α 0.01 × j ,   j =   0 ,   1 ,   2 , ,   100 .
Next, the length of the jth horizontal line is depicted as
l r j = 1 1 z 0.005 × j m u 1 1 + z 0.005 × j m u    d c   ,   j =   0 ,   1 ,   2 , ,   100 .
Obviously, when l r 0 = l r 1 and l r 100   = 0, then the area is defined as
A R = j = 1 100 ( 0.01 ) × l r j 1 + l r j 2 = 0.5 × l r 1 + 0.01 × j = 1 99 l r j .
Moreover, let A S be the area in the curved graph of g c ( x ) but to the right of the vertical line x   = λ ^ L 0 , expressed as
A S = x , α λ ^ L 0 x d c 2 ( α ) , 0 α b ,
where α = b. Then d c 2 ( b )   =   λ ^ L 0 . Similar to A R , let h = [100 ×  b], and then j = 0, 1, 2, …, h for 0   α    b. [100 ×  b] represents the largest integer, less than or equivalent to 100 ×  b. Clearly, b = 0.01 ×  h and h + 1 horizontal lines of α = j × 0.01, (j = 0, 1, 2, …, h) cut A S into h quasi-trapezoid blocks. Therefore, the area of the jth block is written as
A S j = x , α λ ^ L 0 x d c 2 ( 0.01 × j ) , 0.01 × j 1 α 0.01 × j ,   j =   0 ,   1 ,   2 , ,   h .
Then, the length l s j of the jth horizontal line is denoted by
l s j = 1 1 z 0.005 m u   d c   λ ^ L 0 ,   j =   0 ,   1 ,   2 , ,   h .
Obviously, as l s 0 = l s 1 and l s h = 0 , then the area of A R is written as
A S = j = 1 h ( 0.01 ) × l s j 1 + l s j 2 = 0.5 × l s 1 + 0.01 × j = 1 h 1 l s j .
According to some studies, the ratio of A S to A R ( A S / A R ) can be applied to a fuzzy decision, as shown below:
A S A R = 0.5 × l s 1 + 0.01 × j = 1 h 1 l s j 0.5 × l r 1 + 0.01 × j = 1 99 l r j .
It is extremely complicated to calculate the block areas of A S and A R using Equations (47) and (51), respectively. As a result, the calculation of the value of A S / A R is relatively complex (in Equation (52)). In order to facilitate industrial applications, the membership functions of g ( x ) and g c ( x ) in Figure 3 both fall within the asymmetric graphical distributions. Based on the method by Chen and Yu [16], this paper employs the bottom length l S = d c R λ ^ L 0 of block A S to correspond to the area for A s , since the function g c ( x ) is of the asymmetric graphical distribution. Additionally, in Figure 3, this paper uses the bottom length l R = 2( d c R d c ) of block A R to correspond to the area for A R , where l R ≥ ( d c R d c L ). Next, l S l R corresponding to A S A R is utilized as a fuzzy evaluation tool, where l S and l R are defined as follows:
l S = d c R λ ^ L 0 = 1 1 z 0.005 / m u   d c   λ ^ L 0 ,
and
l R = 2 ( d c R d c ) =   2 d c 1 1 z 0.005 m u 1 .
According to Chen and Yu [16], if we let η   = l S l R , then process engineers can determine two values, η 1 and η 2 , as the basis for fuzzy decision-making using their past accumulated experiences or long-term production data analysis cases. As noted by Buckley [29], we let 0 < η 1 < η 2 < 0.5. Applying these two numbers should be taken into consideration for fuzzy testing. Therefore, the rules for fuzzy testing are established below [17]:
(1)
If η   <   η 1 then H 0 is rejected and λ L    e is concluded.
(2)
If η 1     η     η 2 , then a decision on “reject/not reject” is not made.
(3)
If η 2   <   η   <  0.5, then H 0 is rejected and λ L  < e is concluded.

4. Practical Example

This section illustrates the fuzzy testing method demonstrated in Section 3 using a numerical example. According to Chen and Yu [16], the value of the lifetime performance index must be greater than 3. Therefore, the null hypothesis is H 0 : λ L   3, whereas the alternative hypothesis is H 1 : λ L   < 3.
Here are 30 sample data points (t1, t2, …, t30) =
(2.31869876,0.450464354,5.959543721,0.19815643,2.36292942,
1.711931047,1.463690428,2.177143634,1.336923286,0.160569092,
0.117929117,2.045192427,5.126371449,0.718544566,0.279190037,
5.722989698,0.527396462,0.760546501,2.090882944,2.353623477,
1.492238468,0.967454182,1.050979265,1.878164671,2.85121394,
1.06125565,8.876851026,1.889987275,0.185337992,1.2043438),
representing the lifetimes T of electronic components generated from an exponential distribution with a mean of 3.
We work with the relative lifetimes X (=T/L), (x1, x2, …, x30) = (t1, t2, …, t30) with L = 1 of electronic components. Additionally, let the given censoring time be x c = P60. Based on Equation (28), we obtain the number of uncensored data, mu = n × k = 18, and then the censoring time x c = P60. = x ( 18 ) = 1.878164671. Using Equation (2), the type-I observation censored data, (y1, y2, …, y30) =
(0.117929117,0.160569092,0.1853379920.19815643,0.279190037,
0.450464354,0.527396462,0.718544566,0.760546501,0.967454182,
1.06125565,1.2043438,1.336923286,1.463690428,1.492238468,
1.711931047,1.878164671,1.878164671,1.878164671,1.878164671,
1.8781646711.8781646711.8781646711.8781646711.878164671
1.8781646711.8781646711.8781646711.8781646711.878164671).
Using Equation (28), the observed value of the estimate is expressed as follows.
λ ^ L 0 =   i = 1 n y i m u = 38.10309 18 =   2.11684
According to Equations (35) and (36), the values of λ L 0 * * and λ R 0 * * can be calculated as follows:
λ L 0 * * = 1 1 + z 0.005 m u   λ ^ L 0   = 1 1 + 2.57583 18   × 2.11684   = 1.31716
and
λ R 0 * * = 1 1 z 0.005 m u   λ ^ L 0   = 1 1 2.57583 18   × 2.11684   = 5.38813
Applying Equations (56) and (57) to Equation (37), we define the membership function of the fuzzy number λ ˜ L 0 * * as
d ( x ) = 0 ,   i f   x < 1.31716 2 1 F Z 18 2.11684 x 1 ,   i f   1.31716 x <   2.11684 1   ,   i f   x = 2.11684   2 1 F Z 18 1 2.11684 x ,   i f   2.11684 < x 5.38813   0   ,   i f   5.38813 < x
The two hypotheses are denoted by H 0 : λ L 3 versus H 1 : λ L < 3, given the significant level α = 0.05. Using Equation (27), we can obtain the critical value as follows:
d c   =   e z α e m u     e   = 3   = 1.83681 .
According to Equations (41) and (42), we find the values of d c L and d c R below:
d c L = 1 1 + 2.57583 18   ×   1.83681   = 1.14291
and
d c R = 1 1 2.57583 18   ×   1.83681   = 4.67535 .
Applying Equations (60) and (61) to Equation (43), we define the membership function of the fuzzy number d ˜ c as
d c ( x ) = 0 ,   i f   x < 1.14291 2 1 F Z 18 1.83681 x 1 ,   i f   1.14291 x <   1.83681 1 ,   i f   x = 1.83681 2 1 F Z 18 1 1.83681 x ,   i f   1.83681 < x 4.67535   0 ,   i f   4.67535 < x
Subsequently, the curves of d ( x ) and d c ( x ) are depicted in Figure 4.
Using Equations (53) and (54), we obtain the values of l S and l R :
l S = d c R λ ^ L 0 = 1 1 z 0.005 m u   d c   λ ^ L 0 = 1 1 2.57583 18   × 1.83681     2.1168 4 = 2.55851
and
l R =   2 ( d c R d c ) =   2 d c 1 1 z 0.005 m u 1 = 2 × 1.83681 × 1 1 2.57583 18 1    =   5.67708
Therefore, η  =   l S / l R   = 2.55851/5.67708 = 0.45067. Based on statistical inference, λ ^ L 0 =   2.11684 is in the region of “do not reject H 0 ” ( λ ^ L 0 <   1.83681 ; using Equation (59)), the conclusion is “do not reject H 0 , and we conclude that λ L   3”. According to the rules for fuzzy testing, the value of η 1 is 0.2, and the value of η 2 is 0.4 [16,17]. Therefore, we get η 2   <  η < 0.5, indicating that H 0 needs to be rejected and λ L < 3 is concluded. Nevertheless, λ ^ L 0   = 2.11684 is much lower than λ L = 3. Accordingly, based on the fuzzy testing method proposed in this study, H 0 needs to be rejected and λ L < 3 is concluded. According to the above fuzzy testing rules, this result seems reasonable in view of practicability.

5. Conclusions

This paper places great emphasis on developing a model for evaluating a product’s lifetime performance, utilizing the assessment process to enhance its lifetime performance, thereby increasing its value. It can also help attain the goals of green production—saving energy and diminishing waste. In this paper, the electronic component, lifetime (T), was adopted from the exponential distribution that is often used in general studies. In addition, we also employed the electronic component lifetime and the relative electronic component lifetime X (=T/L) of the minimum required unit time L as well as constructed a relative lifetime performance index. The larger the index value, the better the electronic component performance. When collecting electronic component lifetime data, we often fail to collect complete data for several reasons, such as time, cost, and manpower. Consequently, we employed type-I censoring and derived the best estimator (UMVUE) for the lifetime performance index. In the meantime, a fuzzy testing model was proposed and built on the derived confidence intervals. Moreover, this fuzzy testing model prevented misjudgments caused by sampling errors, making it highly suitable for the service and manufacturing industries, which place utmost value on timeliness. Additionally, this model is also applicable to experiments that require destructive testing, higher costs, or longer durations. Since passive components are a cornerstone of the electronics industry and play a vital role in driving peripheral equipment and industries, passive components were utilized as the basis to discover and explain the proposed model in this paper. At the same time, a numerical example was also put forward to illustrate the proposed model.

Author Contributions

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

Funding

This research received no external funding.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Suppipat, S.; Hu, A.H. Achieving sustainable industrial ecosystems by design: A study of the ICT and electronics industry in taiwan. J. Clean. Prod. 2022, 369, 133393. [Google Scholar] [CrossRef]
  2. Wang, C.T.; Chiu, C.S. Competitive strategies for Taiwan’s semiconductor industry in a new world economy. Technol. Soc. 2014, 36, 60–73. [Google Scholar] [CrossRef]
  3. Amindoust, A.; Ahmed, S.; Saghafinia, A.; Bahreininejad, A. Sustainable supplier selection: A ranking model based on fuzzy inference system. Appl. Soft Comput. 2012, 12, 1668–1677. [Google Scholar] [CrossRef]
  4. Kao, F.C.; Huang, S.C.; Lo, H.W. A rough-fermatean DEMATEL approach for sustainable development evaluation for the manufacturing industry. Int. J. Fuzzy Syst. 2022, 24, 3244–3264. [Google Scholar] [CrossRef]
  5. Pearn, W.L.; Cheng, Y.C. Measuring production yield for processes with multiple characteristics. Int. J. Prod. Res. 2010, 48, 4519–4536. [Google Scholar] [CrossRef]
  6. Wu, M.F.; Chen, H.Y.; Chang, T.C.; Wu, C.F. Quality evaluation of internal cylindrical grinding process with multiple quality characteristics for gear products. Int. J. Prod. Res. 2019, 57, 6687–6701. [Google Scholar] [CrossRef]
  7. Yang, C.M. An improved multiple quality characteristic analysis chart for simultaneous monitoring of process mean and variance of steering knuckle pin for green manufacturing. Qual. Eng. 2021, 33, 383–394. [Google Scholar] [CrossRef]
  8. Menon, R.R.; Ravi, V. Using ANP and QFD methodologies to analyze eco-efficiency requirements in an electronic supply chain. Clean. Eng. Technol. 2021, 5, 100350. [Google Scholar] [CrossRef]
  9. Simmou, W.; Govindan, K.; Sameer, I.; Hussainey, K.; Simmou, S. Doing good to be green and live clean!—Linking corporate social responsibility strategy, green innovation, and environmental performance: Evidence from maldivian and moroccan small and medium-sized enterprises. J. Clean. Prod. 2023, 384, 135265. [Google Scholar] [CrossRef]
  10. Yang, Y.; Jiang, Y. Buyer-supplier CSR alignment and firm performance: A contingency theory perspective. J. Bus. Res. 2023, 154, 113340. [Google Scholar] [CrossRef]
  11. Anderson, D.R.; Sweeney, D.J.; Williams, T.A. Statistics for Business and Economics; West Publ. Co.: Saint Paul, MN, USA, 1990. [Google Scholar]
  12. Keller, G.; Warrack, B.; Bartel, H. Statistics for Management and Economics; Duxbury Press: Belmont, CA, USA, 1994. [Google Scholar]
  13. Lee, H.M.; Wu, J.W.; Lei, C.L. Assessing the lifetime performance index of exponential products with step-stress accelerated life-testing data. IEEE Trans. Reliab. 2013, 62, 296–304. [Google Scholar] [CrossRef]
  14. Lee, W.C.; Wu, J.W.; Hong, C.W.; Hong, S.F. Evaluating the lifetime performance index based on the bayesian estimation for the rayleigh lifetime products with the upper record values. J. Appl. Math. 2013, 2013, 547209. [Google Scholar] [CrossRef]
  15. Hsu, B.M.; Shu, M.H.; Chen, B.S. Evaluating lifetime performance for the pareto model with censored and imprecise information. J. Stat. Comput. Simul. 2011, 81, 1817–1833. [Google Scholar] [CrossRef]
  16. Chen, K.S.; Yu, C.M. Lifetime performance evaluation and analysis model of passive component capacitor products. Ann. Oper. Res. 2022, 311, 51–64. [Google Scholar] [CrossRef]
  17. Chiou, K.C. Building up of Fuzzy Evaluation Model of Life Performance Based on Type-II Censored Data. Mathematics 2023, 11, 3686. [Google Scholar] [CrossRef]
  18. Chen, K.S.; Chiou, K.C.; Yu, C.M. Lifetime performance index of electronic products. Microelectron. Reliab. 2020, 113, 113941. [Google Scholar] [CrossRef]
  19. Miller, R.G. Survival Analysis; John Wiley and Sons: New York, NY, USA, 1981. [Google Scholar]
  20. Khorram, E.; Farahani, Z.S.M. Statistical inference of weighted Exponential lifetimes under progressive type-II censoring scheme. Qual. Technol. Quant. Manag. 2014, 11, 433–451. [Google Scholar] [CrossRef]
  21. Wu, S.F.; Hsieh, Y.T. The assessment on the lifetime performance index of products with Gompertz distribution based on the progressive type I interval censored sample. J. Comput. Appl. Math. 2019, 351, 66–76. [Google Scholar] [CrossRef]
  22. Wu, S.F.; Wu, Y.C.; Wu, C.H.; Chang, W.T. Experimental design for the lifetime performance index of weibull products based on the progressive type i interval censored sample. Symmetry 2021, 13, 1691. [Google Scholar] [CrossRef]
  23. Wang, T.C.; Wu, C.W.; Shu, M.H. A variables-type multiple-dependent-state sampling plan based on the lifetime performance index under a weibull distribution. Ann. Oper. Res. 2022, 311, 381–399. [Google Scholar] [CrossRef]
  24. Chen, H.Y.; Lin, K.P. Fuzzy supplier selection model based on lifetime performance index. Expert Syst. Appl. 2022, 208, 118135. [Google Scholar] [CrossRef]
  25. Ahmed, N.; Ali, S.; Shah, I. Control charts for monitoring mean of generalized exponential distribution with type-I censoring. Qual. Reliab. Eng. Int. 2022, 38, 592–614. [Google Scholar] [CrossRef]
  26. Raza, S.M.M.; Riaz, M.; Ali, S. EWMA Control Chart for Poisson-Exponential Lifetime Distribution under Type i Censoring. Qual. Reliab. Eng. Int. 2016, 32, 995–10051. [Google Scholar] [CrossRef]
  27. Wu, S.F.; Jheng, J.J.; Chang, W.T. Sampling design for the lifetime performance index of exponential lifetime distribution under progressive type I interval censoring. Commun. Stat.-Theory Methods 2023, 52, 2766–2782. [Google Scholar] [CrossRef]
  28. Casella, G.; Berger, R.L. Statistical Inference, 2nd ed.; Duxbury Press: Pacific Grove, CA, USA, 2002. [Google Scholar]
  29. Buckley, J.J. Fuzzy statistics: Hypothesis testing. Soft Comput. 2005, 9, 512–518. [Google Scholar] [CrossRef]
Figure 1. Curves for P   k , k = 10(10)100, λ L = 1(1)5, n = 100, α = 0.05.
Figure 1. Curves for P   k , k = 10(10)100, λ L = 1(1)5, n = 100, α = 0.05.
Mathematics 12 01935 g001
Figure 2. Curves for n = 10(10)100, λ L = 1(1)5, P k   (k = 60), α = 0.05.
Figure 2. Curves for n = 10(10)100, λ L = 1(1)5, P k   (k = 60), α = 0.05.
Mathematics 12 01935 g002
Figure 3. The membership functions of g ( x ) and g c ( x ) .
Figure 3. The membership functions of g ( x ) and g c ( x ) .
Mathematics 12 01935 g003
Figure 4. Membership functions d ( x ) and d c ( x ) applied in the numerical example.
Figure 4. Membership functions d ( x ) and d c ( x ) applied in the numerical example.
Mathematics 12 01935 g004
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

Chiou, K.-C.; Huang, T.-H.; Chen, K.-S.; Yu, C.-M. Fuzzy Evaluation Model for Lifetime Performance Using Type-I Censoring Data. Mathematics 2024, 12, 1935. https://doi.org/10.3390/math12131935

AMA Style

Chiou K-C, Huang T-H, Chen K-S, Yu C-M. Fuzzy Evaluation Model for Lifetime Performance Using Type-I Censoring Data. Mathematics. 2024; 12(13):1935. https://doi.org/10.3390/math12131935

Chicago/Turabian Style

Chiou, Kuo-Ching, Tsun-Hung Huang, Kuen-Suan Chen, and Chun-Min Yu. 2024. "Fuzzy Evaluation Model for Lifetime Performance Using Type-I Censoring Data" Mathematics 12, no. 13: 1935. https://doi.org/10.3390/math12131935

APA Style

Chiou, K. -C., Huang, T. -H., Chen, K. -S., & Yu, C. -M. (2024). Fuzzy Evaluation Model for Lifetime Performance Using Type-I Censoring Data. Mathematics, 12(13), 1935. https://doi.org/10.3390/math12131935

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