NOTE: Copyright (c) 1999-2001 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 8.2 (TS2M0) Licensed to GEORGE WASHINGTON UNIVERSITY-CAMPUSWIDE-T/R, Site 0001180006. NOTE: This session is executing on the WIN_ME platform. NOTE: SAS initialization used: real time 13.06 seconds 1 *** S_e from old ST210E; 2 3 options formdlim = '-'; 4 5 filename s_elst 'c:\sasjobs\stat210\jl2003\s_elst.pdf'; 6 ods printer pdf file=s_elst; NOTE: Writing ODS PRINTER output to DISK destination "c:\sasjobs\stat210\jl2003\s_elst.pdf", printer "PDF". 7 8 DATA S210E; 9 INPUT DRUG 1 10 DISEASE 2-3 11 ISBP 4-6 12 MISS 7-9 BSBP 10-13; 13 FSBP=BSBP+ISBP; 14 DRUG2 = 0; IF DRUG=2 THEN DRUG2=1; 15 DRUG3 = 0; IF DRUG=3 THEN DRUG3=1; 16 DRUG4 = 0; IF DRUG=4 THEN DRUG4=1; 17 DISEASE2 = 0; IF DISEASE=2 THEN DISEASE2=1; 18 DISEASE3 = 0; IF DISEASE=3 THEN DISEASE3=1; 19 DRDS22 = DRUG2*DISEASE2; DRDS23 = DRUG2*DISEASE3; 20 DRDS32 = DRUG3*DISEASE2; DRDS33 = DRUG3*DISEASE3; 21 DRDS42 = DRUG4*DISEASE2; DRDS43 = DRUG4*DISEASE3; 22 23 CARDS; NOTE: The data set WORK.S210E has 72 observations and 17 variables. NOTE: DATA statement used: real time 0.87 seconds 96 ; 97 98 DATA one; SET S210E; IF MISS=0; 99 NOTE: There were 72 observations read from the data set WORK.S210E. NOTE: The data set WORK.ONE has 58 observations and 17 variables. NOTE: DATA statement used: real time 0.05 seconds 100 PROC SORT; BY DRUG DISEASE; 101 TITLE1 'FACTORIAL STUDY OF EFFECTS OF 4 DRUGS AND 3 DISEASES'; 102 TITLE2 'ON CHANGE IN BLOOD PRESSURE'; 103 TITLE3 'UNBALANCED DESIGN WITH SOME VALUES RANDOMLY DELETED'; 104 NOTE: There were 58 observations read from the data set WORK.ONE. NOTE: The data set WORK.ONE has 58 observations and 17 variables. NOTE: PROCEDURE SORT used: real time 0.16 seconds 105 PROC ANOVA; CLASS DRUG DISEASE; 106 MODEL ISBP=DRUG DISEASE DRUG*DISEASE; 107 MEANS DRUG DISEASE DRUG*DISEASE; 108 TITLE4 'ANALYSES USING ORDINARY (BIASED) ANOVA'; 109 WARNING: PROC ANOVA has determined that the number of observations in each cell is not equal. PROC GLM may be more appropriate. NOTE: PROCEDURE ANOVA used: real time 1.86 seconds 110 PROC GLM; CLASSES DRUG DISEASE; 111 MODEL ISBP=DRUG DISEASE DRUG*DISEASE / SOLUTION SS1 SS2 SS3 SS4 112 XPX INVERSE; 113 CONTRAST 'DRUG 1 VS REST' 114 DRUG -3 1 1 1; 115 CONTRAST 'DRUG 2 VS 3,4' 116 DRUG 0 -2 1 1; 117 CONTRAST 'DRUG 3 VS 4' 118 DRUG 0 0 -1 1; 119 LSMEANS DRUG DISEASE DRUG*DISEASE / E STDERR PDIFF; 120 OUTPUT OUT=PRD P=PRED R=RES; 121 TITLE4 'ANALYZED USING LEAST SQUARES THROUGH GLM'; 122 TITLE5 'ALL EFFECTS FIXED'; 123 TITLE6 'DRUG EFFECT ADDED FIRST'; 124 NOTE: The data set WORK.PRD has 58 observations and 19 variables. NOTE: PROCEDURE GLM used: real time 2.02 seconds 125 DATA C; SET PRD; NOTE: There were 58 observations read from the data set WORK.PRD. NOTE: The data set WORK.C has 58 observations and 19 variables. NOTE: DATA statement used: real time 0.16 seconds 126 PROC PRINT; 127 TITLE4 'EXAMINATION OF RESIDUALS WITHIN CELLS'; NOTE: There were 58 observations read from the data set WORK.C. NOTE: PROCEDURE PRINT used: real time 0.38 seconds 128 PROC UNIVARIATE FREQ NORMAL PLOT; VAR RES; 129 NOTE: PROCEDURE UNIVARIATE used: real time 0.76 seconds 130 PROC GLM data=one; CLASSES DRUG DISEASE; 131 MODEL ISBP=DISEASE DRUG DRUG*DISEASE / ss1 ss2 ss3; 132 TITLE6 'DISEASE EFFECT ADDED FIRST'; 133 NOTE: PROCEDURE GLM used: real time 0.16 seconds 134 PROC GLM data=one; CLASSES DRUG DISEASE; 135 MODEL ISBP=DISEASE DRUG(DISEASE) / SOLUTION; 136 LSMEANS DISEASE DRUG(DISEASE) / STDERR PDIFF; 137 CONTRAST 'DRUG(DISEASE 1)' 138 DRUG(DISEASE) -3 1 1 1 0 0 0 0 0 0 0 0, 139 DRUG(DISEASE) 0 -2 1 1 0 0 0 0 0 0 0 0, 140 DRUG(DISEASE) 0 0 -1 1 0 0 0 0 0 0 0 0; 141 CONTRAST 'DRUG(DISEASE 2)' 142 DRUG(DISEASE) 0 0 0 0 -3 1 1 1 0 0 0 0, 143 DRUG(DISEASE) 0 0 0 0 0 -2 1 1 0 0 0 0, 144 DRUG(DISEASE) 0 0 0 0 0 0 -1 1 0 0 0 0; 145 CONTRAST 'DRUG(DISEASE 3)' 146 DRUG(DISEASE) 0 0 0 0 0 0 0 0 -3 1 1 1, 147 DRUG(DISEASE) 0 0 0 0 0 0 0 0 0 -2 1 1, 148 DRUG(DISEASE) 0 0 0 0 0 0 0 0 0 0 -1 1; 149 TITLE4 'NESTED EFFECTS, ALL EFFECTS FIXED'; 150 NOTE: PROCEDURE GLM used: real time 0.72 seconds 151 PROC REG data=one; 152 MODEL ISBP = DRUG2 DRUG3 DRUG4 DISEASE2 DISEASE3 153 DRDS22 DRDS23 DRDS32 DRDS33 DRDS42 DRDS43; 154 DRUG: TEST DRUG2=0, DRUG3=0, DRUG4=0; 155 DISEASE: TEST DISEASE2=0, DISEASE3=0; 156 INTERAC: TEST DRDS22=0, DRDS23=0, DRDS32=0, DRDS33=0, 157 DRDS42=0, DRDS43=0; 158 TITLE4 'INTERACTION MODEL FIT THROUGH PROC REG'; 159 NOTE: 58 observations read. NOTE: 58 observations used in computations. NOTE: PROCEDURE REG used: real time 0.71 seconds 160 PROC GLM data=one; CLASSES DRUG DISEASE; 161 MODEL ISBP=DRUG DISEASE 162 / SOLUTION SS1 SS2 SS3; 163 LSMEANS DRUG DISEASE / E STDERR PDIFF; 164 TITLE4 'NO-INTERACTION MODEL'; 165 166 RUN; 167 168 ods printer close; NOTE: ODS PRINTER printed 43 pages to c:\sasjobs\stat210\jl2003\s_elst.pdf. 169 run;