Example 15.41 Specifying Probability Limits

[See SHWXS2 in the SAS/QC Sample Library]This example illustrates how to create and charts with probability limits. The following statements read the kilowatt power output measurements from the data set Turbine (see Creating Charts for Means and Standard Deviations from Raw Data) and create the and charts shown in Output 15.41.1:

ods graphics off;
symbol v=dot h=.8;
title 'Mean and Standard Deviation Charts With Probability Limits';
proc shewhart data=Turbine;
   xschart KWatts*Day  / alpha     = 0.01
                         outlimits = Oillim;
run;

The ALPHA= option specifies the probability () that a subgroup summary statistic is outside the limits. Here, the limits are computed so that the probability that a subgroup mean or standard deviation is less than its lower limit is , and the probability that a subgroup mean or standard deviation is greater than its upper limit is . This assumes that the measurements are normally distributed.

The OUTLIMITS= option names an output data set (Oillim) that saves the probability limits. The data set Oillim is shown in Output 15.41.2.

Output 15.41.1 Probability Limits on and Charts
Probability Limits on X and s Charts

Output 15.41.2 Probability Limit Information
Mean and Standard Deviation Charts with Probability Limits

_VAR_ _SUBGRP_ _TYPE_ _LIMITN_ _ALPHA_ _SIGMAS_ _LCLX_ _MEAN_ _UCLX_ _LCLS_ _S_ _UCLS_ _STDDEV_
KWatts Day ESTIMATE 20 0.01 2.57583 3370.79 3485.41 3600.03 119.432 196.396 283.570 198.996

The variable _ALPHA_ saves the value of . The value of the variable _SIGMAS_ is computed as , where is the inverse standard normal distribution function. Note that, in this case, the probability limits for the mean are equivalent to limits.

Since all the points fall within the probability limits, it can be concluded that the process is in statistical control.