An EWMA-based process mean estimator with dynamic
tuning capability
RUEY-SHAN GUO1 and JIN-JUNG CHEN2
1Department of Industrial Management and Business Administration, National Taiwan University, Taipei, Taiwan 2Department of Mechanical Engineering, National Taiwan University, Taipei, Taiwan
E-mail:rsguo@ccms.ntu.edu.tw or d3502121@ms.cc.ntu.edu.tw Received April 2000 and accepted November 2001
In this paper, we focus on the Exponentially Weighted Moving Average (EWMA) process mean estimator and its application to process adjustment. A novel dynamic-tuning EWMA estimator is proposed that has the capability of adjusting the control parameter dynamically in response to the underlying process random shifts. The current run’s process mean is estimated using the EWMA equation and the newly adjusted control parameter. It is shown that the proposed estimator is very easy to implement and effective under many disturbance situations. Both industrial field data and Monte Carlo simulations are used to validate its performance.
1. Introduction
Manufacturing processes are subject to various special disturbances such as process drifts or shifts. From a quality control point of view, the estimation of the cur-rent process mean provides opportunities for quality monitoring and fault diagnosis. In many cases when the resulting process output deviation can be adjusted to bring the process output to the target value, a good mean estimator will certainly provide a more accurate evalua-tion on how much the adjustment should be made.
In this paper, we focus on the Exponentially Weighted Moving Average (EWMA) process mean estimator and its application to process adjustment. In recent years, the EWMA-based process mean estimator or adjustment has been studied by many researchers including Mac Gregor (1988), Box and Kramer (1992), Ingolfsson and Sachs (1993), Butler and Stefani (1994), Yashchin (1995), Del Castillo (1999) and also Patel and Jenkins (2000). Box and Jenkins (1976) proved that the EWMA-based mean estimation and adjustment has a Minimum Mean Square Error (MMSE) property when the underlying process disturbance follows an IMA(1,1) (first-order integrated moving average) process. In the semiconductor process industry, several applications of the EWMA estimator to process adjustment were proved to be effective by Sachs et al. (1995) and also Del Castillo and Hurwitz (1997) even when the underlying disturbances do not follow the IMA pattern but follow an approximately linear drift pattern (Boning et al., 1996; Smith et al., 1998).
Although effective in many applications, the EWMA-based process mean estimator still has space for en-hancement. To be more specific, the estimation efficiency is highly dependent on the choice of the control para-meter, and in many cases it is pre-determined for processes with small disturbances. A dynamic-tuning capability to accommodate both small and large process disturbances is required.
One such enhancement, the Adaptive Exponentially Weighted (AEW) estimator, was proposed by Yashchin (1995) for the estimation of the current mean of a process subject to abrupt changes. The AEW estimator calculates the mean based on identification of the most recent change point. Another enhancement by Patel and Jenkins (2000) used the recursive scheme for optimizing the gain of an EWMA controller under stability constraints. The objective is to minimize the asymptotic mean square error in the output with minimal a priori information. Both shift and drift disturbances are considered. In this paper, we propose a different enhancement approach. We extend the application of EWMA estimation to processes sub-ject to random shifts of different magnitudes. A novel Dynamic-Tuning EWMA (DT-EWMA) process mean estimator is proposed. It will be shown that the proposed algorithm is very easy to implement and effective under many disturbance situations. Both industrial field data and simulations are used to validate its performance.
The rest of this paper is organized as follows. In tion 2, we introduce the DT-EWMA estimator. In Sec-tion 3, the estimaSec-tion performance among the EWMA, AEW, the adaptive EWMA estimator of Patel and
Jenkins (2000) and the DT-EWMA estimator is analyzed. In Section 4, we discuss the process adjustment applica-tions through simulaapplica-tions. An industrial application of the DT-EWMA estimation is described in Section 5. Conclusions are made in Section 6.
2. Dynamic-tuning EWMA estimator
2.1. Function overview
Before introducing the proposed Dynamic-Tuning EWMA (DT-EWMA) estimator, we first review the EWMA estimation algorithm. Consider a sequence of independent observationsfYtg, where Ythas mean lt and standard deviation r. Then the EWMA estimation of the current process mean of observation t is defined as:
^ l
lt ¼ WYtþ ð1 W Þ^llt1; ¼ WYtþ W ð1 W ÞYt1þ
þ W ð1 W Þt1Y1þ W ð1 W ÞtY0; ð1Þ
where W denotes the EWMA control parameter
(0 W 1). It has been proven that the EWMA esti-mator has optimality properties within the class of linear estimators for estimating the current mean of processes subject to a step change (Yashchin, 1995).
The estimation efficiency of the EWMA estimator, however, is highly dependent on the choice of the control parameter W. A fixed W is not always efficient in mean estimation under different process disturbance condi-tions. An adaptive capability is required to adjust the control parameter dynamically in response to the under-lying process disturbance conditions. Figure 1 shows the function diagram of the proposed DT-EWMA estimator. First the newly collected process output data is plotted using Statistical Process Control (SPC) charts. If there is any evidence of shift occurrence, the DT-EWMA esti-mator will trigger the dynamic-tuning capability and adjust the control parameter. The current process mean is then estimated using the newly adjusted control para-meter. Notice that the SPC charts are used for two pur-poses; first it is used to detect if a shift has occurred; secondly it estimates how many runs ahead (when) the shift occurs.
2.2. Dynamic-tuning algorithm
Suppose that a sudden shift occurs at t; that is, lt ¼ l1 for t < t and lt ¼ l2 for t t. In order to swiftly compensate for the sudden shift, W of the EWMA esti-mator in Equation (1) should be set larger instantly after the change point to capture the shift size. However, the step change occurs only once and the process mean re-mains unchanged after t. If the value of W remains large, the estimator becomes oversensitive to the white noises. To overcome the trade-off between large and small W and
to design a more effective estimator for processes subject to sudden shifts, we propose an estimator with a dy-namically adjusted control parameter:
^ l
lt¼ WtYtþ ð1 WtÞ^llt1: ð2Þ In Equation (2), the values of Wt are dynamically ad-justed over time. The task now is to derive a method for adjusting Wt. Now, suppose the change point is t and we do not have prior knowledge on how the process will shift away. In this case, we rely only on the sample data re-cently obtained from the process to estimate the shift size. In the Appendix, a theorem is provided to determine the optimal Wt that minimizes the Mean Squared Deviation (MSDt) at run t. MSDt is defined as follows:
MSDt¼ E½ð^llt ltÞ 2
: ð3Þ
The following results for t t are obtained: Wtopt¼ 1 t tþ 1: ð4Þ MSDoptt ¼ 1 t tþ 1r 2 : ð5Þ
In Equation (4), Wt decreases over time and diminishes to zero if there is no additional shift. If the change point can be detected exactly, at the detecting point (t¼ t), the optimal initial weight will start with a full weight, i.e., a value of one. Now suppose the process is not only subject
to shifts but also subject to other disturbances such as drifts. Therefore, a minimum value of Wt is required to keep the estimator working for such disturbances. The resulting estimation equation becomes:
^ l
lt¼ WYtþ ð1 W Þ^llt 1 for t t; ð6Þ
W ¼
maxðWt; WbÞ if dynamic-tuning loop is not completed, Wb if dynamic-tuning loop is completed. 8 > < > :
where Wb is the required minimum value of the control parameter.
In reality, the change point t is unknown. As men-tioned earlier, control charts are used to detect the shift occurrence in our approach. Since the signaling of a control chart only takes place some time after the actual change point, the time delay can be roughly estimated by the average run length of the control chart. Suppose now that the control chart detects the shift at the kth run (Fig. 2). Then the change point can be estimated to be:
^tt¼ k d þ 1; ð7Þ
where d is control chart run length. Equations (4) and (6) can be then applied. The complete algorithm of the DT-EWMA estimator is illustrated in Fig. 3.
From Equations (4) and (7), the control parameter of the DT-EWMA estimator has a starting value at the detection point: Wk¼ 1 k ^ttþ 1¼ 1 d: ð8Þ
That is, the estimator will take into consideration d data points prior to the detection point k while estimating lt. Notice that the starting weight value is 1=d, instead of one, to take into account the data points between the detection point and the change point to provide a more accurate estimation of current mean. The control para-meter continues to diminish over time until it reaches Wb and is fixed at that minimum required level afterwards. This control cycle holds on until the control chart signals again.
Notice that the values of d depend on the design of the control chart. Various Average Run Lengths (ARLs) for different shift sizes are reported in the literature by Champ and Woodhall (1987), Crowder (1987) and also Lucas and Saccucci (1990). In this paper, two EWMA control charts were used for the shift detection in the
following simulation case study. Their design parameters and ARLs which were estimated based on the Markovian approximation (Lucas and Saccucci, 1990.) are shown in Table 1. For EWMA 1 control chart, d is roughly equal to two at the detection of a shift size of 3r. For EWMA 2 control chart, d is roughly equal to four at the detection of a shift size of 2r. Figure 4 shows the potential changes of the control parameter over time under these control chart designs.
3. Process mean estimation performance analysis
In this section, we analyze the estimation performance among the EWMA, AEW (Yaschin, 1995), the adaptive EWMA of Patel and Jenkins (2000) and the DT-EWMA estimator using simulations. In these simulations, we consider a process with mean lt¼ l1 and standard de-viation r for t < t, and l
t ¼ l1þ sr for t t (i.e., a sudden shift occurs at t). The estimation performance is then characterized by a performance index, normalized Mean Square Estimation Error (MSEE). The smaller the value, the better the estimation performance. The nor-malized mean square estimation error is denoted as MSEE=r2 and defined as:
MSEEr2¼ 1 n Xn j¼1 ð^lltþj1 ltþj1Þ2 ( ) r2: ð9Þ
Here the MSEE=r2 index is calculated based on the simulation results of 100 runs and 2000 simulation cycles. For each simulation cycle, the shift occurs at the 51st run so t¼ 51 and n ¼ 50.
Figure 5 shows the results of the simulations (MSEE=r2 versus s). Notice that we have chosen each method’s parameters (Wbfor DT-EWMA, r and h for AEW, d and e for Patel and Jenkins) so that they all have roughly the same MSEE=r2when the shift size is zero. As can be seen, the adaptive EWMA estimator of Patel and Jenkins performs the worst under different shift sizes. The other three estimators perform roughly the same when the magnitude of the shift is smaller than 2r. However, both the AEW and DT-EWMA estimators perform better than the EWMA estimator when the magnitude of the shift is larger than 2r. That is, both the AEW and DT-EWMA estimators have a more stable performance un-der different magnitudes of the shifts. This is due to the fact that the AEW estimator has an ‘‘adaptive capability’’ while the DT-EWMA estimator has a ‘‘dynamic-tuning capability’’ once a shift occurs. The adaptive algorithm in the work of Patel and Jenkins (2000) is a very conserva-tive approach in updating the control parameter once a shift occurs. Its adaptive performance is not very effec-tive.
During the simulations, we also found that the algorithm for estimating the change point in the AEW
Fig. 2. The run length for the control chart to detect a process shift.
estimator requires an intensive computational effort. In the DT-EWMA estimator, the location of the change point is
estimated by the ARL of the adopted control charts. Although the search of the change point by the AEW method might be more rigorous than by the DT-EWMA estimator, the performance is roughly the same as indi-cated in Fig. 5. For the computational complexity com-parison, suppose the current process run is the ith run. In the worst case when there is no shift occurrence till the ith run, we have estimated that the required number of computations by the AEW method at the ith run is roughly equal to iði 1Þ=2. In the DT-EWMA estimator, if two control charts are used for shift detection, the number of computations is only two at each run. Suppose there is no shift until run 200, the total number of cal-culations from run 1 to run 200 by the AEW method will
Fig. 3. The complete algorithm of the DT-EWMA estimator.
Table 1. Average run length for two EWMA control charts Shift size ARL EWMA_1 EWMA_2
Wl¼ 0:6 CL¼ 3:285rl Wm¼ 0:33 CL¼ 3:25rm 0r ARL0 1003.44 1007.44 1r ARL1 31.55 15.87 2r ARL1 4.55 3.90 3r ARL1 2.06 2.27 4r ARL1 1.36 1.67
be 1333 300. The number of calculations using the DT-EWMA method is only 400. There is a significant dif-ference. Figure 6 also shows the results for DT-EWMA estimator under different shift magnitudes and Wb. As can be seen, the performance is not very sensitive to the se-lection of Wb.
4. Process adjustment simulations
4.1. Process model and disturbance
In many manufacturing processes, when the process output deviation can be adjusted to bring the process output to the target value, the process mean estimator also provides information on the magnitude of the re-quired adjustment. In this section, we conduct Monte Carlo simulations to study the adjustment performance based on the EWMA and DT-EWMA estimators.
First we need to define the process model and process disturbance. Suppose the underlying process model and disturbance can be approximated as:
Yt¼ a þ bXtþ etþ dt: ð10Þ
et N ð0; r2Þ: ð11Þ
dt ¼ crt þ st: ð12Þ
Here a is the process constant, b is the process gain, X is the input setting, e is the white noise and d is the special disturbance. Notice that in Equation (12) two types of disturbance are simulated: process drift and shift. The process drift follows a linear drift pattern and deviates from the target (T) at the speed of cr. The shift is mod-eled as a random shift with the following pattern:
st¼
st1 with probability 1p, st1þ S with probability p,
ð13Þ
where p is the probability of the shift occurrence. The shift magnitude S is a random variable and S N ðlS;r2Þ. In the Monte Carlo simulations, the values of the above parameters are chosen as: a¼ 3000, b ¼ 10, r¼ 50, c ¼ 0:02, and T ¼ 3500. As for the random shift, it is simulated in four cases and the occurrence proba-bility (p) and magnitude (S) are summarized in Table 2. For example, in the first case the random shift has a low probability of occurrence and a smaller shift magnitude.
Fig. 4. Dynamic and minimum control parameters.
In the fourth case, the random shift has a high probability of occurrence and a larger shift magnitude.
In the process adjustment simulations, once the current run’s mean is estimated, the adjustment for the next run is achieved by tuning the input setting as:
Xtþ1 Xt ¼ T ^llt
b : ð14Þ
As for the adjustment performance characterization, we use the normalized mean square error as the perfor-mance index. The smaller the value, the better the per-formance. The normalized mean square error is denoted as MSE=r2 and is defined as:
MSEr2¼ 1 n Xn t¼1 ðYt T Þ2 ( ) r2: ð15Þ
Here the MSE=r2 index is calculated based on the simulation results of 2000 runs and 200 simulation cycles.
4.2. Monte Carlo simulation results
Simulation results for Case 3 and Case 4 are shown in Figs. 7 and 8. Figure 9 shows the results of the average of Case 1 to Case 4 and represents the situation when the disturbance information is unknown and all types of disturbance are possible.
Based on these results, we have the following conclu-sions:
• When there are only small process drifts or small shifts, as happened in Cases 1 and 3, the process
adjustment based on the EWMA estimator is very effective in compensating for these small distur-bances. As a result, the DT-EWMA estimator has a similar adjustment performance as the EWMA esti-mator.
• When there are large process shifts, as happened in Cases 2 and 4, the adjustment performance based on the DT-EWMA estimator is much better than that of the EWMA estimator with fixed control parameter. The performance of the DT-EWMA estimator is not very dependent on the minimum control parameter. • The DT-EWMA estimator is very effective at com-pensating for disturbances of random shifts and linear drifts.
5. CMP process case study
In this section, we apply both EWMA and DT-EWMA estimators to an industrial process and compare their estimation performance. The selected industrial process, Chemical-Mechanical Polishing (CMP) process (Fig. 10), is considered to be a strategically important technology for IC fabrication. Because the process includes me-chanical abrasion of the surface, the polishing pad wears rapidly, which results in a decreasing and irregular run-to-run process removal rate.
A simple strategy for controlling the CMP process is to estimate the run-to-run process removal rate and then adjust the polishing time based on the estimated removal rate (Boning et al., 1996). The estimation, however, is not
Fig. 6. Process mean estimation performance of the DT-EWMA estimator under different shift magnitudes and Wb.
Table 2. Four cases modelled in our simulations
Case 1 Case 2 Case 3 Case 4
P 0.005 0.005 0.05 0.05
easy for a field process engineer. Figure 11 shows data on the removal rate sampled from a local wafer fab. It can be observed that there are several shifts of different magni-tudes (because of the replacement of the pads) and a decreasing trend of removal rate (because of the wear-out of the pads). The pattern of the field data is actually more complicated than the process model we assume in the previous sections. Nevertheless, they are still valuable to our study.
Since the replacement of the pad usually results in a shift, we have used this prior information in our case study. In applying the DT-EWMA and EWMA to the
CMP data, we reset the estimation right after the pad replacement so that the estimation can be more accurate. Figure 12 shows the results of MSEE versus Wb for EWMA and DT-EWMA estimators. The fact that the minimum MSEE occurs at around Wb¼ 0:8 indicates that the underlying process mean is changing very frequently and a higher value of W is required. Although both the EWMA and DT-EWMA estimators have the same esti-mation performance at around Wb¼ 0:8, the DT-EWMA estimator is less sensitive to the selection of Wb. In real practice, the value of the control parameter is usually se-lected to be less than 0.5. As a result, such a practice will certainly fail using the conventional EWMA estimator. On the other hand, the DT-EWMA estimator appears to be quite robust and the selection of Wb, even for values less than 0.5, will not cause a significant increase in MSEE.
6. Conclusions
In this paper, we present an EWMA-based mean estima-tor with dynamic-tuning capability for processes subject to random shifts of different magnitudes. The DT-EWMA estimator uses the criteria of minimum mean square de-viation to derive the optimal control parameter in re-sponse to the random shift disturbances. Through the simulation and case study, we conclude that:
• The DT-EWMA estimator, similar to the AEW es-timator, has a more stable estimation performance under different magnitudes of random shifts than the other approaches. A possible extension of this ap-proach could be the determination of the change point by a more rigorous method instead of the control chart and the ARL approach if the compu-tation effort is not an issue.
• The dynamic-tuning capability of the DT-EWMA estimator is analytically derived and is very easy to implement.
• The process adjustment that is based on the DT-EWMA estimator is very effective for processes
Fig. 7. Simulation results for case 3.
Fig. 8. Simulation results for case 4.
Fig. 9. Average results of cases 1, 2, 3, and 4.
subject to different magnitudes of random shifts and linear drifts.
Acknowledgement
The authors would like to thank Dr. Argon Chen of the National Taiwan University for his valuable comments during the preparation of this paper. The work was supported by the National Science Council of the Re-public of China under Grant NSC88-2416-H-002-055 and NSC89-2416-H-002-044.
References
Boning, D., Moyne, W., Smith, T., Moyne, J., Hurwitz, A., Shellman, S., Moyne, J., Trefeyan, R. and Taylor, J. (1996) Run by run
process control of chemical – mechanical polishing. IEEE Transactions on Components, Packaging, and Manufacturing Technology – Part C, 19(4), 307–313.
Box, G. and Jenkins, M. (1976) Time Series Analysis – Forecasting and Control, Holden-Day, San Francisco, CA.
Box, G. and Kramer, T. (1992) Statistical process monitoring and feedback adjustment – a discussion. Technometrics, 34(3), 251–267. Butler, S. and Stefani, J. (1994) Supervisory run-to-run control of polysilicon gate etch using in situ ellipsometry. IEEE Transactions on Semiconductor Manufacturing, 7(2), 193–201.
Champ, C.W. and Woodall, W.H. (1987) Exact result for Shewhart control charts with supplementary runs rules. Technometrics, 29(4), 393–399.
Crowder, S.V. (1987) A simple method for studying run-length distri-butions of exponentially weighted moving average charts. Tech-nometrics, 29(4), 401–407.
Del Castillo, E. (1999) Long run and transient analysis of a double EWMA feedback controller. IIE Transactions, 31(12), 1157–1169. Del Castillo, E. and Hurwitz, A.N. (1997) Run-to-run process control: literature review and extensions. Journal of Quality Technology, 29(2), 184–186.
Fig. 11. Removal rate data of the CMP process.
Ingolfsson, A. and Sachs, E. (1993) Stability and sensitivity of an EWMA controller. Journal of Quality Technology, 25(4), 271–287. Lucas, J.M. and Saccucci, M.S. (1990) Exponentially weighted moving average control schemes: properties and enhancements. Techno-metrics, 32(1), 1–29.
MacGregor, J.F. (1988) On line statistical process control. Chemical Engineering Progress, 84, 21–31.
Patel, N.S. and Jenkins, S.T. (2000) Adaptive optimization of run-to-run controllers: the EWMA example. IEEE Transactions on Semiconductor Manufacturing, 13(1), 97–107.
Sachs, E., Hu, A. and Ingolfsson, A. (1995) Run by run process con-trol: combining SPC and feedback control. IEEE Transactions on Semiconductor Manufacturing, 8(1), 26–43.
Smith, T., Boning, D., Stefani, J. and Butler, S. (1998) Run by run advanced process control of metal sputter deposition. IEEE Transactions on Semiconductor Manufacturing, 11(2), 276–284. Yashchin, E. (1995) Estimating the current mean of a process subject to
abrupt changes. Technometrics, 37(3), 311–323.
Appendix
Theorem 1. Given a process model, suppose that ltchanges once somewhere between time t 1 and t; that is, l
t ¼ l1 for t < t and lt¼ l2 for t t but we do not have any prior information on the size of the l2. Then the EWMA estimator with a dynamically adjusted control parameter in Equation(2) minimizes MSDt if
Wt¼ 1 t tþ 1;
for all t t and the observation at the time t is known. Proof. Given the weight
Wt¼ 1 t tþ 1; for any t t < s, ^ l ls¼ WsYsþ ð1 WsÞ^lls1; ¼ WsYsþ ð1 WsÞ½WsYs1þ ð1 WsÞ^lls2 : Substituting Ws1¼ 1 s t; into the above equation, we obtain
^ l ls¼ Wsðl2þ esÞ þ ð1 WsÞ 1 s tðl2þ es1Þ þ s t 1 s t ll^s2 ; ¼ Wsðl2þ esÞ þ ð1 WsÞ 1 s tðl2þ es1Þ þ s t 1 s t ½Ws2Ys2þ ð1 Ws2Þ^lls3 : Again, substituting Ws2¼ 1 s t 1; into the above equation, we obtain
^ l ls¼ Wsðl2þ esÞ þ ð1 WsÞ 1 s tðl2þ es1Þ þ s t 1 s t 1 s t 1ðl2þ es2Þ þs t 2 s t 1^lls3 : ¼ Wsðl2þ esÞ þ ð1 WsÞ 1 s tðl2þ es1Þ þs t 1 s t 1 s t 1ðl2þ es2Þ þs t 2 s t 1ðWs3Ys3þ ð1 Ws3Þ^lls4Þ : ¼ Wsðl2þ esÞ þ ð1 WsÞ 2 s tðl2þ es1Þ þ Wses þð1 WsÞ 1 s t Xs1 i¼s2 ei þ s t 2 s t 1ðWs3Ys3 þð1 Ws3Þ^lls4Þ : By recursively substituting Ws3¼ 1 s t 2; . . . ; Wt ¼ 1; into the equation, we finally have
^ l ls¼ l2þ Wsesþ ð1 WsÞ 1 s t Xs1 i¼t ei: The MSD at s is: MSDs¼ E½ð^lls lsÞ 2 ; ¼ E Wsesþ ð1 WsÞ 1 s t Xs1 i¼t ei !2 2 4 3 5; ¼ Ws2r 2 þ ð1 WsÞ2 1 s tr 2 :
By minimizing MSDs with respect to Ws, we obtain the optimal control parameter at s.
Wsopt¼ 1 s tþ 1; and the optimal MSD
MSDopts ¼ 1 s tþ 1r
2: We have shown that given
Wt ¼ 1 t tþ 1 for t t < s; Wsopt¼ 1 s tþ 1: That is, if
Wtopt ¼
1
t tþ 1¼ 1; then recursively we can obtain
Wtopt¼ 1 t tþ 1 for8 t t:
Now, we show that Wtopt¼ 1. Because
^ l
lt¼ WtYtþ ð1 WtÞ^llt1;
¼ Ytþ ð1 WtÞð^llt1 YtÞ:
Given the output at time t is observed to be yt
yt ¼ ltþ et; the MSD at t is: MSDt ¼ E½ð^llt ltÞ2jYt ¼ yt ; ¼ Ef½Ytþ ð1 WtÞð^llt1 YtÞ ðyt etÞ 2jYt¼ ytg; ¼ Ef½ð1 WtÞð^llt1 YtÞ þ etÞ 2g; ¼ ð1 WtÞ2ð^llt1 ytÞ2þ r2
Minimizing MSDt with respect to Wt, we obtain optimal
weight Wtopt ¼ 1 and MSDoptt ¼ r2.
Therefore, Wtopt¼ 1 t tþ 1: and MSDopttþ1¼ 1 t tþ 1r 2 for 8 t t: j Biographies
Ruey-Shan Guo received his Ph.D. degree in Mechanical Engineer-ing from MIT in 1991 and an MBA degree from San Jose State University in 1994. From 1991 to 1995, he was a Senior Process Engineer at National Semiconductor Fairchild Research Center, Santa Clara, CA. From 1995 till now, he has been an Associate Professor in the Graduate Institute of Business Administration and Graduate Institute of Industrial Engineering of National Taiwan University. During his current position, he has been a principal in-vestigator and consultant to many industry- and government-funded projects. His current research interests include advanced quality control, run-to-run process control, CIM, and supply chain man-agement. He won the excellent teaching award from National Tai-wan University twice and the best annual journal paper from the Chinese Management Academics Association. Dr. Guo is a member of IEEE.
Jin-Jung Chen received a Ph.D. degree in Mechanical Engineering from the National Taiwan University in 2000. He is currently a man-ager in the Product Division of the Mustang Industrial Corporation in Taiwan.
Contributed by the Design of Experiments and Robust Designs Depart-ment