EWMACHART Statement: MACONTROL Procedure

Plotted Points

Each point on the chart indicates the value of the exponentially weighted moving average (EWMA) for that subgroup. The EWMA for the th subgroup () is defined recursively as

where is a weight parameter . Some authors (for example, Hunter 1986 and Crowder 1987a,b) use the symbol instead of for the weight. You can specify the weight with the WEIGHT= option in the EWMACHART statement or with the variable _WEIGHT_ in a LIMITS= data set. If you specify a known value () for , ; otherwise, .

The preceding equation can be rewritten as

which expresses the current EWMA as the previous EWMA plus the weighted error in the prediction of the current mean based on the previous EWMA.

The EWMA for the th subgroup can also be written as

which expresses the EWMA as a weighted average of past subgroup means, where the weights decline exponentially, and the heaviest weight is assigned to the most recent subgroup mean.

Central Line

By default, the central line on an EWMA chart indicates an estimate for , which is computed as

If you specify a known value () for , the central line indicates the value of .

Control Limits

You can compute the limits in the following ways:

• as a specified multiple () of the standard error of above and below the central line. The default limits are computed with (these are referred to as limits).

• as probability limits defined in terms of , a specified probability that exceeds the limits

The following table presents the formulas for the limits:

These formulas assume that the data are normally distributed. If standard values and are available for and , respectively, replace with and with in Table 9.3. Note that the limits vary with both and .

If the subgroup sample sizes are constant (), the formulas for the control limits simplify to

Consequently, when the subgroup sample sizes are constant, the width of the control limits increases monotonically with . For probability limits, replace with in the previous equations. Refer to Roberts (1959) and Montgomery (1996).

As becomes large, the upper and lower control limits approach constant values:

Some authors base the control limits for EWMA charts on the asymptotic expressions in the two previous equations. For asymptotic probability limits, replace with in these equations. You can display asymptotic limits by specifying the ASYMPTOTIC option.

Uniformly weighted moving average charts and exponentially weighted moving average charts have similar properties, and their asymptotic control limits are identical provided that

where is the weight factor for uniformly weighted moving average charts. Refer to Wadsworth, Stephens, and Godfrey (1986) and the American Society for Quality Control (1983).

You can specify parameters for the EWMA limits as follows:

• Specify with the SIGMAS= option or with the variable _SIGMAS_ in a LIMITS= data set.

• Specify with the ALPHA= option or with the variable _ALPHA_ in a LIMITS= data set.

• Specify a constant nominal sample size for the control limits with the LIMITN= option or with the variable _LIMITN_ in a LIMITS= data set.

• Specify with the WEIGHT= option or with the variable _WEIGHT_ in a LIMITS= data set.

• Specify with the MU0= option or with the variable _MEAN_ in a LIMITS= data set.

• Specify with the SIGMA0= option or with the variable _STDDEV_ in a LIMITS= data set.

Choosing the Value of the Weight Parameter

Various approaches have been proposed for choosing the value of .

• Hunter (1986) states that the choice "can be left to the judgment of the quality control analyst" and points out that the smaller the value of , "the greater the influence of the historical data."

• Hunter (1986) also discusses a least squares procedure for estimating from the data, assuming an exponentially weighted moving average model for the data. In this context, the fitted EWMA model provides a forecast of the process that is the basis for dynamic process control. You can use the ARIMA procedure in SAS/ETS^® software to compute the least squares estimate of . (Refer to SAS/ETS User's Guide for information about PROC ARIMA.) Also see Autocorrelation in Process Data.

• A number of authors have studied the design of EWMA control schemes based on average run length (ARL) computations. The ARL is the expected number of points plotted before a shift is detected. Ideally, the ARL should be short when a shift occurs, and it should be long when there is no shift (the process is in control.) The effect of on the ARL was described by Roberts (1959), who used simulation methods. The ARL function was approximated and tabulated by Robinson and Ho (1978), and a more general method for studying run-length distributions of EWMA charts was given by Crowder (1987a, 1987b). Unlike Hunter (1986), these authors assume the data are independent and identically distributed; typically the normal distribution is assumed for the data, although the methods extend to nonnormal distributions. A more detailed discussion of the ARL approach follows.

Average run lengths for two-sided EWMA charts are shown in Table 9.4, which is patterned after Table 1 of Crowder (1987a, 1987b). The ARLs were computed using the EWMAARL DATA step function (see EWMAARL Function for details on the EWMAARL function). Note that Crowder (1987a, 1987b). uses the notation L in place of and the notation in place of .

You can use Table 9.4 to find a combination of and that yields a desired ARL for an in-control process () and for a specified shift of . Note that is assumed to be standardized; in other words, if a shift of is to be detected in the process mean , and if is the process standard deviation, you should select the table entry with

where is the subgroup sample size. Thus, can be regarded as the shift in the sampling distribution of the subgroup mean.

For example, suppose you want to construct an EWMA scheme with an in-control ARL of 90 and an ARL of 9 for detecting a shift of . Table 9.4 shows that the combination and yields an in-control ARL of 91.17 and an ARL of 8.27 for .

Crowder (1987a, 1987b) cautions that setting the in-control ARL at a desired level does not guarantee that the probability of an early false signal is acceptable. For further details concerning the distribution of the ARL, refer to Crowder (1987a, 1987b).

In addition to using Table 9.4 or the EWMAARL DATA step function to choose a EWMA scheme with desired average run length properties, you can use them to evaluate an existing EWMA scheme. For example, the "Getting Started" section of this chapter contains EWMA schemes with and . The following statements use the EWMAARL function to compute the in-control ARL and the ARLs for shifts of and :

data arlewma;
   arlin = ewmaarl(  0,0.3,3.0);
   arl1  = ewmaarl(.25,0.3,3.0);
   arl2  = ewmaarl(.50,0.3,3.0);
run;

The in-control ARL is 465.553, the ARL for is 178.741, and the ARL for is 53.1603. See Example 9.5 for an illustration of how to use the EWMAARL function to compute average run lengths for various EWMA schemes and shifts.