1 The SAS System 13:46 Thursday, April 30, 1998 NOTE: Copyright (c) 1989-1996 by SAS Institute Inc., Cary, NC, USA. NOTE: SAS (r) Proprietary Software Release 6.11 TS040 Licensed to UNIVERSITY OF WISCONSIN, Site 0002176032. This message is contained in the SAS news file, and is presented upon initialization. Edit the files "news" in the "misc/base" directory to display site-specific news and information in the program log. The command line option "-nonews" will prevent this display. NOTE: AUTOEXEC processing beginning; file is /afs/.cs.wisc.edu/s/sas-6.11/alpha_osf20/sas611/autoexec.sas. NOTE: SAS initialization used: real time 9.60 seconds cpu time 0.90 seconds NOTE: AUTOEXEC processing completed. 1 /* Example of Three-Period Crossover MJ 32.3 */ 2 options nocenter linesize=80 pagesize=50; 3 data cross; 4 infile 'cross33.dat' firstobs = 2; 5 input seq period trt $ u1 u2 u3 u4 u5 u6; NOTE: The infile 'cross33.dat' is: File Name=/afs/cs.wisc.edu/p/stat/Data/MJ/cross33.dat, Owner Name=yandell,Group Name=6013, Access Permission=rw-r-----, File Size (bytes)=689 NOTE: 19 records were read from the infile 'cross33.dat'. The minimum record length was 0. The maximum record length was 35. NOTE: SAS went to a new line when INPUT statement reached past the end of a line. NOTE: The data set WORK.CROSS has 18 observations and 9 variables. NOTE: DATA statement used: real time 3.98 seconds cpu time 0.41 seconds 6 proc sort; by seq period; 2 The SAS System 13:46 Thursday, April 30, 1998 NOTE: The data set WORK.CROSS has 18 observations and 9 variables. NOTE: PROCEDURE SORT used: real time 0.51 seconds cpu time 0.06 seconds 7 data b; set cross; 8 if seq = 1 or seq = 2 then do; 9 carry = '0'; LA = 0; LB = 0; LC = 0; 10 if period = 2 then do; carry = 'A'; LA = 1; end; 11 if period = 3 then do; LA = 0; 12 if seq = 1 then do; carry = 'B'; LB = 1; end; 13 if seq = 2 then do; carry = 'C'; LC = 1; end; 14 end; 15 end; 16 if seq = 3 or seq = 4 then do; 17 carry = '0'; LA = 0; LB = 0; LC = 0; 18 if period = 2 then do; carry = 'B'; LB = 1; end; 19 if period = 3 then do; LB = 0; 20 if seq = 3 then do; carry = 'A'; LA = 1; end; 21 if seq = 4 then do; carry = 'C'; LC = 1; end; 22 end; 23 end; 24 if seq = 5 or seq = 6 then do; 25 carry = '0'; LA = 0; LB = 0; LC = 0; 26 if period = 2 then do; carry = 'C'; LC = 1; end; 27 if period = 3 then do; LC = 0; 28 if seq = 5 then do; carry = 'A'; LA = 1; end; 29 if seq = 6 then do; carry = 'B'; LB = 1; end; 30 end; 31 end; 32 wpcarry = 'AB'; 33 if seq = 2 or seq = 5 then wpcarry = 'AC'; 34 if seq = 4 or seq = 6 then wpcarry = 'BC'; 35 y = u1; person = 1; output; 36 y = u2; person = 2; output; 37 y = u3; person = 3; output; 38 y = u4; person = 4; output; 39 y = u5; person = 5; output; 40 y = u6; person = 6; output; 41 drop u1 u2 u3 u4 u5 u6; NOTE: The data set WORK.B has 108 observations and 10 variables. NOTE: DATA statement used: real time 0.60 seconds cpu time 0.26 seconds 3 The SAS System 13:46 Thursday, April 30, 1998 42 proc print; NOTE: The PROCEDURE PRINT printed pages 1-3. NOTE: PROCEDURE PRINT used: real time 1.20 seconds cpu time 0.30 seconds 43 proc glm; 44 class seq wpcarry trt period carry person; 45 model y = person(seq) trt period carry / solution; 46 lsmeans trt period / stderr pdiff; 47 NOTE: The PROCEDURE GLM printed pages 4-7. NOTE: PROCEDURE GLM used: real time 3.50 seconds cpu time 0.65 seconds 48 proc mixed data = b; 49 class trt wpcarry seq period carry person; 50 model y = trt period carry / ddfm = satterth; 51 random person(seq); 52 lsmeans trt carry / pdiff; NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed pages 8-9. NOTE: PROCEDURE MIXED used: real time 5.38 seconds cpu time 0.76 seconds 53 proc mixed data = b; 54 class seq wpcarry trt period person; 55 model y = trt LA LB LC period / ddfm = satterth; 56 random person(seq); 57 lsmeans trt period / pdiff; 58 59 /* now include whole plot carry */ NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed pages 10-11. NOTE: PROCEDURE MIXED used: real time 0.53 seconds cpu time 0.25 seconds 4 The SAS System 13:46 Thursday, April 30, 1998 60 proc glm; 61 class seq wpcarry trt period carry person; 62 model y = wpcarry person(seq*wpcarry) trt period carry / solution; 63 lsmeans trt period / stderr pdiff; NOTE: The PROCEDURE GLM printed pages 12-16. NOTE: PROCEDURE GLM used: real time 1.30 seconds cpu time 0.30 seconds 64 proc mixed data = b; 65 class trt wpcarry seq period carry person; 66 model y = wpcarry trt period carry / ddfm = satterth; 67 random person(seq*wpcarry); 68 lsmeans trt carry / pdiff; NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed pages 17-18. NOTE: PROCEDURE MIXED used: real time 0.61 seconds cpu time 0.25 seconds 69 proc mixed data = b; 70 class seq wpcarry trt period person; 71 model y = wpcarry trt LA LB LC period / ddfm = satterth; 72 random person(seq*wpcarry); 73 lsmeans trt period / pdiff; NOTE: At least one W.D format was too small for the number to be printed. The decimal may be shifted by the "BEST" format. NOTE: The PROCEDURE MIXED printed pages 19-20. NOTE: PROCEDURE MIXED used: real time 1.05 seconds cpu time 0.25 seconds NOTE: The SAS System used: real time 30.53 seconds cpu time 5.05 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414