L.4.2.1 BIOACC.SAS program statements (Cont.)

PROC GLM NOPRINT DATA=A;

CLASS TRT;

MODEL CONC LOGCONC=TRT;

OUTPUT OUT=Z R=RESID RESIDLOG;

PROC UNIVARIATE NORMAL;

VAR RESID RESIDLOG;

TITLE2 'SHAPIRO-WILKS TEST FOR NORMALITY';

/* Conduct Levene's Test for equality of variances of untransformed and */

/* log-transformed data. */

DATA AY;

MERGE A Y; BY TRT;

ABSDEV=ABS(CONC-MEANCONC);

ABSLOG=ABS(LOGCONC-MEANLOG);

LABEL

ABSDEV='ABSOLUTE DEVIATIONS FROM MEAN CONC.'

ABSLOG='ABSOLUTE DEVIATIONS FROM MEAN LOGCONC.';

PROC GLM;

CLASS TRT;

MODEL ABSDEV ABSLOG=TRT;

TITLE2 'LEVENE''S TEST';

/* Perform LSD on untransformed and log-transformed data. */

PROC GLM DATA=A OUTSTAT=W1;

CLASS TRT;

MODEL CONC=TRT;

MEANS TRT/LSD ALPHA=.1;

TITLE2 'LSD TEST (UNTRANSFORMED DATA)';

PROC GLM DATA=A OUTSTAT=W2;

CLASS TRT;

MODEL LOGCONC=TRT;

MEANS TRT/LSD ALPHA=.1;

TITLE2 'LSD TEST (LOG-TRANSFORMED DATA)';

/* Perform t-tests for each dredged sediment-reference sediment comparison */

/* using untransformed and log-transformed data. */

DATA T1;

SET A;

IF TRT>2 THEN DELETE;

PROC TTEST;

CLASS TRT;

VAR CONC LOGCONC;

DATA T2;

SET A;

IF TRT=2 OR TRT=4 THEN DELETE;

PROC TTEST;

CLASS TRT;

VAR CONC LOGCONC;

DATA T3;

SET A;

IF TRT=2 OR TRT=3 THEN DELETE;

PROC TTEST;

CLASS TRT;

VAR CONC LOGCONC;

/* Test normality and equality of variances of rankits. */

PROC GLM NOPRINT DATA=A;

CLASS TRT;

MODEL RANKIT=TRT;

L50

Appendix L Selected Resource Documents