1 The SAS System 20:58 Sunday, April 2, 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.24 seconds cpu time 0.10 seconds NOTE: AUTOEXEC processing completed. 1 options nocenter ps=64 ls=79; 2 data uneq; 3 infile 'weer.dat'; 4 input trt $ y @@; NOTE: The infile 'weer.dat' is: File Name=/usr1/y/yandell/.SAS/weer.dat, Owner Name=yandell,Group Name=calsstat, Access Permission=rw-r--r--, File Size (bytes)=217 NOTE: 4 records were read from the infile 'weer.dat'. The minimum record length was 48. The maximum record length was 62. NOTE: SAS went to a new line when INPUT statement reached past the end of a line. NOTE: The data set WORK.UNEQ has 31 observations and 2 variables. NOTE: DATA statement used: real time 0.27 seconds cpu time 0.08 seconds 5 proc sort; by trt; 6 NOTE: The data set WORK.UNEQ has 31 observations and 2 variables. NOTE: PROCEDURE SORT used: real time 0.18 seconds cpu time 0.01 seconds 7 proc glm; /* usual analysis of variance */ 8 class trt; 9 model y = trt / ss1; 10 means trt; 11 NOTE: The PROCEDURE GLM printed pages 1-3. NOTE: PROCEDURE GLM used: real time 0.69 seconds cpu time 0.04 seconds 12 proc means noprint data=uneq; by trt; /* get mean and SD */ 13 var y; 2 The SAS System 20:58 Sunday, April 2, 1995 14 output out=means mean=my std=sy; NOTE: The data set WORK.MEANS has 4 observations and 5 variables. NOTE: PROCEDURE MEANS used: real time 0.12 seconds cpu time 0.01 seconds 15 data weight; merge uneq means; by trt; 16 wt = 1 / (sy * sy); /* weights = 1 / variance */ 17 levene = abs(y - my); /* Levene's absolute deviations */ 18 NOTE: The data set WORK.WEIGHT has 31 observations and 8 variables. NOTE: DATA statement used: real time 0.17 seconds cpu time 0.04 seconds 19 proc glm; /* use 1/variance as weight */ 20 class trt; 21 model y = trt / ss1; 22 weight wt; 23 means trt; 24 NOTE: The PROCEDURE GLM printed pages 4-6. NOTE: PROCEDURE GLM used: real time 0.66 seconds cpu time 0.03 seconds 25 proc sort data=uneq; by y; /* replace data by ranks */ NOTE: The data set WORK.UNEQ has 31 observations and 2 variables. NOTE: PROCEDURE SORT used: real time 0.16 seconds cpu time 0.03 seconds 26 data rank; set uneq; ry = _N_; /* rank of y */ NOTE: The data set WORK.RANK has 31 observations and 3 variables. NOTE: DATA statement used: real time 0.16 seconds cpu time 0.03 seconds 27 proc glm; /* Kruskal-Wallis Test */ 28 class trt; 29 model ry = trt / ss1; 30 means trt; 31 NOTE: The PROCEDURE GLM printed pages 7-9. NOTE: PROCEDURE GLM used: real time 0.66 seconds cpu time 0.03 seconds 32 proc glm data=weight; /* Levene's test of Unequal Variances */ 33 class trt; 3 The SAS System 20:58 Sunday, April 2, 1995 34 model levene = trt / ss1; 35 means trt; NOTE: The PROCEDURE GLM printed pages 10-12. NOTE: PROCEDURE GLM used: real time 0.67 seconds cpu time 0.03 seconds NOTE: The SAS System used: real time 4.09 seconds cpu time 0.51 seconds NOTE: SAS Institute Inc., SAS Campus Drive, Cary, NC USA 27513-2414