Chapter Contents 
Previous 
Next 
OUTPUT Statement 
See CPCPMK in the SAS/QC Sample Library 
In recent years, a number of process capability indices that have been proposed in the research literature are gradually being introduced in applications. As shown in this example, you can compute such indices in the DATA step after using the OUTPUT statement in the CAPABILITY procedure to save various summary statistics.
Hardness measurements (in scaled units) for 50 titanium samples are saved as values of the variable HARDNESS in the following SAS data set:
data titanium; label hardness = 'Hardness Measurement'; input hardness @@; datalines; 1.38 1.49 1.43 1.60 1.59 1.34 1.44 1.64 1.83 1.57 1.45 1.74 1.61 1.39 1.63 1.73 1.61 1.35 1.51 1.47 1.46 1.41 1.56 1.40 1.58 1.43 1.53 1.53 1.58 1.62 1.58 1.46 1.26 1.57 1.41 1.53 1.36 1.63 1.36 1.66 1.49 1.55 1.67 1.41 1.39 1.75 1.37 1.36 1.86 1.49 ;
The target value for hardness is 1.6, and the lower and upper specification limits are 0.8 and 2.4, respectively. The samples are produced by an incontrol process, and the measurements are assumed to be normally distributed.
The following statements use the OUTPUT statement to save various descriptive statistics and an estimate of the index C_{pm} in a data set named INDICES:
proc capability data=titanium noprint; var hardness; specs lsl=0.8 target=1.6 usl=2.4; output out=indices n = n mean = avg std = std var = var lsl = lsl target = t usl = usl pnormal = pnormal cpm = cpm ; run;
In addition to C_{pm}, you want to report an estimate for the index C_{pmk}, which is defined as follows:
The following statements compute this estimate:
data indices; set indices; d = 0.5*( USL  LSL ); m = 0.5*( USL + LSL ); num = d  abs( avg  m ); den = 3 * sqrt( (n1)*var/n + (avgt)*(avgt) ); cpmk = num/den; title 'Capability Analysis of Titanium Hardness'; proc print data=indices noobs; var n avg std lsl t usl cpm cpmk pnormal; run;
The results are listed in Output 7.1.1.
Output 7.1.1: Computation of C_{pmk}

Note that the pvalue for the KolmogorovSmirnov test of normality is 0.27693, indicating that the assumption of normality is justified.
The following statements also compute an estimate of the index C_{pm} using the SPECIALINDICES option:
proc capability data=titanium specialindices; var hardness; specs lsl=0.8 target=1.6 usl=2.4; run;Output 7.1.2: Computation of C_{pmk} using the SPECIALINDICES option
Chapter Contents 
Previous 
Next 
Top 
Copyright © 1999 by SAS Institute Inc., Cary, NC, USA. All rights reserved.