Median Chart Example 2
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: SHWMEX2 */
/* TITLE: Median Chart Example 2 */
/* PRODUCT: QC */
/* SYSTEM: ALL */
/* KEYS: Shewhart Charts, Median Charts, */
/* PROCS: SHEWHART */
/* DATA: */
/* */
/* REF: SAS/QC Software: Usage and Reference, Version 6, */
/* First Edition, Volume 1 and Volume 2 */
/* */
/****************************************************************/
data Wire;
input Sample Size @;
do i=1 to Size;
input Breakstrength @@;
output;
end;
drop i Size;
label Breakstrength ='Breaking Strength (lb/in)'
Sample ='Sample Index';
datalines;
1 5 60.6 62.3 62.0 60.4 59.9
2 5 61.9 62.1 60.6 58.9 65.3
3 4 57.8 60.5 60.1 57.7
4 5 56.8 62.5 60.1 62.9 58.9
5 5 63.0 60.7 57.2 61.0 53.5
6 7 58.7 60.1 59.7 60.1 59.1 57.3 60.9
7 5 59.3 61.7 59.1 58.1 60.3
8 5 61.3 58.5 57.8 61.0 58.6
9 6 59.5 58.3 57.5 59.4 61.5 59.6
10 5 61.7 60.7 57.2 56.5 61.5
11 3 63.9 61.6 60.9
12 5 58.7 61.4 62.4 57.3 60.5
13 5 56.8 58.5 55.7 63.0 62.7
14 5 62.1 60.6 62.1 58.7 58.3
15 5 59.1 60.4 60.4 59.0 64.1
16 5 59.9 58.8 59.2 63.0 64.9
17 6 58.8 62.4 59.4 57.1 61.2 58.6
18 5 60.3 58.7 60.5 58.6 56.2
19 5 59.2 59.8 59.7 59.3 60.0
20 5 62.3 56.0 57.0 61.8 58.8
21 4 60.5 62.0 61.4 57.7
22 4 59.3 62.4 60.4 60.0
23 5 62.4 61.3 60.5 57.7 60.2
24 5 61.2 55.5 60.2 60.4 62.4
25 5 59.0 66.1 57.7 58.5 58.9
;
title 'Median Chart for Breaking Strengths';
ods graphics on;
proc shewhart data=Wire;
mchart Breakstrength*Sample / odstitle=title;
run;
title 'Estimates of the Process Standard Deviation';
proc shewhart data=Wire;
mchart Breakstrength*Sample / outlimits=Wirelim1
nochart outindex ='NOWEIGHT-Ranges';
mchart Breakstrength*Sample / outlimits=Wirelim2
stddeviations
nochart outindex ='NOWEIGHT-Stds';
mchart Breakstrength*Sample / outlimits=Wirelim3
smethod =mvlue
nochart outindex ='MVLUE -Ranges';
mchart Breakstrength*Sample / outlimits=Wirelim4
stddeviations
smethod =mvlue
nochart outindex ='MVLUE -Stds';
mchart Breakstrength*Sample / outlimits=Wirelim5
stddeviations
smethod =rmsdf
nochart outindex ='RMSDF -Stds';
run;
data Wlimits;
set Wirelim1 Wirelim2 Wirelim3 Wirelim4 Wirelim5;
keep _index_ _stddev_;
run;
title 'The Wlimits Data Set';
proc print data=Wlimits noobs;
run;