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

OUTPUT OUT=Z2 R=RESID;

TITLE2 'BIOACCUMULATION DATA CONVERTED TO RANKITS';

PROC UNIVARIATE NORMAL;

VAR RESID;

TITLE3 'SHAPIRO-WILKS TEST FOR NORMALITY';

PROC MEANS DATA=A NOPRINT;

BY TRT; VAR RANKIT;

OUTPUT OUT=X MEAN=MEAN;

DATA AX;

MERGE A X; BY TRT;

ABSDEV=ABS(RANKIT-MEAN);

PROC GLM;

CLASS TRT;

MODEL ABSDEV=TRT;

TITLE3 'LEVENE''S TEST';

/* Perform LSD on rankits. */

PROC GLM DATA=A;

CLASS TRT;

MODEL RANKIT=TRT;

MEANS TRT/LSD ALPHA=.1;

TITLE3 'LSD TEST';

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

/* using rankits. */

PROC TTEST DATA=T1;

CLASS TRT;

VAR RANKIT;

PROC TTEST DATA=T2;

CLASS TRT;

VAR RANKIT;

PROC TTEST DATA=T3;

CLASS TRT;

VAR RANKIT;

/* Calculate power of an LSD test to detect true population differences */

/* 10, 25, 50, and 100 % above the reference mean contaminant concentration.

DATA C1;

SET W1;

IF _TYPE_^='ERROR' THEN DELETE;

MSE=SS/DF;

MERGEVAR=1;

KEEP MSE DF MERGEVAR;

DATA C2;

SET Y;

IF TRT^=1 THEN DELETE;

DATA C3;

MERGE C1 C2;

TALPHA=TINV(.95,DF);

LABEL

N='NO. OF REPLICATES, N'

MEANCONC='REFERENCE MEAN CONTAMINANT CONCENTRATION'

MSE='MEAN SQUARE ERROR, MSE'

DF='DEGREES OF FREEDOM, DF'

TALPHA='T VALUE FOR (1-ALPHA=0.95,DF)';

TITLE2 'POWER OF LSD TO DETECT A TRUE POPULATION DIFFERENCE (D)';

TITLE3 'ABOVE REFERENCE MEAN CONTAMINANT CONCENTRATION';

PROC PRINT LABEL NOOBS; VAR N MEANCONC MSE DF TALPHA;

DATA C4;

L51

Appendix L Selected Resource Documents