DM 'LOG;CLEAR;OUT;CLEAR;';
OPTIONS LS=76 PS=54 NONUMBER NODATE;

*************************************************;
***  REPEATABILITY AND REPRODUCIBILITY STUDY  ***;
*************************************************;

********************************************;
*** NOTE: WE ARE ASSUMING NO INTERACTION ***;
***     BETWEEN PARTS AND OPERATORS      ***;
********************************************;

DATA RRS;
     DO PART = 1 TO 20;
     DO OPERATOR = 1 TO 3;
     DO MSRMNT = 1 TO 2;
        INPUT RESPONSE @@; OUTPUT;
     END; END; END;
CARDS;
21 20  20 20  19 21    24 23  24 24  23 24
20 21  19 21  20 22    27 27  28 26  27 28
19 18  19 18  18 21    23 21  24 21  23 22
22 21  22 24  22 20    19 17  18 20  19 18
24 23  25 23  24 24    25 23  26 25  24 25
21 20  20 20  21 20    18 19  17 19  18 19
23 25  25 25  25 25    24 24  23 25  24 25
29 30  30 28  21 20    26 26  25 26  25 27
20 20  19 20  20 20    19 21  19 19  21 23
25 26  25 24  25 25    19 19  18 17  19 17
;
PROC GLM DATA=RRS;
   CLASS PART OPERATOR ;
   MODEL RESPONSE = PART OPERATOR / SS3;
   RANDOM PART OPERATOR / TEST;
          OUTPUT OUT=DIAG P=PRED R=RESID;
TITLE 'REPEATABILITY AND REPRODUCIBILITY STUDY';

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

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

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

**********************;
*** REMOVE PART 15 ***;
**********************;

DATA RRS; SET RRS; IF PART=15 THEN DELETE;

PROC GLM DATA=RRS;
   CLASS PART OPERATOR ;
   MODEL RESPONSE = PART OPERATOR / SS3;
   RANDOM PART OPERATOR / TEST;
          OUTPUT OUT=DIAG P=PRED R=RESID;
TITLE 'REPEATABILITY AND REPRODUCIBILITY STUDY';

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

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

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