Economic Design of EWMA Control Charts with Variable Sampling Intervals for Monitoring the Mean and Standard Deviation under Preventive Maintenance and Taguchi’s Loss Functions

,e control chart and the maintenance management need to control process quality and reduce out-of-control cost. ,ey are two key tools in the production process; however, they have usually been analyzed separately in the literature. Moreover, the existing studies integrating these two concepts suffer from three significant drawbacks as follows: (1) using control charts with fixed parameters to monitor the process, so that the small and middle shifts are detected slowly; (2) monitoring the mean and standard deviation separately, whereas, in real condition, the mean and standard deviation should be monitored simultaneously; (3) quality loss function is not usually used to design economic model, which leads to a large social quality loss in the monitoring process of control chart. To eliminate these weaknesses, the economic design of the exponential weighted moving average (EWMA) control chart with variable sampling intervals (VSI) for monitoring the mean and standard deviation under preventive maintenance and Taguchi’s loss functions is proposed. ,e optimal values of the parameters are determined to minimize the loss-per-item in an average cycle process. In addition, a genetic algorithm is used in a numerical example to search for the optimal values of the parameters. According to the sensitivity analysis, the effect of the model parameters on the solution of the economic model is obtained. Finally, the comparison study shows that the VSI EWMA control charts designed by the joint economic model are expected to reduce loss.


Introduction
ere are two key tools for controlling production processes-statistical process control (SPC) and maintenance management (MM). Although these two tools arise from different research areas, their common goal is to improve product quality and reduce cost by eliminating variances in the manufacturing process. SPC is a powerful tool consisting of different types of control charts for monitoring the process, and it can be ensured that the quality characteristics of a product are at the nominal or required levels [1]. However, the control chart always neglects type II error, in which the control chart does not trigger a signal when the process is out-of-control. erefore, it will cause too high and avoidable costs on the process [2]. Maintenance in industrial situations is considered an important factor contributing to reducing production costs and increasing productivity. One way to decrease type II error and the cost for detecting the out-of-control signal is to combine the control charts with the maintenance management technique. is work is closely related to several research fields, including economic design of control charts, adaptive control charts, integration of control chart and maintenance management, and quality loss function. e existing literatures of such researches are reviewed briefly in the following.
To reduce the cost of process control, Duncan [3] first proposed the economic model to determine the three test parameters in the X-bar control chart, so that the average cost can be minimized when a single assignable cause occurs. Montgomery [4] gave a thorough review of the literature in the economic designs domain. Lorenzen and Vance [5] presented a framework of economic design for the control charts. Since then, considerable attention has been devoted to the optimal economic determination of the three parameters of X-bar charts [6]. Subsequently, many scholars have studied the economic design of various control charts. Franco et al. [7] developed an economic model of Shewhart control charts for monitoring autocorrelated data with skip sampling strategies; Naderkhani and Makis [8] investigated the economic design of multivariate Bayesian control chart with two sampling intervals; Costa [9] investigated the Economic statistical design of ARMA control chart through a modified fitness-based self-adaptive differential evolution. Here, we briefly review some researches on the economic design of EWMA control charts. To consider the uncertainty among the cost and process parameters in practice, Amiri et al. [10] developed a process monitoring strategy based on the robust economic and economic-statistical design of the EWMA control chart. Considering measurement error and taking multiple measurement errors, as well as linear and quadratic Taguchi loss function of poor quality products, Saghaei et al. [11] proposed the economic model of EWMA control charts. Chou et al. [12] developed the economic design of the VSI EWMA control chart to determine the values of the six test parameters of the charts, such that the expected total cost is minimized. ey assumed that the measurements in each sampled subgroup are normally distributed, but this assumption is not always true in practice.
In the case of the nonnormal behavior of measurements, Xue and Liu [13] used Burr distribution to approximate various nonnormal distributions and developed an economic model of the VSI EWMA control chart. In addition, many scholars have studied the economic design of nonparametric control chart. Li et al. [14] proposed an economic model based on the Duncan-type cost function for designing a nonparametric sign chart for monitoring the location parameter of a univariate process. When true process location parameter is unknown, Li et al. [15] developed an economically designed nonparametric control chart for monitoring unknown location parameter based on the Wilcoxon rank-sum (hereafter WRS) statistic. Li and Mukherjee [16] presented two distribution-free cost-efficient Shewhart-type schemes for sequentially monitoring process location with restricted false alarm probability, based, respectively, on the Sign and Wilcoxon rank-sum statistics. Li and Mukherjee [17] proposed two economically optimized nonparametric schemes for monitoring process variability based on two popular two-sample rank statistics for differences in scale parameters, known as the Ansari-Bradley statistic and the mood statistic. In addition to the above research, many scholars have carried out the economic design and reliability group of the sampling plan, for example [18][19][20][21][22][23].
With the development of automation and mechanization, many automatic techniques can be applied to production.
us, the status of equipment plays an increasingly important role in controlling quantity, quality, and cost. Although academia and industry have long recognized the close relationship between equipment maintenance and process quality, research on the integration and association between equipment maintenance decisions and process quality control decisions has not attracted widespread attention from relevant scholars until recently [24][25][26][27]. In this regard, Linderman et al. [28] constructed a generalized analytic model that incorporates statistical process control and maintenance policy to minimize total expected cost. Zhou and Zhu [29] developed a joint model that integrates the maintenance activities with the design of control charts. ey claimed that integrated models perform better than the two stand-alone models. Xiang [30] developed an integrated model of statistical process control and preventive maintenance for a manufacturing process that deteriorates according to a discrete-time Markov chain approach. Shrivastava et al. [31] presented an integrated model for joint optimization of preventive maintenance and quality control policy with cumulative sum (CUSUM) control chart parameters. Bouslah et al. [32] researched the joint economic design of production, continuous sampling inspection, and preventive maintenance of a deteriorating production system. Salmasnia et al. [33] developed an integrated model of economic production quantity, statistical process monitoring, and maintenance in the presence of multiple assignable causes. e particle swarm optimization algorithm (PSO) is used to minimize the expected total cost per production cycle, subject to statistical quality constraints. Wan et al. [34] developed an integrated model of MM and quality control policy with an adaptive synthetic chart, and the corrective maintenance and preventive maintenance were considered in that paper.
Although the research on joint optimization of preventive maintenance and control charts has made great progress, most of the previous research is about the invariance control chart, which means that the parameters are not unfixed. However, to improve the performance of detecting small and moderate shifts, different types of adaptive control charts such as variable sampling intervals (VSI), variable sample sizes (VSS), and variable parameters (VP), in which the design parameters are adaptive depending on the position of previous observation on the chart, have developed in recent years [35][36][37][38][39][40][41][42][43]. e research results have shown that the efficiency of the adaptive control chart over the traditional control charts, in which the parameters are constant. On the other hand, the mean and standard deviation usually were monitored by a control chart separately. However, the mean and standard deviation should be monitored using one kind of control chart simultaneously in real condition. In order to improve the efficiency of control chart monitoring small shifts, EWMA control charts for monitoring the process mean and standard deviation jointly with variable sampling intervals were constructed [44].
e concept of quality loss function was proposed by Taguchi firstly. Due to the high measurement cost, the mass quality loss function has been widely used. It focuses on the loss caused by the deviation of the mass output from the target value. If a characteristic measurement is equal to the target value, the loss is zero, but the further far away from its target value, the greater the quality loss.
erefore, the output characteristics should be as close to their target values as possible. Taguchi suggests using the quality loss function to measure the loss caused by the deviation of the output characteristic of the qualified product from the target value.
is concept has also been used in some models of economic design of control charts [45,46]. In addition, Sultana et al. [47] develop an economic statistical design of the (EWMA) chart using variable sampling interval at fixed times (VSIFT) control scheme considering preventive maintenance and Taguchi loss function. ey considered the possibility of an equipment failure in terms of machine breakdown or improper functioning of the equipment, which results in poor product quality and call for maintenance action, to study the linkage between declining performance of a machine and process. In particular, machine failures are divided into two failure modes in their study: (1) failure mode 1 (FM1) leads to immediate breakdown of the machine. (2) Failure mode 2 (FM2) leads to reduction in process quality owing to shifting of the process mean. eir study mainly monitors the shift of process mean and does not consider monitoring the shift of process mean and standard deviation at the same time. Besides, they do not focus on the determination of a warning limit or fixed time interval beyond which planned maintenance will be carried out.
To bridge the existing gaps in the literature, a novelty mathematical optimization model that integrates the concepts of quality loss, maintenance policy, and designing the control chart is proposed in this study. e main contributions of this study can be summarized as follows: (1) this study integrates the concepts of control charts, maintenance management, and quality loss function. It mainly considers preventive maintenance in one cycle but does not consider downtime maintenance; (2) the mean and standard deviation can be monitored simultaneously by EWMA control chart; (3) the variable sampling interval schemes should be applied to improve the monitoring efficiency of control chart. erefore, in this paper, the economic design of the VSI EWMA control chart for monitoring the process mean and standard deviation based on Taguchi's loss function and preventive maintenance is investigated. Six test parameters of the chart (i.e., the sample size, the long sampling interval, the short sampling interval, the control limit coefficient, the warning limit coefficient, and the exponential weight constant) are determined by minimizing the total loss. A numerical example is used to illustrate how the model works and to demonstrate its use.

VSI EWMA Control Charts for Monitoring the Mean and Standard Deviation
Before the economic model of VSI EWMA control charts is established, we first introduce the VSI EWMA control charts for monitoring the mean and standard deviation. e details of the approach are shown as follows.
Assume that the observations Y are followed by the normal distribution with mean μ and standard deviation σ. When the process is in control, μ � μ 0 , σ � σ 0 . Meanwhile, when the process is on out-of control, μ � μ 0 +δ 1 σ 0 , or, σ � δ 2 σ 0 , where δ 1 and δ 2 are constants, and typically δ 2 is larger than 1. Let X i1 , X i2 , . . . , X in be the ith sample, which the sample size is n, and then is the sample variance of the ith sample. en, the statistic of Semicircle (SC) chart is defined as [48] When the process is in control, T * i ∼ χ 2 n , so we can get the mean and variance of T * i as follows: en, the ith sample statistic of EWMA chart for monitoring the mean and standard deviation is defined as [49] where λ is the exponential weight constant, and 0 < λ ≤ 1, and the starting value Q 0 is set to be n [49]. According to equations (2) and (3), we can get the mean and variance of Q i as follows: When t is getting larger, Var(Q i ) will approach the steady-state value. Equation (5) will become Because Q i ≥ 0, the upper control limits UCL and the upper warning limits UWL can be decided.
where k is the control limit coefficient of the VSI EWMA chart, w is the warning limit coefficient of the VSI EWMA chart, and 0 < w < k. e in-control region of the control chart is divided into two parts, including the safe region and the warning region. Let h 1 and h 2 be the long sampling interval and the short Mathematical Problems in Engineering sampling interval, respectively. If the sample statistic falls in the safe region (i.e., 0 ≤ Q i ≤ UWL), then take the next sample using the sampling interval h 1 . If the sample statistic falls in the warning region (i.e., UWL < Q i ≤ UCL), then take the next sample using the sampling interval h 2 . If the sample statistic falls outside the control limits (i.e., Q i > UCL), then we will find the assignable cause.

Development of Cost Function
To develop the cost function, some assumptions should be made as follows: (1) e process is assumed to be in-control (μ � μ 0 , σ � σ 0 ) initially (2) e process will be disturbed by a single assignable cause that causes a fixed shift in the process (μ � μ 0 + δ 1 σ 0 , σ � δ 2 σ 0 ) (3) e time between occurrences of the assignable causes is exponential with a mean of θ occurrences per hour (4) Once the process goes out of control, it remains out of control until being detected and corrected (5) e process shift is instant and cannot automatically be brought back to an in-control state (6) Preventive maintenance will be conducted when the statistic Q i falls in the warning region e expected cycle length is defined as the time span from the beginning of an in-control state to the next incontrol state. In each production cycle, once an assignable cause occurs, it can be detected and identified. After that, the process will be brought back to the in-control process condition. e production cycle length can be divided into four-time intervals: in-control period T 1 ; out-of-control period starting from the occurrence of an assignable cause to the control chart signals, denoted as T 2 ; the time period for identifying and correcting the assignable cause T 3 , and the time period for sampling, inspecting, evaluating, and plotting the subgroup result T 4 .

e Average Cycle Length
(1) e expected length of the in-control period can be written as where t 0 is expected search time when the signal indicates a false alarm, ANSS 0 is the in-control average run length, and s is the expected number of samples taken when the process is in control and can be calculated by [12] s � e − θh 0 c 1 � 1, if production continues during searches, 0, if production ceases during searches. Note that h 0 can be expressed as where ATS 0 is the in-control average time to signal. (2) e out-of-control period before the last sample of the detecting subgroup is taken can be expressed as where ATS 1 is the out-of-control average time to signal, and τ is the expected time of occurrence of the assignable cause between two samples when the process is in control. It can be shown that [12] τ (3) e time period for conducting preventive maintenance is equal to (4) e time period for identifying and correcting the assignable cause is equal to where t 1 is the expected time to identify the assignable cause, and t 2 is the expected time to correct the assignable cause (T P < t 2 ). (5) Define g as the average sampling, inspecting, evaluating, and plotting time for each sample. e time period for sampling, inspecting, evaluating, and plotting the subgroup result is According to equations (8), (11), (13)- (15), therefore, the average cycle length can be obtained:

e Cost Function.
To reduce the quality loss of the production process, the concept of Taguchi's loss function can be applied to the economic design model of the VSI EWMA control chart for monitoring the mean and standard deviation in this section. en, the expected loss cost generated in a cycle length includes several relevant costs, namely, (1) the loss cost of false detection of signal, sampling, inspecting, evaluating, plotting, and actually finding, repairing the assignable cause, (2) the preventive maintenance loss cost, (3) the social loss cost in control state, and finally (4) the social loss cost in the out of control state. We denote d as the average cost per false alarm, W as the cost to locate and repair an assignable cause, a as the fixed cost per sample, b as the cost per unit sampled, and d as the cost due to loss of a false detection of signal. If production continues during the correction process, c 2 � 1, and if production ceases during the correction process, c 2 � 0.
en, the cost of process sampling, inspecting, evaluating, plotting, and actually finding, repairing the assignable cause, can be shown as where h 0 is the expected sampling interval when the process is in-control, and h 0 / is the expected sampling interval when the process is out-of-control. h 0 / can be expressed as h 0 / � (ATS 1 /ANSS 1 ).
If we let C pm be the average cost per preventive maintenance action, then the preventive maintenance cost can be expressed as where p 01 is the probability that samples fall into the warning region when the process is in-control, and p 11 is the probability that samples fall into the warning region when the process is out-of-control, respectively. Denote the probability that the samples fall into the safe region when the process is in-control by p 00 , and the probability that the samples fall into the safe region when the process is out-ofcontrol by p 10 , respectively. en p 01 and p 11 can be expressed as Suppose that the specification limit of a quality characteristic is m ± Δ, where m is the target value. It will cost M dollars to repair the item and bring the shifted process back to in-control state. en, the value for the coefficient of loss function can be obtained by (M/Δ 2 ). Suppose that the process is still in-control. Let dv be the difference between the average and target, and let σ be the standard deviation of the process. en, the in-control average social loss can be expressed as where y is the production quantity per unit time.
If an assignable cause occurs, it will cause loss M due to a greater percentage of items being outside the control limits.
Hence, from equations (17), (18), (20), and (21), the loss in an average cycle process is given by From equations (16) and (22), the loss-per-item in an average cycle process is given by e economic design of the VSI EWMA control chart for monitoring the mean and standard deviation under preventive maintenance and Taguchi's loss functions can be determined by the optimal values of the eight test parameters (n, h 1 , h 2 , k, w, λ) such that the loss-per-item in an average cycle process ETL is minimized.

An Example and Solution Procedure
From the examination of the probability components in equation (23), it can be seen that determining the economically optimal values of the six test parameters for the VSI EWMA control chart for monitoring the mean and standard deviation is not straightforward. To illustrate the nature of the solutions obtained from economic design of VSI EWMA control chart, a particular numerical example is provided.
A bearing factory produces various shapes of casting, so that the castings manufacturing process should be monitored by the VSI EWMA control chart. It is known that the diameter of casting approximately follows a normal distribution with a mean μ and standard variation σ. When the process is in control, μ � μ 0 � 2.0, σ � σ 0 � 1. When the process is out of control, μ � μ 1 � μ 0 + δ 1 σ 0 , σ � δ 2 σ 0 � 1.5. e cost and model parameters in this example are given as follows:

Mathematical Problems in Engineering
e optimal design parameters n, h 1 , h 2 , k, w, λ can be derived by minimizing the ETL. Because the expression of ETL is very complex, the general optimization method cannot be used. Although there are several kinds of numerical optimization methods, Genetic Algorithm (GA) has advantages in the following aspects: (1) it can be applied for different types of optimization problems; (2) it can avoid traps of the local optimum; (3) it can give the global optimum solution efficiently. Hence, the optimization problem can be solved using GA to get the optimal solution.
To compare the proposed control chart with the existing control chart, we get the optimal design parameters n, h, k, λ of fixed sampling intervals EWMA control chart that considers quality loss and preventive maintenance by minimizing the ETL. e constraint conditions for each test parameter are set as follows: 1 ≤ n ≤ 30, 1 < h < 3, 0.001 < k < 4, 0.001 < λ < 1. After running the toolbox, we obtain the optimal values of n, h 1 , h 2 , k, w, λ as follows: n � 1, h � 0.681, k � 2.764, r � 0.305, ETL � 18.852. e result shows that the VSI EWMA control charts designed by the joint economic model have a less expected loss and are superior to the FSI EWMA control charts designed by the joint economic model.

Sensitivity Analysis
Using orthogonal-array experimental design and multiple regression method, a sensitivity analysis of the economic model of the VSI EWMA chart for monitoring the mean and standard deviation is conducted to study the effect of model parameters on the solution. e model parameters (a, b, θ, d, g, t 1 , t 2 , e, C pm , W, M) are independent variables, and the six design parameters (n, h 1 , h 2 , k, w, λ) and the lossper-item in an average cycle process ETL are dependent variables. Independent variables (model parameters) and their corresponding level planning are shown in Table 1. is is a trial with eleven factors and two levels, so the L 16 orthogonal array is employed, and there are 16 trials. For each trial, the genetic algorithm is applied to give the optimal solution of the economic design, with the following model parameters fixed: c 1 � c 2 � 1, t 0 � 1. e 16 trials selected according to the orthogonal Table L 16 (2 15 ) are shown in Table 2, and the output of the genetic algorithm for each trial is recorded in Table 3.
To study the effect of model parameters on the solution of economic design of the VSI EWMA chart for monitoring the mean and standard deviation, the statistical software SPSS is used to run the regression analysis for each dependent variable with the significance level equal to 0.1. From the output of SPSS with design parameters and lossper-item in an average cycle process ETL, we can obtain the model parameters that significantly impact the solution of the economic design of the VSI EWMA chart.
e SPSS output for sample size n is presented in Table 4. From the ANOVA in Table 4 (1), if the significance level is set to be 0.1, there is one parameter affecting the sample size n significantly at least. From the regression coefficients in Table 4 (2), we find that the cost per unit sampled b, the expected time to identify correct the assignable cause t 2 , the average sampling, inspecting, evaluating, and plotting time for each sample g and the average cost per preventive maintenance action C pm affect the sample size significantly. Since the coefficients are negative, the larger the four parameters, the larger the sample size. e output of SPSS for the long sampling interval h 1 is recorded in Table 5. From the ANOVA in Table 5 (1), there is one parameter affecting the long sampling interval h 1 significantly at least. It can be seen from the regression coefficients in Table 5 (2) that five model parameters significantly affect the long sampling interval. An increase in the fixed cost per sample a and the deviation of process standard deviation δ 2 lengthen the long sampling interval, respectively. An increase in the average sampling, inspecting, evaluating, and plotting time for each sample g, the cost to locate and repair an assignable cause W and the loss incurred by the nonconformity product M negatively affect the long sampling interval. Table 6 is the SPSS output for the short sampling intervals h 2 . It can be seen from Table 6 (2) that five model parameters significantly affect the short sampling intervals. An increase in the average cost per preventive maintenance action C pm , the expected time to identify the assignable cause t 1 and the expected time to correct the assignable cause t 2 lengthen the short sampling intervals, respectively. e cost per unit sampled b and the cost to locate and repair an assignable cause W negatively affect the short sampling interval. Table 7 is the SPSS output for the control limit coefficient k. From the regression coefficients in Table 7 (2), we find that the loss incurred by the nonconformity product M significantly and positively affects the control limit coefficient. Table 8 is the SPSS output for the warning limit coefficient w. It may be seen from Table 8 (2) that five model parameters significantly affect the warning limit coefficient. e larger the average cost per false alarm d, the larger the warning limit coefficient. e occurrences frequency of the assignable cause θ, the loss incurred by the nonconformity product M, the deviation of process standard deviation δ 2 and the expected time to identify correct the assignable cause t 2 , all negatively affect the warning limit coefficient. e output of SPSS for the exponential weight constant λ is recorded in Table 9. From the regression coefficients in Table 9 (2), we find that the process mean shift δ 1 and the deviation of process standard deviation δ 2 both significantly and positively affect the warning limit coefficient. Table 10 is the SPSS output for loss-per-item in an average cycle process ETL. It may be seen from Table 10 (2) that the occurrences frequency of the assignable cause θ, the loss incurred by the nonconformity product M and the deviation of process standard deviation δ 2 significantly and positively affect the loss-per-item in an average cycle process.
In order to intuitively illustrate the influence of the model parameters on the design parameters of the control chart and the loss cost function per unit time, we present the main effect analysis chart obtained by Minitab software in Figure 1.

Comparison Study
Assume that the distribution of the observations Y from a process is normally distributed and has the mean μ and standard deviation σ. When the process is in control μ � μ 0 , σ � σ 0 , and when the process is out of control μ � μ 0 + δ 1 σ 0 , σ � δ 2 σ 0 , where δ 1 and δ 2 are two constants, and typically δ 2 is larger than 1.
is process can be monitored by the VSI EWMA control charts. e optimality of the joint economic model of the VSI EWMA control chart for monitoring the process mean and standard deviation is verified by comparing the following three control charts: (1) VSI EWMA control chart that is designed based on the economic model (23); (2) FSI EWMA control chart that considers quality loss and preventive maintenance; (3) VSI EWMA control chart, which is designed by using the average time to signal (ATS) as the evaluation criterion. e parameters of the economic model (a, b, θ, d, g, t 1 , t 2 , δ, C pm , W, M) are shown in Table 1. e expected loss of the control chart designed by method (1) is denoted by ETL, the expected loss of the control chart designed by method (2) is denoted by ETL 1 , and the expected loss of the control chart designed by method (3) in cases where λ � 0.1, ATS � 300, 400, 500 and λ � 0.2, ATS � 300, 400, 500 is denoted as ETL 2 . e parameters in method     (1) ANOVA  Table 4: SPSS output for sample size n.
First, the economic advantages of the dynamic control chart designed by method (1) and the static control chart designed by method (2) are investigated. As can be seen from Table 11, in all the trials considered here, ETL are less than ETL 1 . For example, in the third test, the expected loss of VSI EWMA control chart designed by method (1) is 38.6429, and that of static EWMA control chart designed by method (2) is 40.1353, which is larger than that by method (1). e average values of 16 experiments were calculated and recorded on the table. e average value of ETL is 39.5824, which is less than the average value of ETL 1 41.8791. at is to say, VSI EWMA control charts designed by the joint economic model (23) are better than FSI EWMA control charts considering quality loss and preventive maintenance. en, the economic efficiency of the dynamic control chart designed by method (1) and method (3) is compared. Table 4 shows that, for the 16 tests considered here, ETL are less than ETL 2 ; for example, in the second test, the corresponding ETL of the VSI EWMA control chart designed by the joint economic model (23) is 38.8650, while ETL 2 of the VSI EWMA control chart that uses the average time to signal

Conclusion
To determine the values of the six design parameters of the VSI EWMA control charts, such that loss-per-item in an average cycle process is minimized, we developed an economic design of the VSI EWMA control chart for monitoring the mean and standard deviation under preventive maintenance and Taguchi's loss function. e method is illustrated using an example, and GA is employed to search for the optimal parameters of the  Figure 1: e main effect analysis chart of the design parameter (n, h 1 , h 2 , k, w, λ) and the loss-per-item in an average cycle process ETL. economic design. e effect of model parameters on the optimal parameters can be obtained from a sensitivity analysis of the economic model. Finally, the comparison study results show that the VSI EWMA control charts designed by the joint economic model is superior to the FSI EWMA control charts designed by the joint economic model and the VSI EWMA control chart that uses ATS as evaluation criteria. In this paper, we only consider cases where observations are univariate. In many applications, multivariate performance variables are involved. us, joint economic design of control charts for monitoring multivariate quality characteristics and preventive maintenance will need to be explored in future research.
In the traditional tests under classical statistics, it is assumed that all observations are crisp in the population or the sample. But, the data obtained from the complex system may not be determined, exact, and certain. In this case, classical statistics are no longer suitable, which are replaced by Neutrosophic statistics. Neutrosophic statistics is used when the data is obtained from the complex process. Almarashi and Aslam [50] proposed the control chart for monitoring the Gamma distributed product under neutrosophic statistics using resampling scheme. In consequence of the existing Anderson-Darling test that cannot be applied for testing the assumption of the Weibull distribution, Aslam [51] presented the Anderson-Darling test under neutrosophic statistics. Khan et al. [52] presented new attribute control charts for monitoring the blood components under the neutrosophic statistics. e applications of these control charts demonstrate that the proposed control charts are quite effective, adequate, flexible, and informative for monitoring the blood components under uncertain environment. Aslam [53] designed a control chart for neutrosophic exponentially weighted moving average (NEWMA) using repetitive sampling, and the application of the proposed NEWMA chart is given to monitoring road traffic crashes (RTC). From these literatures, it is observed that the control chart designed by Neutrosophic statistics is an efficient addition in the tool kit of the quality control personnel; thence, when the data is obtained from the complex process, the proposed chart in this paper can be extended for neutrosophic statistics in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Conflicts of Interest
e authors declare that they have no conflicts of interest.