Next Article in Journal
Multiple Control Policy in Unreliable Two-Phase Bulk Queueing System with Active Bernoulli Feedback and Vacation
Previous Article in Journal
Path Planning of Logistic Robot Using Method of Vector Marks Tree Generation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Residual Control Chart Based on a Convolutional Neural Network and Support Vector Regression for Type-I Censored Data with the Weibull Model

Department of Information Management, Chaoyang University of Technology, Taichung 413310, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(1), 74; https://doi.org/10.3390/math12010074
Submission received: 9 October 2023 / Revised: 20 December 2023 / Accepted: 21 December 2023 / Published: 25 December 2023
(This article belongs to the Section Probability and Statistics)

Abstract

:
Control charts with conditional expected value (CEV) can be used with novel statistical techniques to monitor the means of moderately and lowly censored data. In recent years, machine learning and deep learning have been successfully combined with quality technology to solve many process control problems. This paper proposes a residual control chart combining a convolutional neural network (CNN) and support vector regression (SVR) for type-I censored data with the Weibull model. The CEV and exponentially weighted moving average (EWMA) statistics are used to generate training data for the CNN and SVR. The average run length shows that the proposed chart approach outperforms the traditional EWMA CEV chart approach in various shift sizes and censored rates. The proposed chart approach is suitable to be used in detecting small shift size for highly censored data. An illustrative example presents the application of the proposed method in an electronics industry.

1. Introduction

Product service life has always been a quality highly valued by customers. To ensure customers enjoy good reliability, many manufacturers conduct life tests on their products before they leave the factory. As production technology has advanced, products have achieved longer service life. Practitioners thus need to spend increasing time and cost to collect product life data. To save time and cost in testing, practitioners may set a test termination time at which to halt the test, whether or not there has been a failure. As some testing units have not failed by the termination time, practitioners can only obtain incomplete lifetime records. These incomplete lifetime records are termed type-I right-censored data.
Control charts are used to monitor and track process data in order to quickly identify potential variations, which assists practitioners in understanding whether a production process is in an in-control state. As traditional Shewhart-type and exponentially weighted moving average (EWMA)-type control charts cannot be used with right-censored data, novel control charts have been proposed for this purpose [1,2,3]. Steiner and Mackay made multiple proposals, including integrating the X ¯ control chart with a conditional expected value (CEV) to monitor highly right-censored data and setting a lower control limit (LCL) to monitor lifetime decreases under normally distributed processes, developing Shewhart-type CEV control charts for non-normally distributed right-censored data, and applying an EWMA CEV control chart to increase detection ability for highly right-censored data [4,5,6]. Lee applied an economic design approach to determine the design parameters of the CEV X ¯ control chart [7].
To detect decreases and increases in average lifetimes under Weibull distributions, Zhang and Chen developed two single-sided EWMA CEV control charts [8]. Subsequently, many researchers proposed EWMA-type CEV and conditional median control charts to control non-normally distributed censored data, achieving superior statistical performance [9,10,11]. Double EWMA (DEWMA) control charts for censored data were developed and it was found that they detect shifts of the scale parameter under gamma distributions more efficiently than EWMA charts [12,13]. Many studies for further applications of DEWMA CEV control charts have also been published [12,14,15]. The EWMA CEV chart has also been modified to monitor various censored data. For example, EWMA CEV control charts were extended to detect average lifetimes for multiple censored data and an EWMA chart was implemented with a modified CEV to monitor Weibull window-censored data [16,17]. In the above-mentioned literature, most of the research focuses on increasing the detection efficiency of moderate and low censored rates and moderate shift sizes [10,12,13,14,15]. Few of the literature have proposed to improve the detection efficiency of small average lifetime shifts and high censored rates using larger sample sizes [3]. Since reliability life testing is a destructive test method, increasing the sample size will lead to increased testing costs. If control charts can improve the detection ability of high censored rates and small shift sizes without increasing sample size, this will help practitioners reduce quality control costs.
With the recent rise of artificial intelligence, many studies have introduced machine-learning or deep-learning methods to solve quality problems. Many researchers have all proposed control charts based on machine or deep learning to improve detection ability [18,19,20,21,22,23,24,25]. Deep-learning methods have achieved excellent performance in failure detection in high-dimensional or complex processes [3,26,27,28]. Moreover, some researchers recognized abnormal control chart patterns using machine- or deep-learning methods [29,30,31,32,33,34]. However, to the authors’ knowledge, there have been no relevant studies that use deep-learning or machine-learning methods to monitor right-censored data.
Convolutional neural networks (CNN) and support vector regression (SVR) are types of deep-learning and machine-learning methods, respectively. The literature indicates that both methods have successfully improved the performance of control charts and enhanced failure detection [28,35,36]. To the authors’ knowledge, however, there has been no study that combines these two methods for process control.
Most existing works have explored various statistical methods to obtain good data features to improve the accuracy of SVR. Finding feature-extraction methods is typically time-consuming, and CNNs have a good ability to extract data features. Thus, combining a CNN and SVR may save time in finding suitable data feature extraction methods and improve prediction accuracy. In addition, the Weibull distribution can present skewed and approximately symmetrical distributions using differently shaped parameter values and is widely used in lifetime analysis. Thus, this study will construct a control chart with a combination of CNN and SVR to improve the detection ability in high censored rates and small shift sizes for Weibull type-I right-censored data.
The remainder of this study is organized as follows. Section 2 presents the basic CNN and SVR methodologies. In Section 3, the proposed method is described, and in Section 4, its performance is investigated. Section 5 provides a real-world case to illustrate the implementation of the proposed chart. Finally, Section 6 summarizes key findings and outlines future development directions.

2. Combined CNN and SVR

2.1. CNN

CNNs are feed-forward neural networks that use multiple convolutional and pooling layers to extract higher-level spatial information from input data. The approach has broad applications, particularly in speech recognition and computer vision [37]. A CNN model generally consists of an input layer, an output layer, and multiple hidden layers, which include convolution layers, pooling layers, and fully connected layers. The sample points plotted on a control chart are similar to one-dimensional sequence data; a control chart is suitable for training the sample points of the control map with a one-dimensional CNN model [28].
In a convolution layer, the output feature map is generated through the convolution of multiple input feature maps and multiple convolutional kernels. The output of a convolution layer is expressed by the following:
ζ j l = f R L b j l + i = 1 D l 1 ζ i l 1 × ω i j l 1 ,
where f R L is a rectified linear unit (ReLU) activation function, l is the l -th layer of the CNN network, ω i j l 1 is the one-dimensional filter kernel from the i -th neuron at layer l 1 to the j -th neuron at layer l , b j l is the bias of the j -th neuron at layer l , and D l 1 is the number of the input at layer l 1 .
The pooling layer serves to decrease the feature map’s dimensions and streamline computations. The max pooling method is a commonly implemented pooling operation in this layer. The pooling output at layer l is given by following:
M p j l = m a x ζ j l .
A fully connected layer is the last layer of a CNN and follows several convolution and pooling layers. The operation of this layer is similar to a multilayer neural network and the neurons in this layer connect all the activations in the previous layer. The output of a fully connected layer can consider a categorical variable for multi-class classification problems or a regression response variable for relevance prediction problems.

2.2. SVR

SVR can effectively address regression problems involving high-dimensional features by finding a function that approximates the relationship between the input vector and the output values. Let u 1 , y 1 , , u m , y m be the data set for regression analysis. SVR is then expressed as follows:
min 1 2 w 2 + c s i = 1 m ξ i + ξ ^ i s . t .   y i + w T u i e ε + ξ i w T u i + e y i ε + ξ i ξ i ,   ξ ^ i 0 ,
where is the Euclidian norm, w is the weight vector, e is the bias term, ε is the width of the insensitive zone, ξ i and ξ ^ i are slack variables for relaxing the training errors (where ξ i , ξ ^ i 0 ), and c s is a positive constant that determines the trade-off between minimizing the error and minimizing the slack variables [9,38].
Next, the dual model of Equation (1) can be established using the Lagrangian method as follows:
max i = 1 m y i ν i ν ^ i ε i = 1 m ν i + ν ^ i 1 2 i = 1 m j = 1 m ν i ν ^ i ν j ν ^ j u i T u j s . t .   i = 1 m ν ^ i ν i = 0 ν ^ i ,   ν i 0 , c s ,
where ν ^ and ν i are the Lagrange multipliers and must be non-negative and real. By solving the above, the SVR predictor becomes
f s v r u , ν ^ i ,   ν i = i = 1 m ν ^ i ν i κ u i T , u j + e ,
where κ u i T , u j is a kernel function. The radial basis function (RBF) is a kernel function commonly used in SVR predictors [39].

2.3. Combination of CNN and SVR

In this section, a CNNSVR network model was created that combines CNN and SVR for process control. The next section presents the setup and procedure for a new control chart with the proposed CNNSVR network. In the CNNSVR network, a CNN is trained to generate the feature maps and output the activation values of the last layer, and these activation values are used as the input for the SVR to generate the predicted values, as shown in Figure 1. Because the sample points are one-dimensional sequence data, the CNNSVR network adopts a one-dimensional sequence input layer. A combinational layer is defined as one convolutional layer followed by one pooling layer. Additionally, the ReLU activation function is used in each combinational layer to calculate the output vector. The CNNSVR network consists of N combinational layers for the activation value. This activation value is inputted to the SVR to obtain the predicted value.
In Figure 1, green, brown, and black dashed lines represent the convolution (Co.layer), pooling (Po.layer), and fully connected (FC.layer) layers, respectively. The dark green and orange squares are the kernels of the Co.layer and Po.layer, respectively. Each yellow rectangle and each red rectangle represent a feature map (F.map) of the Co.layer and Po.layer, respectively, the height of the rectangle means the size of the feature map, and Ki is the number of feature maps of i-th layer. The number and size of the convolution kernel (or filter) affects the number and size of the feature map of the convolution layer. The kernel size of the pooling layer determines the size of the feature map of the pooling layer. When using the CNNSVR network, the number of combinational layers and convolution kernels and the size of the convolution and pooling kernels must be determined first to achieve good training and testing results. This paper refers to the literature and uses the trail-error method to optimize these parameter values [18,23].

3. CNNSVR-Based Residual Chart

Deep-learning-based residual control charts have been presented to increase the performance of process control [18,20]. The current paper extends their method to propose a CNNSVR-based residual chart for Weibull type-I right-censored data. The sample point of a control chart in an in-control state is used to train the CNNSVR network to predict the sample point (or statistic) of the next period. In other words, the sample point at period t − 1 is used to predict the sample point at period t. If the predicting sample point is similar to the actual sample point, then this indicates the process is in-control; otherwise, the process mean has shifted. A threshold value is determined by the residual values of the trained CNNSVR network, and it plays a role as the LCL of the traditional control chart to indicate the process state. If the error value is below the threshold value, it indicates a significant difference between the predicted sample point and the actual sample point; therefore, the process is considered to be in an out-of-control state.
The MATLAB R2023a software provides various deep-learning, statistics, and machine-learning toolboxes, which can be easily employed for the simulation processes to investigate the control chart performance.

3.1. EWMA CEV Statistic for a Weibull Distribution

Let S = s 1 , s 2 , , s k be a Weibull random variable and a 0 and b 0 be the scale and shape parameters of an in-control process, respectively. The in-control mean ( M 0 ) and variance ( V 0 ) are as follows:
M 0 = a 0 Γ 1 + 1 / b 0   and V 0 = a 0 2 Γ 1 + 2 / b 0 Γ 1 + 1 / b 0 2 .
The censoring rate for Weibull lifetimes can be expressed as P c = e x p c T / a b , where c T is the censoring time. The CEV of the Weibull distribution is as follows:
C e v = E S S c T = a 0 Υ D w , 1 + 1 / b 0 e x p D w ,
where D w = c T / a 0 b 0 and Υ θ 1 , θ 2 = 0 θ 1 r θ 2 1 e x p r d r [14].
Practitioners take n samples and measure their lifetime values using the reliability lifetime testing method. Let s i be the lifetime of the i-th testing sample. Then, the CEV value x i of the i-th testing sample is as follows:
x i = s i f o r   s i   c T C e v f o r   s i > c T .
The sample mean of size n is X ¯ = i = 1 n x i / n . In practice, n = 5 is always used for sampling and plotting control charts. Let λ be the smoothing parameter of the EWMA chart. Zhang and Chen [8] showed that the EWMA statistic at period t for monitoring the mean decrease is as follows:
E t = m i n M 0 , λ X t ¯ + 1 λ E t 1
For an EWMA CEV chart used to monitor a process, an LCL is established to detect the mean disease, because practitioners always focus on the detection of average lifetime reduction. The most commonly used criterion to assess control charts’ performance is average run length (ARL) [8,14,16,18,29]. A larger ARL value in an in-control process means a lower false alarm, whereas a smaller ARL value in an out-of-control state represents the faster detection of shifts.

3.2. Setup of the Proposed Control Chart

The main requirement to set up the control chart is to find a reasonable threshold value for the CNNSVR-based residual chart in order to maintain its false alarm rate. To generate large numbers of the EWMA statistic E t for CNNSVR network training, the following Monte Carlo simulation procedure is applied:
(1)
Generate the lifetimes s i of n samples from a Weibull distribution with parameters a 0 and b 0 .
(2)
If s i > c T , then x i = C E V ; otherwise, x i = s i .
(3)
Calculate the sample mean X ¯ .
(4)
Calculate the EWMA statistic E t using Equation (2).
(5)
Repeat steps (1) to (4) m times to obtain m EWMA statistics E t .
The CNNSVR network is similar to a univariate autoregressive model of order 1 constructed by a neural network, with a record in the CNNSVR training data consisting of an input value and an output value. Let E 1 , E 2 , , E m 1 , E m be the training data set for the CNNSVR network, where E t 1 , E t represents a record in the training data, and E t 1 and E t are the input value and output value, respectively, of the CNNSVR network for this record. The input and output vectors have dimensions 1 × m 1 . Increasing m can improve the CNNSVR network’s prediction ability, but this requires longer training time. This study uses m =10,000.
After training, the CNNSVR network outputs the predictive values. The residual value can be obtained using the actual value (the output value of the training data set) minus the predicted value. Next, a reasonable threshold value η is determined from these residual values. The threshold value η is similar to the LCL of the control charts and affects the ARL value. With consideration for the practitioners’ interest in the state of the reduced product lifetime, only the training residual values E t < E ^ t are used to determine the threshold value η . Since η is related to the false alarm rate, the practitioner can choose a specific in-control ARL value based on their experience to determine the optimal threshold value η with a simulation method. The determination procedure for η is shown in Figure 2.
This study used the ‘wblrnd’ function to generate complete Weibull lifetime data for training, and applied the ‘sequenceInputLayer’, ‘convolution1dLayer’, ‘maxPooling1dLayer’, ‘reluLayer’, ‘fullyConnectedLayer’, ‘regressionLayer’, ‘trainNetwork’, ‘activations’, ‘predict’, and ‘fitrsvm’ functions to construct and train the CNNSVR network.

3.3. Performance Simulation for the Proposed Control Chart

In the control phase, the sample statistic at the (t − 1)-st period E t 1 is inputted into a trained CNNSVR network to predict the sample statistic E ^ t . At the t-th period, the practitioner takes n samples and computes the sample statistic E t . If the prediction error E t E ^ t is less than the threshold η , then the process is identified as being in an out-of-control state; otherwise, the process is in an in-control state.
The occurrence of an assignable cause results in the reduction of the process mean as being in an out-of-control state. Let M 1 = δ × M 0 be the mean of an out-of-control state where the process variance is unchanged and δ is a mean shift size that can be estimated by δ = M 1 / M 0 . The Weibull parameters in an out-of-control state can be obtained by solving the following system of simultaneous equations:
M 1 = a 1 Γ 1 + 1 / b 1 V 0 = a 1 2 Γ 1 + 2 / b 1 Γ 1 + 1 / b 1 2 ,
where a 1 and b 1 are the Weibull scale and shape parameters in an out-of-control state, respectively. Out-of-control lifetime data are generated using the Weibull distribution with parameters a 1 and b 1 .
Control charts have different detection efficiencies at different shift sizes. Therefore, the out-of-control ARL value must be evaluated for various possible shift sizes before this control chart is implemented in order to monitor the process. An out-of-control ARL simulation procedure of the CNNSVR-based residual control chart is proposed, as shown in Figure 3. The same MATLAB functions as in Section 3.2 are applied.

4. Performance Comparison

This section compares the ARL performance of the CNNSVR-based residual chart and the EWMA CEV chart. The parameter choice for the Weibull distribution considers the cases of b 0 = 0.8 , 1 , 2.5 ,   and   4 and a 0 = 1 for comparison. Table 1 shows the skewness values (Sk) of the chosen parameter values. The smaller shape parameter values resulted in a greater skewness of the Weibull distribution. For b 0 = 4 , the Weibull distribution approximates a symmetric probability distribution. In the comparison, the sample size n is fixed at 5, the smoothing parameter λ of the EWMA statistic is 0.1 or 0.2, and P c is 0.2 , 0.4, and 0.6 for lower and moderate censoring rates [8,13,15]. This paper builds the CNNSVR network with three combinational layers. The maximum epoch of network training is set to 1000 and the size and number of kernels for each combinational layer are as shown in Table 2. The complete CNNSVR network architecture is shown in Figure 4.
Table 3 shows the values of LCL for the EWMA CEV chart (EWMA chart) and η for the CNNSVR-based residual chart (CNNSVR chart), where these values fix their in-control ARL values at 200.
Table 4 and Table 5 show the out-of-control ARL values for shift sizes δ = 0.9 , 0.8 , 0.7 , 0.6 , 0.5 ,   and   0.2 , where # indicates when the CNNSVR chart is less efficient than the EWMA chart. Recall that the out-of-control ARL values of the CNNSVR chart can be obtained via the simulation procedure of Figure 3. For the case of λ = 0.1 in Table 4, as the skewness values become small, the out-of-control ARL values for P c = 0.2 and 0.4 also decrease. Therefore, both charts with an approximately symmetric distribution exhibit better detection efficiency for the lower censored rate. In the case of P c = 0.6 , as the skewness increases, both charts in all shift sizes become slow. Few out-of-control ARL values of the CNNSVR chart at a shift size of δ = 0.9 are larger than the values of the EWMA chart but, at other shift sizes, the CNNSVR chart outperforms the EWMA chart. The ARL values of both charts increase as the censoring rate increases. For most shift sizes, the CNNSVR chart demonstrates a superior ability to detect mean reduction.
As seen in Table 5, both charts produce the same results for λ = 0.1   and   0.2 . However, when comparing Table 4 and Table 5, the detection efficiency of λ = 0.2 outperforms that of λ = 0.1 for most shift sizes. Thus, the choice of a larger λ value significantly improves the detection efficiency.
For the cases of λ = 0.1   and   0.2 , the mean reduction of P c = 0.8 cannot be effectively detected, because the ARL values of both charts at small shift sizes ( δ = 0.9 , 0.8 ,   and   0.7 ) almost exceed 1000. This study increases λ to 0.8 for the censored data of P c = 0.8 to simulate the ARL values. Table 6 shows the ARL values, LCL, and η for the censored data of P c = 0.8 . The CNNSVR still exhibits a better detection efficiency than the EWMA chart at most shift sizes for higher censored data. As the skewness values become smaller, CNNSVR is clearly better than EWMA statistic in terms of detection.

5. Illustrative Example for Implementing a CNNSVR-Based Residual Chart

For a liquid crystal display module (LCM), the reliability test is the primary method to determine its lifetime and it is conducted at a temperature of 70 °C and humidity of 80%. LCM lifetimes follow a Weibull distribution of a 0 = 3.1 and b 0 = 2.2 . According to customer requirements, five units of each batch of LCMs must be randomly sampled to test lifetime values with a 60% censored rate (censored time is 2.284 h, and CEV is 0.997) before shipment. The process mean of an in-control state is 2.745. Practitioners referred to the parameters of the EWMA chart and developed the CNNSVR-based residual chart using λ = 0.2 and the in-control ARL = 200.
Practitioners may refer to Section 3.1 to set up the CNNSVR-based residual chart, including generating 10,000 EWMA statistics for the training data, such that the training data consist of 9999 records. Table 7 shows the calculation of the first three EWMA statistics. The s i values are five random numbers generated by the Weibull distribution of a 0 = 3.1 and b 0 = 2.2 . If s i > 2.284 , then x i = 0.997 ; otherwise, x i = s i . Then, X t ¯ is the average of the five x i values. The EWMA statistic E t is calculated using Equation (2). The first two records in this training dataset are (input, output) = 2.385 , 2.099 and 2.099 , 1.863 . Following the procedure in Figure 2, the training dataset is inputted into the CNNSVR network and the optimal threshold value η is calculated to be 0.124 .
Next, the CNNSVR-based residual chart is implemented to monitor the lifetime of LCM. Practitioners sample five units from each batch for life testing under the conditions of 70 °C and 80% humidity and halt testing when the test time reaches 1.568 h. Data from five testing units of each batch are used to calculate the EWMA statistic. The EWMA statistic is inputted into the CNNSVR network to predict the statistic of the next period. As shown in Table 8, the EWMA statistic of the 101st batch is 0.29, so the practitioners input 0.29 into the CNNSVR network and the CNNSVR network outputs the predicted value of 2.71 for the 102nd batch. After the 102nd batch has been produced and the lifetime tests of five units are completed, the actual statistic of this batch is calculated to be 2.72 and the error value of this batch is 0.01 (=2.72 − 2.71). Practitioners plot the error value of 0.01 on the CNNSVR-based residual chart. Figure 5 shows the CNNSVR-based residual chart for the 29 error values shown in Table 8. The error values of batches 120–130 are below the η value, and so this chart signals a variation at the 120th batch.

6. Conclusions

In this paper, a residual chart combining a convolutional neural network (CNN) and support vector regression (SVR) with an exponentially weighted moving average (EWMA) and conditional expected value (CEV) statistics was proposed to detect the mean decrease for Weibull type-I right-censored data. As demonstrated by a comparison of the average run lengths, the CNNSVR-based residual chart generally outperforms the EWMA CEV chart at low, moderate, and high censoring rates. The CNNSVR-based residual chart for the Weibull parameters with approximate symmetry achieves better detection efficiency. As the skew coefficient value of the Weibull distribution increases, the detection efficiency decreases. In a heavily skewed Weibull distribution or higher censored rate, the use of a larger smoothing parameter of the EWMA statistic increases the performance of the CNNSVR-based residual chart. Future work could be extended to use double EWMA CEV statistics to monitor other types of censored data with normal or non-normal distributions. In addition, combining multiple deep-learning methods to build control charts is an interesting future study.

Author Contributions

P.-H.L.: Conceptualization, methodology, writing, supervision, and project administration; S.-L.L.: software, validation, and investigation. 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 for this study come from computer simulations and industry data provided by a private company. The raw data cannot be made public due to privacy or ethical restrictions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khan, N.; Aslam, M.; Raza, S.M.M.; Jun, C.H. A new variable control chart under failure-censored reliability tests for Weibull distribution. Qual. Reliab. Eng. Int. 2019, 35, 572–581. [Google Scholar] [CrossRef]
  2. Xu, S.; Jeske, D.R. Weighted EWMA charts for monitoring type I censored Weibull lifetimes. J. Qual. Technol. 2018, 50, 220–230. [Google Scholar] [CrossRef]
  3. Yu, M.; Zhao, W.; Zhou, Y.; Wu, C. Robust online detection on highly censored data using a semi-parametric EWMA chart. J. Stat. Comput. Simul. 2023, 93, 1403–1419. [Google Scholar] [CrossRef]
  4. Steiner, S.H.; MacKay, R.J. Frontiers in Statistical Quality Control, 6th ed.; Springer: Berlin, Germany, 2001; pp. 275–289. [Google Scholar] [CrossRef]
  5. Steiner, S.H.; MacKay, R.J. Monitoring processes with data censored owing to competing risks by using exponentially weighted moving average control charts. J. R. Stat. Soc. Ser. C (Appl. Stat.) 2001, 50, 293–302. [Google Scholar] [CrossRef]
  6. Steiner, S.H.; MacKay, R.J. Monitoring processes with highly censored data. J. Qual. Technol. 2000, 32, 199–208. [Google Scholar] [CrossRef]
  7. Lee, P.-H. Economic design of a CEV X ¯ control chart for determining optimal right-censored times. Qual. Technol. Quant. Manag. 2021, 18, 418–431. [Google Scholar] [CrossRef]
  8. Zhang, L.; Chen, G. EWMA charts for monitoring the mean of censored Weibull lifetimes. J. Qual. Technol. 2004, 36, 321–328. [Google Scholar] [CrossRef]
  9. Bizuneh, B.; Wang, F.K. Comparison of different control charts for a Weibull process with type-I censoring. Commun. Stat.-Simul. Comput. 2019, 48, 1088–1100. [Google Scholar] [CrossRef]
  10. Raza, S.M.M.; Siddiqi, A.F. EWMA and DEWMA control charts for poissonExponential distribution: Conditional median approach for censored data. Qual. Reliab. Eng. Int. 2017, 33, 387–399. [Google Scholar] [CrossRef]
  11. Tsai, T.R.; Lin, C.C. The design of EWMA control chart for average with type-I censored data. Int. J. Qual. Reliab. Manag. 2009, 26, 397–405. [Google Scholar] [CrossRef]
  12. 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–1005. [Google Scholar] [CrossRef]
  13. Raza, S.M.; Riaz, M.; Ali, S. On the performance of EWMA and DEWMA control charts for censored data. J. Chin. Inst. Eng. 2015, 38, 714–722. [Google Scholar] [CrossRef]
  14. Ali, S.; Raza, S.M.; Aslam, M.; Butt, M.M. CEV-Hybrid Dewma charts for censored data using Weibull distribution. Commun. Stat.-Simul. Comput. 2021, 50, 446–461. [Google Scholar] [CrossRef]
  15. Raza, S.M.; Ali, S.; Butt, M.M. DEWMA control charts for censored data using Rayleigh lifetimes. Qual. Reliab. Eng. Int. 2018, 34, 1675–1684. [Google Scholar] [CrossRef]
  16. Lee, P.-H.; Torng, C.-C.; Kao, P.-H.; Chou, C.-Y. Performance of two-sided EWMA CEV control charts with multiple censored data. Int. J. Ind. Eng. Theory Appl. Pract. 2022, 29, 487–498. [Google Scholar] [CrossRef]
  17. Zhao, W.; Wu, C. Monitoring the alternating renewal processes with Weibull window-censored data. Qual. Technol. Quant. Manag. 2023, 20, 468–484. [Google Scholar] [CrossRef]
  18. Chen, S.; Yu, J. Deep recurrent neural network-based residual control chart for autocorrelated processes. Qual. Reliab. Eng. Int. 2019, 35, 2687–2708. [Google Scholar] [CrossRef]
  19. Farokhnia, M.; Niaki, S.T.A. Principal component analysis-based control charts using support vector machines for multivariate non-normal distributions. Commun. Stat.-Simul. Comput. 2020, 49, 1815–1838. [Google Scholar] [CrossRef]
  20. Kim, J.-M.; Ha, I.D. Deep learning-based residual control chart for count data. Qual. Eng. 2022, 34, 370–381. [Google Scholar] [CrossRef]
  21. Lee, S.; Kim, S.B. Time-adaptive support vector data description for nonstationary process monitoring. Eng. Appl. Artif. Intell. 2018, 68, 18–31. [Google Scholar] [CrossRef]
  22. Maboudou-Tchao, E.M. Change detection using least squares one-class classification control chart. Qual. Technol. Quant. Manag. 2020, 17, 609–626. [Google Scholar] [CrossRef]
  23. Wang, F.-K.; Mamo, T. Hybrid approach for remaining useful life prediction of ball bearings. Qual. Reliab. Eng. Int. 2019, 35, 2494–2505. [Google Scholar] [CrossRef]
  24. Yeganeh, A.; Abbasi, S.A.; Shongwe, S.C.; Malela-Majika, J.-C.; Shadman, A.R. Evolutionary support vector regression for monitoring Poisson profiles. Soft Comput. 2023. [Google Scholar] [CrossRef]
  25. Yeganeh, A.; Chukhrova, N.; Johannssen, A.; Fotuhi, H. A network surveillance approach using machine learning based control charts. Expert Syst. Appl. 2023, 219, 119660. [Google Scholar] [CrossRef]
  26. Maged, A.; Lui, C.F.; Haridy, S.; Xie, M. Variational AutoEncoders-LSTM based fault detection of time-dependent high dimensional processes. Int. J. Prod. Res. 2023. [Google Scholar] [CrossRef]
  27. Yu, J.; Li, S.; Liu, X.; Gao, Y.; Wang, S.; Liu, C. Dynamic convolutional gated recurrent unit attention auto-encoder for feature learning and fault detection in dynamic industrial processes. Int. J. Prod. Res. 2023, 61, 7434–7452. [Google Scholar] [CrossRef]
  28. Zhang, C.; Yu, J.; Wang, S. Fault detection and recognition of multivariate process based on feature learning of one-dimensional convolutional neural network and stacked denoised autoencoder. Int. J. Prod. Res. 2021, 59, 2426–2449. [Google Scholar] [CrossRef]
  29. Lee, P.-H.; Torng, C.-C.; Lin, C.-H.; Chou, C.-Y. Control chart pattern recognition using spectral clustering technique and support vector machine under gamma distribution. Comput. Ind. Eng. 2022, 171, 108437. [Google Scholar] [CrossRef]
  30. Lu, Z.; Wang, M.; Dai, W. A condition monitoring approach for machining process based on control chart pattern recognition with dynamically-sized observation windows. Comput. Ind. Eng. 2020, 142, 106360. [Google Scholar] [CrossRef]
  31. Xu, J.; Lv, H.; Zhuang, Z.; Lu, Z.; Zou, D.; Qin, W. Control Chart Pattern Recognition Method Based on Improved One-dimensional Convolutional Neural Network. IFAC-Pap. 2019, 52, 1537–1542. [Google Scholar] [CrossRef]
  32. Xue, L.; Wu, H.; Zheng, H.; He, Z. Control chart pattern recognition for imbalanced data based on multi-feature fusion using convolutional neural network. Comput. Ind. Eng. 2023, 182, 109410. [Google Scholar] [CrossRef]
  33. Yu, Y.; Zhang, M. Control chart recognition based on the parallel model of CNN and LSTM with GA optimization. Expert Syst. Appl. 2021, 185, 115689. [Google Scholar] [CrossRef]
  34. Zan, T.; Liu, Z.; Wang, H.; Wang, M.; Gao, X. Control chart pattern recognition using the convolutional neural network. J. Intell. Manuf. 2020, 31, 703–716. [Google Scholar] [CrossRef]
  35. Sabahno, H.; Amiri, A. New statistical and machine learning based control charts with variable parameters for monitoring generalized linear model profiles. Comput. Ind. Eng. 2023, 184, 109562. [Google Scholar] [CrossRef]
  36. Wang, F.-K.; Bizuneh, B.; Cheng, X.-B. One-sided control chart based on support vector machines with differential evolution algorithm. Qual. Reliab. Eng. Int. 2019, 35, 1634–1645. [Google Scholar] [CrossRef]
  37. Yu, J.; Liu, X. One-dimensional residual convolutional auto encoder for fault detection in complex industrial processes. Int. J. Prod. Res. 2022, 60, 5655–5674. [Google Scholar] [CrossRef]
  38. Beniwal, M.; Singh, A.; Kumar, N. Forecasting long-term stock prices of global indices: A forward-validating Genetic Algorithm optimization approach for Support Vector Regression. Appl. Soft Comput. 2023, 145, 110566. [Google Scholar] [CrossRef]
  39. Huang, Y.; Lei, Y.; Luo, X.; Fu, C. Prediction of compressive strength of rice husk ash concrete: A comparison of different metaheuristic algorithms for optimizing support vector regression. Case Stud. Constr. Mater. 2023, 18, e02201. [Google Scholar] [CrossRef]
Figure 1. CNNSVR network.
Figure 1. CNNSVR network.
Mathematics 12 00074 g001
Figure 2. Procedure for determining the optimal threshold value η .
Figure 2. Procedure for determining the optimal threshold value η .
Mathematics 12 00074 g002
Figure 3. ARL simulation procedure.
Figure 3. ARL simulation procedure.
Mathematics 12 00074 g003
Figure 4. CNNSVR network architecture for Weibull censored data, KS: Kernel size; KN: Number of kernels; AF: Activation function.
Figure 4. CNNSVR network architecture for Weibull censored data, KS: Kernel size; KN: Number of kernels; AF: Activation function.
Mathematics 12 00074 g004
Figure 5. CNNSVR-based residual control chart for LCM process control.
Figure 5. CNNSVR-based residual control chart for LCM process control.
Mathematics 12 00074 g005
Table 1. Skewness of the Weibull distribution.
Table 1. Skewness of the Weibull distribution.
The Parameters for
Weibull   Distribution   ( a 0 ,   b 0 )
Process Mean
M 0
Process Standard
Deviation   V 0
Skewness
Coefficient Sk
(1, 0.8)1.131.432.81
(1, 1)112
(1, 2.5)0.890.380.36
(1, 4)0.910.25−0.09
Table 2. Sizes and numbers of kernels for the CNNSVR network.
Table 2. Sizes and numbers of kernels for the CNNSVR network.
LayerConvolution LayerPooling Layer
Combinational
Layer
#1Kernel size51
Number of kernels200-
Activation functionReLU
#2Kernel size31
Number of kernels70-
Activation functionReLU
#3Kernel size31
Number of kernels50-
Activation functionReLU
Table 3. Design parameter values of CNNSVR and EWMA charts.
Table 3. Design parameter values of CNNSVR and EWMA charts.
The Parameters for
Weibull   Distribution   ( a 0 ,   b 0 )
Censored
Rate   P c
Smoothing   Parameter   λ = 0.1 Smoothing   Parameter   λ = 0.2
EWMACNNSVREWMACNNSVR
LCL η LCL η
(1, 0.8)0.20.741−0.1090.642−0.123
0.40.352−0.0480.321−0.072
0.60.122−0.0180.112−0.028
(1, 1)0.20.776−0.1100.675−0.130
0.40.419−0.0450.387−0.077
0.60.168−0.0210.156−0.031
(1, 2.5)0.20.823−0.0590.773−0.088
0.40.712−0.0540.679−0.078
0.60.391−0.0300.374−0.046
(1, 4)0.20.863−0.0400.828−0.068
0.40.837−0.0490.806−0.069
0.60.497−0.0310.480−0.046
Table 4. ARL values of CNNSVR and EWMA charts for λ = 0.1 .
Table 4. ARL values of CNNSVR and EWMA charts for λ = 0.1 .
Censored
Rate   P c
Charts Shift   Size   δ Shift   Size   δ
0.90.80.70.60.50.20.90.80.70.60.50.2
( a 0 , b 0 ) = (1, 0.8)( a 0 , b 0 ) = (1, 1)
0.2EWMA45.821.413.810.38.05.438.519.09.57.25.43.4
CNNSVR71.1#19.013.17.94.42.360.5#18.79.25.33.91.4
0.4EWMA57.935.126.721.918.813.451.230.022.017.614.810.2
CNNSVR73.6#34.623.017.814.39.450.124.515.211.79.46.3
0.6EWMA126.179.655.543.936.725.5110.967.548.037.630.920.7
CNNSVR120.976.449.539.030.720.6132.0#66.244.832.025.616.1
( a 0 , b 0 ) = (1, 2.5)( a 0 , b 0 ) = (1, 4)
0.2EWMA27.39.34.63.22.51.316.24.82.62.01.51.0
CNNSVR27.29.33.42.21.61.016.14.31.71.21.11.0
0.4EWMA25.212.58.26.14.93.112.15.33.42.62.11.3
CNNSVR29.1#11.76.44.33.21.710.33.82.11.31.11.0
0.6EWMA225.699.845.527.819.710.41729.1922.180.727.317.08.1
CNNSVR256.3#89.742.122.215.07.0847.0534.879.923.812.95.1
Table 5. ARL values of CNNSVR and EWMA charts for λ = 0.2 .
Table 5. ARL values of CNNSVR and EWMA charts for λ = 0.2 .
Censored
Rate   P c
Charts Shift   Size   δ Shift   Size   δ
0.90.80.70.60.50.20.90.80.70.60.50.2
( a 0 , b 0 ) = (1, 0.8) ( a 0 , b 0 ) = (1, 1)
0.2EWMA48.619.110.97.45.83.543.416.38.76.04.52.6
CNNSVR44.818.29.64.32.61.442.915.06.94.52.81.3
0.4EWMA49.523.016.012.510.37.145.519.613.410.18.35.6
CNNSVR61.0#21.712.57.45.93.550.4#17.110.36.34.62.5
0.6EWMA114.558.434.924.820.012.9100.849.730.221.316.510.5
CNNSVR114.153.430.118.914.38.194.240.425.715.511.56.5
( a 0 , b 0 ) = (1, 2.5)( a 0 , b 0 ) = (1, 4)
0.2EWMA28.39.44.43.02.21.215.74.52.51.81.31.1
CNNSVR27.29.23.72.11.81.015.04.21.91.21.11.0
0.4EWMA20.98.75.54.03.22.011.94.42.72.11.71.0
CNNSVR20.27.13.02.01.41.013.4#3.41.61.21.01.0
0.6EWMA194.181.130.115.910.75.5852.2443.556.016.39.34.2
CNNSVR171.066.324.910.25.92.3634.3343.150.812.86.11.7
Table 6. ARL values of CNNSVR and EWMA charts for P c = 0.8 .
Table 6. ARL values of CNNSVR and EWMA charts for P c = 0.8 .
The Parameters for
Weibull   Distribution   ( a 0 , b 0 )
Charts LCL   or   η *
Shift   Size   δ
0.90.80.70.60.50.2
(1, 0.8)EWMA0.014133.088.058.033.120.15.2
CNNSVR−0.01692.486.656.431.828.4#5.0
(1, 1)EWMA0.023310.1187.6105.051.026.54.7
CNNSVR−0.030107.674.767.650.535.2#4.6
(1, 2.5)EWMA0.096642.9675.2715.1427.2138.73.4
CNNSVR−0.109112.383.392.4118.4101.73.2
(1, 4)EWMA0.1365080.111,220.012,193.08056.03636.03.5
CNNSVR−0.078317.91201.185.265.350.13.3
Table 7. The generation of training data for LCM.
Table 7. The generation of training data for LCM.
tSimulation Data
s i
Conditional Expected Value
x i
Sample Mean
X t ¯
EWMA
Statistic   E t
10.7341.5812.1151.6651.8510.7340.9970.9970.9970.9970.9442.385
20.8390.9362.3121.9831.6590.8390.9360.9970.9970.9970.9532.099
30.9360.6601.8462.5322.6150.9360.6600.9970.9970.9970.9171.863
Table 8. Predicted and error values of the CNNSVR network for the LCM lifetime test.
Table 8. Predicted and error values of the CNNSVR network for the LCM lifetime test.
No. of BatchActual Value
E t
Predicted Value
E t ^
Error Value
E t E t ^
No. of BatchActual Value
E t
Predicted Value
E t ^
Error Value
E t E t ^
1010.29--1162.692.470.22
1022.722.710.101172.682.74−0.06
1032.692.580.111182.672.490.18
1042.742.030.711192.752.740.01
1050.830.84−0.011202.332.60−0.27
1061.361.350.011211.722.04−0.32
1072.701.920.781222.072.41−0.34
1082.752.87−0.121232.352.59−0.24
1092.752.630.121240.120.45−0.33
1102.752.370.381251.772.20−0.43
1112.722.240.481261.852.41−0.56
1121.041.020.021270.181.14−0.96
1131.011.08−0.071281.202.07−0.87
1140.060.050.011291.042.10−1.06
1152.842.290.551300.241.74−1.50
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

Lee, P.-H.; Liao, S.-L. Residual Control Chart Based on a Convolutional Neural Network and Support Vector Regression for Type-I Censored Data with the Weibull Model. Mathematics 2024, 12, 74. https://doi.org/10.3390/math12010074

AMA Style

Lee P-H, Liao S-L. Residual Control Chart Based on a Convolutional Neural Network and Support Vector Regression for Type-I Censored Data with the Weibull Model. Mathematics. 2024; 12(1):74. https://doi.org/10.3390/math12010074

Chicago/Turabian Style

Lee, Pei-Hsi, and Shih-Lung Liao. 2024. "Residual Control Chart Based on a Convolutional Neural Network and Support Vector Regression for Type-I Censored Data with the Weibull Model" Mathematics 12, no. 1: 74. https://doi.org/10.3390/math12010074

APA Style

Lee, P. -H., & Liao, S. -L. (2024). Residual Control Chart Based on a Convolutional Neural Network and Support Vector Regression for Type-I Censored Data with the Weibull Model. Mathematics, 12(1), 74. https://doi.org/10.3390/math12010074

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