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
[
11,
12,
13]. Thus, the probability density function (
p.d.f.) is
for
t > 0 and the cumulative function (
c.d.f.) of
T is
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
, expressed as
X~
. The second parameter
=
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
. Its probability density function (
p.d.f.) is
for x > 0 and the cumulative distribution function (
c.d.f.) is
, for x > 0. According to Chen et al. [
18], the relative failure rate is denoted by
, and the reliability function of the relative lifetime X is
for x > 0.
Based on Chen et al. [
18], product reliability (
) represents the probability showing that the relative lifetime
X is greater than 1, such that
, indicating its direct mathematical relationship with the relative lifetime performance index
. Obviously, when the value of
(=
) is larger than that of
L, then the relative lifetime performance index is larger. Additionally,
increases as the value of the indicator
rises. Therefore, the relative lifetime performance index
can reveal the lifetime performance of the electronic component in a reasonable manner.
If the relative lifetime performance index
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
. Concerning the censored observation data (
), 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
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 (
)% confidence intervals using type-I right censoring sample data. Subsequently, in
Section 3, these 100 (
)% 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
, where
denotes the relative lifetime performance index for an electronic component as follows:
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
is called type-I right censoring, then
can be rewritten as follows:
where
i = 1, 2, …,
n. Let
be an indicator function, expressed as follows:
Next, let
mu be the number of uncensored data, defined as follows:
In fact,
mu follows a Binomial distribution, denoted by
mu~
B(
n,
p), where
p =
P(
X ≤
xc). According to the probability density function and the cumulative distribution function of X, we can find the distribution of estimate
using type-I right censoring. Let the likelihood function of pair
be as follows:
where
y1,
y2, …,
yn are the type-I censored observation data for
Y1,
Y2, …,
Yn. Let
. Then, the likelihood function of type-I right censoring is expressed as follows:
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.
Setting Equation (7) to 0, we can find the estimator of
as follows:
Using Equation (8), the Fisher information of
is defined as follows:
Based on the large-sample theorem, the asymptotic variance of the estimate
can be expressed as:
Then, the asymptotic distribution of estimate
is presented as follows:
When
xc → ∞ and
mu ≈
n, the estimate
is defined as
Given that the censoring time is xc and (Y1, Y2, …, Yn) is a random sample, the number of uncensored data in Equation (4) is unknown, and it follows a Binomial distribution (n, p) with p = P(X ≤ xc) = . Next, we can get the expected value E(mu) = np = .
As noted above, this study needs to estimate the lifetime performance index using the sample data since its real value is unknown. Let the electronic component lifetime T have an exponential distribution with mean . Since the relative lifetime X = T/L, X is distributed as an exponential distribution with mean . 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
as follows:
Let the lifetime performance index
be the function
h(
),
h(
) =
. Using Equation (12) and the delta method (Casella and Berger [
28]), we can get the asymptotic normal distribution for the estimator of
=
h(
) as follows:
where the asymptotic variance of
(=
h(
)) is
=
. Thus, we have the following equation:
Given
the estimator of
is unbiased and consistent with
, displayed as follows:
According to the Lehmann–Scheffé theorem,
represents the uniformly minimum variance unbiased estimator (
UMVUE) of
. Using Equation (16), we can get the asymptotic standard normal distribution. The random variable
Z, as noted above, is displayed in the following equation:
Thus, the (1 − α) × 100% confidence interval of
is
, where the functions
and
are written as follows:
and
where
represents the upper
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
be
Then the expected value of
is written as follows:
The electronic component lifetime
(T) has an exponential distribution with mean
, while the relative lifetime
X (=
T/
L) has an exponential distribution with mean
. Let
be the k
th percentile of the relative lifetime
X,
k = 1, 2, …,99, such that
P(
X ≤
) =
k%. Then, the
kth percentile (
) is expressed as follows:
Given the confidence interval level (1 − α) × 100% = 95%, sample size
n = 100, relative lifetime performance index
= 1(1)5 and censoring time
xc =
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 =
, then the value of
E(
mu) is received as follows:
When the number of samples,
n = 100, is fixed, based on Equation (22), we can get the distribution for
, as shown in
Figure 1. As
is fixed, the larger the censoring time
xc (=
) is, the smaller the average confidence interval length of
is, which means the estimate of the index
is better. That is, when the number of uncensored data
mu is higher, the estimated indicator
is better.
When the
k value of quantile
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
gets smaller, indicating the estimate of the index
is better, as shown in
Figure 2.
Based on Equations (23) and (25), the expected value
of the confidence interval length for the relative lifetime performance index
is related to the number of uncensored observations
mu. When
mu is larger, then
is smaller, that is, the estimate of index
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
. 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
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 : e versus alternative hypothesis : e—to evaluate if the lifetime performance index has reached the desired level. Before performing this method, we examine the testing rules of statistical hypotheses as follows:
- (1)
Given that , do not reject and conclude that e.
- (2)
Given that , reject and conclude that e.
In the rules, the critical value
is obtained using Equation (16). Then, the significant level α is received and expressed as
Next, we obtain the critical value
as follows:
Using Equations (2), (9), and (14), the observed value of estimate
is calculated using the following equation:
Based on Buckley’s approach [
29], Equations (20) and (21), the
of the triangular fuzzy number
is denoted by the following equation [
16,
17]:
where
and
Clearly, the value of
is not equivalent to that of
given that
α < 1. As
1 and −
=
, we obtain the result
=
, where
=
. Therefore, let
, such that the
of the triangular fuzzy number
is defined as
where
and
Obviously, the value of
is equal to that of
with
1, and the newly transformed triangular fuzzy number is represented as
, where
,
and
Furthermore, the membership function of the fuzzy number
is denoted by
where
is a cumulative distribution function of the standard normal distribution. Similar to the fuzzy number
, the
of triangular-shaped fuzzy critical value number
is expressed as
where
and
Clearly, the value of
is equivalent to that of
given that
Then, a newly transformed triangular fuzzy number is denoted by
, where
,
and
Furthermore, the membership function of the fuzzy number
is represented as
where
is a cumulative distribution function of the standard normal distribution. Subsequently, the curves of
and
are depicted below:
Based on Buckley [
29] and Chiou [
17], we let
be the area in the curved graph of
and be defined as
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 Riemann sum is performed in the study to find the area of the block. The procedures are described as follows:
(1) We divide
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
is cut into 100 trapezoidal blocks by these 101 horizontal lines. Therefore, the area of the
jth block is denoted by
Next, the length of the
jth horizontal line is depicted as
Obviously, when
and
0, then the area is defined as
Moreover, let
be the area in the curved graph of
but to the right of the vertical line
, expressed as
where
α = b. Then
. Similar to
, 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
into
h quasi-trapezoid blocks. Therefore, the area of the
jth block is written as
Then, the length
of the
jth horizontal line is denoted by
Obviously, as
and
, then the area of
is written as
According to some studies, the ratio of
to
(
/
) can be applied to a fuzzy decision, as shown below:
It is extremely complicated to calculate the block areas of
and
using Equations (47) and (51), respectively. As a result, the calculation of the value of
/
is relatively complex (in Equation (52)). In order to facilitate industrial applications, the membership functions of
and
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
=
−
of block
to correspond to the area for
, since the function
is of the asymmetric graphical distribution. Additionally, in
Figure 3, this paper uses the bottom length
= 2(
) of block
to correspond to the area for
, where
≥ (
−
). Next,
corresponding to
is utilized as a fuzzy evaluation tool, where
and
are defined as follows:
and
According to Chen and Yu [
16], if we let
, then process engineers can determine two values,
and
, 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 <
<
< 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 then is rejected and e is concluded.
- (2)
If , then a decision on “reject/not reject” is not made.
- (3)
If 0.5, then is rejected and < 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
:
3, whereas the alternative hypothesis is
:
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
=
P60. Based on Equation (28), we obtain the number of uncensored data,
mu =
n ×
k = 18, and then the censoring time
=
P60. =
= 1.878164671. Using Equation (2), the type-I observation censored data, (
y1,
y2, …,
y30) =
(0.117929117, | 0.160569092, | 0.185337992 | 0.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.878164671 | 1.878164671 | 1.878164671 | 1.878164671 | 1.878164671 |
1.878164671 | 1.878164671 | 1.878164671 | 1.878164671 | 1.878164671). |
Using Equation (28), the observed value of the estimate is expressed as follows.
According to Equations (35) and (36), the values of
and
can be calculated as follows:
and
Applying Equations (56) and (57) to Equation (37), we define the membership function of the fuzzy number
as
The two hypotheses are denoted by
:
3 versus
:
3, given the significant level
α = 0.05. Using Equation (27), we can obtain the critical value as follows:
According to Equations (41) and (42), we find the values of
and
below:
and
Applying Equations (60) and (61) to Equation (43), we define the membership function of the fuzzy number
as
Subsequently, the curves of
and
are depicted in
Figure 4.
Using Equations (53) and (54), we obtain the values of
and
:
and
Therefore,
η /
2.55851/5.67708
0.45067. Based on statistical inference,
2.11684 is in the region of “do not reject
” (
; using Equation (59)), the conclusion is “do not reject
, and we conclude that
3”. According to the rules for fuzzy testing, the value of
is 0.2, and the value of
is 0.4 [
16,
17]. Therefore, we get
η 0.5, indicating that
needs to be rejected and
< 3 is concluded. Nevertheless,
= 2.11684 is much lower than
= 3. Accordingly, based on the fuzzy testing method proposed in this study,
needs to be rejected and
< 3 is concluded. According to the above fuzzy testing rules, this result seems reasonable in view of practicability.