/****************************************************************/
/*               E S E G   L I B R A R Y                        */
/*                                                              */
/*    NAME: Anvar                                               */
/*   TITLE: Analysis of variance of true defect rate            */
/*          of PBR experiment                                   */
/* PRODUCT: SAS/STAT                                            */
/*  SYSTEM:                                                     */
/*    KEYS: ANOVA GLM REPEATED                                  */
/*   PROCS: PRINT MEANS GLM                                     */
/*    DATA: AllGrSub.dat, AllNASub.dat                          */
/*                                                              */
/*   PURPOSE OF PROGRAM: Perform four analyses of variance      */
/*     (one for each combination of exp. run and problem domain)*/
/*      considering a two-factor experiment with repeated       */
/*      measures in blocks of size two                          */
/*      - grp is the block variable                             */
/*      - teamt and doc are repeated measure factors            */
/*    Subjects are nested within grp                            */
/*                                                              */
/****************************************************************/
options ls=78 ps=60 nodate nonumber;
data allgrsu;
   infile 'AllGrSub.dat';
   input subj $ run $ grp $ prsp $ yexp 
         usutdrt pbrtdrt atmtdrt pgtdrt;

/*Pilot Study in the Generic Problem Domain*/
data allgrpi;
   set allgrsu;
   if run = 'PILOT';
proc print data=allgrpi;
   var subj run grp prsp usutdrt pbrtdrt atmtdrt pgtdrt;
run;
proc means data=allgrpi;
   var usutdrt pbrtdrt atmtdrt pgtdrt;
   run;
proc glm data=allgrpi;
   classes grp;
   model usutdrt pbrtdrt = grp / nouni;
   repeated rtech;
   means grp;
  run;
proc glm data=allgrpi;
   classes grp;
   model atmtdrt pgtdrt = grp / nouni;
   repeated doc;
   means grp;
  run;

/*1995 Run in the Generic problem Domain*/
data allgr95;
   set allgrsu;
   if run = '1995';
proc print data=allgr95;
   var subj run grp prsp usutdrt pbrtdrt atmtdrt pgtdrt;
run;
proc means data=allgr95;
   var usutdrt pbrtdrt atmtdrt pgtdrt;
   run;
proc glm data=allgr95;
   classes grp;
   model usutdrt pbrtdrt = grp / nouni;
   repeated rtech;
   lsmeans grp / stderr;
  run;
proc glm data=allgr95;
   classes grp;
   model atmtdrt pgtdrt = grp / nouni;
   repeated doc;
   lsmeans grp / stderr;
  run;

data allnasu;
   infile 'AllNASub.dat';
   input subj $ run $ grp $ prsp $ yexp 
         usutdrt pbrtdrt atdrt btdrt;

/*Pilot Study in the NASA Problem Domain*/
data allnapi;
   set allnasu;
   if run = 'PILOT';
proc print data=allnapi;
     var subj run grp prsp usutdrt pbrtdrt atdrt btdrt;
run;
proc means data=allnapi;
   var usutdrt pbrtdrt atdrt btdrt;
   run;
proc glm data=allnapi;
   classes grp;
   model usutdrt pbrtdrt = grp / nouni;
   repeated rtech;
   means grp;
  run;
proc glm data=allnapi;
   classes grp;
   model atdrt btdrt = grp / nouni;
   repeated doc;
   means grp;
  run;
 
/*1995 Run in the NASA Domain*/
data allna95;
   set allnasu;
   if run = '1995';
proc print data=allna95;
     var subj run grp prsp usutdrt pbrtdrt atdrt btdrt;
run;
proc means data=allna95;
   var usutdrt pbrtdrt atdrt btdrt;
   run;
proc glm data=allna95;
   classes grp;
   model usutdrt pbrtdrt = grp / nouni;
   repeated rtech;
   lsmeans grp / stderr;
  run;
proc glm data=allna95;
   classes grp;
   model atdrt btdrt = grp / nouni;
   repeated doc;
   lsmeans grp / stderr;
  run;