Chapter Contents
Chapter Contents
Previous
Previous
Next
Next
Specialized Control Charts

Analyzing the Difference from Nominal

The following example* is adapted from an application in aircraft component manufacturing. A metal extrusion process is used to make three slightly different models of the same component. The three product types (labeled M1, M2, and M3) are produced in small quantities because the process is expensive and time-consuming.

Figure 49.14 shows the structure of a SAS data set named OLD, which contains the diameter measurements for various short runs. Samples 1 to 30 are to be used to estimate the process standard deviation \sigma for the differences from nominal.

 
sample prodtype diameter
1 M3 13.99
2 M3 14.69
3 M3 13.86
4 M3 14.32
5 M3 13.23
6 M1 17.55
7 M1 14.26
8 M1 14.62
9 M1 12.97
10 M2 16.18
11 M2 15.29
12 M2 16.20
13 M3 13.89
14 M3 12.71
15 M3 14.32
16 M3 15.35
17 M2 15.08
18 M2 14.72
19 M2 14.79
20 M2 15.27
21 M2 15.95
22 M1 14.78
23 M1 15.19
24 M1 15.41
25 M1 16.26
26 M3 16.68
27 M3 15.60
28 M3 14.86
29 M3 16.67
30 M3 14.35
Figure 49.14: Diameter Measurements in the Data Set OLD

In short run applications involving many product types, it is common practice to maintain a database for the nominal values for the product types. Here, the nominal values are saved in a SAS data set named NOMVAL, which is listed in Figure 49.15.

 
prodtype nominal
M1 15.0
M2 15.5
M3 14.8
M4 15.2
Figure 49.15: Nominal Values for Product Types in the Data Set NOMVAL

To compute the differences from nominal, you must merge the data with the nominal values. You can do this with the following SAS statements. Note that an IN= variable is used in the MERGE statement to allow for the fact that NOMVAL includes nominal values for product types that are not represented in OLD. Figure 49.16 lists the merged data set OLD.

   proc sort data=old;
      by prodtype;
   data old;
      format diff 5.2 ;
      merge nomval old(in = a);
         by prodtype;
      if a;
      diff = diameter - nominal;
   proc sort data=old;
      by sample;
   run;

 
sample prodtype diameter nominal diff
1 M3 13.99 14.8 -0.81
2 M3 14.69 14.8 -0.11
3 M3 13.86 14.8 -0.94
4 M3 14.32 14.8 -0.48
5 M3 13.23 14.8 -1.57
6 M1 17.55 15.0 2.55
7 M1 14.26 15.0 -0.74
8 M1 14.62 15.0 -0.38
9 M1 12.97 15.0 -2.03
10 M2 16.18 15.5 0.68
11 M2 15.29 15.5 -0.21
12 M2 16.20 15.5 0.70
13 M3 13.89 14.8 -0.91
14 M3 12.71 14.8 -2.09
15 M3 14.32 14.8 -0.48
16 M3 15.35 14.8 0.55
17 M2 15.08 15.5 -0.42
18 M2 14.72 15.5 -0.78
19 M2 14.79 15.5 -0.71
20 M2 15.27 15.5 -0.23
21 M2 15.95 15.5 0.45
22 M1 14.78 15.0 -0.22
23 M1 15.19 15.0 0.19
24 M1 15.41 15.0 0.41
25 M1 16.26 15.0 1.26
26 M3 16.68 14.8 1.88
27 M3 15.60 14.8 0.80
28 M3 14.86 14.8 0.06
29 M3 16.67 14.8 1.87
30 M3 14.35 14.8 -0.45
Figure 49.16: Data Merged with Nominal Values

Assume that the variability in the process is constant across product types. To estimate the common process standard deviation \sigma, you first estimate \sigma for each product type based on the average of the moving ranges of the differences from nominal. You can do this in several steps, the first of which is to sort the data and compute the average moving range with the SHEWHART procedure.

   proc sort data=old;
      by prodtype;

   proc shewhart data=old;
      irchart diff*sample /
         nochart
         outlimits=baselim;
      by prodtype;
   run;
The purpose of this procedure step is simply to save the average moving range for each product type in the OUTLIMITS= data set BASELIM, which is listed in Figure 49.17 (note that PRODTYPE is specified as a BY variable).

 
Control Limits By Product Type

prodtype _VAR_ _SUBGRP_ _TYPE_ _LIMITN_ _ALPHA_ _SIGMAS_ _LCLI_ _MEAN_ _UCLI_ _LCLR_ _R_ _UCLR_ _STDDEV_
M1 diff sample ESTIMATE 2 .002699796 3 -3.13258 0.13000 3.39258 0 1.22714 4.00850 1.08753
M2 diff sample ESTIMATE 2 .002699796 3 -1.77795 -0.06500 1.64795 0 0.64429 2.10458 0.57098
M3 diff sample ESTIMATE 2 .002699796 3 -3.22641 -0.19143 2.84356 0 1.14154 3.72887 1.01166
Figure 49.17: Values of \bar{R} by Product Type

To obtain a combined estimate of \sigma, you can use the MEANS procedure to average the average ranges in BASELIM and then divide by the unbiasing constant d2.

   proc means data=baselim noprint;
      var _r_;
      output out=difflim (keep=_r_) mean=_r_;

   data difflim;
      set difflim;
      drop _r_;
      length _var_ _subgrp_ $ 8;
      _var_    = 'diff';
      _subgrp_ = 'sample';
      _mean_   = 0.0;
      _stddev_ = _r_ / d2(2);
      _limitn_ = 2;
      _sigmas_ = 3;
   run;
The data set DIFFLIM is structured for subsequent use by the SHEWHART procedure as an input LIMITS= data set. The variables in a LIMITS= data set provide pre-computed control limits or -as in this case - the parameters from which control limits are to be computed. These variables have reserved names that begin and end with the underscore character. Here, the variable _STDDEV_ saves the estimate of \sigma, and the variable _MEAN_ saves the mean of the differences from nominal. Recall that this mean is zero, since the nominal values are assumed to represent the process mean for each product type. The identifier variables _VAR_ and _SUBGRP_ record the names of the process and subgroup variables (these variables are critical in applications involving many product types). The variable _LIMITN_ is assigned a value of 2 to specify moving ranges of two consecutive measurements, and the variable _SIGMAS_ is assigned a value of 3 to specify 3\sigma limits. The data set DIFFLIM is listed in Figure 49.18.

 
Control Limit Parameters For Differences

_var_ _subgrp_ _mean_ _stddev_ _limitn_ _sigmas_
diff sample 0 0.89006 2 3
Figure 49.18: Estimates of Mean and Standard Deviation

Now that the control limit parameters are saved in DIFFLIM, diameters for an additional 30 parts (samples 31 to 60) are measured and saved in a SAS data set named NEW. You can construct short run control charts for this data by merging the measurements in NEW with the corresponding nominal values in NOMVAL, computing the differences from nominal, and then contructing the short run individual measurements and moving range charts.

   proc sort data=new;
      by prodtype;

   data new;
      format diff 5.2 ;
      merge nomval new(in = a);
         by prodtype;
      if a;
      diff = diameter - nominal;
      label sample   = 'Sample Number'
            prodtype = 'Model';

   proc sort data=new;
      by sample;
   symbol1 v=dot    c=yellow;
   symbol2 v=plus   c=red;
   symbol3 v=circle c=green;
   title 'Chart for Difference from Nominal';
   proc shewhart data=new limits=difflim;
      irchart diff*sample=prodtype / split    = '/'
                                     cframe   = vigb
                                     cinfill  = vlib
                                     cconnect = yellow;
      label diff = 'Difference/Moving Range';
   run;
The chart is displayed in Figure 49.19. Note that the product types are identified with symbol markers as requested by specifying PRODTYPE as a symbol-variable.

sshort6.gif (7895 bytes)

Figure 49.19: Short Run Control Chart

You can also identify the product types with a legend by specifying PRODTYPE as a _PHASE_ variable.

   symbol v=dot c=yellow;
   title 'Chart for Difference from Nominal';
   proc shewhart data=new (rename=(prodtype=_phase_))
      limits=difflim;
      irchart diff*sample /
         readphases = all
         phaseref
         phasebreak
         phaselegend
         split      = '/'
         cframe     = vigb
         cconnect   = yellow
         cinfill    = vlib;
      label diff = 'Difference/Moving Range';
   run;
The display is shown in Figure 49.20. Note that the PHASEBREAK option is used to suppress the connection of adjacent points in different phases (product types).

sshort7.gif (8048 bytes)

Figure 49.20: Identification of Product Types

In some applications, it may be useful to replace the moving range chart with a plot of the nominal values. You can do this with the TRENDVAR= option in the XCHART statement* provided that you reset the value of _LIMITN_ to 1 to specify a subgroup sample of size one.

   data difflim;
      set difflim;
      _var_    = 'diameter';
      _limitn_ = 1;
   symbol v=dot c=yellow;
   title 'Differences and Nominal Values';
   proc shewhart data=new limits=difflim;
      xchart diameter*sample (prodtype) /
         nolimitslegend
         nolegend
         split         = '/'
         blockpos      = 3
         blocklabtype  = scaled
         blocklabelpos = left
         xsymbol       = xbar
         trendvar      = nominal
         cframe        = vigb
         cinfill       = vlib
         cconect       = yellow;
      label diameter   = 'Difference/Nominal'
            prodtype   = 'Product';
   run;

The display is shown in Figure 49.22. Note that you identify the product types by specifying PRODTYPE as a block variable enclosed in parentheses after the subgroup variable SAMPLE. The BLOCKLABTYPE= option specifies that values of the block variable are to be scaled (if necessary) to fit the space available in the block legend. The BLOCKLABELPOS= option specifies that the label of the block variable is to be displayed to the left of the block legend.

sshort8.gif (7911 bytes)

Figure 49.21: Short Run Control Chart with Nominal Values

Chapter Contents
Chapter Contents
Previous
Previous
Next
Next
Top
Top

Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.