OPTIONS LINESIZE=80;
TITLE1 'Factorial Design of Metal Part Finish ';
TITLE2 'Factor A (DEPTH): 4 Levels of Depth of Cut (.15, .18, .2, .25 )';
TITLE3 'Factor B (FEED) : 3 Levels of Feed Rate (0.2, 0.25, 0.3 in/min)';
TITLE4 '* FSYM is the alphabetical symbol for Feed Rate values';
TITLE5 ' (0.02 = O, 0.25 = P, 0.30 = Q)';
DATA SURFNS;
/* @@ means a loop in reading input variable sequence */
INPUT DEPTH FEED FSYM $ FINISH @@;
CARDS;
0.15 0.20 O 74 0.15 0.20 O 64 0.15 0.20 O 60
0.15 0.25 P 92 0.15 0.25 P 86 0.15 0.25 P 88
0.15 0.30 Q 99 0.15 0.30 Q 98 0.15 0.30 Q 102
0.18 0.20 O 79 0.18 0.20 O 68 0.18 0.20 O 73
0.18 0.25 P 98 0.18 0.25 P 104 0.18 0.25 P 88
0.18 0.30 Q 104 0.18 0.30 Q 99 0.18 0.30 Q 95
0.20 0.20 O 82 0.20 0.20 O 88 0.20 0.20 O 92
0.20 0.25 P 99 0.20 0.25 P 108 0.20 0.25 P 95
0.20 0.30 Q 108 0.20 0.30 Q 110 0.20 0.30 Q 99
0.25 0.20 O 99 0.25 0.20 O 104 0.25 0.20 O 96
0.25 0.25 P 104 0.25 0.25 P 110 0.25 0.25 P 99
0.25 0.30 Q 114 0.25 0.30 Q 111 0.25 0.30 Q 107
;
RUN;
/* Print the original data set */
PROC PRINT;
RUN;
/* Create Table & tabulate input data set */
/* Nicety, but not absolutely necessary for */
/* the analysis - you may skip */
PROC TABULATE;
CLASS DEPTH FEED ;
VAR FINISH;
TABLE FEED, DEPTH*FINISH/CONDENSE;
KEYLABEL SUM='Y(IJ)';
RUN;
/* General Linear Model (GLM) */
PROC GLM;
CLASS DEPTH FEED FSYM;
/* Dependent = Effects */
MODEL FINISH= DEPTH FEED DEPTH*FEED;
OUTPUT OUT = A2 P = YHAT R = RESID;
RUN;
PROC PLOT DATA = A2;
PLOT YHAT*DEPTH = FSYM/HAXIS = 0.1 TO 0.3 BY 0.05;
PLOT RESID*YHAT;
RUN;
/* Compute Normal Scores of input data, 'A2' */
/* using cum. normal function/Residual */
/* BLOM -> yi = þ-1(ri -3/8)/(n+1/4) */
/* ri = rank of the ith obs */
PROC RANK DATA = A2 NORMAL = BLOM OUT = NPLOT2;
VAR RESID;
RANKS NSCORE;
RUN;
PROC PLOT DATA = NPLOT2;
PLOT NSCORE*RESID;
RUN;
|