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 1.03 seconds 1 *** S_C.sas from old ST210C; 2 3 options formdlim = '-'; 4 5 filename s_clst 'c:\sasjobs\stat210\jl2003\s_clst.pdf'; 6 ods printer pdf file=s_clst; NOTE: Writing ODS PRINTER output to DISK destination "c:\sasjobs\stat210\jl2003\s_clst.pdf", printer "PDF". 7 8 DATA ST210C; 9 INPUT DRUG 1 10 DISEASE 2-3 11 ISBP 4-6 12 MISS 7-9; 13 I=DRUG; A=4; 14 J=DISEASE; B=3; 15 Z=J+((I-1)*B); 16 CARDS; NOTE: The data set WORK.ST210C has 72 observations and 9 variables. NOTE: DATA statement used: real time 0.04 seconds 89 ; 90 91 PROC PRINT; var drug disease isbp z; 92 TITLE1 'FACTORIAL STUDY OF EFFECTS OF 4 DRUGS AND 3 DISEASES'; 93 TITLE2 'ON CHANGE IN BLOOD PRESSURE'; 94 NOTE: There were 72 observations read from the data set WORK.ST210C. NOTE: PROCEDURE PRINT used: real time 0.10 seconds 95 PROC ANOVA; CLASS DRUG DISEASE; 96 MODEL ISBP=DRUG DISEASE DRUG*DISEASE; 97 MEANS DRUG DISEASE DRUG*DISEASE; 98 TITLE3 'COMPLETELY BALANCED ANOVA'; 99 NOTE: PROCEDURE ANOVA used: real time 0.21 seconds 100 PROC GLM; CLASS DRUG DISEASE; 101 MODEL ISBP=DRUG DISEASE DRUG*DISEASE / SOLUTION; 102 LSMEANS DRUG DISEASE DRUG*DISEASE / STDERR PDIFF; 103 OUTPUT OUT=PRD P=PRED R=RES; 104 TITLE3 'IMPLEMENTED USING LEAST SQUARES IN GLM'; 105 NOTE: The data set WORK.PRD has 72 observations and 11 variables. NOTE: PROCEDURE GLM used: real time 0.49 seconds 106 DATA C; SET PRD; 107 IF DRUG > . AND DISEASE > .; NOTE: There were 72 observations read from the data set WORK.PRD. NOTE: The data set WORK.C has 72 observations and 11 variables. NOTE: DATA statement used: real time 0.05 seconds 108 PROC PRINT; var drug disease isbp pred res; 109 TITLE4 'EXAMINATION OF RESIDUALS WITHIN CELLS'; NOTE: There were 72 observations read from the data set WORK.C. NOTE: PROCEDURE PRINT used: real time 0.05 seconds 110 PROC UNIVARIATE NORMAL PLOT; VAR RES; 111 NOTE: PROCEDURE UNIVARIATE used: real time 0.05 seconds 112 PROC glm; class z; 113 model isbp=z; 114 means Z / HOVTEST=BARTLETT; 115 TITLE3 'TEST OF HOMOSCEDASTICITY ASSUMPTION OVER CELLS'; 116 NOTE: PROCEDURE GLM used: real time 0.50 seconds 117 PROC GLM; CLASS DRUG DISEASE; 118 MODEL ISBP=DISEASE DRUG(DISEASE) / SOLUTION; 119 LSMEANS DISEASE DRUG(DISEASE) / STDERR PDIFF; 120 CONTRAST 'DRUG(DISEASE 1)' 121 DRUG(DISEASE) -3 1 1 1 0 0 0 0 0 0 0 0, 122 DRUG(DISEASE) 0 -2 1 1 0 0 0 0 0 0 0 0, 123 DRUG(DISEASE) 0 0 -1 1 0 0 0 0 0 0 0 0; 124 CONTRAST 'DRUG(DISEASE 2)' 125 DRUG(DISEASE) 0 0 0 0 -3 1 1 1 0 0 0 0, 126 DRUG(DISEASE) 0 0 0 0 0 -2 1 1 0 0 0 0, 127 DRUG(DISEASE) 0 0 0 0 0 0 -1 1 0 0 0 0; 128 CONTRAST 'DRUG(DISEASE 3)' 129 DRUG(DISEASE) 0 0 0 0 0 0 0 0 -3 1 1 1, 130 DRUG(DISEASE) 0 0 0 0 0 0 0 0 0 -2 1 1, 131 DRUG(DISEASE) 0 0 0 0 0 0 0 0 0 0 -1 1; 132 TITLE4 'NESTED EFFECTS FOR DRUG GROUP WITHIN DISEASE CATEGORY'; 133 NOTE: PROCEDURE GLM used: real time 0.37 seconds 134 PROC GLM; CLASS DRUG DISEASE; 135 MODEL ISBP=DRUG DISEASE / SOLUTION; 136 LSMEANS DRUG DISEASE / STDERR PDIFF; 137 MEANS DRUG DISEASE / SIDAK TUKEY; 138 TITLE3 'MAIN EFFECTS ONLY MODEL'; 139 140 RUN; NOTE: Means from the MEANS statement are not adjusted for other terms in the model. For adjusted means, use the LSMEANS statement. 141 142 ods printer close; NOTE: ODS PRINTER printed 32 pages to c:\sasjobs\stat210\jl2003\s_clst.pdf. 143 run;