1 The SAS System 17:43 Wednesday, March 8, 1995 NOTE: Copyright(c) 1989 by SAS Institute Inc., Cary, NC USA. NOTE: SAS (r) Proprietary Software Release 6.09 TS027 Licensed to UNIVERSITY OF WISCONSIN, Site 0002176014. NOTE: AUTOEXEC processing beginning; file is /usr/sas609/autoexec.sas. NOTE: SAS initialization used: real time 0.26 seconds cpu time 0.10 seconds NOTE: AUTOEXEC processing completed. 1 options ls=79 ps=50; 2 3 /* read data */ 4 data dat; 5 infile 'leprosy.dat' firstobs=2; 6 input drug $ before after; 7 NOTE: The infile 'leprosy.dat' is: File Name=/usr1/y/yandell/.SAS/leprosy.dat, Owner Name=yandell,Group Name=calsstat, Access Permission=rw-r--r--, File Size (bytes)=496 NOTE: 30 records were read from the infile 'leprosy.dat'. The minimum record length was 15. The maximum record length was 15. NOTE: The data set WORK.DAT has 30 observations and 3 variables. NOTE: DATA statement used: real time 0.28 seconds cpu time 0.08 seconds 8 proc glm; /* before and after by drug */ 9 class drug; 10 model before after = drug; 11 means drug; 12 NOTE: The PROCEDURE GLM printed pages 1-4. NOTE: PROCEDURE GLM used: 2 The SAS System 17:43 Wednesday, March 8, 1995 real time 0.75 seconds cpu time 0.04 seconds 13 proc glm data=dat; /* parallel lines */ 14 class drug; 15 model after = before drug; 16 lsmeans drug / stderr pdiff; 17 output out=lepfit p=pafter r=rafter; NOTE: The data set WORK.LEPFIT has 30 observations and 5 variables. NOTE: The PROCEDURE GLM printed pages 5-7. NOTE: PROCEDURE GLM used: real time 0.82 seconds cpu time 0.04 seconds 18 proc plot; /* scatter plot with lines */ 19 plot after*before=drug pafter*before='*' / overlay; 20 plot rafter*pafter=drug; /* diagnostic plot of residuals */ 21 22 NOTE: The PROCEDURE PLOT printed pages 8-9. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.01 seconds 23 proc glm data=dat; /* partial plot adjusted for drug */ 24 class drug; 25 model after = drug; 26 output out=drug r=rafter; NOTE: The data set WORK.DRUG has 30 observations and 4 variables. NOTE: The PROCEDURE GLM printed pages 10-11. NOTE: PROCEDURE GLM used: real time 0.81 seconds cpu time 0.04 seconds 27 proc glm; 28 model rafter = before; 29 output out=parta p=prafter; NOTE: The data set WORK.PARTA has 30 observations and 5 variables. NOTE: The PROCEDURE GLM printed pages 12-13. NOTE: PROCEDURE GLM used: 3 The SAS System 17:43 Wednesday, March 8, 1995 real time 0.78 seconds cpu time 0.02 seconds 30 proc plot data=parta; 31 plot rafter*before=drug prafter*before='*' / overlay; 32 NOTE: The PROCEDURE PLOT printed page 14. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.01 seconds 33 proc glm data=dat; /* partial plot adjusted for before */ 34 class drug; 35 model after = before; 36 output out=before r=rafter p=pafter; NOTE: The data set WORK.BEFORE has 30 observations and 5 variables. NOTE: The PROCEDURE GLM printed pages 15-16. NOTE: PROCEDURE GLM used: real time 0.83 seconds cpu time 0.03 seconds 37 proc glm; 38 class drug pafter; 39 model rafter = drug; 40 output out=partb p=mrafter; NOTE: The data set WORK.PARTB has 30 observations and 6 variables. NOTE: The PROCEDURE GLM printed pages 17-18. NOTE: PROCEDURE GLM used: real time 0.81 seconds cpu time 0.03 seconds 41 proc plot; 42 plot rafter*pafter=drug mrafter*pafter='*' / overlay; 43 NOTE: The PROCEDURE PLOT printed page 19. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.00 seconds 4 The SAS System 17:43 Wednesday, March 8, 1995 44 proc glm data=dat; /* sequential tests */ 45 class drug; 46 model after = drug before / ss1; NOTE: The PROCEDURE GLM printed pages 20-21. NOTE: PROCEDURE GLM used: real time 0.63 seconds cpu time 0.03 seconds 47 proc glm data=dat; 48 class drug; 49 model after = before drug / ss1; 50 NOTE: The PROCEDURE GLM printed pages 22-23. NOTE: PROCEDURE GLM used: real time 0.64 seconds cpu time 0.02 seconds 51 proc glm data=dat; /* adjusted tests */ 52 class drug; 53 model after = before drug / ss3; 54 NOTE: The PROCEDURE GLM printed pages 24-25. NOTE: PROCEDURE GLM used: real time 0.67 seconds cpu time 0.05 seconds 55 proc glm data=dat; /* different slopes by factor */ 56 class drug; 57 model after = before * drug; 58 output out=lepifit p=pafter r=rafter; NOTE: The data set WORK.LEPIFIT has 30 observations and 5 variables. NOTE: The PROCEDURE GLM printed pages 26-27. NOTE: PROCEDURE GLM used: real time 0.83 seconds cpu time 0.08 seconds 59 proc plot; 60 plot after*before=drug pafter*before='*'; 61 5 The SAS System 17:43 Wednesday, March 8, 1995 NOTE: The PROCEDURE PLOT printed pages 28-29. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.01 seconds 62 proc glm data=lepfit; /* check after parallel slope fit */ 63 class drug; 64 model rafter = pafter * drug; 65 output out=leprfit r=rrafter p=prafter; NOTE: The data set WORK.LEPRFIT has 30 observations and 7 variables. NOTE: The PROCEDURE GLM printed pages 30-31. NOTE: PROCEDURE GLM used: real time 0.84 seconds cpu time 0.03 seconds 66 proc plot; 67 plot rafter*pafter=drug prafter*pafter='*' / overlay; NOTE: The PROCEDURE PLOT printed page 32. NOTE: PROCEDURE PLOT used: real time 0.01 seconds cpu time 0.01 seconds NOTE: The SAS System used: real time 9.12 seconds cpu time 0.72 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414