OPTIONS NODATE NONUMBER PAGESIZE=66 LINESIZE=72;

********************************;
***  GAUGE CAPABILITY STUDY  ***;
********************************;

DATA RRS;
     DO PART = 1 TO 20;
     DO MSRMNT = 1 TO 2;
        INPUT RESPONSE @@; OUTPUT;
     END; END;
CARDS;
21 20    24 23    20 21    27 27  
19 18    23 21    22 21    19 17  
24 23    25 23    21 20    18 19  
23 25    24 24    29 30    26 26  
20 20    19 21    25 26    19 19  
;
PROC GLM DATA=RRS;
   CLASS PART ;
   MODEL RESPONSE = PART / SS3;
         OUTPUT OUT=DIAG P=PRED R=RESID;
   ID MSRMNT;
TITLE 'GAUGE CAPABILITY STUDY';
TITLE2 'EXAMPLE 9-7 (MONTGOMERY PAGES 455-456)';

PROC PLOT DATA=DIAG VPERCENT=50;
     PLOT RESID*PART=MSRMNT / VREF=0;

************************************;
*** VARIANCE COMPONENTS ANALYSIS ***;
************************************;

PROC VARCOMP DATA=RRS METHOD=REML;
     CLASS PART ;
     MODEL RESPONSE = PART ;

RUN;